On the importance of blood rheology for bulk flow in hemodynamic models of the carotid bifurcation

Modeling Hemodynamic Effects of the Closure Procedure of Carotid Artery Endarterectomy

Currently, it is unclear why some patients experience restenosis after carotid endarterectomy (CEA) and whether the closing procedure is linked to greater rates of restenosis. Here, the morphology and hemodynamics are compared for the carotid bulb of two patients post-CEA. One carotid bulb was closed with a patch which later suffered a restenosis, while the other patient's bulb was treated using primary closure and did not. Contrast-enhanced magnetic resonance angiography (CE-MRA) was segmented to provide the domain for computational fluid dynamics (CFD). Flowrate waveforms measured with phase-contrast MR were provided for the common carotid artery (CCA) and internal carotid artery (ICA), while only the mean flow rate was provided for the external carotid artery (ECA), requiring the ECA waveform to be calculated. A Womersley profile was applied to the CCA inlet and ECA outlet, with a traction-free boundary condition applied to the ICA outlet. The patch patient who restenosed exhibited a nonphysiological hemodynamic environment that differed from the flow environment observed in the healthy, contralateral bulb. In contrast, the hemodynamics of the primary closure patient who underwent a successful CEA showed more favorable levels and trends of WSS as well as healthy mixing from vortices that were both present in the healthy, contralateral bulb. Changing model parameters such as flow rate, wall compliance, and flow waveforms did not alter these conclusions. Therefore, the geometry of the carotid bulb, as opposed to flow characteristics, seems responsible for the observed differences between these two cases in hemodynamic environments and subsequent outcomes.

Elucidating the hemodynamic impact of residual stenosis post-carotid artery stenting: A numerical study

BACKGROUND

Residual stenosis (RS) and hemodynamics demonstrate a significant correlation with postoperative in-stent restenosis/thrombosis following carotid artery stenting (CAS).

PURPOSE

This study endeavors to elucidate the potential associations between RS and adverse postoperative hemodynamic factors.

METHODS

This study utilized 46 patient-specific carotid artery models post-stenting, which were categorized into two groups based on the presence of RS: the normal group (N, n = 23) and the RS group (RS, n = 23). A comparative analysis was conducted to evaluate the discrepancies in geometry and adverse hemodynamic parameters, alongside investigating the potential correlation between hemodynamic and geometric parameters.

RESULTS

The results reveal that a higher reflux flow volume is discernible in the RS group during low-velocity phases of the cardiac cycle, concomitant with an augmented extent of areas exposed to oscillatory shear stress and extended particle residence time. Moreover, the adverse hemodynamic parameters exhibit a positive correlation with the degree of stent expansion, stent length in the common carotid artery (CCA), and the distal slope of the RS.

CONCLUSION

The distal slope and tortuosity of RS significantly influence the development of adverse hemodynamic conditions post-stenting, exacerbating the hemodynamic environment near the stenosis. Interestingly, while an extended stent length in the internal carotid artery (ICA) region improves hemodynamics by reducing flow disturbance, a longer stent in the CCA significantly worsens these conditions. Hence, it is prudent to analyze the characteristics of the local lesion regions to optimize the strategy for stent implantation.

Mathematical modeling in assessing the risk of restenosis after carotid endarterectomy

Carotid endarterectomy is the main way to combat atherosclerosis of the carotid arteries, which disrupts cerebral circulation. The generally accepted marker of atherogenesis risk are hemodynamic indices associated with near-wall shear stress. The purpose of the work is to conduct a comparative analysis of hemodynamic indices in various carotid bifurcation models. The influence of a virtual change in the geometric shape of the model in order to optimize hemodynamic indices is also being studied. On the basis of computed angiography data, carotid bifurcation models are constructed, in which critical zones of hemodynamic indices are built using computational fluid dynamics. A comparative analysis of the critical zones for different classes of models is carried out. Comparison of averaged indices for critical zones between 'normal' and post-operative groups gave more than 5-x worse results for the latter. The same results for the near-bifurcation parts of the zones give a 25% better result for postoperative models. Virtual 'removal' of insignificant plaques leads to a deterioration of the indices of up to 40% in the places of the plaque's former location. The described method makes it possible to build the indices critical zones and compare them for various types of models. A technique for virtual changing the shape of a vessel (virtual surgery) is proposed. The novelty of the approach lies in the use for comparative analysis both real vessel models and hypothetical 'improved' virtual ones, as well in the proposed division of post-operative model's critical zones into subzones of different genesis.

Prediction of Hemodynamic-Related Hemolysis in Carotid Stenosis and Aiding in Treatment Planning and Risk Stratification Using Computational Fluid Dynamics

Atherosclerosis affects human health in many ways, leading to disability or premature death due to ischemic heart disease, stroke, or limb ischemia. Poststenotic blood flow disruption may also play an essential role in artery wall impairment linked with hemolysis related to shear stress. The maximum shear stress in the atherosclerotic plaque area is the main parameter determining hemolysis risk. In our work, a 3D internal carotid artery model was built from CT scans performed on patients qualified for percutaneous angioplasty due to its symptomatic stenosis. The obtained stenosis geometries were used to conduct a series of computer simulations to identify critical parameters corresponding to the increase in shear stress in the arteries. Stenosis shape parameters responsible for the increase in shear stress were determined. The effect of changes in the carotid artery size, length, and degree of narrowing on the change in maximum shear stress was demonstrated. Then, a correlation for the quick initial diagnosis of atherosclerotic stenoses regarding the risk of hemolysis was developed. The developed relationship for rapid hemolysis risk assessment uses information from typical non-invasive tests for treated patients. Practical guidelines have been developed regarding which stenosis shape parameters pose a risk of hemolysis, which may be adapted in medical practice.

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?".

Hemolysis of red blood cells in blood vessels modeled via computational fluid dynamics

The research aims to verify the universal relationship between vessel shape and the risk of hemolysis using a rheological model of blood reflecting the physiological processes related to blood for any blood vessel. Blood is a multi-component fluid, the rheology of which depends on many factors, such as the concentration of red blood cells and local shear stress, which significantly affect the process of hemolysis. Blood rheology models used so far cannot be used for all flows and geometries. Therefore, a new rheology model has been introduced suitable for modeling hemolytic flows observed in arteries with atherosclerotic lesions in the in vivo environment. The previously presented model also has advantages in modeling local viscosity in stenosis. Geometries of the blood vessels from computed tomography scans and simplified models of the actual arteries observed during medical procedures were used in the calculations. Population Balance Based Rheology model predicts the concentration of single, deagglomerated red blood cells and the concentration and size of red blood cell agglomerates, which affect blood rheology and hemolysis. Based on the simulations carried out, a correlation was found between the shape of the vessel cavity and the risk of hemolysis. Presented results can be used in the future to create a correlation between the shape of the atherosclerotic lesions and the risk of hemolysis in the blood to make an initial risk assessment for a given patient.

