Primary stenosis progression versus secondary stenosis formation in the left coronary bifurcation: A mechanical point of view

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%.

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.

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.

Long-term prognostic impact of paravalvular leakage on coronary artery disease requires patient-specific quantification of hemodynamics

Transcatheter aortic valve replacement (TAVR) is a frequently used minimally invasive intervention for patient with aortic stenosis across a broad risk spectrum. While coronary artery disease (CAD) is present in approximately half of TAVR candidates, correlation of post-TAVR complications such as paravalvular leakage (PVL) or misalignment with CAD are not fully understood. For this purpose, we developed a multiscale computational framework based on a patient-specific lumped-parameter algorithm and a 3-D strongly-coupled fluid-structure interaction model to quantify metrics of global circulatory function, metrics of global cardiac function and local cardiac fluid dynamics in 6 patients. Based on our findings, PVL limits the benefits of TAVR and restricts coronary perfusion due to the lack of sufficient coronary blood flow during diastole phase (e.g., maximum coronary flow rate reduced by 21.73%, 21.43% and 21.43% in the left anterior descending (LAD), left circumflex (LCX) and right coronary artery (RCA) respectively (N = 6)). Moreover, PVL may increase the LV load (e.g., LV load increased by 17.57% (N = 6)) and decrease the coronary wall shear stress (e.g., maximum wall shear stress reduced by 20.62%, 21.92%, 22.28% and 25.66% in the left main coronary artery (LMCA), left anterior descending (LAD), left circumflex (LCX) and right coronary artery (RCA) respectively (N = 6)), which could promote atherosclerosis development through loss of the physiological flow-oriented alignment of endothelial cells. This study demonstrated that a rigorously developed personalized image-based computational framework can provide vital insights into underlying mechanics of TAVR and CAD interactions and assist in treatment planning and patient risk stratification in patients.

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.

Hemodynamic Characteristics of Patients With Suspected Coronary Heart Disease at Their Initial Visit

Purpose

It is difficult for doctors to decide whether patients with suspected coronary heart disease classified as Coronary Artery Disease Reporting and Data System (CAD-RADS) < 3 should be administered preventive treatment, or whether non-atherosclerotic chest pain should be considered. The aim of the current study was to investigate coronary hemodynamic characteristics in such patients, which may provide more information on their stenosis and be helpful for initial diagnoses.

Methods

Two patient-specific models were reconstructed based on the coronary computed tomographic angiography underwent in 2012. Patient 1 was classified as CAD-RADS 0, and was readmitted to hospital due to coronary artery disease within 5 years. Patient 2 was classified as CAD-RADS 2, and has experienced no adverse events to date. Computational fluid dynamics (CFD) analysis was used to obtain hemodynamic parameters including flow rate waveform, flow streamlines, time-average wall shear stress (TAWSS), and oscillatory shear index (OSI).

Results

Patient 1 exhibited no physiological characteristics of right coronary artery flow waveform, large areas of low TAWSS, and slow blood flow in the proximal and middle segments of the left anterior descending branch. Patient 2 exhibited reduced coronary supply, small and separate areas of abnormal TAWSS, and a higher left anterior descending branch OSI than patient 1.

Conclusion

Hemodynamic abnormalities may play an important role in the prognosis of patients with coronary stenosis, and patient-specific hemodynamic characteristics may facilitate more accurate initial diagnosis, and better management. Overall hemodynamics (along the whole vessel) warranted attention at the time of the initial visit in patients classified as CAD-RADS < 3.

Correlation of Computational Instantaneous Wave-Free Ratio With Fractional Flow Reserve for Intermediate Multivessel Coronary Disease

This study computationally assesses the accuracy of an instantaneous wave-free ratio (iFR) threshold range compared to standard modalities such as fractional flow reserve (FFR) and coronary flow reserve (CFR) for multiple intermediate lesions near the left main (LM) coronary bifurcation. iFR is an adenosine-independent index encouraged for assessment of coronary artery disease (CAD), but different thresholds are debated. This becomes particularly challenging in cases of multivessel disease when sensitivity to downstream lesions is unclear. Idealized LM coronary arteries with 34 different intermediate stenoses were created and categorized (Medina) as single and multiple lesion groups. Computational fluid dynamics modeling was performed with physiologic boundary conditions using an open-source software (simvascular1) to solve the time-dependent Navier-Stokes equations. A strong linear relationship between iFR and FFR was observed among studied models, indicating computational iFR values of 0.92 and 0.93 are statistically equivalent to an FFR of 0.80 in single and multiple lesion groups, respectively. At the clinical FFR value (i.e., 0.8), a triple-lesion group had smaller CFR compared to the single and double lesion groups (e.g., triple = 3.077 versus single = 3.133 and double = 3.132). In general, the effect of additional intermediate downstream lesions (minimum lumen area > 3 mm2) was not statistically significant for iFR and CFR. A computational iFR of 0.92 best predicts an FFR of 0.80 and may be recommended as threshold criteria for computational assessment of LM stenosis following additional validation using patient-specific models.

