Predicting Lipid-Rich Plaque Progression in Coronary Arteries Using Multimodal Imaging and Wall Shear Stress Signatures.

BACKGROUND: Plaque composition and wall shear stress (WSS) magnitude act as well-established players in coronary plaque progression. However, WSS magnitude per se does not completely capture the mechanical stimulus to which the endothelium is subjected, since endothelial cells experience changes in the WSS spatiotemporal configuration on the luminal surface. This study explores WSS profile and lipid content signatures of plaque progression to identify novel biomarkers of coronary atherosclerosis. METHODS: Thirty-seven patients with acute coronary syndrome underwent coronary computed tomography angiography, near-infrared spectroscopy intravascular ultrasound, and optical coherence tomography of at least 1 nonculprit vessel at baseline and 1-year follow-up. Baseline coronary artery geometries were reconstructed from intravascular ultrasound and coronary computed tomography angiography and combined with flow information to perform computational fluid dynamics simulations to assess the time-averaged WSS magnitude (TAWSS) and the variability in the contraction/expansion action exerted by WSS on the endothelium, quantifiable in terms of topological shear variation index (TSVI). Plaque progression was measured as intravascular ultrasound-derived percentage plaque atheroma volume change at 1-year follow-up. Plaque composition information was extracted from near-infrared spectroscopy and optical coherence tomography. RESULTS: Exposure to high TSVI and low TAWSS was associated with higher plaque progression (4.00±0.69% and 3.60±0.62%, respectively). Plaque composition acted synergistically with TSVI or TAWSS, resulting in the highest plaque progression (≥5.90%) at locations where lipid-rich plaque is exposed to high TSVI or low TAWSS. CONCLUSIONS: Luminal exposure to high TSVI, solely or combined with a lipid-rich plaque phenotype, is associated with enhanced plaque progression at 1-year follow-up. Where plaque progression occurred, low TAWSS was also observed. These findings suggest TSVI, in addition to low TAWSS, as a potential biomechanical predictor for plaque progression, showing promise for clinical translation to improve patient prognosis.

Current and Future Applications of Computational Fluid Dynamics in Coronary Artery Disease.

Hemodynamics interacts with the cellular components of human vessels, influencing function and healthy status. Locally acting hemodynamic forces have been associated-by a steadily increasing amount of scientific evidence-with nucleation and evolution of atherosclerotic plaques in several vascular regions, resulting in the formulation of the 'hemodynamic risk hypothesis' of the atherogenesis. At the level of coronary arteries, however, the complexity of both anatomy and physiology made the study of this vascular region particularly difficult for researchers. Developments in computational fluid dynamics (CFD) have recently allowed an accurate modelling of the intracoronary hemodynamics, thus offering physicians a unique tool for the investigation of this crucial human system by means of advanced mathematical simulations. The present review of CFD applications in coronary artery disease was set to concisely offer the medical reader the theoretical foundations of quantitative intravascular hemodynamics-reasoned schematically in the text in its basic (i.e., pressure and velocity) and derived quantities (e.g., fractional flow reserve, wall shear stress and helicity)-along with its current implications in clinical research. Moreover, attention was paid in classifying computational modelling derived from invasive and non-invasive imaging modalities with unbiased remarks on the advantages and limitations of each procedure. Finally, an extensive description-aided by explanatory figures and cross references to recent clinical findings-was presented on the role of near-wall hemodynamics, in terms of shear stress, and of intravascular flow complexity, in terms of helical flow.

Coronary Artery Stenting Affects Wall Shear Stress Topological Skeleton.

