Fluid-structure interaction analysis of the left coronary artery with variable angulation

Patient-Specific Numerical Simulations of Coronary Artery Hemodynamics and Biomechanics: A Pathway to Clinical Use

PURPOSE

Numerical models that simulate the behaviors of the coronary arteries have been greatly improved by the addition of fluid-structure interaction (FSI) methods. Although computationally demanding, FSI models account for the movement of the arterial wall and more adequately describe the biomechanical conditions at and within the arterial wall. This offers greater physiological relevance over Computational Fluid Dynamics (CFD) models, which assume the walls do not move or deform. Numerical simulations of patient-specific cases have been greatly bolstered by the use of imaging modalities such as Computed Tomography Angiography (CTA), Magnetic Resonance Imaging (MRI), Optical Coherence Tomography (OCT), and Intravascular Ultrasound (IVUS) to reconstruct accurate 2D and 3D representations of artery geometries. The goal of this study was to conduct a comprehensive review on CFD and FSI models on coronary arteries, and evaluate their translational potential.

METHODS

This paper reviewed recent work on patient-specific numerical simulations of coronary arteries that describe the biomechanical conditions associated with atherosclerosis using CFD and FSI models. Imaging modality for geometry collection and clinical applications were also discussed.

RESULTS

Numerical models using CFD and FSI approaches are commonly used to study biomechanics within the vasculature. At high temporal and spatial resolution (compared to most cardiac imaging modalities), these numerical models can generate large amount of biomechanics data.

CONCLUSIONS

Physiologically relevant FSI models can more accurately describe atherosclerosis pathogenesis, and help to translate biomechanical assessment to clinical evaluation.

Bio-inspired microfluidics: A review

Biomicrofluidics, a subdomain of microfluidics, has been inspired by several ideas from nature. However, while the basic inspiration for the same may be drawn from the living world, the translation of all relevant essential functionalities to an artificially engineered framework does not remain trivial. Here, we review the recent progress in bio-inspired microfluidic systems via harnessing the integration of experimental and simulation tools delving into the interface of engineering and biology. Development of "on-chip" technologies as well as their multifarious applications is subsequently discussed, accompanying the relevant advancements in materials and fabrication technology. Pointers toward new directions in research, including an amalgamated fusion of data-driven modeling (such as artificial intelligence and machine learning) and physics-based paradigm, to come up with a human physiological replica on a synthetic bio-chip with due accounting of personalized features, are suggested. These are likely to facilitate physiologically replicating disease modeling on an artificially engineered biochip as well as advance drug development and screening in an expedited route with the minimization of animal and human trials.

The role of retinal vessel geometry as an indicator of systemic arterial stiffness assessed by cardio-ankle vascular index

Objective

To determine whether retinal vessel geometry is associated with systemic arterial stiffness, as determined by the cardio-ankle vascular index (CAVI).

Methods

This single-center retrospective cross-sectional study included 407 eyes of 407 subjects who underwent routine health exams, including CAVI and fundus photography. Retinal vessel geometry was measured using a computer-assisted program (Singapore "I" Vessel Assessment). Subjects were classified into two groups based on CAVI values: high CAVI (≥9) or low CAVI (<9). The main outcome measures included the association of retinal vessel geometry and CAVI value evaluated using multivariable logistic regression models.

Results

Three hundred forty-three subjects (343, 84.3%) were in the low CAVI group, and 64 (15.7%) subjects were in the high CAVI group. Multivariable logistic linear regression analyses adjusted for age, sex, body mass index, smoking status, mean arterial pressure, and the presence of hypertension, diabetes mellitus, and dyslipidemia showed a significant association between high CAVI values and the following retinal vessel geometry parameters: central retinal arteriolar equivalent caliber (CRAE; adjusted odds ratio [AOR], 0.95; 95% confidence interval [CI], 0.89-1.00; P = 0.043), fractal dimension of arteriolar network (FDa; AOR, 4.21 × 10^-4; 95% CI, 2.32 × 10^-7-0.77; P = 0.042), and arteriolar branching angle (BAa; AOR, 0.96; 95% CI, 0.93-0.99; P = 0.007).

Conclusions

Increased systemic arterial stiffness had a significant association with retinal vessel geometry related to arterial narrowing (CRAE), less branching complexity of the arterial tree (FDa), and acute arteriolar bifurcation (BAa).

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.

Time-resolved simulation of blood flow through left anterior descending coronary artery: effect of varying extent of stenosis on hemodynamics

BACKGROUND AND OBJECTIVES