Application of Patient-Specific Computational Fluid Dynamics in Anomalous Aortic Origin of Coronary Artery: A Systematic Review

Anomalous aortic origin of a coronary artery (AAOCA) is a rare congenital heart condition with fixed and dynamic stenotic elements, potentially causing ischemia. Invasive coronary angiography under stress is the established method for assessing hemodynamics in AAOCA, yet it is costly, technically intricate, and uncomfortable. Computational fluid dynamics (CFD) simulations offer a noninvasive alternative for patient-specific hemodynamic analysis in AAOCA. This systematic review examines the role of CFD simulations in AAOCA, encompassing patient-specific modeling, noninvasive imaging-based boundary conditions, and flow characteristics. Screening articles using AAOCA and CFD-related terms prior to February 2023 yielded 19 publications, covering 370 patients. Over the past four years, 12 (63%) publications (259 patients) employed dedicated CFD models, whereas 7 (37%) publications (111 patients) used general-purpose CFD models. Dedicated CFD models were validated for fixed stenosis but lacked dynamic component representation. General-purpose CFD models exhibited variability and limitations, with fluid-solid interaction models showing promise. Interest in CFD modeling of AAOCA has surged recently, mainly utilizing dedicated models. However, these models inadequately replicate hemodynamics, necessitating novel CFD approaches to accurately simulate pathophysiological changes in AAOCA under stress conditions.

Vortical Structures Promote Atheroprotective Wall Shear Stress Distributions in a Carotid Artery Bifurcation Model

Carotid artery diseases, such as atherosclerosis, are a major cause of death in the United States. Wall shear stresses are known to prompt plaque formation, but there is limited understanding of the complex flow structures underlying these stresses and how they differ in a pre-disposed high-risk patient cohort. A 'healthy' and a novel 'pre-disposed' carotid artery bifurcation model was determined based on patient-averaged clinical data, where the 'pre-disposed' model represents a pathological anatomy. Computational fluid dynamic simulations were performed using a physiological flow based on healthy human subjects. A main hairpin vortical structure in the internal carotid artery sinus was observed, which locally increased instantaneous wall shear stress. In the pre-disposed geometry, this vortical structure starts at an earlier instance in the cardiac flow cycle and persists over a much shorter period, where the second half of the cardiac cycle is dominated by perturbed secondary flow structures and vortices. This coincides with weaker favorable axial pressure gradient peaks over the sinus for the 'pre-disposed' geometry. The findings reveal a strong correlation between vortical structures and wall shear stress and imply that an intact internal carotid artery sinus hairpin vortical structure has a physiologically beneficial role by increasing local wall shear stresses. The deterioration of this beneficial vortical structure is expected to play a significant role in atherosclerotic plaque formation.

Fluid flow to mimic organ function in 3D in vitro models

Many different strategies can be found in the literature to model organ physiology, tissue functionality, and disease in vitro; however, most of these models lack the physiological fluid dynamics present in vivo. Here, we highlight the importance of fluid flow for tissue homeostasis, specifically in vessels, other lumen structures, and interstitium, to point out the need of perfusion in current 3D in vitro models. Importantly, the advantages and limitations of the different current experimental fluid-flow setups are discussed. Finally, we shed light on current challenges and future focus of fluid flow models applied to the newest bioengineering state-of-the-art platforms, such as organoids and organ-on-a-chip, as the most sophisticated and physiological preclinical platforms.

Computational Fluid Dynamics for the Prediction of Endograft Thrombosis in the Superficial Femoral Artery

Purpose:

Contemporary diagnostic modalities, including contrast-enhanced computed tomography (CTA) and duplex ultrasound, have been insufficiently able to predict endograft thrombosis. This study introduces an implementation of image-based computational fluid dynamics (CFD), by exemplification with 4 patients treated with an endograft for occlusive disease of the superficial femoral artery (SFA). The potential of personalized CFD for predicting endograft thrombosis is investigated.

Materials and Methods:

Four patients treated with endografts for an occluded SFA were retrospectively included. CFD simulations, based on CTA and duplex ultrasound, were compared for patients with and without endograft thrombosis to investigate potential flow-related causes of endograft thrombosis. Time-averaged wall shear stress (TAWSS) was computed, which highlights areas of prolonged residence times of coagulation factors in the graft.

Results:

CFD simulations demonstrated normal TAWSS (>0.4 Pa) in the SFA for cases 1 and 2, but low levels of TAWSS (<0.4 Pa) in cases 3 and 4, respectively. Primary patency was achieved in cases 1 and 2 for over 2 year follow-up. Cases 3 and 4 were complicated by recurrent endograft thrombosis.

Conclusion:

The presence of a low TAWSS was associated with recurrent endograft thrombosis in subjects with otherwise normal anatomic and ultrasound assessment and a good distal run-off.

Fluid-structure interaction study of hemodynamics and its biomechanical influence on carotid artery atherosclerotic plaque deposits

Atherosclerotic plaque deposits are common causes of blood flow disruption in the carotid artery bifurcation and the associated fluid mechanics has been extensively studied using Computational Fluid Dynamics (CFD) and Fluid Structure Interaction (FSI). However, the elastic responses of plaques to hemodynamics in the carotid artery bifurcation has not been deeply studied using either of the above-mentioned numerical techniques. In this study, a two-way FSI study was coupled with CFD technique, using Arbitrary-Lagrangian-Eulerian method, to study the biomechanics of blood flow on nonlinear and hyperelastic calcified plaque deposits in a realistic geometry of the carotid sinus. FSI parameters such as total mesh displacement and von Misses stress on the plaque, as well as flow velocity and blood pressure around the plaques, were analyzed and compared to variables such as velocity streamline, pressure and wall shear stress obtained from CFD simulation in a healthy model. The blood flow simulations reveal complete reversed blood flow behavior in the internal carotid artery, ICAs and external carotid artery, ECAs for both cases. In particular, this study suggests that plaques, irrespective of the masses, possess a high yielding response to hemodynamic forces at the attaching edges, while the surfaces are vulnerable to rupture.

Shear stress metrics associated with pro-atherogenic high-risk anatomical features in a carotid artery bifurcation model

BACKGROUND

Diseases associated with atherosclerotic plaques in the carotid artery are a major cause of deaths in the United States. Blood-flow-induced shear-stresses are known to trigger plaque formation. Prior literature suggests that the internal carotid artery sinus is prone to atherosclerosis, but there is limited understanding of why only certain patients are predisposed towards plaque formation.

METHODS

We computationally investigate the effect of vessel geometry on wall-shear-stress distribution by comparing flowfields and wall-shear-stress-metrics between a low-risk and a novel predisposed high-risk carotid artery bifurcation anatomy. Both models were developed based on clinical risk estimations and patient-averaged anatomical features. The high-risk geometry has a larger internal carotid artery branching angle and a lower internal-to-carotid-artery-diameter-ratio. A patient-averaged physiological carotid artery inflow waveform is used.

FINDINGS

