A comprehensive fluid-structure interaction model of the left coronary artery

Non-invasive fractional flow reserve estimation in coronary arteries using angiographic images

Coronary artery disease is the leading global cause of mortality and Fractional Flow Reserve (FFR) is widely regarded as the gold standard for assessing coronary artery stenosis severity. However, due to the limitations of invasive FFR measurements, there is a pressing need for a highly accurate virtual FFR calculation framework. Additionally, it's essential to consider local haemodynamic factors such as time-averaged wall shear stress (TAWSS), which play a critical role in advancement of atherosclerosis. This study introduces an innovative FFR computation method that involves creating five patient-specific geometries from two-dimensional coronary angiography images and conducting numerical simulations using computational fluid dynamics with a three-element Windkessel model boundary condition at the outlet to predict haemodynamic distribution. Furthermore, four distinct boundary condition methodologies are applied to each geometry for comprehensive analysis. Several haemodynamic features, including velocity, pressure, TAWSS, and oscillatory shear index are investigated and compared for each case. Results show that models with average boundary conditions can predict FFR values accurately and observed errors between invasive FFR and virtual FFR are found to be less than 5%.

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.

A hemodynamic model of artery bypass graft considering microcirculation function

BACKGROUND

The incidence of arterial stenosis is increasing year by year. In order to better diagnose and treat arterial stenosis, numerical simulation technology has become a popular method.

OBJECTIVE

A novel model is constructed to investigate the influence of microcirculation on the hemodynamics of artery bypass graft.

METHODS

In this paper, a severely narrow artery bypass graft model is considered. The geometric shape includes a narrow artery tube and a bypass graft of the same diameter with a 45° suture angle. The fluid-structure interaction model is considered by finite element numerical calculation, and the flow is simulated with microcirculation as the outlet boundary condition. The changes of blood flow velocity, pressure and wall shear stress are analyzed.

RESULTS

The results show that blood almost entirely flows into the graft tube and there is no recirculation area at the anastomosis.

CONCLUSION

The artery bypass graft model considering microcirculation function could simulate the physiological characteristics of blood flow more reasonably, and it provide helps for clinicians to diagnose and treat arterial stenosis.

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.

Numerical analysis of the hemodynamics of rat aorta based on magnetic resonance imaging and fluid-structure interaction

Murine models have been widely used to investigate the mechanobiology of aortic atherosclerosis and dissections, which develop preferably at different anatomic locations of aorta. Based MRI and finite element analysis with fluid-structure interaction, we numerically investigated factors that may affect the blood flow and structural mechanics of rat aorta. The results indicated that aortic root motion greatly increases time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), displacement of the aorta, and enhances helical flow pattern but has limited influence on effective stress, which is highly modulated by blood pressure. Moreover, the influence of the motion component on these indicators is different with axial motion more obvious than planar motion. Surrounding fixation of the intercostal arteries and the branch vessels on aortic arch would reduce the influence of aortic root motion. The compliance of the aorta has different influences at different regions, leading to decrease in TAWSS and helical flow, increase in OSI, RRT at the aortic arch, but has reversed effects on the branch vessels. When compared with the steady flow, the pulsatile blood flow would obviously increase the WSS, the displacement, and the effective stress in most regions. In conclusion, to accurately quantify the blood flow and structural mechanics of rat aorta, the motion of the aortic root, the compliance of aortic wall, and the pulsation of blood flow should be considered. However, when only focusing on the effective stress in rat aorta, the motion of the aortic root may be neglected.

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.

Platelet adhesion potential estimation in a normal and diseased coronary artery model: effects of shear stress magnitude versus shear stress history