Real-time blood flow variation is crucial for understanding the dynamic development of coronary atherosclerosis. The main objective of this study is to investigate the effect of varying extent of stenosis on the hemodynamic features in left anterior descending coronary artery.

METHODS

Various Computational fluid dynamics (CFD) models were constructed with patient-specific CT image data, using actual fractional flow reserve (FFR) as boundary conditions to provide a real-time quantitative description of hemodynamic properties. The hemodynamic parameters, such as the local and instantaneous wall shear stress (WSS), oscillating shear index (OSI) and relative residence time (RRT), blood flow velocity and pressure drop during various phases of cardiac cycle were provided in detail.

RESULTS

There was no evident variation in hemodynamic parameters in the cases of less than 50% stenosis while there were abrupt and dramatic changes in hemodynamics when the stenosis aggravated from 60 to 70%. Furthermore, when the stenosis was beyond 70%, there existed substantial pressure difference, WSS, and blood flow velocity in the center of the stenosis. Although OSI and RRT increased along with the aggravation of stenosis, they appeared with obvious abnormalities across all cases, even in mild stenosis.

CONCLUSION

The simulation could present a dynamic and comprehensive profile of how hemodynamic parameters vary in accordance with divergent severities of stenosis, which could serve as an effective reference for the clinicians to have a deeper insight into the pathological mechanism of coronary atherosclerosis and stenosis.

A spatiotemporal analysis of the left coronary artery biomechanics using fluid-structure interaction models

Biomechanics plays a critical role in coronary artery disease development. FSI simulation is commonly used to understand the hemodynamics and mechanical environment associated with atherosclerosis pathology. To provide a comprehensive characterization of patient-specific coronary biomechanics, an analysis of FSI simulation in the spatial and temporal domains was performed. In the current study, a three-dimensional FSI model of the LAD coronary artery was built based on a patient-specific geometry using COMSOL Multiphysics. The effect of myocardial bridging was simulated. Wall shear stress and its derivatives including time-averaged wall shear stress, wall shear stress gradient, and OSI were calculated across the cardiac cycle in multiple locations. Arterial wall strain (radial, circumferential, and longitudinal) and von Mises stress were calculated. To assess perfusion, vFFR was calculated. The results demonstrated the FSI model could identify regional and transient differences in biomechanical parameters within the coronary artery. The addition of myocardial bridging caused a notable change in von Mises stress and an increase in arterial strain during systole. The analysis performed in this manner takes greater advantage of the information provided in the space and time domains and can potentially assist clinical evaluation.

Influence of Rigid-Elastic Artery Wall of Carotid and Coronary Stenosis on Hemodynamics

Cardiovascular system abnormalities can result in serious health complications. By using the fluid-structure interaction (FSI) procedure, a comprehensive realistic approach can be employed to accurately investigate blood flow coupled with arterial wall response. The hemodynamics was investigated in both the coronary and carotid arteries based on the arterial wall response. The hemodynamics was estimated based on the numerical simulation of a comprehensive three-dimensional non-Newtonian blood flow model in elastic and rigid arteries. For stenotic right coronary artery (RCA), it was found that the maximum value of wall shear stress (WSS) for the FSI case is higher than the rigid wall. On the other hand, for the stenotic carotid artery (CA), it was found that the maximum value of WSS for the FSI case is lower than the rigid wall. Moreover, at the peak systole of the cardiac cycle (0.38 s), the maximum percentage of arterial wall deformation was found to be 1.9%. On the other hand, for the stenotic carotid artery, the maximum percentage of arterial wall deformation was found to be 0.46%. A comparison between FSI results and those obtained by rigid wall arteries is carried out. Findings indicate slight differences in results for large-diameter arteries such as the carotid artery. Accordingly, the rigid wall assumption is plausible in flow modeling for relatively large diameters such as the carotid artery. Additionally, the FSI approach is essential in flow modeling in small diameters.

Angular difference in human coronary artery governs endothelial cell structure and function

Blood vessel branch points exhibiting oscillatory/turbulent flow and lower wall shear stress (WSS) are the primary sites of atherosclerosis development. Vascular endothelial functions are essentially dependent on these tangible biomechanical forces including WSS. Herein, we explored the influence of blood vessel bifurcation angles on hemodynamic alterations and associated changes in endothelial function. We generated computer-aided design of a branched human coronary artery followed by 3D printing such designs with different bifurcation angles. Through computational fluid dynamics analysis, we observed that a larger branching angle generated more complex turbulent/oscillatory hemodynamics to impart minimum WSS at branching points. Through the detection of biochemical markers, we recorded significant alteration in eNOS, ICAM1, and monocyte attachment in EC grown in microchannel having 60^o vessel branching angle which correlated with the lower WSS. The present study highlights the importance of blood vessel branching angle as one of the crucial determining factors in governing atherogenic-endothelial dysfunction.