The high-risk geometry experiences stronger flow separation in the sinus. Furthermore, it experiences a more equal flow split at the bifurcation, thereby reducing internal carotid artery flowrate and increasing atherosclerosis-prone low-velocity areas. Lowest time-averaged-wall-shear-stresses are present at the sinus outer wall, where plaques are often found, for both geometries. The high-risk geometry has significantly high, unfavorable oscillatory-shear-index values not found in the low-risk geometry. High oscillatory-shear-index areas are located at the vessels outside walls distal to the bifurcation and on the sinus wall.

INTERPRETATION

These results highlight the effectiveness of oscillatory-shear-index, to augment classical time-averaged-wall-shear-stress, in evaluating pro-atherogenic geometry features. Furthermore, the flow split at the bifurcation is a promising clinical indicator for atherosclerosis risk as it can be directly accessed using clinical imaging, whereas shear-stress-metrics cannot.

MRI-based mechanical analysis of carotid atherosclerotic plaque using a material-property-mapping approach: A material-property-mapping method for plaque stress analysis

BACKGROUND AND OBJECTIVE

Atherosclerosis is a major underlying cause of cardiovascular conditions. In order to understand the biomechanics involved in the generation and rupture of atherosclerotic plaques, numerical analysis methods have been widely used. However, several factors limit the practical use of this information in a clinical setting. One of the key challenges in finite element analysis (FEA) is the reconstruction of the structure and the generation of a mesh. The complexity of the shapes associated with carotid plaques, including multiple components, makes the generation of meshes for biomechanical computation a difficult and in some cases, an impossible task. To address these challenges, in this study, we propose a novel material-property-mapping method for carotid atherosclerotic plaque stress analysis that aims to simplify the process.

METHODS

The different carotid plaque components were identified and segmented using magnetic resonance imaging (MRI). For the mapping method, this information was used in conjunction with an in-house code, which provided the coordinates for each pixel/voxel and tissue type within a predetermined region of interest. These coordinates were utilized to assign specific material properties to each element in the volume mesh which provides a region of transition. The proposed method was subsequently compared to the traditional method, which involves creating a composed mesh for the arterial wall and plaque components, based on its location and size.

RESULTS

The comparison between the proposed material-property-mapping method and the traditional method was performed in 2D, 3D structural-only, and fluid-structure interaction (FSI) simulations in terms of stress, wall shear stress (WSS), time-averaged WSS (TAWSS), and oscillatory shear index (OSI). The stress contours from both methods were found to be similar, although the proposed method tended to produce lower local maximum stress values. The WSS contours were also in agreement between the two methods. The velocity contours generated by the proposed method were verified against phase-contrast magnetic resonance imaging (MRI) measurements, for a higher level of confidence.

CONCLUSION

This study shows that a material-property-mapping method can effectively be used for analyzing the biomechanics of carotid plaques in a patient-specific manner. This approach has the potential to streamline the process of creating volume meshes for complex biological structures, such as carotid plaques, and to provide a more efficient and less labor-intensive method.

Classification of Blood Rheological Models through an Idealized Symmetrical Bifurcation

The assumed rheological behavior of blood influences the hemodynamic characteristics of numerical blood flow simulations. Until now, alternative rheological specifications have been utilized, with uncertain implications for the results obtained. This work aims to group sixteen blood rheological models in homogeneous clusters, by exploiting data generated from numerical simulations on an idealized symmetrical arterial bifurcation. Blood flow is assumed to be pulsatile and is simulated using a commercial finite volume solver. An appropriate mesh convergence study is performed, and all results are collected at three different time instants throughout the cardiac cycle: at peak systole, early diastole, and late diastole. Six hemodynamic variables are computed: the time average wall shear stress, oscillatory shear index, relative residence time, global and local non-Newtonian importance factor, and non-Newtonian effect factor. The resulting data are analyzed using hierarchical agglomerative clustering algorithms, which constitute typical unsupervised classification methods. Interestingly, the rheological models can be partitioned into three homogeneous groups, whereas three specifications appear as outliers which do not belong in any partition. Our findings suggest that models which are defined in a similar manner from a mathematical perspective may behave substantially differently in terms of the data they produce. On the other hand, models characterized by different mathematical formulations may belong to the same statistical group (cluster) and can thus be considered interchangeably.

Computer Modeling of Carotid Endarterectomy With the Different Shape Patches and Prediction of the Atherosclerotic Plaque Formation Zones

The article describes a method for constructing geometric models of the carotid bifurcation and computer simulation of endarterectomy surgery with the patches of various configurations. The purpose of this work is to identify the areas of the greatest risk of restenosis in the constructed models and to conduct a comparative analysis of risk factors when using the patches of different widths and shapes. The method is demonstrated on a reconstructed model of a healthy vessel. Its building is based on a preoperative computed tomography study of a particular patient's affected vessel. The flow in the vessel is simulated by computational fluid dynamics using data from the patient's ultrasound Doppler velocimetry. Risk factors are assessed through the hemodynamic indices on the vessel wall associated with Wall Shear Stress. The distribution of risk zones in the healthy vessel, presumably leading to its observed lesion (plaque), is analyzed. Comparative evaluation of 10 various patches implantation results is carried out and the optimal variant is determined. The proposed method can be used to predict the hemodynamic results of surgery using patches of various sizes and shapes.

Modelling coronary flows: impact of differently measured inflow boundary conditions on vessel-specific computational hemodynamic profiles

BACKGROUND AND OBJECTIVES

The translation of hemodynamic quantities based on wall shear stress (WSS) or intravascular helical flow into clinical biomarkers of coronary atherosclerotic disease is still hampered by the assumptions/idealizations required by the computational fluid dynamics (CFD) simulations of the coronary hemodynamics. In the resulting budget of uncertainty, inflow boundary conditions (BCs) play a primary role. Accordingly, in this study we investigated the impact of the approach adopted for in vivo coronary artery blood flow rate assessment on personalized CFD simulations where blood flow rate is used as inflow BC.

METHODS

CFD simulations were carried out on coronary angiograms by applying personalized inflow BCs derived from four different techniques assessing in vivo surrogates of flow rate: continuous thermodilution, intravascular Doppler, frame count-based 3D contrast velocity, and diameter-based scaling law. The impact of inflow BCs on coronary hemodynamics was evaluated in terms of WSS- and helicity-based quantities.

RESULTS

As main findings, we report that: (i) coronary flow rate values may differ based on the applied flow derivation technique, as continuous thermodilution provided higher flow rate values than intravascular Doppler and diameter-based scaling law (p = 0.0014 and p = 0.0023, respectively); (ii) such intrasubject differences in flow rate values lead to different surface-averaged values of WSS magnitude and helical blood flow intensity (p<0.0020); (iii) luminal surface areas exposed to low WSS and helical flow topological features showed robustness to the flow rate values.

CONCLUSIONS

