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

Modelling lower-limb peripheral arterial disease using clinically available datasets: impact of inflow boundary conditions on hemodynamic indices for restenosis prediction

BACKGROUND AND OBJECTIVES

The integration of hemodynamic markers as risk factors in restenosis prediction models for lower-limb peripheral arteries is hindered by fragmented clinical datasets. Computed tomography (CT) scans enable vessel geometry reconstruction and can be obtained at different times than the Doppler ultrasound (DUS) images, which provide information on blood flow velocity. Computational fluid dynamics (CFD) simulations allow the computation of near-wall hemodynamic indices, whose accuracy depends on the prescribed inlet boundary condition (BC), derived from the DUS images. This study aims to: (i) investigate the impact of different DUS-derived velocity waveforms on CFD results; (ii) test whether the same vessel areas, subjected to altered hemodynamics, can be detected independently of the applied inlet BC; (iii) suggest suitable DUS images to obtain reliable CFD results.

METHODS

CFD simulations were conducted on three patients treated with bypass surgery, using patient-specific DUS-derived inlet BCs recorded at either the same or different time points than the CT scan. The impact of the chosen inflow condition on bypass hemodynamics was assessed in terms of wall shear stress (WSS)-derived quantities. Patient-specific critical thresholds for the hemodynamic indices were applied to identify critical luminal areas and compare the results with a reference obtained with a DUS image acquired in close temporal proximity to the CT scan.

RESULTS

The main findings indicate that: (i) DUS-derived inlet velocity waveforms acquired at different time points than the CT scan led to statistically significantly different CFD results (p<0.001); (ii) the same luminal surface areas, exposed to low time-averaged WSS, could be identified independently of the applied inlet BCs; (iii) similar outcomes were observed for the other hemodynamic indices if the prescribed inlet velocity waveform had the same shape and comparable systolic acceleration time to the one recorded in close temporal proximity to the CT scan.

CONCLUSIONS

Despite a lack of standardised data collection for diseased lower-limb peripheral arteries, an accurate estimation of luminal areas subjected to altered near-wall hemodynamics is possible independently of the applied inlet BC. This holds if the applied inlet waveform shares some characteristics - derivable from the DUS report - as one matching the acquisition time of the CT scan.

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.

Influence of intracoronary hemodynamic forces on atherosclerotic plaque phenotypes

BACKGROUND AND AIMS

Coronary hemodynamics impact coronary plaque progression and destabilization. The aim of the present study was to establish the association between focal vs. diffuse intracoronary pressure gradients and wall shear stress (WSS) patterns with atherosclerotic plaque composition.

METHODS

Prospective, international, single-arm study of patients with chronic coronary syndromes and hemodynamic significant lesions (fractional flow reserve [FFR] ≤ 0.80). Motorized FFR pullback pressure gradient (PPG), optical coherence tomography (OCT), and time-average WSS (TAWSS) and topological shear variation index (TSVI) derived from three-dimensional angiography were obtained.

RESULTS

One hundred five vessels (median FFR 0.70 [Interquartile range (IQR) 0.56-0.77]) had combined PPG and WSS analyses. TSVI was correlated with PPG (r = 0.47, [95% Confidence Interval (95% CI) 0.30-0.65], p < 0.001). Vessels with a focal CAD (PPG above the median value of 0.67) had significantly higher TAWSS (14.8 [IQR 8.6-24.3] vs. 7.03 [4.8-11.7] Pa, p < 0.001) and TSVI (163.9 [117.6-249.2] vs. 76.8 [23.1-140.9] m-1, p < 0.001). In the 51 vessels with baseline OCT, TSVI was associated with plaque rupture (OR 1.01 [1.00-1.02], p = 0.024), PPG with the extension of lipids (OR 7.78 [6.19-9.77], p = 0.003), with the presence of thin-cap fibroatheroma (OR 2.85 [1.11-7.83], p = 0.024) and plaque rupture (OR 4.94 [1.82 to 13.47], p = 0.002).

CONCLUSIONS

Focal and diffuse coronary artery disease, defined using coronary physiology, are associated with differential WSS profiles. Pullback pressure gradients and WSS profiles are associated with atherosclerotic plaque phenotypes. Focal disease (as identified by high PPG) and high TSVI are associated with high-risk plaque features.