In silico and in vitro investigations reveal angular differences in the blood vessel branching points differentially alter the hemodynamics to impact endothelial structure and function.

Numerical investigation of different viscosity models on pulsatile blood flow of thoracic aortic aneurysm (TAA) in a patient-specific model

Aortic aneurysm is one of the most common aortic diseases that can lead to unfortunate consequences. Numerical simulations have an important role in the prediction of the aftereffects of vascular diseases including aneurysm. In this research, numerical simulation of pulsatile blood flow is performed for a 3-dimensional patient-specific model of a thoracic aortic aneurysm (TAA). Since the choice of blood viscosity model may have a significant impact on the simulation results, the effects of four non-Newtonian models of blood viscosity namely Carreau, Casson, Herschel-Bulkley, power low, and the Newtonian model on the wall shear stress (WSS) distribution, shear rate, and oscillatory shear index (OSI) have been analyzed. Simulation results showed that all the non-Newtonian and Newtonian models generally, predict similar patterns for blood flow and shear rate. At high flow rates in the cardiac cycle, the WSS value for all the models are similar to each other except for the power-law model due to the shear thinning behavior. All models predict high values of OSI on the inner wall of the ascending aorta and broad areas of the inner wall of the aneurysm sac. However, the Newtonian model predicts the OSI less than the non-Newtonian models in some areas of the aneurysm sac. Results indicated that the Newtonian model generally can predict the hemodynamic parameters of the blood flow similar to the non-Newtonian but for more precise analysis and to predict the regions prone to rupture and atherosclerosis, choosing a proper non-Newtonian model is recommended.

A new and automated risk prediction of coronary artery disease using clinical endpoints and medical imaging-derived patient-specific insights: protocol for the retrospective GeoCAD cohort study

INTRODUCTION

Coronary artery disease (CAD) is the leading cause of death worldwide. More than a quarter of cardiovascular events are unexplained by current absolute cardiovascular disease risk calculators, and individuals without clinical risk factors have been shown to have worse outcomes. The 'anatomy of risk' hypothesis recognises that adverse anatomical features of coronary arteries enhance atherogenic haemodynamics, which in turn mediate the localisation and progression of plaques. We propose a new risk prediction method predicated on CT coronary angiography (CTCA) data and state-of-the-art machine learning methods based on a better understanding of anatomical risk for CAD. This may open new pathways in the early implementation of personalised preventive therapies in susceptible individuals as a potential key in addressing the growing burden of CAD.

METHODS AND ANALYSIS

GeoCAD is a retrospective cohort study in 1000 adult patients who have undergone CTCA for investigation of suspected CAD. It is a proof-of-concept study to test the hypothesis that advanced image-derived patient-specific data can accurately predict long-term cardiovascular events. The objectives are to (1) profile CTCA images with respect to variations in anatomical shape and associated haemodynamic risk expressing, at least in part, an individual's CAD risk, (2) develop a machine-learning algorithm for the rapid assessment of anatomical risk directly from unprocessed CTCA images and (3) to build a novel CAD risk model combining traditional risk factors with these novel anatomical biomarkers to provide a higher accuracy CAD risk prediction tool.

ETHICS AND DISSEMINATION

The study protocol has been approved by the St Vincent's Hospital Human Research Ethics Committee, Sydney-2020/ETH02127 and the NSW Population and Health Service Research Ethics Committee-2021/ETH00990. The project outcomes will be published in peer-reviewed and biomedical journals, scientific conferences and as a higher degree research thesis.

A Study of the Fluid–Structure Interaction of the Plaque Circumferential Distribution in the Left Coronary Artery

Atherosclerotic plaques within the coronary arteries can prevent blood from flowing to downstream tissues, causing coronary heart disease and a myocardial infarction over time. The degree of stenosis is an important reference point during percutaneous coronary intervention (PCI). However, clinically, patients with the same degree of stenosis exhibit different degrees of disease severity. To investigate the connection between this phenomenon and the plaque circumferential distribution, in this paper, four models with different plaque circumferential locations were made based on the CT data. The blood in the coronary arteries was simulated using the fluid–structure interaction method in ANSYS Workbench software. The results showed that the risk of plaque rupture was less affected by the circumferential distribution of plaque, and the distribution of blood in each branch was affected by the circumferential distribution of plaque. Low TAWSS areas were found posterior to the plaque, and the TAWSS < 0.4 Pa area was ranked from highest to lowest in each model species: plaque on the side away from the left circumflex branch, plaque on the side away from the heart; plaque on the side close to the heart; and plaque on the side close to the left circumflex branch. The same trend was also found in the OSI. It was concluded that the circumferential distribution of plaques affects their further development. This finding will be useful for clinical treatment.