Although the absence of a clinically applicable gold standard approach prevents a general recommendation for one coronary blood flow rate derivation technique, our findings indicate that the inflow BC may impact computational hemodynamic results, suggesting that a standardization would be desirable to provide comparable results among personalized CFD simulations of the coronary hemodynamics.

Computational Fluid Dynamics for the Prediction of Endograft Thrombosis in the Superficial Femoral Artery

PURPOSE

Contemporary diagnostic modalities, including contrast-enhanced computed tomography (CTA) and duplex ultrasound, have been insufficiently able to predict endograft thrombosis. This study introduces an implementation of image-based computational fluid dynamics (CFD), by exemplification with 4 patients treated with an endograft for occlusive disease of the superficial femoral artery (SFA). The potential of personalized CFD for predicting endograft thrombosis is investigated.

MATERIALS AND METHODS

Four patients treated with endografts for an occluded SFA were retrospectively included. CFD simulations, based on CTA and duplex ultrasound, were compared for patients with and without endograft thrombosis to investigate potential flow-related causes of endograft thrombosis. Time-averaged wall shear stress (TAWSS) was computed, which highlights areas of prolonged residence times of coagulation factors in the graft.

RESULTS

CFD simulations demonstrated normal TAWSS (>0.4 Pa) in the SFA for cases 1 and 2, but low levels of TAWSS (<0.4 Pa) in cases 3 and 4, respectively. Primary patency was achieved in cases 1 and 2 for over 2 year follow-up. Cases 3 and 4 were complicated by recurrent endograft thrombosis.

CONCLUSION

The presence of a low TAWSS was associated with recurrent endograft thrombosis in subjects with otherwise normal anatomic and ultrasound assessment and a good distal run-off.

Simulation of Mechanical Heart Valve Dysfunction and the Non-Newtonian Blood Model Approach

The mechanical heart valve (MHV) is commonly used for the treatment of cardiovascular diseases. Nonphysiological hemodynamic in the MHV may cause hemolysis, platelet activation, and an increased risk of thromboembolism. Thromboembolism may cause severe complications and valve dysfunction. This paper thoroughly reviewed the simulation of physical quantities (velocity distribution, vortex formation, and shear stress) in healthy and dysfunctional MHV and reviewed the non-Newtonian blood flow characteristics in MHV. In the MHV numerical study, the dysfunction will affect the simulation results, increase the pressure gradient and shear stress, and change the blood flow patterns, increasing the risks of hemolysis and platelet activation. The blood flow passes downstream and has obvious recirculation and stagnation region with the increased dysfunction severity. Due to the complex structure of the MHV, the non-Newtonian shear-thinning viscosity blood characteristics become apparent in MHV simulations. The comparative study between Newtonian and non-Newtonian always shows the difference. The shear-thinning blood viscosity model is the basics to build the blood, also the blood exhibiting viscoelastic properties. More details are needed to establish a complete and more realistic simulation.

Non-Newtonian Endothelial Shear Stress Simulation: Does It Matter?

Patient-specific coronary endothelial shear stress (ESS) calculations using Newtonian and non-Newtonian rheological models were performed to assess whether the common assumption of Newtonian blood behavior offers similar results to a more realistic but computationally expensive non-Newtonian model. 16 coronary arteries (from 16 patients) were reconstructed from optical coherence tomographic (OCT) imaging. Pulsatile CFD simulations using Newtonian and the Quemada non-Newtonian model were performed. Endothelial shear stress (ESS) and other indices were compared. Exploratory indices including local blood viscosity (LBV) were calculated from non-Newtonian simulation data. Compared to the Newtonian results, the non-Newtonian model estimates significantly higher time-averaged ESS (1.69 (IQR 1.36)Pa versus 1.28 (1.16)Pa, p < 0.001) and ESS gradient (0.90 (1.20)Pa/mm versus 0.74 (1.03)Pa/mm, p < 0.001) throughout the cardiac cycle, under-estimating the low ESS (<1Pa) area (37.20 ± 13.57% versus 50.43 ± 14.16%, 95% CI 11.28-15.18, p < 0.001). Similar results were also found in the idealized artery simulations with non-Newtonian median ESS being higher than the Newtonian median ESS (healthy segments: 0.8238Pa versus 0.6618Pa, p < 0.001 proximal; 0.8179Pa versus 0.6610Pa, p < 0.001 distal; stenotic segments: 0.8196Pa versus 0.6611Pa, p < 0.001 proximal; 0.2546Pa versus 0.2245Pa, p < 0.001 distal) On average, the non-Newtonian model has a LBV of 1.45 times above the Newtonian model with an average peak LBV of 40-fold. Non-Newtonian blood model estimates higher quantitative ESS values than the Newtonian model. Incorporation of non-Newtonian blood behavior may improve the accuracy of ESS measurements. The non-Newtonian model also allows calculation of exploratory viscosity-based hemodynamic indices, such as local blood viscosity, which may offer additional information to detect underlying atherosclerosis.

Hemodynamic Impact of Stenting on Carotid Bifurcation: A Potential Role of the Stented Segment and External Carotid Artery

Carotid stenting near the bifurcation carina is associated with adverse events, especially in-stent restenosis, thrombosis, and side branch occlusion in clinical data. This study is aimed at determining the potential biomechanical mechanisms for these adverse events after carotid stenting. The patient-specific carotid models were constructed with different stenting scenarios to study the flow distribution and hemodynamic parameters, such as wall shear stress (WSS), flow velocity, relative residence time (RRT), and oscillating shear index (OSI) in the carotid bifurcation. The results suggested that the existing stents surely reduced blood flow to the external carotid artery (ECA) but enhanced local flow disturbance both in ECA and stented internal carotid artery (ICA), and the inner posterior wall of the stented ICA and the outer posterior wall of ECA might endure a relatively low level of WSS and remarkably elevated OSI and RRT. In addition, the implanted stent leads to more ECA adverse flow than ICA after stenting. While disturbed flow near the strut increased as stent length increased, blood flow and areas of local flow disturbance in ECA slightly decreased as stent length increased. In conclusion, the results revealed that ECA might be in relatively high levels of abnormal local hemodynamics after stenting, followed by stented ICA, leading to potential adverse events after intervention.

Comparison of Non-Newtonian Models of One-Dimensional Hemodynamics

The paper is devoted to the comparison of different one-dimensional models of blood flow. In such models, the non-Newtonian property of blood is considered. It is demonstrated that for the large arteries, the small parameter is observed in the models, and the perturbation method can be used for the analytical solution. In the paper, the simplified nonlinear problem for the semi-infinite vessel with constant properties is solved analytically, and the solutions for different models are compared. The effects of the flattening of the velocity profile and hematocrit value on the deviation from the Newtonian model are investigated.

Comparative evaluation of ballet-type and conventional stent graft configurations for endovascular aneurysm repair: A CFD analysis

PURPOSE:

