Study of Coronary Atherosclerosis Using Blood Residence Time

Computational Modeling Approach to Profile Hemodynamical Behavior in a Healthy Aorta

Cardiovascular diseases (CVD) remain the leading cause of mortality among older adults. Early detection is critical as the prognosis for advanced-stage CVD is often poor. Consequently, non-invasive diagnostic tools that can assess hemodynamic function, particularly of the aorta, are essential. Computational fluid dynamics (CFD) has emerged as a promising method for simulating cardiovascular dynamics efficiently and cost-effectively, using increasingly accessible computational resources. This study developed a CFD model to assess the aorta geometry using tetrahedral and polyhedral meshes. A healthy aorta was modeled with mesh sizes ranging from 0.2 to 1 mm. Key hemodynamic parameters, including blood pressure waveform, pressure difference, wall shear stress (WSS), and associated wall parameters like relative residence time (RRT), oscillatory shear index (OSI), and endothelial cell activation potential (ECAP) were evaluated. The performance of the CFD simulations, focusing on accuracy and processing time, was assessed to determine clinical viability. The CFD model demonstrated clinically acceptable results, achieving over 95% accuracy while reducing simulation time by up to 54%. The entire simulation process, from image construction to the post-processing of results, was completed in under 120 min. Both mesh types (tetrahedral and polyhedral) provided reliable outputs for hemodynamic analysis. This study provides a novel demonstration of the impact of mesh type in obtaining accurate hemodynamic data, quickly and efficiently, using CFD simulations for non-invasive aortic assessments. The method is particularly beneficial for routine check-ups, offering improved diagnostics for populations with limited healthcare access or higher cardiovascular disease risk.

Investigation of the size and shape of nano-microcarriers for targeted drug delivery to atherosclerotic plaque in ischemic stroke prevention

Recognizing the significance of drug carriers in the treatment of atherosclerotic plaque is crucial in light of the worldwide repercussions of ischemic stroke. Conservative methodologies, specifically targeted drug delivery, present encouraging substitutes that mitigate the hazards linked to invasive procedures. With the intention of illuminating their considerable significance and prospective benefits, this study examines the impact of the geometry and dimensions of drug-loaded nano-microcarriers on atherosclerotic plaque. The research utilizes a finite element approach to simulate the motion and fluid dynamics of nano-microcarriers loaded with drugs within the carotid arteries. Carriers are available in a variety of shapes and sizes to accommodate patient-specific geometries, pulsatile fluid flow, and non-Newtonian blood properties. Optimization of drug delivery is achieved through the examination of carrier interaction with the inner wall. The results demonstrated that the interaction data between particles and the inner wall of atherosclerotic plaques exhibits micro- and nanoscale patterns that are distinct. Symmetric plaques demonstrate that nanoparticles with a 0.4 shape factor and diameters below 200 nm show the highest interaction rate. Conversely, larger particles (200 and 500 nm) with shape factors of 1 demonstrate comparatively elevated interaction rates. The optimal shape factor for drug-loaded microparticles has been determined to be one, and the number of interactions increases as the diameter of the nanoparticles increases, with a significant increase observed at a shape factor of one. Asymmetric plaques exhibit the maximum interaction rates among particles that have a shape factor of 0.4 and have diameters smaller than 500 µm. The findings establish a foundation for novel therapeutic strategies, establishing nano-microparticles as auspicious contenders for accurate and efficacious drug delivery systems that inhibit plaque proliferation.

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

Effect of PCI on ophthalmic artery hemodynamics in patients with acute coronary syndrome

Purpose

We aimed to explore the effects of percutaneous coronary intervention (PCI) on the ophthalmic artery (OA) hemodynamics in patients with acute coronary syndrome (ACS).

Methods

A total of 73 participants (Group0: healthy controls, Group1: Patients with ACS underwent PCI < 3 months, Group2: Patients with ACS underwent PCI ≥ 3 months) were enrolled. Computed tomographic angiography images were used to construct three-dimensional models of participants' OAs. Numerical simulations based on computational fluid dynamics were used to acquire hemodynamic parameters.