Despite the important advancements in the stent technology for the treatment of diseased coronary arteries, major complications still affect the postoperative long-term outcome. The stent-induced flow disturbances, and especially the altered wall shear stress (WSS) profile at the strut level, play an important role in the pathophysiological mechanisms leading to stent thrombosis (ST) and in-stent restenosis (ISR). In this context, the analysis of the WSS topological skeleton is gaining more and more interest by extending the current understanding of the association between local hemodynamics and vascular diseases. This study aims to analyze the impact that a deployed coronary stent has on the WSS topological skeleton. Computational fluid dynamics (CFD) simulations were performed in three stented human coronary artery geometries reconstructed from clinical images. The selected cases presented stents with different designs (i.e., two contemporary drug-eluting stents and one bioresorbable scaffold) and included regions with stent malapposition or overlapping. A recently proposed Eulerian-based approach was applied to analyze the WSS topological skeleton features. The results highlighted that the presence of single or multiple stents within a coronary artery markedly impacts the WSS topological skeleton. In particular, repetitive patterns of WSS divergence were observed at the luminal surface, highlighting a WSS contraction action exerted proximal to the stent struts and a WSS expansion action distal to the stent struts. This WSS action pattern was independent from the stent design. In conclusion, these findings could contribute to a deeper understanding of the hemodynamics-driven processes underlying ST and ISR.

Smartphone-based particle tracking velocimetry for the in vitro assessment of coronary flows.

The present study adopts a smartphone-based approach for the experimental characterization of coronary flows. Technically, Particle Tracking Velocimetry (PTV) measurements were performed using a smartphone camera and a low-power continuous wave laser in realistic healthy and stenosed phantoms of left anterior descending artery with inflow Reynolds numbers approximately ranging from 20 to 200. A Lagrangian-Eulerian mapping was performed to convert Lagrangian PTV velocity data to a Eulerian grid. Eulerian velocity and vorticity data obtained from smartphone-based PTV measurements were compared with Particle Image Velocimetry (PIV) measurements performed with a smartphone-based setup and with a conventional setup based on a high-power double-pulsed laser and a CMOS camera. Smartphone-based PTV and PIV velocity flow fields substantially agreed with conventional PIV measurements, with the former characterized by lower average percentage differences than the latter. Discrepancies emerged at high flow regimes, especially at the stenosis throat, due to particle image blur generated by smartphone camera shutter speed and image acquisition frequency. In conclusion, the present findings demonstrate the feasibility of PTV measurements using a smartphone camera and a low-power light source for the in vitro characterization of cardiovascular flows for research, industrial and educational purposes, with advantages in terms of costs, safety and usability.

Modelling blood flow in coronary arteries: Newtonian or shear-thinning non-Newtonian rheology?

BACKGROUND: The combination of medical imaging and computational hemodynamics is a promising technology to diagnose/prognose coronary artery disease (CAD). However, the clinical translation of in silico hemodynamic models is still hampered by assumptions/idealizations that must be introduced in model-based strategies and that necessarily imply uncertainty. This study aims to provide a definite answer to the open question of how to properly model blood rheological properties in computational fluid dynamics (CFD) simulations of coronary hemodynamics. METHODS: The geometry of the right coronary artery (RCA) of 144 hemodynamically stable patients with different stenosis degree were reconstructed from angiography. On them, unsteady-state CFD simulations were carried out. On each reconstructed RCA two different simulation strategies were applied to account for blood rheological properties, implementing (i) a Newtonian (N) and (ii) a shear-thinning non-Newtonian (non-N) rheological model. Their impact was evaluated in terms of wall shear stress (WSS magnitude, multidirectionality, topological skeleton) and helical flow (strength, topology) profiles. Additionally, luminal surface areas (SAs) exposed to shear disturbances were identified and the co-localization of paired N and non-N SAs was quantified in terms of similarity index (SI). RESULTS: The comparison between paired N vs. shear-thinning non-N simulations revealed remarkably similar profiles of WSS-based and helicity-based quantities, independent of the adopted blood rheology model and of the degree of stenosis of the vessel. Statistically, for each paired N and non-N hemodynamic quantity emerged negligible bias from Bland-Altman plots, and strong positive linear correlation (r > 0.94 for almost all the WSS-based quantities, r > 0.99 for helicity-based quantities). Moreover, a remarkable co-localization of N vs. non-N luminal SAs exposed to disturbed shear clearly emerged (SI distribution 0.95 [0.93, 0.97]). Helical flow topology resulted to be unaffected by blood rheological properties. CONCLUSIONS: This study, performed on 288 angio-based CFD simulations on 144 RCA models presenting with different degrees of stenosis, suggests that the assumptions on blood rheology have negligible impact both on WSS and helical flow profiles associated with CAD, thus definitively answering to the question "is Newtonian assumption for blood rheology adequate in coronary hemodynamics simulations?".