Cross limb stent graft (SG) configuration technique for endovascular aneurysm repair (EVAR) is employed for splayed aortic bifurcations to avoid device kinking and smoothen cannulation. The present study investigates three types of stent graft (SG) configurations for endovascular aneurysm repair (EVAR) in abdominal aortic aneurysm. A computational fluid dynamic analysis was performed on the pulsatile non-Newtonian flow characteristics in three ideally modeled geometries of abdominal aortic (AA) SG configurations.

METHODS:

The three planar and crosslimb SG configurations were ideally modeled, namely, top-down nonballet-type, top-down ballet-type, and bottom-up nonballet-type configurations. In top-down SG configuration, most of the device is deployed in the main body in the vicinity of renal artery and the limbs are extended to the iliac artery. While in the bottom-up configuration, some of the SG device is deployed in the main body, the limbs are deployed in aortic bifurcation, and the extra stent graft of the main body is extended to the proximal aorta until the below of the renal artery. The effects of non-Newtonian pulsatile flow on the wall stresses and flow patterns of the three models were investigated and compared. Moreover, the average wall shear stress (AWSS), oscillatory shear stress index (OSI), absolute helicity, pressure distribution, graft displacement and flow visualization plots were analyzed.

RESULTS:

The top-down ballet-type showed less branch blockage effect than the top-down nonballet-type models. Furthermore, the top-down ballet-type configuration showed an increased tendency to sustain high WSS and higher helicity characteristics than that of the bottom-up and top-down non-ballet type configurations. However, displacement forces of the top-down ballet-type configuration were 40% and 9.6% higher than those of the bottom-up and top-down nonballet-type configurations, respectively.

CONCLUSIONS:

Some complications such as graft tearing, thrombus formation, limb disconnection during long term follow up periods might be relevant to hemodynamic characteristics according to the configurations of EVAR. Hence, the reported data required to be validated with the clinical results.

Evaluating the Impact of Calcification on Plaque Vulnerability from the Aspect of Mechanical Interaction Between Blood Flow and Artery Based on MRI

Acute cerebral ischemic events and thrombosis are associated with the rupture/erosion of carotid atherosclerotic plaques. The aim of the present study was to determine the impact of calcification deposition on the wall shear stress (WSS) and stresses within the plaques using 3D fluid-structure interaction (FSI) models. Six patients with calcified carotid atherosclerosis underwent multisequence magnetic resonance imaging (MRI) and were divided into three groups according to the calcification volume. To evaluate the role of the calcification deposition on the stresses, the calcification content was replaced by lipids and arterial tissue, respectively. By comparing the results from the simulation with calcification, and when changing it to lipids there was a significant increment in the stresses at the fibrous cap (p = 0.004). Instead, by changing it to arterial tissue, there was no significant difference (p = 0.07). The calcification shapes that presented the highest stresses were thin concave arc-shaped (AS1) and thin convex arc-shaped (AS3), with mean stress values of 107 ± 54.2 and 99.6 ± 23.4 kPa, respectively. It was also observed that, the calcification shape has more influence on the level of stress than its distance to the lumen. Higher WSS values were associated with the presence of calcification. Calcification shape plays an important role in producing high stresses in the plaque. This work further clarifies the impact of calcification on plaque vulnerability.

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 (R² = 0.51, p < 0.05), and (2) high fixed point residence times, weighted by WSS contraction/expansion strength (R² = 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.

Blood flow simulations in patient-specific geometries of the carotid artery: A systematic review

Computational Fluid Dynamics (CFD) and Fluid-Structure Interaction (FSI) are currently widely applied in the study of blood flow parameters and their alterations under pathological conditions, which are important indicators for diagnosis of atherosclerosis. In this manuscript, a systematic review of the published literature was conducted, according to the guidelines of the Preferred Reporting Items for Systematic Reviews and Meta-Analyses, on the simulation studies of blood flow in patient-specific geometries of the carotid artery bifurcation. Scopus, PubMed and ScienceDirect databases were used in the literature search, which was completed on the 3rd of August 2020. Forty-nine articles were included after the selection process and were organized in two distinct categories: the CFD studies (36/49 articles), which comprise only the fluid analysis and the FSI studies (13/49 articles), which includes both fluid and Fluid-Structure domain in the analysis. The data of the research works was structured in different categories (Geometry, Viscosity models, Type of Flow, Boundary Conditions, Flow Parameters, Type of Solver and Validation). The aim of this systematic review is to demonstrate the methodology in the modelling, simulation and analysis of carotid blood flow and also identify potential gaps and challenges in this research field.

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.

Three dimensional numerical analysis of hemodynamic of stenosed artery considering realistic outlet boundary conditions

BACKGROUND AND OBJECTIVE

Mortality rate increases globally among which one third is due to diseased blood vessels. Due to late diagnoses of the disease in vessels (severe stenoses), qualitative and rapid assessment becomes difficult. Earlier assessment of stenoses can lead to formulation of effective treatment protocol. It is often found that proliferation of secondary stenoses at downstream of a stenosed vessel depends on the degree of severity of primary stenoses. Numerical investigation of flow dynamics of blood in such condition helps in prediction of distributed field of secondarystenoses. This investigation also requires consideration of rigorous boundary conditions at inlet and outlet of defined flow domain.

METHODS

Patient-specific geometry of aortic arch with stenoses in descending aorta was considered for numerical estimation of biofluid dynamics. Boundary conditionsat inlet and outlet were extracted from time-resolved pulsed Doppler Ultrasound imaging at appropriate sections of the vessel. Womersley inlet flux was considered. Flow parameters like wall shear stress, oscillatory shear index, etc. were evaluated at upper and lower aortic arch of the vessel at different combinations of boundary conditions at inlet and four outlets respectively.

RESULTS

Effect of outlet boundary conditions were acknowledged for the progression of secondary stenoses. Severity of primary stenoses was found influencing the progression of secondary stenoses. It was found that the outlets Left Subclavian Artery and Left Common Carotid Artery greatly influence the flow dynamic structure within the stenosed aortic arch. Simultaneously, lower wall of aortic-arch had shown more affinity for secondary stenoses progression.

CONCLUSION

Aortic arch is a vital anatomical region of circulatory system which is vulnerable to progression of secondary stenoses in presence of primary stenoses in ascending or descending aorta. It also drives the author to speculate the influence of anurysm in descending aorta on this landmark of aortic arch.

THE EFFECT OF PATIENT-SPECIFIC NON-NEWTONIAN BLOOD VISCOSITY ON ARTERIAL HEMODYNAMICS PREDICTIONS

Blood flow simulations can identify arterial regions that are vulnerable to atherosclerotic or thrombotic evolution. To accurately define vulnerable arterial regions, hemodynamic parameters such as arterial geometry, blood flow velocity and blood viscosity (BV) must be measured individually. However, previous numerical studies have largely employed either a single representative value or simply used a nonspecific curve of non-Newtonian characteristics of BV. This study aimed to evaluate whether various BV models could produce similar arterial wall shear stress (WSS) results. We performed a blood flow simulation in carotid arteries obtained from time-of-flight magnetic resonance (TOF MR) angiography using the hemodynamic characteristics of subjects via carotid duplex ultrasonography. The BV models were categorized into the following five types: patient-specific non-Newtonian BV (model 1), representative non-Newtonian BV based on the Carreau model (BV model 2), Newtonian BV measured at a specific shear rate of 300[Formula: see text][Formula: see text] (BV model 3), Newtonian BV obtained from a hematocrit-based equation (BV model 4) and a representative Newtonian BV of 4[Formula: see text]cP (BV model 5). In total, 20 carotid arteries from 20 healthy volunteers (mean age [Formula: see text] SD of [Formula: see text] years; 50% women) were examined. Compared with the mean values of carotid WSS in BV model 1 as the reference model, the other four models showed significant differences in both genders for all carotid segments. To obtain reliable physiological WSS results, patient-specific non-Newtonian BV should be carefully employed.

Mechanotransmission of haemodynamic forces by the endothelial glycocalyx in a full-scale arterial model

The glycocalyx has been identified as a key mechano-sensor of the shear forces exerted by streaming blood onto the vascular endothelial lining. Although the biochemical reaction to the blood flow has been extensively studied, the mechanism of transmission of the haemodynamic shear forces to the endothelial transmembrane anchoring structures and, consequently, to the subcellular elements in the cytoskeleton, is still not fully understood. Here we apply a multiscale approach to elucidate how haemodynamic shear forces are transmitted to the transmembrane anchors of endothelial cells. Wall shear stress time histories, as obtained from image-based computational haemodynamics models of a carotid bifurcation, are used as a load and a continuum model is applied to obtain the mechanical response of the glycocalyx all along the cardiac cycle. The main findings of this in silico study are that: (1) the forces transmitted to the transmembrane anchors are in the range of 1-10 pN, which is in the order of magnitude reported for the different conformational states of transmembrane mechanotranductors; (2) locally, the forces transmitted to the anchors of the glycocalyx structure can be markedly different from the near-wall haemodynamic shear forces both in amplitude and frequency content. The findings of this in silico approach warrant future studies focusing on the actual forces transmitted to the transmembrane mechanotransductors, which might outperform haemodynamic descriptors of disturbed shear as localizing factors of vascular disease.

HEMODYNAMIC EFFECT OF OBSTRUCTION TO RENAL ARTERIES CAUSED BY STENT GRAFTS IN PATIENTS WITH ABDOMINAL AORTIC ANEURYSMS

To investigate the hemodynamic effects of partial obstruction to the renal orifice caused by inappropriate stent-graft location. Pre-operative and deployment models of a stent graft with various degrees of obstruction to the renal orifice are constructed based on medical images of abdominal aortic aneurysm. Hemodynamics, including flow pattern, time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), as well as relative residence time (RRT) are analyzed using numerical simulations. Flow rate distributions are assessed and verified by in vitro experiments. Results show that partial blockage to the renal branch orifice leads to flow recirculation and vortices with low wall shear stress around the renal ostia, whereas OSI and RRT on the renal arteries increase with the degree of obstruction. While the in vitro experiment indicates a decreasing flow rate to the bilateral renal arteries as renal artery ostia are obstructed. In conclusion, obstruction to the renal arteries induced by an inappropriate stent graft location causes stenosis in the renal artery in the long term. This study reveals a possible pathological mechanism of renal complications due to the implantation of a stent graft.