Results

The angle between the OA and internal carotid artery in Group2 was significantly larger compared with Group0 and Group1 (P = 0.003 and P = 0.044). Hemodynamic simulation showed a significantly slower OA blood velocity in Group1 than in the control (P < 0.001) and Group2 (P = 0.033). Lower wall shear stress was found in Group1 than that in control (P = 0.040). Patients after PCI had a higher wall pressure than healthy controls (P = 0.012 and P = 0.004). Mass flow ratios were decreased in Group1 and Group2 (P = 0.021 and P = 0.002). The hemodynamic parameters of OA were correlated with several clinical indicators.

Conclusions

The OA blood flow velocity of patients with ACS after PCI initially slowed down, which increased the risk of plaque formation, and then showed an increasing trend. There was a correlation between OA hemodynamic parameters and clinical indexes related to cardiac stress. Ischemia-reperfusion injury and changes in blood flow status after PCI may affect OA morphology and hemodynamics, leading to ocular lesions.

Trial registration

ChiCTR2100050428.

Longitudinal wall shear stress evaluation using centerline projection approach in the numerical simulations of the patient-based carotid artery

In this numerical study, areas of the carotid bifurcation and of a distal stenosis in the internal carotid artery are closely observed to evaluate the patient's current risks of ischemic stroke. An indicator for the vessel wall defects is the stress exerted by blood on the vessel tissue, typically expressed by the amplitude of the wall shear stress vector (WSS) and its oscillatory shear index. To detect negative shear stresses corresponding with reversal flow, we perform orientation-based shear evaluation. We investigate the longitudinal component of the wall shear vector, where tangential vectors aligned longitudinally with the vessel are necessary. However, resulting from imaging segmentation resolution of patients' computed tomography angiography scans and stenotic regions, the geometry model's mesh is non-smooth on its surface areas and the automatically generated tangential vector field is discontinuous and multi-directional, making an interpretation of our orientation-based risk indicators unreliable. We improve the evaluation of longitudinal shear stress by applying the projection of the vessel's centerline to the surface to construct smooth tangential field aligned longitudinally with the vessel. We validate our approach for the longitudinal WSS component and the corresponding oscillatory index by comparing them to results obtained using automatically generated tangents in both rigid and elastic vessel modeling and to amplitude-based indicators. We present the major benefit of our longitudinal WSS evaluation based on its directionality for the cardiovascular risk assessment, which is the detection of negative WSS indicating persistent reversal or transverse flow. This is impossible in the case of the amplitude-based WSS.

Biomimicking Atherosclerotic Vessels: A Relevant and (Yet) Sub-Explored Topic

Atherosclerosis represents the etiologic source of several cardiovascular events, including myocardial infarction, cerebrovascular accidents, and peripheral artery disease, which remain the leading cause of mortality in the world. Numerous strategies are being delineated to revert the non-optimal projections of the World Health Organization, by both designing new diagnostic and therapeutic approaches or improving the interventional procedures performed by physicians. Deeply understanding the pathological process of atherosclerosis is, therefore, mandatory to accomplish improved results in these trials. Due to their availability, reproducibility, low expensiveness, and rapid production, biomimicking physical models are preferred over animal experimentation because they can overcome some limitations, mainly related to replicability and ethical issues. Their capability to represent any atherosclerotic stage and/or plaque type makes them valuable tools to investigate hemodynamical, pharmacodynamical, and biomechanical behaviors, as well as to optimize imaging systems and, thus, obtain meaningful prospects to improve the efficacy and effectiveness of treatment on a patient-specific basis. However, the broadness of possible applications in which these biomodels can be used is associated with a wide range of tissue-mimicking materials that are selected depending on the final purpose of the model and, consequently, prioritizing some materials' properties over others. This review aims to summarize the progress in fabricating biomimicking atherosclerotic models, mainly focusing on using materials according to the intended application.