Impact of wall displacements on the large-scale flow coherence in ascending aorta.

In the context of aortic hemodynamics, uncertainties affecting blood flow simulations hamper their translational potential as supportive technology in clinics. Computational fluid dynamics (CFD) simulations under rigid-walls assumption are largely adopted, even though the aorta contributes markedly to the systemic compliance and is characterized by a complex motion. To account for personalized wall displacements in aortic hemodynamics simulations, the moving-boundary method (MBM) has been recently proposed as a computationally convenient strategy, although its implementation requires dynamic imaging acquisitions not always available in clinics. In this study we aim to clarify the real need for introducing aortic wall displacements in CFD simulations to accurately capture the large-scale flow structures in the healthy human ascending aorta (AAo). To do that, the impact of wall displacements is analyzed using subject-specific models where two CFD simulations are performed imposing (1) rigid walls, and (2) personalized wall displacements adopting a MBM, integrating dynamic CT imaging and a mesh morphing technique based on radial basis functions. The impact of wall displacements on AAo hemodynamics is analyzed in terms of large-scale flow patterns of physiological significance, namely axial blood flow coherence (quantified applying the Complex Networks theory), secondary flows, helical flow and wall shear stress (WSS). From the comparison with rigid-wall simulations, it emerges that wall displacements have a minor impact on the AAo large-scale axial flow, but they can affect secondary flows and WSS directional changes. Overall, helical flow topology is moderately affected by aortic wall displacements, whereas helicity intensity remains almost unchanged. We conclude that CFD simulations with rigid-wall assumption can be a valid approach to study large-scale aortic flows of physiological significance.

Wall shear stress-related plaque growth of lipid-rich plaques in human coronary arteries: an near-infrared spectroscopy and optical coherence tomography study.

AIMS: Low wall shear stress (WSS) is acknowledged to play a role in plaque development through its influence on local endothelial function. Also, lipid-rich plaques (LRPs) are associated with endothelial dysfunction. However, little is known about the interplay between WSS and the presence of lipids with respect to plaque progression. Therefore, we aimed to study the differences in WSS-related plaque progression between LRPs, non-LRPs, or plaque-free regions in human coronary arteries. METHODS AND RESULTS: In the present single-centre, prospective study, 40 patients who presented with an acute coronary syndrome successfully underwent near-infrared spectroscopy intravascular ultrasound (NIRS-IVUS) and optical coherence tomography (OCT) of at least one non-culprit vessel at baseline and completed a 1-year follow-up. WSS was computed applying computational fluid dynamics to a three-dimensional reconstruction of the coronary artery based on the fusion of the IVUS-segmented lumen with a CT-derived centreline, using invasive flow measurements as boundary conditions. For data analysis, each artery was divided into 1.5 mm/45° sectors. Plaque growth based on IVUS-derived percentage atheroma volume change was compared between LRPs, non-LRPs, and plaque-free wall segments, as assessed by both OCT and NIRS. Both NIRS- and OCT-detected lipid-rich sectors showed a significantly higher plaque progression than non-LRPs or plaque-free regions. Exposure to low WSS was associated with a higher plaque progression than exposure to mid or high WSS, even in the regions classified as a plaque-free wall. Furthermore, low WSS and the presence of lipids had a synergistic effect on plaque growth, resulting in the highest plaque progression in lipid-rich regions exposed to low shear stress. CONCLUSION: This study demonstrates that NIRS- and OCT-detected lipid-rich regions exposed to low WSS are subject to enhanced plaque growth over a 1-year follow-up. The presence of lipids and low WSS proves to have a synergistic effect on plaque growth.