CLINICAL TRIAL REGISTRATION

https://clinicaltrials,gov/ct2/show/NCT03782688.

Blood Flow Energy Identifies Coronary Lesions Culprit of Future Myocardial Infarction

The present study establishes a link between blood flow energy transformations in coronary atherosclerotic lesions and clinical outcomes. The predictive capacity for future myocardial infarction (MI) was compared with that of established quantitative coronary angiography (QCA)-derived predictors. Angiography-based computational fluid dynamics (CFD) simulations were performed on 80 human coronary lesions culprit of MI within 5 years and 108 non-culprit lesions for future MI. Blood flow energy transformations were assessed in the converging flow segment of the lesion as ratios of kinetic and rotational energy values (KER and RER, respectively) at the QCA-identified minimum lumen area and proximal lesion sections. The anatomical and functional lesion severity were evaluated with QCA to derive percentage area stenosis (%AS), vessel fractional flow reserve (vFFR), and translesional vFFR (ΔvFFR). Wall shear stress profiles were investigated in terms of topological shear variation index (TSVI). KER and RER predicted MI at 5 years (AUC = 0.73, 95% CI 0.65-0.80, and AUC = 0.76, 95% CI 0.70-0.83, respectively; p < 0.0001 for both). The predictive capacity for future MI of KER and RER was significantly stronger than vFFR (p = 0.0391 and p = 0.0045, respectively). RER predictive capacity was significantly stronger than %AS and ΔvFFR (p = 0.0041 and p = 0.0059, respectively). The predictive capacity for future MI of KER and RER did not differ significantly from TSVI. Blood flow kinetic and rotational energy transformations were significant predictors for MI at 5 years (p < 0.0001). The findings of this study support the hypothesis of a biomechanical contribution to the process of plaque destabilization/rupture leading to MI.

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

Toward a Mesoscopic Modeling Approach of Magnetohydrodynamic Blood Flow in Pathological Vessels: A Comprehensive Review

The investigation of magnetohydrodynamic (MHD) blood flow within configurations that are pertinent to the human anatomy holds significant importance in the realm of scientific inquiry because of its practical implications within the medical field. This article presents an exhaustive appraisal of the diverse applications of magnetohydrodynamics and their computational modeling in biological contexts. These applications are classified into two categories: simple flow and pulsatile flow. An alternative approach of traditional CFD methods called Lattice Boltzmann Method (LBM), a mesoscopic method based on kinetic theory, is introduced to solve complex problems, such as hemodynamics. The results show that the flow velocity reduces considerably by increasing the magnetic field intensity, and the flow separation area is minimized by the increase of magnetic field strength. The LBM with BGK collision model has shown good results in terms of precision. Finally, this literature review has revealed a number of potential avenues for further research. Suggestions for future works are proposed accordingly.

Impact of the Spatial Velocity Inlet Distribution on the Hemodynamics of the Thoracic Aorta

The impact of the distribution in space of the inlet velocity in the numerical simulations of the hemodynamics in the thoracic aorta is systematically investigated. A real healthy aorta geometry, for which in-vivo measurements are available, is considered. The distribution is modeled through a truncated cone shape, which is a suitable approximation of the real one downstream of a trileaflet aortic valve during the systolic part of the cardiac cycle. The ratio between the upper and the lower base of the truncated cone and the position of the center of the upper base are selected as uncertain parameters. A stochastic approach is chosen, based on the generalized Polynomial Chaos expansion, to obtain accurate response surfaces of the quantities of interest in the parameter space. The selected parameters influence the velocity distribution in the ascending aorta. Consequently, effects on the wall shear stress are observed, confirming the need to use patient-specific inlet conditions if interested in the hemodynamics of this region. The surface base ratio is globally the most important parameter. Conversely, the impact on the velocity and wall shear stress in the aortic arch and descending aorta is almost negligible.

Is spontaneous coronary artery dissection (SCAD) related to local anatomy and hemodynamics? An exploratory study

AIMS

Spontaneous coronary artery dissection (SCAD) is an increasingly diagnosed cause of myocardial infarction with unclear pathophysiology. The aim of the study was to test if vascular segments site of SCAD present distinctive local anatomy and hemodynamic profiles.