Numerical and Experimental Investigation of Novel Blended Bifurcated Stent Grafts with Taper to Improve Hemodynamic Performance

The typical helical flow within the human arterial system is widely used when designing cardiovascular devices, as this helical flow can be generated using the “crossed limbs” strategy of the bifurcated stent graft (BSG) and enhanced by the tapered structure of arteries. Here, we propose the use of a deflected blended bifurcated stent graft (BBSG) with various tapers, using conventional blended BSGs with the same degree of taper as a comparison. Hemodynamic performances, including helical strength and wall shear stress- (WSS-) based indicators, were assessed. Displacement forces that may induce stent-graft migration were assessed using numerical simulations and in vitro experiments. The results showed that as the taper increased, the displacement force, helicity strength, and time-averaged wall shear stress (TAWSS) within the iliac grafts increased, whereas the oscillating shear index (OSI) and relative residence time (RRT) gradually decreased for both types of BBSGs. With identical tapers, deflected BBSGs, compared to conventional BBSGs, exhibited a wider helical structure and lower RRT on the iliac graft and lower displacement force; however, there were no differences in hemodynamic indicators. In summary, the presence of tapering facilitated helical flow and produced better hemodynamic performance but posed a higher risk of graft migration. Conventional and deflected BBSGs with taper might be the two optimal configurations for endovascular aneurysm repair, given the helical flow. The deflected BBSG provides a better configuration, compared to the conventional BBSG, when considering the reduction of migration risk.

Accounting for residence-time in blood rheology models: do we really need non-Newtonian blood flow modelling in large arteries?

Patient-specific computational fluid dynamics (CFD) is a promising tool that provides highly resolved haemodynamics information. The choice of blood rheology is an assumption in CFD models that has been subject to extensive debate. Blood is known to exhibit shear-thinning behaviour, and non-Newtonian modelling has been recommended for aneurysmal flows. Current non-Newtonian models ignore rouleaux formation, which is the key player in blood's shear-thinning behaviour. Experimental data suggest that red blood cell aggregation and rouleaux formation require notable red blood cell residence-time (RT) in a low shear rate regime. This study proposes a novel hybrid Newtonian and non-Newtonian rheology model where the shear-thinning behaviour is activated in high RT regions based on experimental data. Image-based abdominal aortic and cerebral aneurysm models are considered and highly resolved CFD simulations are performed using a minimally dissipative solver. Lagrangian particle tracking is used to define a backward particle RT measure and detect stagnant regions with increased rouleaux formation likelihood. Our novel RT-based non-Newtonian model shows a significant reduction in shear-thinning effects and provides haemodynamic results qualitatively identical and quantitatively close to the Newtonian model. Our results have important implications in patient-specific CFD modelling and suggest that non-Newtonian models should be revisited in large artery flows.

Numerical and Experimental Investigation of the Hemodynamic Performance of Bifurcated Stent Grafts with Various Torsion Angles

The “crossed limbs” strategy for bifurcated stent grafts (BSGs) is widely employed when abdominal aortic aneurysm (AAA) patients have unfavorable neck or highly splayed iliac arteries. Helical flow is regarded as a typical flow pattern within the human arterial system and is believed to have the positive physiological effects of inhibiting thrombosis formation and atherosclerosis. The “crossed limbs” strategy may induce helical flow and improve the stent graft outcome. To verify the performance of this strategy by considering hemodynamics, we constructed a series of idealized BSGs with various torsion angles and evaluated the hemodynamic performance, including the helical strength, time-averaged wall shear stress (TAWSS), oscillatory shear index, relative resident time (RRT), and displacement force. Our numerical results indicate that an increased torsion angle enhances the helicity strength at the iliac outlets. However, with increasing torsion angle, the TAWSS in the iliac graft decreases and the RRT increases. In addition, our numerical simulations and in vitro experiments reveal that the displacement force increases gradually with increasing torsion angle. In summary, the “crossed limbs” strategy may have benefits for AAA treatment in terms of helical flow, but because of the unfavorable hemodynamic performance verified by analyzing the hemodynamic indicators, the risk of stent graft migration increases with increasing torsion angle. Therefore, the “crossed limbs” strategy should be carefully employed in surgical AAA treatment.