Role of secondary flows in coronary artery bifurcations before and after stenting: What is known so far?

Coronary arteries are uniformly exposed to traditional cardiovascular risk factors. However, atherosclerotic lesions occur in preferential regions of the coronary tree, especially in areas with disturbed local blood flow, such as coronary bifurcations. Over the latest years, secondary flows have been linked to the inception and progression of atherosclerosis. Most of these novel findings have been obtained in the field of computational fluid dynamic (CFD) analysis and biomechanics but remain poorly understood by cardiovascular interventionalists, despite the important impact that they may have in clinical practice. We aimed to summarize the current available data regarding the pathophysiological role of secondary flows in coronary artery bifurcation, providing an interpretation of these findings from an interventional perspective.

Numerical computation of pulsatile hemodynamics and diagnostic concern of coronary bifurcated artery flow for Newtonian and non-Newtonian fluid

Atherosclerotic with the high occurrence of plaque formation due to stenosis has attracted wide attention among researchers. The left coronary artery has been studied in two-dimensional and in three-dimensional (3D) bifurcation as the models of blood flow through Newtonian and non-Newtonian fluids to better understand the physical mechanism. The computational Fluid Dynamics (CFD) technique is incorporated in COMSOL Multiphysics and then it is justified by satisfactory validation. It is found that the Newtonian model shows larger recirculation zones than non-Newtonian does. The present study also focuses on the evaluations of the lesion of diagnostic and the coefficient of pressure drop assessments on the basis of the diagnostic parameter's critical values affected by the rheological model. Nevertheless, the leading concentration of the subsisting investigation works is confined to the change of importance factor (IFc) affected by arterial blockage. But the IFc of non-Newtonian fluid for 3D left coronary artery bifurcation model decreases with increasing bifurcation angle and the time-averaged inlet pressure is the least for smaller bifurcation angles. The current research further concentrates that the flow separation length reduces with developing bifurcation angle in bifurcated geometry. It is significant to mention that non-Newtonian blood flow model incorporating hemodynamic and diagnostic parameters has great impacts on instantaneous flow systems.

Conduit design with expanding diameter for enhanced flow

Conduits are commonly used for treating lesions in arteries and veins. The conventional stents are cylindrical in shape, which increases flow resistance with length. This study presents a design of stents and conduits where the conduit caliber expands gradually to reduce resistance while avoiding flow separation. Inflow was provided from a header tank at two different pressures (i.e., 10 and 25 mm Hg pressure) into a cylindrical or expanding conduit. The initial conduit calibers were 2-, 3-, 4-, and 5-mm and 160-, 310-, and 620-mm lengths in each case. The flow rates of expanding caliber conduits (at a rate of r4-6/cm where r is the initial conduit radius) were compared to traditional cylindrical conduits of constant radius. The expanded caliber yields a significantly increased flow of 16-55% for R4/L expansion, 9-44% for R5/L expansion, and 1-28% for R6/L expansion. Simulated flow models using computational fluid dynamics (CFD) were used to validate and expand the experimental findings. Flow separation was detected for certain simulations by flow pathlines and wall shear stress (WSS) calculations. The results showed that a caliber expansion rate of r6/cm is the optimal rate of expansion for most potential applications with minimum flow separation, lower resistance, and increased flow.

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

BACKGROUND

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

METHODS

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

FINDINGS

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

INTERPRETATION

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

Non-invasive diagnostics of blockage growth in the descending aorta-computational approach