Divergence of the normalized wall shear stress as an effective computational template of low-density lipoprotein polarization at the arterial blood-vessel wall interface.

BACKGROUND AND OBJECTIVE: Near-wall transport of low-density lipoproteins (LDL) in arteries plays a relevant role in the initiation of atherosclerosis. Although it can be modelled in silico by coupling the Navier-Stokes equations with the 3D advection-diffusion (AD) equation, the associated computational cost is high. As wall shear stress (WSS) represents a first-order approximation of the near-wall velocity in arteries, we aimed at identifying computationally convenient WSS-based quantities to infer LDL near-wall transport based on the underlying near-wall hemodynamics in five models of three human arterial districts (aorta, carotid bifurcations, coronary arteries). The simulated LDL transport and its WSS-based surrogates were qualitatively compared with in vivo longitudinal measurements of wall thickness growth on the coronary artery models. METHODS: Numerical simulations of blood flow coupled with AD equations for LDL transport and blood-wall transfer were performed. The co-localization of the simulated LDL concentration polarization patterns with luminal surface areas characterized by low cycle-average WSS, near-wall flow stagnation and WSS attracting patterns was quantitatively assessed by the similarity index (SI). In detail, the latter two represent features of the WSS topological skeleton, obtained respectively through the Lagrangian tracking of surface-born particles, and the Eulerian analysis of the divergence of the normalized cycle-average WSS vector field. RESULTS: Convergence of the solution of the AD problem required the simulation of 3 (coronary artery) to 10 (aorta) additional cardiac cycles with respect to the Navier-Stokes problem. Co-localization results underlined that WSS topological skeleton features indicating near-wall flow stagnation and WSS attracting patterns identified LDL concentration polarization profiles more effectively than low WSS, as indicated by higher SI values (SI range: 0.17-0.50 for low WSS; 0.24-0.57 for WSS topological skeleton features). Moreover, the correspondence between the simulated LDL uptake and WSS-based quantities profiles with the in vivo measured wall thickness growth in coronary arteries appears promising. CONCLUSIONS: The recently introduced Eulerian approach for identifying WSS attracting patterns from the divergence of normalized WSS provides a computationally affordable template of the LDL polarization at the arterial blood-wall interface without simulating the AD problem. It thus candidates as an effective biomechanical tool for elucidating the mechanistic link amongst LDL transfer at the arterial blood-wall interface, WSS and atherogenesis.

Early Atherosclerotic Changes in Coronary Arteries are Associated with Endothelium Shear Stress Contraction/Expansion Variability.

Although unphysiological wall shear stress (WSS) has become the consensus hemodynamic mechanism for coronary atherosclerosis, the complex biomechanical stimulus affecting atherosclerosis evolution is still undetermined. This has motivated the interest on the contraction/expansion action exerted by WSS on the endothelium, obtained through the WSS topological skeleton analysis. This study tests the ability of this WSS feature, alone or combined with WSS magnitude, to predict coronary wall thickness (WT) longitudinal changes. Nine coronary arteries of hypercholesterolemic minipigs underwent imaging with local WT measurement at three time points: baseline (T1), after 5.6 ± 0.9 (T2), and 7.6 ± 2.5 (T3) months. Individualized computational hemodynamic simulations were performed at T1 and T2. The variability of the WSS contraction/expansion action along the cardiac cycle was quantified using the WSS topological shear variation index (TSVI). Alone or combined, high TSVI and low WSS significantly co-localized with high WT at the same time points and were significant predictors of thickening at later time points. TSVI and WSS magnitude values in a physiological range appeared to play an atheroprotective role. Both the variability of the WSS contraction/expansion action and WSS magnitude, accounting for different hemodynamic effects on the endothelium, (1) are linked to WT changes and (2) concur to identify WSS features leading to coronary atherosclerosis.