Index of microcirculatory resistance: state-of-the-art and potential applications in computational simulation of coronary artery disease

The dysfunction of coronary microcirculation is an important cause of coronary artery disease (CAD). The index of microcirculatory resistance (IMR) is a quantitative evaluation of coronary microcirculatory function, which provides a significant reference for the prediction, diagnosis, treatment, and prognosis of CAD. IMR also plays a key role in investigating the interaction between epicardial and microcirculatory dysfunctions, and is closely associated with coronary hemodynamic parameters such as flow rate, distal coronary pressure, and aortic pressure, which have been widely applied in computational studies of CAD. However, there is currently a lack of consensus across studies on the normal and pathological ranges of IMR. The relationships between IMR and coronary hemodynamic parameters have not been accurately quantified, which limits the application of IMR in computational CAD studies. In this paper, we discuss the research gaps between IMR and its potential applications in the computational simulation of CAD. Computational simulation based on the combination of IMR and other hemodynamic parameters is a promising technology to improve the diagnosis and guide clinical trials of CAD.

The effects of plaque morphological characteristics on the post-stenotic flow in left main coronary artery bifurcation

Local post-stenotic hemodynamics has critical influence in the atherosclerotic plaque progression occurring in susceptible arterial sites, in particular the left main coronary artery (LMCA) bifurcation. Understanding the effects of plaque morphological characteristics: stenosis severity (SS), eccentricity index (EI) and lesion length (LL) on the post-stenotic flow behavior can significantly improve treatment planning. In order to investigate these effects, we have employed computational fluid dynamics (CFD) simulations in twenty computer-generated and five patient-specific LMCA models and the hemodynamic parameters: velocity, pressure (P), wall pressure gradient (WPG), wall shear stress (WSS), time averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT) and helicity intensity (h₂) were analyzed. Our results revealed that the effect of stenosis eccentricity varied significantly for different values of stenosis severity and lesion length. Regions with low WSS, low TAWSS and high RRT were more prominent in models having higher stenosis severity. For smaller lesion length, at low and moderate stenosis severity, surface area with low TAWSS and high RRT decreased with increasing eccentricity index, whereas for high stenosis severity models, low TAWSS region and average RRT values increased with eccentricity. However, for models with longer lesion length, regions with high OSI and RRT overall increased gradually with eccentricity. The helicity intensity (h₂) of all models remained very low except at the most eccentric model with longer lesion length. The presence of very high helical flow in this model suggests the possibility of atheroprotective flow. It can be concluded that all plaque morphological characteristics covered under this investigation play an important role in plaque progression.

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.

Assessment with clinical data of a coupled bio-hemodynamics numerical model to predict leukocyte adhesion in coronary arteries

Numerical simulations of coupled hemodynamics and leukocyte transport and adhesion inside coronary arteries have been performed. Realistic artery geometries have been obtained for a set of four patients from intravascular ultrasound and angiography images. The numerical model computes unsteady three-dimensional blood hemodynamics and leukocyte concentration in the blood. Wall-shear stress dependent leukocyte adhesion is also computed through agent-based modeling rules, fully coupled to the hemodynamics and leukocyte transport. Numerical results have a good correlation with clinical data. Regions where high adhesion is predicted by the simulations coincide to a good approximation with artery segments presenting plaque increase, as documented by clinical data from baseline and six-month follow-up exam of the same artery. In addition, it is observed that the artery geometry and, in particular, the tortuosity of the centerline are a primary factor in determining the spatial distribution of wall-shear stress, and of the resulting leukocyte adhesion patterns. Although further work is required to overcome the limitations of the present model and ultimately quantify plaque growth in the simulations, these results are encouraging towards establishing a predictive methodology for atherosclerosis progress.

Effect of stenosis and dilatation on the hemodynamic parameters associated with left coronary artery

BACKGROUND AND OBJECTIVE

The main objective of the work is to examine the curvature effects of stenosis/dilatation region pertaining to left coronary artery. The hemodymamic features during the cardiac cycle is thoroughly examined.

METHODS

A numerical fluid structure interaction model incorporating multi- layered elastic artery wall, non-Newtonian blood viscosity and pulsating boundary conditions is developed. The composite arterial wall consists of a thin layer tunica intima, atheroma and a thick wall. Higher stiffness of atheroma is captured by using higher Young's modulus. The CFD and FSI models are validated with available experimental and analytical data. Computations are done with five different non-Newtonian models and arterial wall with various elasticity levels. The local and time averaged WSS, velocity contours downstream of stenosis, wall pressure and pressure drop during various phases of cardiac cycle are provided in detail.