Platelets play a salient role in the pathogenesis of coronary diseases; primarily through their adhesion to other platelets, endothelial cells and plasma proteins. It is necessary for platelets to activate in order for them to adhere to these different substrates. One of the key regulatory mechanical factors in platelet activation is shear stress, which has been shown to alter multiple platelet functions through the activation of mechanoreceptors. Our goal was to investigate how different numerical shear stress tracking techniques affect platelet adhesion estimates within physiologically relevant computational models. Previously, we developed a physiological coronary artery computational fluid dynamics model. Shear stress waveforms, obtained from these models, were used to monitor in vitro platelet and endothelial cell adhesion marker expression. In this work, the adhesion marker expression data was regressed to obtain numerical functions for receptor expression predictions. These functions were input into a customized adhesion model utilizing different shear stress tracking techniques. For the normal vascular conditions and minimal pathological disease models, shear stress tracking did not significantly affect the adhesion estimates. However, for the severe pathological model, the two shear stress tracking methods had vastly different estimates. Therefore, shear stress tracking methods must be chosen accurately to predict platelet adhesion potentials for accurate modeling techniques.

Transient numerical simulation of the right coronary artery originating from the left sinus and the effect of its acute take-off angle on hemodynamics

Background

An anomalous origin of the right coronary artery from the left coronary artery sinus is usually characterized by an acute take-off angle. Most affected patients have no clinical symptoms; however, some patients have decreased blood flow into the right coronary artery during exercise, which can lead to symptoms such as myocardial ischemia. Most researchers who have studied an anomalous origin of the right coronary artery from the left coronary artery sinus have done so through clinical cases. In this study, we used numerical simulation to evaluate the hemodynamics of this condition and the effect of an acute take-off angle on hemodynamic parameters. We expect that the results of this study will help in further understanding the clinical symptoms of this anomaly and the hemodynamic impact of an acute take-off angle.

Methods

Three-dimensional models were reconstructed based on the computed tomography images from 16 patients with a normal right coronary artery and 26 patients with an anomalous origin of the right coronary artery from the left coronary artery sinus. A numerical simulation of a two-way fluid-structure interaction was executed with ANSYS Workbench software. The blood was assumed to be an incompressible Newtonian fluid, and the vessel was assumed to be an isotropic, linear elastic material. Hemodynamic parameters and the effect of an acute take-off angle were statistically analyzed.

Results

During the systolic period, the wall pressure in the right coronary artery was significantly reduced in patients with an anomalous origin of the right coronary artery (t =1.32 s, P=0.0001; t =1.34-1.46 s, P<0.0001). The wall shear stress in the abnormal group was higher at the beginning of the systolic period (t =1.24 s, P=0.0473; t =1.26 s, P=0.0193; t =1.28 s, P=0.0441). The acute take-off angle was smaller in patients with clinical symptoms (27.81°±4.406°) than in patients without clinical symptoms (31.86°±2.789°; P=0.017). In the symptomatic group, pressure was negatively correlated with the acute take-off angle (P=0.0185-0.0341, r=-0.459 to -0.4167).

Conclusions

This study shows that an anomalous origin of the right coronary artery from the left coronary artery sinus causes changes in hemodynamic parameters, and that an acute take-off angle in patients with this anomaly is associated with terminal ischemia of the right coronary artery.

Effect of Wall Elasticity on Hemodynamics and Wall Shear Stress in Patient-Specific Simulations in the Coronary Arteries

Wall shear stress (WSS) has been shown to be associated with myocardial infarction (MI) and progression of atherosclerosis. Wall elasticity is an important feature of hemodynamic modeling affecting WSS calculations. The objective of this study was to investigate the role of wall elasticity on WSS, and justify use of either rigid or elastic models in future studies. Digital anatomic models of the aorta and coronaries were created based on coronary computed tomography angiography (CCTA) in four patients. Hemodynamics was computed in rigid and elastic models using a finite element flow solver. WSS in five timepoints in the cardiac cycle and time averaged wall shear stress (TAWSS) were compared between the models at each 3 mm subsegment and 4 arcs in cross sections along the centerlines of coronaries. In the left main (LM), proximal left anterior descending (LAD), left circumflex (LCX), and proximal right coronary artery (RCA) of the elastic model, the mean percent radial increase 5.95 ± 1.25, 4.02 ± 0.97, 4.08 ± 0.94, and 4.84 ± 1.05%, respectively. WSS at each timepoint in the cardiac cycle had slightly different values; however, when averaged over the cardiac cycle, there were negligible differences between the models. In both the subsegments (n = 704) and subarc analysis, TAWSS in the two models were highly correlated (r = 0.99). In investigation on the effect of coronary wall elasticity on WSS in CCTA-based models, the results of this study show no significant differences in TAWSS justifying using rigid wall models for future larger studies.