Hemodynamic performance within crossed stent grafts: computational and experimental study on the effect of cross position and angle

Background and aims

The crossed limbs stent graft technique is regularly employed to treat abdominal aortic aneurysm patients with unfavorable aneurysm necks or widely splayed common iliac arteries. This article numerically evaluates the hemodynamic performance of the crossed limbs strategy by analyzing numerical simulations and conducting experiments using two series of idealized bifurcated stent grafts with different cross angles and cross positions.

Results

Results demonstrated that the absolute helicity at outlets decreased with increased cross angles and increased with decreased cross positions. The time-averaged wall shear stress remained approximately unchanged, whereas the oscillating shear index and relative resident time decreased slightly when the cross angle increased and cross position decreased in iliac grafts. Additionally, both numerical and in vitro experimental results indicate the displacement force acting on the stent graft gradually increased as cross angles increased and cross positions decreased. Results further indicated that strip areas with a high oscillating shear index and high relative resident time, which may be vulnerable to thrombosis formation, exist along the outer surface of the iliac artery grafts.

Conclusion

Given this information, the optimal crossed limbs configuration may contain a small cross angle and low cross position; however, low cross positions may increase the risk of migration.

Numerical study on the hemodynamics of patient-specific carotid bifurcation using a new mesh approach

The definition of a suitable mesh to simulate blood flow in the human carotid bifurcation has been investigated. In this research, a novel mesh generation method is proposed: hexahedral cells at the center of the vessel and a fine grid of tetrahedral cells near the artery wall, in order to correctly simulate the large blood velocity gradients associated with specific locations. The selected numerical examples to show the pertinence of the novel generation method are supported by carotid ultrasound image data of a patient-specific case. Doppler systolic blood velocities measured during ultrasound examination are compared with simulated velocities using 4 different combinations of hexahedral and tetrahedral meshes and different fluid dynamic simulators. The Lin's test was applied to show the concordance of the results. Wall shear stress-based descriptors and localized normalized helicity descriptor emphasize the performance of the new method. Another feature is the reduced computation time required by the developed methodology. With the accurate combined mesh, different flow rate partitions, between the internal carotid artery and external carotid artery, were studied. The overall effect of the partitions is mainly in the blood flow patterns and in the hot-spot modulation of atherosclerosis-susceptible regions, rather than in their distribution along the bifurcation.

Comparison between a generalized Newtonian model and a network-type multiscale model for hemodynamic behavior in the aortic arch: Validation with 4D MRI data for a case study

Blood is a complex fluid in which the presence of the various constituents leads to significant changes in its rheological properties. Thus, an appropriate non-Newtonian model is advisable; and we choose a Modified version of the rheological model of Phan-Thien and Tanner (MPTT). The different parameters of this model, derived from the rheology of polymers, allow characterization of the non-Newtonian nature of blood, taking into account the behavior of red blood cells in plasma. Using the MPTT model that we implemented in the open access software OpenFOAM, numerical simulations have been performed on blood flow in the thoracic aorta for a healthy patient. We started from a patient-specific model which was constructed from medical images. Exiting flow boundary conditions have been developped, based on a 3-element Windkessel model to approximate physiological conditions. The parameters of the Windkessel model were calibrated with in vivo measurements of flow rate and pressure. The influence of the selected viscosity of red blood cells on the flow and wall shear stress (WSS) was investigated. Results obtained from this model were compared to those of the Newtonian model, and to those of a generalized Newtonian model, as well as to in vivo dynamic data from 4D MRI during a cardiac cycle. Upon evaluating the results, the MPTT model shows better agreement with the MRI data during the systolic and diastolic phases than the Newtonian or generalized Newtonian model, which confirms our interest in using a complex viscoelastic model.

WALL SHEAR STRESS AND OSCILLATORY SHEAR INDEX DISTRIBUTION IN CAROTID ARTERY WITH VARYING DEGREE OF STENOSIS: A HEMODYNAMIC STUDY

A significant proportion of cerebral stroke is a consequence of the arterial stenotic plaque rupture causing local thrombosis or distal embolization. The formation and subsequent rupture of the plaque depends on wall shear stress (WSS) and oscillatory shear index (OSI). The purpose of the present study was to understand the effect of hemodynamics on the spatial and temporal variations of WSS and OSI using realistic models with varying degree of carotid artery stenosis (DOS). Multiple CT volumes were obtained from subjects in the carotid bifurcation zone and the 3D models were generated. A finite volume-based computational fluid dynamics (CFD) method was utilized to understand the hemodynamics in pulsatile flow conditions. It was observed that high stenosis models occupied a large value of normalized WSS in the internal carotid artery (ICA) whereas they had smaller values of normalized WSS in the common carotid artery (CCA). For clinical use, the authors recommend using the spatial average value of oscillatory shear rather than the maximum value for an accurate knowledge about the severity of stenosis. The resultant vorticity changes the direction of spin after the bifurcation zone. Additionally, we propose the use of limiting streamlines as a novel and convenient method to identify the disturbed flow regions that are prone to atherogenesis.

Assessment of Influences of Stenoses in Right Carotid Artery on Left Carotid Artery Using Wall Stress Marker

Purpose. Atherosclerosis is a diseased condition of blood vessel. It causes partial blockage in lumen of vessel and affects hemodynamic of localized flowing blood. Complex geometries like region of bifurcation also affects hemodynamic to a larger extent. Complexity further increases in presence of stenoses at region of bifurcation. Such morphological change in vessel largely affects parent as well as corresponding sister and daughter vessels. In this paper, complexity in hemodynamic of blood in pair of carotid arteries (left and right carotid arteries) is evaluated in presence of stenoses at basilar segment of right artery in three-dimensional domain using reconstructed tomographic images of patient. Methods. Transient information of blood flow is obtained using four-dimensional phase-contrast MRI technique. Haematocrit component of blood at diseased condition is considered using Power Law and Quemada model. Numerical techniques are used to solve pressure-coupled governing equations of flowing blood. Results. Dysfunctions of endothelial cells near the wall are characterised by evaluating shear stress markers. Wall shear stress and its gradient based and harmonic based descriptors are calculated over complete geometry during one cardiac cycle. Conclusion. Internal branch of left carotid artery and external branch of right carotid artery are found prone to secondary stenoses in presence of primary stenoses at basilar segment of right carotid artery.

Computational fluid dynamic simulation of human carotid artery bifurcation based on anatomy and volumetric blood flow rate measured with magnetic resonance imaging