RESULTS

The influence of non-Newtonian effects of blood viscosity is found to be significant especially at stenosis regions. The flexible wall caused wall deformation and the associated flow and pressure wave propagation affecting WSS and pressure drop compared to the rigid wall. Flow recirculation is noticed at stenosis downstream locations and its strength increases with increased severity of the stenosis. A stenosis is characterised by a sudden drop in wall pressure and a slower two stage recovery during peak velocity periods of the cardiac cycle.

CONCLUSIONS

The pressure drop, local WSS at stenosis centre, and radial velocity increase are significantly higher for stenosis cases and the effect is severe during peak diastole. The variation in hemodynamic parameters is found to be less significant for dilatation. Significantly lower WSS is noticed for the recirculation regions downstream of stenosis which can enhance the tendency for monocytes to attach to the endothelium. The radius of curvature of the stenosis is found to be the most sensitive parameter affecting the hemodynamic characteristics rather than the detailed geometry of the stenosis. The main effect of variation of artery wall stiffness is noted at recirculation regions present downstream of stenosis. The results from the study may be useful for predicting wall shear stress signatures associated with stenosis/dilatation changes and the management of specific cases.

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.

Numerical study of hemodynamics in a complete coronary bypass with venous and arterial grafts and different degrees of stenosis

Cardiovascular diseases are among the leading causes of death in the world. The coronary blockage is one of most common types of these diseases that in the majority of cases has been treated by bypass surgery. In the bypass surgery, a graft is implemented to alter the blocked coronary and allow the blood supply process. The hemodynamic characteristics of the bypass strongly depend on the geometry and mechanical properties of the graft. In the present study, the fluid-structure interaction (FSI) analysis is conducted to investigate the bypass performance for a thoracic artery as well as a saphenous vein graft. Blood flow introduces a pressure on the walls of the graft which behaves as a hyperelastic material. A complete coronary bypass with stenosis degrees of 70% and 100% is modeled. To consider the nonlinear stress-strain behavior of the grafts, a five parameter Mooney-Rivlin hyperplastic model is implemented for the structural analysis and blood is assumed to behave as a Newtonian fluid. The simulations are performed for a structured grid to solve the governing equations using finite element method (FEM). The results show that wall shear stress (WSS) for saphenous vein is larger than that of thoracic artery while the total deformation of the thoracic artery is larger compared to the saphenous vein. Also, for the venous grafts or lower stenosis degree, the oscillatory shear index (OSI) is higher at both left and right anastomoses meaning that venous grafts as well as lower degree of stenosis are more critical in terms of restenosis.

Hemodynamic impacts of hematocrit level by two-way coupled FSI in the left coronary bifurcation

Cardiovascular disease is now under the influence of several factors that encourage researchers to investigate the flow of these vessels. Oscillation influences the blood circulation in the volume of red blood cells (RBC) strongly. Therefore, in this study, its effects have been considered on hemodynamic parameters in the elastic wall and coronary bifurcation. In this study, a 3D geometry of non-Newtonian and pulsatile blood circulation is considered in the left coronary artery bifurcation. The Casson model with various hematocrits is analyzed in elastic and rigid walls. The wall shear stress (WSS) cannot show the stenosis artery alone, therefore, the oscillatory shear index (OSI) is represented as a hemodynamic parameter of WSS individually of time. The results are determined using two-way fluid-structure interaction (FSI) coupling method using an arbitrary Lagrangian-Eulerian method. The most prominent difference in velocity happened in the bifurcation and at hematocrit 30 with yield stress 6.59E-04 Pa. The backflow and vortex flow in the LCx branch grown with increasing shear rates. The likelihood of plaque generation at the ending of the LM branch is observed in hematocrits 10 and 20, while the WSS magnitude is normal in the hematocrit 60 with the greatest yield stress in the bifurcation. The shear stress among the rigid and elastic models is the highest at the ending of the LM branch. The wall shear stress magnitude among the models decreased at most of 24.49% by dividing the flow. Time-independent results for models showed that there is the highest value of OSI at the bifurcation, which then quickly dropped.

Coronary Smooth Muscle Cell Calcium Dynamics: Effects of Bifurcation Angle on Atheroprone Conditions