A Hemodynamic Comparison of Myocardial Bridging and Coronary Atherosclerotic Stenosis: A Computational Model With Experimental Evaluation

Myocardial bridging (MB) and coronary atherosclerotic stenosis can impair coronary blood flow and may cause myocardial ischemia or even heart attack. It remains unclear how MB and stenosis are similar or different regarding their impacts on coronary hemodynamics. The purpose of this study was to compare the hemodynamic effects of coronary stenosis and MB using experimental and computational fluid dynamics (CFD) approaches. For CFD modeling, three MB patients with different levels of lumen obstruction, mild, moderate, and severe were selected. Patient-specific left anterior descending (LAD) coronary artery models were reconstructed from biplane angiograms. For each MB patient, the virtually healthy and stenotic models were also simulated for comparison. In addition, an in vitro flow-loop was developed, and the pressure drop was measured for comparison. The CFD simulations results demonstrated that the difference between MB and stenosis increased with increasing MB/stenosis severity and flowrate. Experimental results showed that increasing the MB length (by 140%) only had significant impact on the pressure drop in the severe MB (39% increase at the exercise), but increasing the stenosis length dramatically increased the pressure drop in both moderate and severe stenoses at all flow rates (31% and 93% increase at the exercise, respectively). Both CFD and experimental results confirmed that the MB had a higher maximum and a lower mean pressure drop in comparison with the stenosis, regardless of the degree of lumen obstruction. A better understanding of MB and atherosclerotic stenosis may improve the therapeutic strategies in coronary disease patients and prevent acute coronary syndromes.

Computational fluid dynamics simulation of time-resolved blood flow in Budd-Chiari syndrome with inferior vena cava stenosis and its implication for postoperative efficacy assessment

BACKGROUND

This study aimed to adopt computational fluid dynamics to simulate the blood flow dynamics in inferior vena cava stenosis based on time-dependent patient-specific models of Budd-Chiari syndrome as well as a normal model. It could offer valuable references for a retrospective insight into the underlying mechanisms of Budd-Chiari syndrome pathogenesis as well as more accurate evaluation of postoperative efficacy.

METHODS

Three-dimensional inferior vena cava models of Budd-Chiari syndrome patient-specific (preoperative and postoperative) and normal morphology model were reconstructed as per magnetic resonance images using Simpleware. Moreover, computational fluid dynamics of time-resolved inferior vena cava blood flow were simulated using actual patient-specific measurements to reflect time-dependent flow rates.

FINDINGS

The assessment of the preoperative model revealed the dramatic variations of hemodynamic parameters of the stenotic inferior vena cava. Moreover, the comparison of the preoperative and postoperative models with the normal model as benchmark showed that postoperative hemodynamic parameters were markedly ameliorated via stenting, with the attenuation of overall velocity and wall shear stress, and the increase of pressure. However, the comparative analysis of the patient-specific simulations revealed that some postoperative hemodynamic profiles still bore some resemblance to the preoperative ones, indicating potential risks of restenosis.

INTERPRETATION

Computational fluid dynamics simulation of time-resolved blood flow could reveal the tight correlation between the hemodynamic characteristics and the pathological mechanisms of inferior vena cava stenosis. Furthermore, such time-resolved hemodynamic profiles could provide a quantitative approach to diagnosis, operative regimen and postoperative evaluation of Budd-Chiari syndrome with inferior vena cava stenosis.

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.

Numerical computation of blood hemodynamic through constricted human left coronary artery: Pulsatile simulations

BACKGROUND AND OBJECTIVE

The accumulation of plaque in the coronary artery of the human heart restricts the path of blood flow in that region and leads to Coronary Artery Disease. This study's goal is to present the pulsatile blood flow conduct through four different levels of constrictions, i.e., healthy, 25%, 50%, and 75% in human left coronary arteries.

METHODS

Using CT scan data of a healthy person, the two-dimensional coronary model is constructed. A non-Newtonian Carreau model is used to study the maximum flow velocity, streamline effect, and maximum Wall Shear Stress at the respective constricted areas over the entire cardiac cycle. Finite Volume Method is executed for solving the governing equations. The fluctuating Wall Shear Stress (WSS) at different levels was assessed using Computational Fluid Dynamics (CFD).