Blood flow patterns and local hemodynamic parameters have been widely associated with the onset and progression of atherosclerosis in the carotid artery. Assessment of these parameters can be performed noninvasively using cine phase-contrast (PC) magnetic resonance imaging (MRI). In addition, in the last two decades, computational fluid dynamics (CFD) simulation in three dimensional models derived from anatomic medical images has been employed to investigate the blood flow in the carotid artery. This study developed a workflow of a subject-specific CFD analysis using MRI to enhance estimating hemodynamics of the carotid artery. Time-of-flight (TOF) MRI scans were used to construct three-dimensional computational models. PC-MRI measurements were utilized to impose the boundary condition at the inlet and a 0-dimensional lumped parameter model was employed for the outflow boundary condition. The choice of different viscosity models of blood flow as a source of uncertainty was studied, by means of the axial velocity, wall shear stress, and oscillatory shear index. The sequence of workflow in CFD analysis was optimized for a healthy subject using PC-MRI. Then, a patient with carotid artery stenosis and its hemodynamic parameters were examined. The simulations indicated that the lumped parameter model used at the outlet gives physiologically reasonable values of hemodynamic parameters. Moreover, the dependence of hemodynamics parameters on the viscosity models was observed to vary for different geometries. Other factors, however, may be required for a more accurate CFD analysis, such as the segmentation and smoothness of the geometrical model, mechanical properties of the artery's wall, and the prescribed velocity profile at the inlet.

Computational simulation of carotid stenosis and flow dynamics based on patient ultrasound data - A new tool for risk assessment and surgical planning

PURPOSE

There is nowadays extensive experimental and computational investigation on the pathophysiology of atherosclerosis, searching correlations between its focal nature and local hemodynamic environment. The goal of this work is to present a methodology for patient-specific hemodynamics study of the carotid artery bifurcation based on the use of ultrasound (US) morphological and blood flow velocity patient data.

MATERIALS/METHODS

Subject-specific studies were performed for two patients, using a developed finite element code. Geometrical models were obtained from the acquisition of longitudinal and sequential cross-sectional ultrasound images and boundary conditions from Doppler velocity measurements at the common carotid artery.

RESULTS

There was a good agreement between ultrasound imaging data and computational simulated results. For a normal and a stenosed carotid bifurcation the velocity, wall shear stress (WSS) and WSS descriptors analysis illustrated the extremely complex hemodynamic behavior along the cardiac cycle. Different patterns were found, associated with morphology and hemodynamic patient-specific conditions. High values of time-averaged WSS (TAWSS) were found at stenosis site and for both patients TAWSS fields presented low values within areas of high oscillating shear index and relative residence time values, corresponding to recirculation zones.

CONCLUSION

Simulated hemodynamic parameters were able to capture the disturbed flow conditions in a normal and a stenosed carotid artery bifurcation, which play an important role in the development of local atherosclerotic plaques. Computational simulations based on clinic US might help improving diagnostic and treatment management of carotid atherosclerosis.

Numerical study of wall shear stress-based descriptors in the human left coronary artery

The present work is about the application of wall shear stress descriptors - time averaged wall shear stress (TAWSS), oscillating shear index (OSI) and relative residence time (RRT) - to the study of blood flow in the left coronary artery (LCA). These descriptors aid the prediction of disturbed flow conditions in the vessels and play a significant role in the detection of potential zones of atherosclerosis development. Hemodynamic descriptors data were obtained, numerically, through ANSYS® software, for the LCA of a patient-specific geometry and for a 3D idealized model. Comparing both cases, the results are coherent, in terms of location and magnitude. Low TAWSS, high OSI and high RRT values are observed in the bifurcation - potential zone of atherosclerosis appearance. The dissimilarities observed in the TAWSS values, considering blood as a Newtonian or non-Newtonian fluid, releases the importance of the correct blood rheologic caracterization. Moreover, for a higher Reynolds number, the TAWSS values decrease in the bifurcation and along the LAD branch, increasing the probability of plaques deposition. Furthermore, for a stenotic LCA model, very low TAWSS and high RRT values in front and behind the stenosis are observed, indicating the probable extension, in the flow direction, of the lesion.

Insights into the co-localization of magnitude-based versus direction-based indicators of disturbed shear at the carotid bifurcation

The observed co-localization of disturbed flow and lesion prevalence at predisposed districts such as the carotid bifurcation has led to the identification of the wall shear stress (WSS) as biomechanical localizing factor of vascular dysfunction. In particular, a proatherogenic role is attributed to low and oscillatory WSS. However, the endothelial cells (ECs) are exposed to a complex hemodynamic milieu that can be only partially described by low/oscillatory WSS. Recently, in the attempt to close this gap, descriptors of the complex multidirectional nature of WSS have been proposed, i.e., the axial component of WSS (aligned with the vessel׳s centerline, to quantify flow reversal), and the transverse WSS (transWSS, quantifying the WSS component orthogonal to the cycle averaged WSS direction). Here we explore the relationship between recently-proposed indicators quantifying WSS multidirectionality and "established" WSS-based hemodynamic descriptors of low/oscillatory WSS, in a representative sample (N=46) of subject-specific computational hemodynamics models of ostensibly normal carotid bifurcations. To do it, we quantitatively assess the co-localization of those descriptors at the luminal surface, aiming at providing connections among the peculiar hemodynamic features captured by the different descriptors. According to our findings: (1) regions of flow reversal are moderately co-localized with low WSS; (2) high WSS oscillations (quantified by the oscillatory shear index, OSI) at the carotid bulb are prevalently aligned with the main flow, where flow reversal is predominant; (3) regions where transWSS is high do not co-localize with the other descriptors. We suggest that the investigated WSS-based descriptors might represent different hemodynamic disturbances with different impact on ECs homeostasis, potentially being part of WSS phenotypes more effective in localizing the map of vascular atherosclerotic lesions.

Role of biomechanical forces in the natural history of coronary atherosclerosis

Atherosclerosis remains a major cause of morbidity and mortality worldwide, and a thorough understanding of the underlying pathophysiological mechanisms is crucial for the development of new therapeutic strategies. Although atherosclerosis is a systemic inflammatory disease, coronary atherosclerotic plaques are not uniformly distributed in the vascular tree. Experimental and clinical data highlight that biomechanical forces, including wall shear stress (WSS) and plaque structural stress (PSS), have an important role in the natural history of coronary atherosclerosis. Endothelial cell function is heavily influenced by changes in WSS, and longitudinal animal and human studies have shown that coronary regions with low WSS undergo increased plaque growth compared with high WSS regions. Local alterations in WSS might also promote transformation of stable to unstable plaque subtypes. Plaque rupture is determined by the balance between PSS and material strength, with plaque composition having a profound effect on PSS. Prospective clinical studies are required to ascertain whether integrating mechanical parameters with medical imaging can improve our ability to identify patients at highest risk of rapid disease progression or sudden cardiac events.