This work investigates the effect of arterial bifurcation angulation on atherosclerosis development through in-silico simulations of coupled cell dynamics. The computational model presented here combines cellular pathways, fluid dynamics, and physiologically-realistic vessel geometries as observed in the human vasculature. The coupled cells model includes endothelial cells (ECs) and smooth muscle cells (SMCs) with ion dynamics, hetero and homotypic coupling, as well as electro-diffusive coupling. Three arterial bifurcation surface models were used in the coupled cells simulations. All three simulations showed propagating waves of Ca²⁺ in both the SMC and EC layers, following the introduction of a luminal agonist, in this case ATP. Immediately following the introduction of ATP concentration Ca²⁺ waves propagate from the area of high ATP toward the areas of low ATP concentration, forming complex patterns where waves interact with eachother, collide and fade. These dynamic phenomena are repeated with a series of waves of slower velocity. The underlying motivation of this research was to examine the macro-scale phenomena, given that the characteristic length scales of atherosclerotic plaques are much larger than a single cell. The micro-scale dynamics were modeled on macro-scale arterial bifurcation surfaces containing over one million cells. The results of the simulations presented here suggest that susceptibility to atherosclerosis development depends on the bifurcation angulation. In conjunction with findings reported in the literature, the simulation results demonstrate that arterial bifurcations containing wider angles have a more prominent influence on the coupled cells pathways associated with the development of atherosclerosis, by means of disturbed flow and lower SMC Ca²⁺ concentrations. The discussion of the results considers the findings of this research within the context of the potential link between information transport through frequency encoding of Ca²⁺ wave dynamics and development of atheroprone conditions.

Art care: A multi-modality coronary 3D reconstruction and hemodynamic status assessment software

BACKGROUND

Due to the incremental increase of clinical interest in the development of software that allows the 3-dimensional (3D) reconstruction and the functional assessment of the coronary vasculature, several software packages have been developed and are available today.

OBJECTIVE

Taking this into consideration, we have developed an innovative suite of software modules that perform 3D reconstruction of coronary arterial segments using different coronary imaging modalities such as IntraVascular UltraSound (IVUS) and invasive coronary angiography images (ICA), Optical Coherence Tomography (OCT) and ICA images, or plain ICA images and can safely and accurately assess the hemodynamic status of the artery of interest.

METHODS

The user can perform automated or manual segmentation of the IVUS or OCT images, visualize in 3D the reconstructed vessel and export it to formats, which are compatible with other Computer Aided Design (CAD) software systems. We employ finite elements to provide the capability to assess the hemodynamic functionality of the reconstructed vessels by calculating the virtual functional assessment index (vFAI), an index that corresponds and has been shown to correlate well to the actual fractional flow reserve (FFR) value.

RESULTS

All the modules of the proposed system have been thoroughly validated. In brief, the 3D-QCA module, compared to a successful commercial software of the same genre, presented very good correlation using several validation metrics, with a Pearson's correlation coefficient (R) for the calculated volumes, vFAI, length and minimum lumen diameter of 0.99, 0.99, 0.99 and 0.88, respectively. Moreover, the automatic lumen detection modules for IVUS and OCT presented very high accuracy compared to the annotations by medical experts with the Pearson's correlation coefficient reaching the values of 0.94 and 0.99, respectively.

CONCLUSIONS

In this study, we have presented a user-friendly software for the 3D reconstruction of coronary arterial segments and the accurate hemodynamic assessment of the severity of existing stenosis.

Atherosclerosis at arterial bifurcations: evidence for the role of haemodynamics and geometry

Atherosclerotic plaques are found at distinct locations in the arterial system, despite the exposure to systemic risk factors of the entire vascular tree. From the study of arterial bifurcation regions, emerges ample evidence that haemodynamics are involved in the local onset and progression of the atherosclerotic disease. This observed co-localisation of disturbed flow regions and lesion prevalence at geometrically predisposed districts such as arterial bifurcations has led to the formulation of a ‘haemodynamic hypothesis’, that in this review is grounded to the most current research concerning localising factors of vascular disease. In particular, this review focuses on carotid and coronary bifurcations because of their primary relevance to stroke and heart attack. We highlight reported relationships between atherosclerotic plaque location, progression and composition, and fluid forces at vessel’s wall, in particular shear stress and its ‘easier-tomeasure’ surrogates, i.e. vascular geometric attributes (because geometry shapes the flow) and intravascular flow features (because they mediate disturbed shear stress), in order to give more insight in plaque initiation and destabilisation. Analogous to Virchow’s triad for thrombosis, atherosclerosis must be thought of as subject to a triad of, and especially interactions among, haemodynamic forces, systemic risk factors, and the biological response of the wall.

Determination of hemodynamic risk for vascular disease in planar artery bifurcations