RESULTS

The comparative study of the diseased arteries showcases that at the systolic phase, the 75% blocked artery attains the maximum velocity of 0.14 m/s and 0.53 m/s at t=0.005 s and t=0.115 s, respectively. While the maximum velocity takes a significant drop at t=0.23 s and t=0.345 s, this marks the diastolic phase. The streamline contour showcased the blood flow conduct at different phases of the cardiac cycle. At the peak systolic phase, a dense flow separation was observed near the blocked regions. It highlights the disturbed flow in that particular region. The most severely diseased artery acquires the maximum WSS of 18.81 Pa at the peak systolic phase, i.e., at t=0.115 s.

CONCLUSIONS

The computational study of the hemodynamic parameters can aid in the early anticipation of the degree of the severity of the diseased arteries. This study, in a way, could benefit doctors/surgeons to plan an early treatment/surgery on the grounds of the severity of the disease. Thus, a before time prognosis could restrain the number of deaths caused due to Coronary Artery Disease.

The association between genetic variants in the genes for cytochrome P450 B1 and ATP-binding cassette transporter genes and breast cancer risk

Breast cancer is among the most common malignancies in women. Recent studies have shown that polymorphisms in genes involved in the metabolism and transport of anticancer drugs are associated with outcomes of several malignancies, e.g., breast cancer. In this study we evaluate whether CYP1B1/rs1056836 and ABCB1/rs2032582 gene variants are associated with breast cancer. Eighty eight cases and 200 controls, were genotyped for polymorphisms of the CYP1B1 and ABCB1 genes using Taqman®-based methods. Logistic regression was also used to test the associations between breast cancer risk and the various genotypes involved. The GG genotype of rs2032582 locus had a frequency of 43.5% with 0.38 MAF; while the GT and TT genotypes in the control group were 40% and 16.5%, respectively. The GG, GT and TT genotype frequencies in the patients with breast cancer were 45.5%, 12.5% and 26.1%, respectively. An association was observed between the TT genotype of ABCB1/rs2032582 locus and a larger breast cancer tumor size (P < 0.05). However, neither the relationship between the CYP1B1 polymorphism and breast cancer type nor the risk of breast cancer were statistically significant. Our data suggest a potential association of the ABCB1 genetic variant with breast cancer tumor size, however further investigation in a larger population is necessary to show its value as a risk stratification biomarker.

Numerical study of hemodynamic and diagnostic parameters affected by stenosis in bifurcated artery

Stenosis in the bifurcated coronary artery has attracted wide attention among the researchers. Many investigations have been carried out by means of Computational Fluid Dynamics (CFD) to better understand the physical mechanism inside the stenotic bifurcated artery. However, the main focus of the existing publications is limited to the variation of hemodynamic parameters affected by the stenosis and bifurcation structure. The present study aims to make further evaluations of stenosis development and diagnostic lesion assessments based on the critical values of hemodynamic and practical diagnostic parameters. The bifurcated coronary artery with initial stenosis source has been studied in 2 D unsteady model. Different locations of initial stenosis source have been found to greatly affect the orientation of the further stenosis development. In addition, different stenosis severities (diameter stenosis: 30%, 50% and 70%), different stenosis lengths (3 mm, 6 mm and 9 mm) and different pulse rates (75 bpm, 100 bpm and 120 bpm) as controlling parameters have been investigated.

Numerical study to identify the effect of fluid presence on the mechanical behavior of the stents during coronary stent expansion

In this study, structural analysis and one-way fluid-structure interaction (FSI) analysis were performed to identify the effect of fluid presence on the mechanical behavior of the stents during stent expansion. An idealized vessel model with stenosis was used for simulation, and stents made of metal and polymer were assumed, respectively. The bilinear model was applied to the stents, and the Mooney-Rivlin model was applied to the arterial wall and plaque. The blood used in the FSI analysis was assumed to be a non-Newtonian fluid. As a result of all numerical simulations, the von Mises stress, the first principal stress and the displacement were calculated as the mechanical behaviors. Through the comparison of the results of the structural analysis with those of the one-way FSI analysis, our results indicated the fluid had no significant influence on the expansion of the metal stent. However, it was found that the expansion of the polymer stent affected by the presence of fluid. These findings meant the one-way FSI technique was suggested to achieve an accurate analysis when targeting a polymer stent for numerical simulation.