Endothelial Cell Biomechanical Responses are Dependent on Both Fluid Shear Stress and Tensile Strain

Introduction

The goal of this study was to investigate how concurrent shear stress and tensile strain affect endothelial cell biomechanical responses.

Methods

Human coronary artery endothelial cells were exposed to concurrent pulsatile shear stress and cyclic tensile strain in a programmable shearing and stretching device. Three shear stress-tensile strain conditions were used: (1) pulsatile shear stress at 1 Pa and cyclic tensile strain at 7%, simulating normal stress/strain conditions in a healthy coronary artery; (2) shear stress at 3.7 Pa and tensile strain at 3%, simulating pathological stress/strain conditions near a stenosis; (3) shear stress at 0.7 Pa and tensile strain at 5%, simulating pathological stress/strain conditions in a recirculation zone. Cell morphology was quantified using immunofluorescence microscopy. Cell surface PECAM-1 phosphorylation, ICAM-1 expression, ERK1/2 and NF-κB activation were measured using ELISA or Western blot.

Results

Simultaneous stimulation from pulsatile shear stress and cyclic tensile strain induced a significant increase in cell area, compared to that induced by shear stress or tensile strain alone. The combined stimulation caused significant increases in PECAM-1 phosphorylation. The combined stimulation also significantly enhanced EC surface ICAM-1 expression (compared to that under shear stress alone) and transcriptional factor NF-κB activation (compared to that under control conditions).

Conclusion

Pulsatile shear stress and cyclic tensile strain could induce increased but not synergistic effect on endothelial cell morphology or activation. The combined mechanical stimulation can be relayed from cell membrane to nucleus. Therefore, to better understand how mechanical conditions affect endothelial cell mechanotransduction and cardiovascular disease development, both shear stress and tensile strain need to be considered.

Three-dimensional evaluation of the spatial morphology of stented coronary artery segments in relation to restenosis

To investigate the correlations between the three-dimensional (3D) parameters of target coronary artery segments and restenosis after stent implantation. Sixty-four patients after single, cobalt chromium platform stent (27 BM stents and 37 DES) implantation were investigated retrospectively 12 ± 6 months after the index procedure. 3D coronary artery reconstruction was performed before and after the stent implantation using appropriate projections by a dedicated reconstruction software. The curve of the target segment was characterized by the ratio of the vessel length measured at midline (arc: A) and the distance between the edge points of the stent (chord: C): A/C ratio (ACr). Age, diabetes and hyperlipidaemia were taken into account for the statistical evaluation. 22 patients were diagnosed with ISR, while 42 patients without any restenosis served as controls. The two groups did not differ regarding major cardiovascular risk factors, proportion of the treated vessels or the type of stents. Higher initial ACr values were associated with greater straightening of the vessel curvature in all groups (p < 0.001). Significant negative correlations were found in cases of proximal or distal edge bending angles (p < 0.001). Pre-stent edge bending angles < 7° often showed an increase after the stent implantation, while in case of higher initial values, the bending angles generally decreased. Using multivariate logistic regression modelling we found that the pre-stent ACr was an independent predictor of in-stent restenosis (odds ratio for 1% increase of the ACr: 1.08; p = 0.012). Changes of angles at the stent edges following stent implantation correlate with the initial local bending angles. The ACr predispose to chronic shear stress in the vessel wall, which may contribute to the pathological intimal proliferation.