Comparison of Swine and Human Computational Hemodynamics Models for the Study of Coronary Atherosclerosis.

Coronary atherosclerosis is a leading cause of illness and death in Western World and its mechanisms are still non completely understood. Several animal models have been used to 1) study coronary atherosclerosis natural history and 2) propose predictive tools for this disease, that is asymptomatic for a long time, aiming for a direct translation of their findings to human coronary arteries. Among them, swine models are largely used due to the observed anatomical and pathophysiological similarities to humans. However, a direct comparison between swine and human models in terms of coronary hemodynamics, known to influence atherosclerotic onset/development, is still lacking. In this context, we performed a detailed comparative analysis between swine- and human-specific computational hemodynamic models of coronary arteries. The analysis involved several near-wall and intravascular flow descriptors, previously emerged as markers of coronary atherosclerosis initiation/progression, as well as anatomical features. To do that, non-culprit coronary arteries (18 right-RCA, 18 left anterior descending-LAD, 13 left circumflex-LCX coronary artery) from patients presenting with acute coronary syndrome were imaged by intravascular ultrasound and coronary computed tomography angiography. Similarly, the three main coronary arteries of ten adult mini-pigs were also imaged (10 RCA, 10 LAD, 10 LCX). The geometries of the imaged coronary arteries were reconstructed (49 human, 30 swine), and computational fluid dynamic simulations were performed by imposing individualized boundary conditions. Overall, no relevant differences in 1) wall shear stress-based quantities, 2) intravascular hemodynamics (in terms of helical flow features), and 3) anatomical features emerged between human- and swine-specific models. The findings of this study strongly support the use of swine-specific computational models to study and characterize the hemodynamic features linked to coronary atherosclerosis, sustaining the reliability of their translation to human vascular disease.

The definition of low wall shear stress and its effect on plaque progression estimation in human coronary arteries.

Wall shear stress (WSS), the frictional force of the blood on the vessel wall, plays a crucial role in atherosclerotic plaque development. Low WSS has been associated with plaque growth, however previous research used different approaches to define low WSS to investigate its effect on plaque progression. In this study, we used four methodologies to allocate low, mid and high WSS in one dataset of human coronary arteries and investigated the predictive power of low WSS for plaque progression. Coronary reconstructions were based on multimodality imaging, using intravascular ultrasound and CT-imaging. Vessel-specific flow was measured using Doppler wire and computational fluid dynamics was performed to calculate WSS. The absolute WSS range varied greatly between the coronary arteries. On the population level, the established pattern of most plaque progression at low WSS was apparent in all methodologies defining the WSS categories. However, for the individual patient, when using measured flow to determine WSS, the absolute WSS values range so widely, that the use of absolute thresholds to determine low WSS was not appropriate to identify regions at high risk for plaque progression.

Lipid-rich Plaques Detected by Near-infrared Spectroscopy Are More Frequently Exposed to High Shear Stress.

High wall shear stress (WSS) and near-infrared spectroscopy (NIRS) detected lipid-rich plaque (LRP) are both known to be associated with plaque destabilization and future adverse cardiovascular events. However, knowledge of spatial co-localization of LRP and high WSS is lacking. This study investigated the co-localization of LRP based on NIRS and high WSS. Fifty-three patients presenting acute coronary syndrome underwent NIRS-intravascular-ultrasound (NIRS-IVUS) imaging of a non-culprit coronary artery. WSS was obtained using WSS profiling in 3D-reconstructions of the coronary arteries based on fusion of IVUS-segmented lumen and CT-derived 3D-centerline. Thirty-eight vessels were available for final analysis and divided into 0.5 mm/45° sectors. LRP sectors, as identified by NIRS, were more often colocalized with high WSS than sectors without LRP. Moreover, there was a dose-dependent relationship between lipid content and high WSS exposure. This study is a first step in understanding the evolution of LRPs to vulnerable plaques. Graphical Abstract.