Understanding hemodynamics in blood circulation is crucial in order to unveil the mechanisms underlying the formation of stenosis and atherosclerosis. In fact, there are experimental evidences pointing out to the existence of some given vessel configurations that are more likely to develop the above mentioned pathologies. Along this manuscript, we performed an exhaustive investigation in a simplified model aiming to characterize by means of physical quantities those regions and configurations in vessel bifurcations that are more likely to develop such pathologies. The two-fold analysis is based, on the one hand, on numerical simulations (via CFD) and, on the other hand, on experiments realized in an ad-hoc designed polydimethylsiloxane (PDMS) channel with the appropriate parameters and appropriate fluid flows. The results obtained demonstrate that low velocity regions and low shear stress zones are located in the outer walls of bifurcations. In fact, we found that there is a critical range of bifurcation angles that is more likely to vascular disease than the others in correspondence with some experimental evidence. The effect of the inflow velocity on this critical range is also analyzed.

An investigation of correlation between left coronary bifurcation angle and hemodynamic changes in coronary stenosis by coronary computed tomography angiography-derived computational fluid dynamics

Background

To investigate the correlation between left coronary bifurcation angle and coronary stenosis as assessed by coronary computed tomography angiography (CCTA)-generated computational fluid dynamics (CFD) analysis when compared to the CCTA analysis of coronary lumen stenosis and plaque lesion length with invasive coronary angiography (ICA) as the reference method.

Methods

Thirty patients (22 males, mean age: 59±6.9 years) with calcified plaques at the left coronary artery were included in the study with all patients undergoing CCTA and ICA examinations. CFD simulation was performed to analyze hemodynamic changes to the left coronary artery models in terms of wall shear stress, wall pressure and flow velocity, with findings correlated to the coronary stenosis and degree of bifurcation angle. Calcified plaque length was measured in the left coronary artery with diagnostic value compared to that from coronary lumen and bifurcation angle assessments.

Results

Of 26 significant stenosis at left anterior descending (LAD) and 13 at left circumflex (LCx) on CCTA, only 14 and 5 of them were confirmed to be >50% stenosis at LAD and LCx respectively on ICA, resulting in sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV) of 100%, 52%, 49% and 100%. The mean plaque length was measured 5.3±3.6 and 4.4±1.9 mm at LAD and LCx, respectively, with diagnostic sensitivity, specificity, PPV and NPV being 92.8%, 46.7%, 61.9% and 87.5% for extensively calcified plaques. The mean bifurcation angle was measured 83.9±13.6º and 83.8±13.3º on CCTA and ICA, respectively, with no significant difference (P=0.98). The corresponding sensitivity, specificity, PPV and NPV were 100%, 78.6%, 84.2% and 100% based on bifurcation angle measurement on CCTA, 100%, 73.3%, 78.9% and 100% based on bifurcation angle measurements on ICA, respectively. Wall shear stress was noted to increase in the LAD and LCx models with significant stenosis and wider angulation (>80º), but demonstrated little or no change in most of the coronary models with no significant stenosis and narrower angulation (<80º).

Conclusions

This study further clarifies the relationship between left coronary bifurcation angle and significant stenosis, with angulation measurement serving as a more accurate approach than coronary lumen assessment or plaque lesion length for determining significant coronary stenosis. Left coronary bifurcation angle is suggested to be incorporated into coronary artery disease (CAD) assessment when diagnosing significant CAD.

Healthy and diseased coronary bifurcation geometries influence near-wall and intravascular flow: A computational exploration of the hemodynamic risk

Local hemodynamics has been identified as one main determinant in the onset and progression of atherosclerotic lesions at coronary bifurcations. Starting from the observation that atherosensitive hemodynamic conditions in arterial bifurcation are majorly determined by the underlying anatomy, the aim of the present study is to investigate how peculiar coronary bifurcation anatomical features influence near-wall and intravascular flow patterns. Different bifurcation angles and cardiac curvatures were varied in population-based, idealized models of both stenosed and unstenosed bifurcations, representing the left anterior descending (LAD) coronary artery with its diagonal branch. Local hemodynamics was analyzed in terms of helical flow and exposure to low/oscillatory shear stress by performing computational fluid dynamics simulations. Results show that bifurcation angle impacts lowly hemodynamics in both stenosed and unstenosed cases. Instead, curvature radius influences the generation and transport of helical flow structures, with smaller cardiac curvature radius associated to higher helicity intensity. Stenosed bifurcation models exhibit helicity intensity values one order of magnitude higher than the corresponding unstenosed cases. Cardiac curvature radius moderately affects near-wall hemodynamics of the stenosed cases, with smaller curvature radius leading to higher exposure to low shear stress and lower exposure to oscillatory shear stress. In conclusion, the proposed controlled benchmark allows investigating the effect of various geometrical features on local hemodynamics at the LAD/diagonal bifurcation, highlighting that cardiac curvature influences near wall and intravascular hemodynamics, while bifurcation angle has a minor effect.