METHODS

Coronary arteries with spontaneously healed SCAD (confirmed by follow-up angiography) underwent three-dimensional reconstruction, morphometric analysis with definition of vessel local curvature and torsion, and computational fluid dynamics (CFD) simulations with derivation of time-averaged wall shear stress (TAWSS) and topological shear variation index (TSVI). The (reconstructed) healed proximal SCAD segment was visually inspected for co-localization with curvature, torsion, and CFD-derived quantities hot spots.

RESULTS

Thirteen vessels with healed SCAD underwent the morpho-functional analysis. Median time between baseline and follow-up coronary angiograms was 57 (interquartile range [IQR] 45-95) days. In seven cases (53.9%), SCAD was classified as type 2b and occurred in the left anterior descending artery or near a bifurcation. In all cases (100%), at least one hot spot co-localized within the healed proximal SCAD segment, in 9 cases (69.2%) ≥3 hot spots were identified. Healed SCAD in proximity of a coronary bifurcation presented lower TAWSS peak values (6.65 [IQR 6.20-13.2] vs. 3.81 [2.53-5.17] Pa, p = 0.008) and hosted less frequently TSVI hot spots (100% vs. 57.1%, p = 0.034).

CONCLUSION

Vascular segments of healed SCAD were characterized by high curvature/torsion and WSS profiles reflecting increased local flow disturbances. Hence, a pathophysiological role of the interaction between vessel anatomy and shear forces in SCAD is hypothesized.

On the nonlinear relationship between wall shear stress topology and multi-directionality in coronary atherosclerosis

BACKGROUND AND OBJECTIVE

In this paper we investigate twelve multi-directional/topological wall shear stress (WSS) derived metrics and their relationships with the formation of coronary plaques in both computational fluid dynamics (CFD) and dynamic fluid-structure interaction (FSI) frameworks. While low WSS is one of the most established biomechanical markers associated with coronary atherosclerosis progression, alone it is limited. Multi-directional and topological WSS derived metrics have been shown to be important in atherosclerosis related mechanotransduction and near-wall transport processes. However, the relationships between these twelve WSS metrics and the influence of both FSI simulations and coronary dynamics is understudied.

METHODS

We first investigate the relationships between these twelve WSS derived metrics, stenosis percentage and lesion length through a parametric, transient CFD study. Secondly, we extend the parametric study to FSI, both with and without the addition of coronary dynamics, and assess their correlations. Finally, we present the case of a patient who underwent invasive coronary angiography and optical coherence tomography imaging at two time points 18 months apart. Associations between each of the twelve WSS derived metrics in CFD, static FSI and dynamic FSI simulations were assessed against areas of positive/negative vessel remodelling, and changes in plaque morphology.

RESULTS

22-32% stenosis was the threshold beyond which adverse multi-directional/topological WSS results. Each metric produced a different relationship with changing stenoses and lesion length. Transient haemodynamics was impacted by coronary dynamics, with the topological shear variation index suppressed by up to 94%. These changes appear more critical at smaller stenosis levels, suggesting coronary dynamics could play a role in the earlier stages of atherosclerosis development. In the patient case, both dynamics and FSI vs CFD changes altered associations with measured changes in plaque morphology. An appendix of the linear fits between the various FSI- and CFD-based simulations is provided to assist in scaling CFD-based results to resemble the compliant walled characteristics of FSI more accurately.

CONCLUSIONS

These results highlight the potential for coronary dynamics to alter multi-directional/topological WSS metrics which could impact associations with changes in coronary atherosclerosis over time. These results warrant further investigation in a wider range of morphological settings and longitudinal cohort studies in the future.

Smartphone-based particle image velocimetry for cardiovascular flows applications: A focus on coronary arteries