Finite Element Modeling Application in Forensic Practice: A Periprosthetic Femoral Fracture Case Study.

The incidence of periprosthetic fractures has rapidly increased in the last two decades and has been the cause of a large number of revision surgeries and permanent physical disability for many patients, as well as a significant socioeconomic burden for many nations. This research deals with a periprosthetic femur fracture real event, occurred following a total hip arthroplasty and treated with one of the most widespread internal fixation methods: the implant of a periprosthetic femur plate system. A Finite Element analysis was performed to investigate the implanted femur plate break after a short follow-up and to understand the plate break causes. Such events are currently object of forensic debate as more and more often hospitals, surgeons, and medical device manufacturers are denounced by patients to whom similar events occur. In this work, different load situations acting on the femur during daily and incidental activities were simulated, in order to validate the correct behavior of the plate, according to the intended use recommended by the manufacturer. The analysis demonstrates that the plate failure can occur in situations of unconventional loading such as that caused by stumbling and in presence of incomplete bone healing.

Wall Shear Stress Topological Skeleton Independently Predicts Long-Term Restenosis After Carotid Bifurcation Endarterectomy.

Wall Shear Stress (WSS) topological skeleton, composed by fixed points and the manifolds linking them, reflects the presence of blood flow features associated to adverse vascular response. However, the influence of WSS topological skeleton on vascular pathophysiology is still underexplored. This study aimed to identify direct associations between the WSS topological skeleton and markers of vascular disease from real-world clinical longitudinal data of long-term restenosis after carotid endarterectomy (CEA). Personalized computational hemodynamic simulations were performed on a cohort of 13 carotid models pre-CEA and at 1 month after CEA. At 60 months after CEA, intima-media thickness (IMT) was measured to detect long-term restenosis. The analysis of the WSS topological skeleton was carried out by applying a Eulerian method based on the WSS vector field divergence. To provide objective thresholds for WSS topological skeleton quantitative analysis, a computational hemodynamic dataset of 46 ostensibly healthy carotid bifurcation models was considered. CEA interventions did not completely restore physiological WSS topological skeleton features. Significant associations emerged between IMT at 60 months follow-up and the exposure to (1) high temporal variation of WSS contraction/expansion (R2 = 0.51, p < 0.05), and (2) high fixed point residence times, weighted by WSS contraction/expansion strength (R2 = 0.53, p < 0.05). These WSS topological skeleton features were statistically independent from the exposure to low WSS, a previously reported predictor of long-term restenosis, therefore representing different hemodynamic stimuli and potentially impacting differently the vascular response. This study confirms the direct association between WSS topological skeleton and markers of vascular disease, contributing to elucidate the mechanistic link between flow disturbances and clinical observations of vascular lesions.

Deciphering ascending thoracic aortic aneurysm hemodynamics in relation to biomechanical properties.

The degeneration of the arterial wall at the basis of the ascending thoracic aortic aneurysm (ATAA) is a complex multifactorial process, which may lead to clinical complications and, ultimately, death. Individual genetic, biological or hemodynamic factors are inadequate to explain the heterogeneity of ATAA development/progression mechanisms, thus stimulating the analysis of their complex interplay. Here the disruption of the hemodynamic environment in the ATAA is investigated integrating patient-specific computational hemodynamics, CT-based in vivo estimation of local aortic stiffness and advanced fluid mechanics methods of analysis. The final aims are (1) deciphering the ATAA spatiotemporal hemodynamic complexity and its link to near-wall topological features, and (2) identifying the existing links between arterial wall degeneration and hemodynamic insult. Technically, two methodologies are applied to computational hemodynamics data, the wall shear stress (WSS) topological skeleton analysis, and the Complex Networks theory. The same analysis was extended to the healthy aorta. As main findings of the study, we report that: (1) different spatiotemporal heterogeneity characterizes the ATAA and healthy hemodynamics, that markedly reflect on their WSS topological skeleton features; (2) a link (stronger than canonical WSS-based descriptors) emerges between the variation of contraction/expansion action exerted by WSS on the endothelium along the cardiac cycle, and ATAA wall stiffness. The findings of the study suggest the use of advanced methods for a deeper understanding of the hemodynamics disruption in ATAA, and candidate WSS topological skeleton features as promising indicators of local wall degeneration.