Atherosclerosis causes blockages to the main arteries such as the aorta preventing blood flow from delivering oxygen to the organs. Non-invasive diagnosis of these blockages is difficult, particularly in primary healthcare. In this paper, the effect of arterial blockage development and growth is investigated at the descending aorta on some possible non-invasive assessment parameters including the blood pressure waveform, wall shear stress (WSS), time-average WSS (TAWSS) and the oscillation shear index (OSI). Blockage severity growth is introduced in a simulation model as 25%, 35%, 50% and 65% stenosis at the descending aorta based on specific healthy control aorta data clinically obtained. A 3D aorta model with invasive pulsatile waveforms (blood flow and pressure) is used in the CFD simulation. Blockage severity is assessed by using blood pressure measurements at the left subclavian artery. An arterial blockage growth more than 35% of the lumen diameter significantly affects the pressure. A strong correlation is also observed between the ascending aorta pressure values, pressure at the left subclavian artery and the relative residence time (RRT). An increase of RRT downstream from the stenosis indicates a 35% stenosis at the descending aorta which results in high systolic and diastolic pressure readings. The findings of this study could be further extended by transferring the waveform reading from the left subclavian artery to the brachial artery.

Graphical abstract

Structure (Epicardial Stenosis) and Function (Microvascular Dysfunction) That Influence Coronary Fractional Flow Reserve Estimation

Background. The treatment of coronary stenosis is decided by performing high risk invasive surgery to generate the fractional flow reserve diagnostics index, a ratio of distal to proximal pressures in respect of coronary atherosclerotic plaques. Non-invasive methods are a need of the times that necessitate the use of mathematical models of coronary hemodynamic physiology. This study proposes an extensible mathematical description of the coronary vasculature that provides an estimate of coronary fractional flow reserve. Methods. By adapting an existing computational model of human coronary blood flow, the effects of large vessel stenosis and microvascular disease on fractional flow reserve were quantified. Several simulations generated flow and pressure information, which was used to compute fractional flow reserve under several conditions including focal stenosis, diffuse stenosis, and microvascular disease. Sensitivity analysis was used to uncover the influence of model parameters on fractional flow reserve. The model was simulated as coupled non-linear ordinary differential equations and numerically solved using our implicit higher order method. Results. Large vessel stenosis affected fractional flow reserve. The model predicts that the presence, rather than severity, of microvascular disease affects coronary flow deleteriously. Conclusions. The model provides a computationally inexpensive instrument for future in silico coronary blood flow investigations as well as clinical-imaging decision making. A combination of focal and diffuse stenosis appears to be essential to limit coronary flow. In addition to pressure measurements in the large epicardial vessels, diagnosis of microvascular disease is essential. The independence of the index with respect to heart rate suggests that computationally inexpensive steady state simulations may provide sufficient information to reliably compute the index.

Numerical simulation in the abdominal aorta and the visceral arteries with or without stenosis based on 2D PCMRI

It is important to obtain accurate boundary conditions (BCs) in hemodynamic simulations. This article aimed to improve the accuracy of BCs in computational fluid dynamics (CFD) simulation and analyze the differences in hemodynamics between healthy volunteers and patients with visceral arterial stenosis (VAS). The geometric models of seven cases were reconstructed using the magnetic resonance angiogram (MRA) or computed tomography angiogram (CTA) imaging data. The physiological flow waveforms obtained from 2D Phase Contrast Magnetic Resonance Imaging (PCMRI) were imposed on the aortic inlet and the visceral arteries' outlets. The individualized RCR values of the three-element Windkessel model were imposed on the aortic outlet. CFD simulations were run in the open-source software: svSolver. Two specific time points were selected to compare the hemodynamics of healthy volunteers and patients with VAS. The results suggested that blood in the stenotic visceral arteries flowed at high speed throughout the cardiac cycle. The low pressure is distributed at stenotic lesions. The wall shear stress (WSS) reached 4 Pa near stenotic locations. The low time average wall shear stress (TAWSS), high oscillatory shear index (OSI), and high relative residence time (RRT) concentrated in the abdominal aorta. Besides, the ratios of the areas with low TAWSS, high OSI, and high RRT to the computational domain were higher in patients with VAS than which in the healthy volunteers. The individualized BCs were used for hemodynamic simulations and results suggest that patients with stenosis have a higher risk of blood retention and atherosclerosis formation in the abdominal aorta.