An experimental set-up is presented for the in vitro characterization of the fluid dynamics in personalized phantoms of healthy and stenosed coronary arteries. The proposed set-up was fine-tuned with the aim of obtaining a compact, flexible, low-cost test-bench for biomedical applications. Technically, velocity vector fields were measured adopting a so-called smart-PIV approach, consisting of a smartphone camera and a low-power continuous laser (30 mW). Experiments were conducted in realistic healthy and stenosed 3D-printed phantoms of left anterior descending coronary artery reconstructed from angiographic images. Time resolved image acquisition was made possible by the combination of the image acquisition frame rate of last generation commercial smartphones and the flow regimes characterizing coronary hemodynamics (velocities in the order of 10 cm/s). Different flow regimes (Reynolds numbers ranging from 20 to 200) were analyzed. The smart-PIV approach was able to provide both qualitative flow visualizations and quantitative results. A comparison between smart-PIV and conventional PIV (i.e., the gold-standard experimental technique for bioflows characterization) measurements showed a good agreement in the measured velocity vector fields for both the healthy and the stenosed coronary phantoms. Displacement errors and uncertainties, estimated by applying the particle disparity method, confirmed the soundness of the proposed smart-PIV approach, as their values fell within the same range for both smart and conventional PIV measured data (≈5% for the normalized estimated displacement error and below 1.2 pixels for displacement uncertainty). In conclusion, smart-PIV represents an easy-to-implement, low-cost methodology for obtaining an adequately robust experimental characterization of cardiovascular flows. The proposed approach, to be intended as a proof of concept, candidates to become an easy-to-handle test bench suitable for use also outside of research labs, e.g., for educational or industrial purposes, or as first-line investigation to direct and guide subsequent conventional PIV measurements.

Wall Shear Stress Differences Between Arterial and Venous Coronary Artery Bypass Grafts One Month After Surgery

Although coronary artery bypass graft (CABG) surgery is a well-established intervention, graft failure can occur, and the underlying mechanisms remain incompletely understood. The purpose of this prospective study is to utilize computational fluid dynamics (CFD) to investigate how graft hemodynamics one month post surgery may vary among graft types, which have different long-term patency rates. Twenty-four grafts from 10 participants (64.6 ± 8.5 years, 9 men) were scanned with coronary CT angiography and 4D flow MRI one month after CABG surgery. Grafts included 10 left internal mammary arteries (LIMA), 3 radial arteries (RA), and 11 saphenous vein grafts (SVG). Image-guided CFD was used to quantify blood flow rate and wall area exposed to abnormal wall shear stress (WSS). Arterial grafts had a lower abnormal WSS area than venous grafts (17.9% vs. 70.1%; p = 0.001), and a similar trend was observed for LIMA vs. SVG (13.8% vs. 70.1%; p = 0.001). Abnormal WSS area correlated positively to lumen diameter (p < 0.001) and negatively to flow rate (p = 0.001). This CFD study is the first of its kind to prospectively reveal differences in abnormal WSS area 1 month post surgery among CABG types, suggesting that WSS may influence the differential long-term graft failure rates observed among these groups.

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.

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 post-operative 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. The present study aims to analyze the impact that a deployed coronary stent has on the WSS topological skeleton. Computational fluid dynamics 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 proximal to the struts and a WSS expansion action distal to the struts. This WSS action pattern was independent from the stent design. In conclusions, these findings could contribute to a deeper understanding of the hemodynamic-driven processes underlying ST and ISR.

Assessing mechanical vibration-altered wall shear stress in digital arteries

The aim of this study is to implement and validate a method for assessing acute vibration-altered Wall Shear Stress (WSS) in the proper volar digital artery of the non-exposed left forefinger when subjecting the right hand to mechanical vibration. These changes of WSS may be involved in Vibration White Finger. Hence, an experimental device was set-up to link a vibration shaker and an ultra-high frequency ultrasound scanner. The Womersley-based WSS was computed by picking up the maximum velocity from pulse Wave Doppler measurements and extracting the artery diameter from B-mode images through an in-house image processing technique. The parameters of the former method were optimised on numerical ultrasound phantoms of cylindrical and lifelike arteries. These phantoms were computed with the FIELD II and FOCUS platforms which mimicked our true ultrasound device. The Womersley-based WSS were compared to full Fluid Structure Interaction (FSI) and rigid wall models built from resonance magnetic images of a volunteer-specific forefinger artery. Our FSI model took into account the artery's surrounding tissues. The diameter computing procedure led to a bias of 4%. The Womersley-based WSS resulted in misestimating the FSI model by roughly 10% to 20%. No difference was found between the rigid wall computational model and FSI simulations. Regarding the WSS measured on a group of 20 volunteers, the group-averaged basal value was 3 Pa, while the vibration-altered WSS was reduced to 1 Pa, possibly triggering intimal hyperplasia mechanisms and leading to the arterial stenoses encountered in patients suffering from vibration-induced Raynaud's syndrome.