Aortic valve leaflet wall shear stress characterization revisited: impact of coronary flow

Computational characterizations of aortic valve hemodynamics have typically discarded the effects of coronary flow. The objective of this study was to complement our previous fluid-structure interaction aortic valve model with a physiologic coronary circulation model to quantify the impact of coronary flow on aortic sinus hemodynamics and leaflet wall shear stress (WSS). Coronary flow suppressed vortex development in the two coronary sinuses and altered WSS magnitude and directionality on the three leaflets, with the most substantial differences occurring in the belly and tip regions.

Intravascular Ultrasound Classification of Plaque in Angiographic True Bifurcation Lesions of the Left Main Coronary Artery

Background:

Accurately, characterizing plaques is critical for selecting the optimal intervention strategy for the left main coronary artery (LMCA) bifurcation. Coronary angiography cannot precisely assess the location or nature of plaques in bifurcation lesions. Few intravascular ultrasound (IVUS) classification scheme has been reported for angiographic imaging of true bifurcation lesions of the unprotected LMCA thus far. In addition, the plaque composition at the bifurcation has not been elucidated. This study aimed to detect plaque composition at LMCA bifurcation lesions by IVUS.

Methods:

Fifty-eight patients were recruited. The location, concentricity or eccentricity, site of maximum thickness, and composition of plaques of the distal LMCA, ostial left anterior descending (LAD) coronary artery and, left circumflex (LCX) coronary artery were assessed using IVUS and described using illustrative diagrams.

Results:

True bifurcation lesions of the unprotected LMCA were classified into four types: Type A, with continuous involvement from the distal LMCA to the ostial LAD and the ostial LCX with eccentric plaques; Type B, with concentric plaques at the distal LMCA, eccentric plaques at the ostial LAD, and no plaques at the LCX; Type C, with continuous involvement from the distal LMCA to the ostial LCX, with eccentric plaques, and to the ostial LAD, with eccentric plaques; and Type D, with continuous involvement from the distal LMCA to the ostial LAD, with eccentric plaques, and to the ostial LCX, with concentric plaques. The carina was involved in only 3.5% of the plaques. A total of 51.7% of the plaques at the ostium of the LAD were soft, while 44.8% and 44.6% were fibrous in the distal LMCA and in the ostial LCX, respectively.

Conclusions:

We classified LMCA true bifurcation lesions into four types. The carina was always free from disease. Plaques at the ostial LAD tended to be soft, whereas those at the ostial LCX and the distal LMCA tended to be fibrous.

Numerical analysis of the effect of turbulence transition on the hemodynamic parameters in human coronary arteries

BACKGROUND

Local hemodynamics plays an important role in atherogenesis and the progression of coronary atherosclerosis disease (CAD). The primary biological effect due to blood turbulence is the change in wall shear stress (WSS) on the endothelial cell membrane, while the local oscillatory nature of the blood flow affects the physiological changes in the coronary artery. In coronary arteries, the blood flow Reynolds number ranges from few tens to several hundreds and hence it is generally assumed to be laminar while calculating the WSS calculations. However, the pulsatile blood flow through coronary arteries under stenotic condition could result in transition from laminar to turbulent flow condition.

METHODS

In the present work, the onset of turbulent transition during pulsatile flow through coronary arteries for varying degree of stenosis (i.e., 0%, 30%, 50% and 70%) is quantitatively analyzed by calculating the turbulent parameters distal to the stenosis. Also, the effect of turbulence transition on hemodynamic parameters such as WSS and oscillatory shear index (OSI) for varying degree of stenosis is quantified. The validated transitional shear stress transport (SST) k-ω model used in the present investigation is the best suited Reynolds averaged Navier-Stokes turbulence model to capture the turbulent transition. The arterial wall is assumed to be rigid and the dynamic curvature effect due to myocardial contraction on the blood flow has been neglected.

RESULTS

Our observations shows that for stenosis 50% and above, the WSSavg, WSSmax and OSI calculated using turbulence model deviates from laminar by more than 10% and the flow disturbances seems to significantly increase only after 70% stenosis. Our model shows reliability and completely validated.

CONCLUSIONS

Blood flow through stenosed coronary arteries seems to be turbulent in nature for area stenosis above 70% and the transition to turbulent flow begins from 50% stenosis.