Exploring wall shear stress spatiotemporal heterogeneity in coronary arteries combining correlation-based analysis and complex networks with computational hemodynamics.

Atherosclerosis at the early stage in coronary arteries has been associated with low cycle-average wall shear stress magnitude. However, parallel to the identification of an established active role for low wall shear stress in the onset/progression of the atherosclerotic disease, a weak association between lesions localization and low/oscillatory wall shear stress has been observed. In the attempt to fully identify the wall shear stress phenotype triggering early atherosclerosis in coronary arteries, this exploratory study aims at enriching the characterization of wall shear stress emerging features combining correlation-based analysis and complex networks theory with computational hemodynamics. The final goal is the characterization of the spatiotemporal and topological heterogeneity of wall shear stress waveforms along the cardiac cycle. In detail, here time-histories of wall shear stress magnitude and wall shear stress projection along the main flow direction and orthogonal to it (a measure of wall shear stress multidirectionality) are analyzed in a representative dataset of 10 left anterior descending pig coronary artery computational hemodynamics models. Among the main findings, we report that the proposed analysis quantitatively demonstrates that the model-specific inlet flow-rate shapes wall shear stress time-histories. Moreover, it emerges that a combined effect of low wall shear stress magnitude and of the shape of the wall shear stress-based descriptors time-histories could trigger atherosclerosis at its earliest stage. The findings of this work suggest for new experiments to provide a clearer determination of the wall shear stress phenotype which is at the basis of the so-called arterial hemodynamic risk hypothesis in coronary arteries.

Hemodialysis arterio-venous graft design reducing the hemodynamic risk of vascular access dysfunction.

Although arterio-venous grafts (AVGs) represent the second choice as permanent vascular access for hemodialysis, this solution is still affected by a relevant failure rate due to graft thrombosis, and development of neointimal hyperplasia (IH) at the distal vein. As a key role in these processes has been attributed to the abnormal hemodynamics establishing in the distal vein, the optimization of AVGs design aimed at minimizing flow disturbances would reduce AVG hemodynamic-related risks. In this study we used computational fluid dynamics to investigate the impact of alternative AVG designs on the reduction of IH and thrombosis risk at the distal venous anastomosis. The performance of the newly designed AVGs was compared to that of commercially available devices. In detail, a total of eight AVG models in closed-loop configuration were constructed: two models resemble the commercially available straight conventional and helical-shaped AVGs; six models are characterized by the insertion of a flow divider (FD), straight or helical shaped, differently positioned inside the graft. Unfavorable hemodynamic conditions were analyzed by assessing the exposure to disturbed shear at the distal vein. Bulk flow was investigated in terms of helical blood flow features, potential thrombosis risk, and pressure drop over the graft. Findings from this study clearly show that using a helically-shaped FD located at the venous side of the graft could induce beneficial helical flow patterns that, minimizing flow disturbances, reduce the IH-related risk of failure at the distal vein, with a clinically irrelevant increase in thrombosis risk and pressure drop over the graft.

A Eulerian method to analyze wall shear stress fixed points and manifolds in cardiovascular flows.