In Vivo Based Fluid-Structure Interaction Biomechanics of the Left Anterior Descending Coronary Artery

A fluid-structure interaction-based biomechanical model of the entire left anterior descending coronary artery is developed from in vivo imaging via the finite element method in this paper. Included in this investigation is ventricle contraction, three-dimensional motion, all angiographically visible side branches, hyper/viscoelastic artery layers, non-Newtonian and pulsatile blood flow, and the out-of-phase nature of blood velocity and pressure. The fluid-structure interaction model is based on in vivo angiography of an elite athlete's entire left anterior descending coronary artery where the influence of including all alternating side branches and the dynamical contraction of the ventricle is investigated for the first time. Results show the omission of side branches result in a 350% increase in peak wall shear stress and a 54% decrease in von Mises stress. Peak von Mises stress is underestimated by up to 80% when excluding ventricle contraction and further alterations in oscillatory shear indices are seen, which provide an indication of flow reversal and has been linked to atherosclerosis localization. Animations of key results are also provided within a video abstract. We anticipate that this model and results can be used as a basis for our understanding of the interaction between coronary and myocardium biomechanics. It is hoped that further investigations could include the passive and active components of the myocardium to further replicate in vivo mechanics and lead to an understanding of the influence of cardiac abnormalities, such as arrythmia, on coronary biomechanical responses.

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.

Role of the left coronary artery geometry configuration in atherosusceptibility: CFD simulations considering sPTT model for blood

The achievement of clinically viable methodologies to simulate the hemodynamics in patient-specific coronary arteries is still a major challenge. Therefore, the novelty of this work is attained by the introduction of the viscoelastic property of blood in the numerical simulations, to study the role of the left coronary artery (LCA) geometry configuration in the atherosusceptibility. Apparently healthy patients were used and four different methodologies were tested. The methodology giving the most accurate results at the same time of having the lowest computational time is the one considering the viscoelastic property of blood and computational fluid dynamics. A Pearson correlation analysis was used to highlight relationships between geometric configuration and hemodynamic descriptors based on the simulated wall shear stress (WSS). The left main stem (LMS) has the greatest atherosusceptibility followed by the left anterior descending artery (LAD) since the relative residence time (RRT) average values are 3.81 and 3.70 Pa^-1, respectively. The geometric parameters with relevant contribution to directional flow change are the cross-sectional areas, especially the one of LMS segment (A_LMS), and the curvature of LMS segment. For LMS and LAD segments, when A_LMS increases, blood flow disturbance (r = 0.81 in LMS and r = 0.74 in LAD) and atherosusceptibility (r = 0.84 in LMS and r = 0.85 in LAD) increases. When the LMS curvature decreases, the WSS magnitude (r = 0.80 in LMS and r = 0.83 in LAD) decreases, and disturbance (r=-0.80 in LMS and r=-0.91 in LAD) and atherosusceptibility (r=-0.74 in LMS and r=-0.74 in LAD) increases.

Blood Flow Modeling in Coronary Arteries: A Review

Atherosclerosis is one of the main causes of cardiovascular events, namely, myocardium infarction and cerebral stroke, responsible for a great number of deaths every year worldwide. This pathology is caused by the progressive accumulation of low-density lipoproteins, cholesterol, and other substances on the arterial wall, narrowing its lumen. To date, many hemodynamic studies have been conducted experimentally and/or numerically; however, this disease is not yet fully understood. For this reason, the research of this pathology is still ongoing, mainly, resorting to computational methods. These have been increasingly used in biomedical research of atherosclerosis because of their high-performance hardware and software. Taking into account the attempts that have been made in computational techniques to simulate realistic conditions of blood flow in both diseased and healthy arteries, the present review aims to give an overview of the most recent numerical studies focused on coronary arteries, by addressing the blood viscosity models, and applied physiological flow conditions. In general, regardless of the boundary conditions, numerical studies have been contributed to a better understanding of the development of this disease, its diagnosis, and its treatment.