Based upon dynamical systems theory, a fixed point of a vector field such as the wall shear stress (WSS) at the luminal surface of a vessel is a point where the vector field vanishes. Unstable/stable manifolds identify contraction/expansion regions linking fixed points. The significance of such WSS topological features lies in their strong link with "disturbed" flow features like flow stagnation, separation and reversal, deemed responsible for vascular dysfunction initiation and progression. Here, we present a Eulerian method to analyze WSS topological skeleton through the identification and classification of WSS fixed points and manifolds in complex vascular geometries. The method rests on the volume contraction theory and analyzes the WSS topological skeleton through the WSS vector field divergence and Poincar[Formula: see text] index. The method is here applied to computational hemodynamics models of carotid bifurcation and intracranial aneurysm. An in-depth analysis of the time dependence of the WSS topological skeleton along the cardiac cycle is provided, enriching the information obtained from cycle-average WSS. Among the main findings, it emerges that on the carotid bifurcation, instantaneous WSS fixed points co-localize with cycle-average WSS fixed points for a fraction of the cardiac cycle ranging from 0 to [Formula: see text]; a persistent instantaneous WSS fixed point confined on the aneurysm dome does not co-localize with the cycle-average low-WSS region. In conclusion, the here presented approach shows the potential to speed up studies on the physiological significance of WSS topological skeleton in cardiovascular flows, ultimately increasing the chance of finding mechanistic explanations to clinical observations.

The impact of helical flow on coronary atherosclerotic plaque development.

BACKGROUND AND AIMS: Atherosclerosis has been associated with near-wall hemodynamics and wall shear stress (WSS). However, the role of coronary intravascular hemodynamics, in particular of the helical flow (HF) patterns that physiologically develop in those arteries, is rarely considered. The purpose of this study was to assess how HF affects coronary plaque initiation and progression, definitively demonstrating its atheroprotective nature. METHODS: The three main coronary arteries of five adult hypercholesterolemic mini-pigs on a high fat diet were imaged by computed coronary tomography angiography (CCTA) and intravascular ultrasound (IVUS) at 3 (T1, baseline) and 9.4 ± 1.9 (T2) months follow-up. The baseline geometries of imaged coronary arteries (n = 15) were reconstructed, and combined with pig-specific boundary conditions (based on in vivo Doppler blood flow measurements) to perform computational fluid dynamic simulations. Local wall thickness (WT) was measured on IVUS images at T1 and T2, and its temporal changes were assessed. Descriptors of HF and WSS nature were computed for each model, and statistically compared to WT data. RESULTS: HF intensity was strongly positively associated with WSS magnitude (p < 0.001). Overall, coronary segments exposed to high baseline levels of HF intensity exhibited a significantly lower WT growth (p < 0.05), compared to regions with either mid or low HF intensity. CONCLUSIONS: This study confirms the physiological significance of HF in coronary arteries, revealing its protective role against atherosclerotic WT growth and its potential in predicting regions undergoing WT development. These findings support future in vivo measurement of coronary HF as atherosclerotic risk marker, overcoming current limitations of in vivo WSS assessment.

Does the inflow velocity profile influence physiologically relevant flow patterns in computational hemodynamic models of left anterior descending coronary artery?

Patient-specific computational fluid dynamics is a powerful tool for investigating the hemodynamic risk in coronary arteries. Proper setting of flow boundary conditions in computational hemodynamic models of coronary arteries is one of the sources of uncertainty weakening the findings of in silico experiments, in consequence of the challenging task of obtaining in vivo 3D flow measurements within the clinical framework. Accordingly, in this study we evaluated the influence of assumptions on inflow velocity profile shape on coronary artery hemodynamics. To do that, (1) ten left anterior descending coronary artery (LAD) geometries were reconstructed from clinical angiography, and (2) eleven velocity profiles with realistic 3D features such as eccentricity and differently shaped (single- and double-vortex) secondary flows were generated analytically and imposed as inflow boundary conditions. Wall shear stress and helicity-based descriptors obtained prescribing the commonly used parabolic velocity profile were compared with those obtained with the other velocity profiles. Our findings indicated that the imposition of idealized velocity profiles as inflow boundary condition is acceptable as long the results of the proximal vessel segment are not considered, in LAD coronary arteries. As a pragmatic rule of thumb, a conservative estimation of the length of influence of the shape of the inflow velocity profile on LAD local hemodynamics can be given by the theoretical entrance length for cylindrical conduits in laminar flow conditions.