Functional assessment of the stenotic carotid artery by CFD-based pressure gradient evaluation

Hemodynamic effects of bifurcation and stenosis geometry on carotid arteries with different degrees of stenosis

Objective.Carotid artery stenosis (CAS) is a key factor in pathological conditions, such as thrombosis, which is closely linked to hemodynamic parameters. Existing research often focuses on analyzing the influence of geometric characteristics at the stenosis site, making it difficult to predict the effects of overall vascular geometry on hemodynamic parameters. The objective of this study is to comprehensively examine the influence of geometric morphology at different degrees of CAS and at bifurcation sites on hemodynamic parameters.Approach.A three-dimensional model is established using computed tomography angiography images, and eight geometric parameters of each patient are measured by MIMICS. Then, computational fluid dynamics is utilized to investigate 60 patients with varying degrees of stenosis (10%-95%). Time and grid tests are conducted to optimize settings, and results are validated through comparison with reference calculations. Subsequently, correlation analysis using SPSS is performed to examine the relationship between the eight geometric parameters and four hemodynamic parameters. In MATLAB, prediction models for the four hemodynamic parameters are developed using back propagation neural networks (BPNN) and multiple linear regression.Main results.The BPNN model significantly outperforms the multiple linear regression model, reducing mean absolute error, mean squared error, and root mean squared error by 91.7%, 93.9%, and 75.5%, respectively, and increasingR2from 19.0% to 88.0%. This greatly improves fitting accuracy and reduces errors. This study elucidates the correlation and patterns of geometric parameters of vascular stenosis and bifurcation in evaluating hemodynamic parameters of CAS.Significance.This study opens up new avenues for improving the diagnosis, treatment, and clinical management strategies of CAS.

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

Turbulent blood dynamics in the left heart in the presence of mitral regurgitation: a computational study based on multi-series cine-MRI

In this work, we performed a computational image-based study of blood dynamics in the whole left heart, both in a healthy subject and in a patient with mitral valve regurgitation. We elaborated multi-series cine-MRI with the aim of reconstructing the geometry and the corresponding motion of left ventricle, left atrium, mitral and aortic valves, and aortic root of the subjects. This allowed us to prescribe such motion to computational blood dynamics simulations where, for the first time, the whole left heart motion of the subject is considered, allowing us to obtain reliable subject-specific information. The final aim is to investigate and compare between the subjects the occurrence of turbulence and the risk of hemolysis and of thrombi formation. In particular, we modeled blood with the Navier-Stokes equations in the arbitrary Lagrangian-Eulerian framework, with a large eddy simulation model to describe the transition to turbulence and a resistive method to manage the valve dynamics, and we used a finite element discretization implemented in an in-house code for the numerical solution.

The relationship between geometry and hemodynamics of the stenotic carotid artery based on computational fluid dynamics

OBJECTIVE

The purpose of this work was to investigate the relationship between the geometric factors and the hemodynamics of the stenotic carotid artery.

METHODS

We retrospectively reviewed data of patients with carotid stenosis (40%-95%). The Navier-Stokes equations were solved using ANSYS CFX 18.0. Correlation analysis was based on Spearman's test. Geometric variables (p < 0.1 in the univariate analysis) were entered into the logistical regression. A receiver-operating characteristics analysis was used to detect hemodynamically significant lesions.

RESULTS

81 patients (96 arteries) were evaluated. The logistic regression analysis revealed that the translesional pressure ratio was significantly correlated with the stenosis degree (OR = 1.147, p < 0.001) and the angle between internal carotid artery and external carotid artery (angle γ) (OR = 0.933, p = 0.01). The translesional wall shear stress ratio was significantly correlated with stenosis degree (OR = 1.094, p < 0.001), lesion length (OR = 0.873, p = 0.01), lumen area of internal carotid artery (OR = 0.867, p = 0.002), and lumen area of common carotid artery (OR = 1.058, p = 0.01). For predicting low translesional pressure ratio, the AUC was 0.71 (p < 0.001) for angle γ, and was 0.87 (p < 0.001) for stenosis degree. For predicting high translesional wall shear stress ratio, the AUC was 0.62 (p = 0.04) for lumen area of internal carotid artery, and was 0.77 (p < 0.001) for stenosis degree.

CONCLUSIONS

Apart from stenosis degree, other geometric characteristics of lesions may also have an influence on hemodynamics of the stenotic carotid artery.

Assessment of extracranial carotid artery disease using digital twins - A pilot study

To improve risk stratification in extracranial internal carotid artery disease (CAD), patients who would benefit maximally from revascularization must be identified. In cardiology, the fractional flow reserve (FFR) has become a reference standard for evaluating the functional severity of coronary artery stenosis, and noninvasive surrogates thereof relying on computational fluid dynamics (CFD) have been developed. Here, we present a CFD-based workflow using digital twins of patients' carotid bifurcations derived from computed tomography angiography for the noninvasive functional assessment of CAD. We reconstructed patient-specific digital twins of 37 carotid bifurcations. We implemented a CFD model using common carotid artery peak systolic velocity (PSV) acquired with Doppler ultrasound (DUS) as inlet boundary condition and a two-element Windkessel model as oulet boundary condition. The agreement between CFD and DUS on the PSV in the internal carotid artery (ICA) was then compared. The relative error for the agreement between DUS and CFD was 9% ± 20% and the intraclass correlation coefficient was 0.88. Furthermore, hyperemic simulations in a physiological range were feasible and unmasked markedly different pressure drops along two ICA stenoses with similar degree of narrowing under comparable ICA blood flow. Hereby, we lay the foundation for prospective studies on noninvasive CFD-based derivation of metrics similar to the FFR for the assessment of CAD.

Correlation between Carotid Stenosis Degree and Blood Pressure Variability in Patients with Carotid Stenosis

Background

Carotid artery stenosis is one of the most serious diseases that endanger human health in contemporary times. It is a frequently occurring and common disease of the middle-aged and elderly people. Its incidence is increasing year by year, bringing a heavy economic burden to society and families. Whether there is a relationship between the degree of carotid artery stenosis and blood pressure variability is less studied.

Aims

To investigate the correlation between the degree of carotid stenosis and blood pressure variability in patients with carotid stenosis.

Materials and Methods

A total of 200 patients with carotid artery stenosis who were treated in our hospital from January 2017 to January 2020 were selected as the subjects of prospective study and were divided into mild stenosis according to the degree of carotid stenosis (carotid artery stenosis rate was 0-50%), moderate stenosis (carotid artery stenosis rate was between 50% and 70%), severe stenosis (carotid artery stenosis rate≧70%), and the control group with 50 cases each. The correlations between the hemodynamics, the degree of carotid artery stenosis, and blood pressure variability in patients with carotid artery stenosis were analyzed.

Results

The levels of 24hSSD, 24hDSD, dSSD, dDSD, and nSSD in the mild stenosis group and moderate stenosis group were significantly higher than those in the control group. In the stenosis group, the levels of 24hSSD, 24hDSD, dSSD, dDSD, and nSSD in the severe stenosis group were significantly higher than those in the moderate stenosis group, with statistical significance (P < 0.05). The levels of PSV, EDV, and MV in the mild stenosis group and moderate stenosis group were lower than those in the control group, while the PI and RI indexes were higher than those in the control group. PI and RI levels were significantly higher than those in the mild stenosis group and moderate stenosis group (P < 0.05). Logistic analysis showed that EDV (P = 0.001, OR = 2.245, 95%CI = 1.638 ~ 3.078), SSD (P = 0.014, OR = 0.725, 95%CI = 0.528 ~ 0.996), and PSV (P = 0.001, OR = 1.970, 95%CI = 1.300 ~ 2.990) were closely related to the degree of carotid artery stenosis.

Conclusion

Hemodynamics and blood pressure variability are related to the severity of carotid stenosis, which provides a reference and basis for clinical treatment of carotid stenosis.

Optimal Design and Biomechanical Analysis of a Biomimetic Lightweight Design Plate for Distal Tibial Fractures: A Finite Element Analysis

The treatment of fractures of the distal tibia can be problematic due to the insubstantial soft-tissue covering this part of the anatomy. This study investigates a novel strategy for minimally invasive plate osteosynthesis of distal tibia fractures called bionic lightweight design plating. Following the structure of the animal trabecular bone, we utilized topological mathematical methods to redesign the material layout of the internal fixation device to fulfill the desired lightweight design within given boundary conditions. The results showed that this method can maintain the same stability of the construct as the original plate after a reduction in the original volume by 30%, and the differences in strain energy of plates and maximum node displacement of constructs between the constructs [RP construct vs. LP construct] were not statistically significant (p > 0.05). In the safety assessment of the constructs, the peak stress of plates between constructs was found to not be statistically significantly different under a doubled physiological load (p > 0.05). The average stress of the plates’ elements exceeding the allowable stress was analyzed, and no statistically significant differences were found between the two constructs under axial compression stress conditions (p > 0.05). The average stress of the plates’ elements in the redesigned plating construct under torsional stress conditions was 3.08% less than that of the locked plating construct (p < 0.05). Under the double physiological load condition, 89% of the elements of the plate in the redesigned plating construct and 85% of the elements of the plate in the locked plating construct were lower than the maximum safe stress of the plate, which was 410 MPa (secondary allowable stresses). That reminds us the topology optimization offer a possible way to improve the capacity of soft tissue protection while ensuring the safety of the RP construct by reducing the volume of the implants.

Underlying mechanism of hemodynamics and intracranial aneurysm

In modern society, subarachnoid hemorrhage, mostly caused by intracranial aneurysm rupture, is accompanied by high disability and mortality rate, which has become a major threat to human health. Till now, the etiology of intracranial aneurysm has not been entirely clarified. In recent years, more and more studies focus on the relationship between hemodynamics and intracranial aneurysm. Under the physiological condition, the mechanical force produced by the stable blood flow in the blood vessels keeps balance with the structure of the blood vessels. When the blood vessels are stimulated by the continuous abnormal blood flow, the functional structure of the blood vessels changes, which becomes the pathophysiological basis of the inflammation and atherosclerosis of the blood vessels and further promotes the occurrence and development of the intracranial aneurysm. This review will focus on the relationship between hemodynamics and intracranial aneurysms, will discuss the mechanism of occurrence and development of intracranial aneurysms, and will provide a new perspective for the research and treatment of intracranial aneurysms.

Validation and Diagnostic Performance of a CFD-Based Non-invasive Method for the Diagnosis of Aortic Coarctation

Purpose: The clinical diagnosis of aorta coarctation (CoA) constitutes a challenge, which is usually tackled by applying the peak systolic pressure gradient (PSPG) method. Recent advances in computational fluid dynamics (CFD) have suggested that multi-detector computed tomography angiography (MDCTA)-based CFD can serve as a non-invasive PSPG measurement. The aim of this study was to validate a new CFD method that does not require any medical examination data other than MDCTA images for the diagnosis of CoA.

Materials and methods: Our study included 65 pediatric patients (38 with CoA, and 27 without CoA). All patients underwent cardiac catheterization to confirm if they were suffering from CoA or any other congenital heart disease (CHD). A series of boundary conditions were specified and the simulated results were combined to obtain a stenosis pressure-flow curve. Subsequently, we built a prediction model and evaluated its predictive performance by considering the AUC of the ROC by 5-fold cross-validation.

Results: The proposed MDCTA-based CFD method exhibited a good predictive performance in both the training and test sets (average AUC: 0.948 vs. 0.958; average accuracies: 0.881 vs. 0.877). It also had a higher predictive accuracy compared with the non-invasive criteria presented in the European Society of Cardiology (ESC) guidelines (average accuracies: 0.877 vs. 0.539).

Conclusion: The new non-invasive CFD-based method presented in this work is a promising approach for the accurate diagnosis of CoA, and will likely benefit clinical decision-making.

Geometric determinants of local hemodynamics in severe carotid artery stenosis

In cases of severe carotid artery stenosis (CAS), carotid endarterectomy (CEA) is performed to recover lumen patency and alleviate stroke risk. Under current guidelines, the decision to surgically intervene relies primarily on the percent loss of native arterial lumen diameter within the stenotic region (i.e. the degree of stenosis). An underlying premise is that the degree of stenosis modulates flow-induced wall shear stress elevations at the lesion site, and thus indicates plaque rupture potential and stroke risk. Here, we conduct a retrospective study on pre-CEA computed tomography angiography (CTA) images from 50 patients with severe internal CAS (>60% stenosis) to better understand the influence of plaque and local vessel geometry on local hemodynamics, with geometrical descriptors that extend beyond the degree of stenosis. We first processed CTA images to define a set of multipoint geometric metrics characterizing the stenosed region, and next performed computational fluid dynamics simulations to quantify local wall shear stress and associated hemodynamic metrics. Correlation and regression analyses were used to relate obtained geometric and hemodynamic metrics, with inclusion of patient sub-classification based on the degree of stenosis. Our results suggest that in the context of severe CAS, prediction of shear stress-based metrics can be enhanced by consideration of readily available, multipoint geometric metrics in addition to the degree of stenosis.

Intelligent Diabetes Detection System based on Tongue Datasets

Background:

Scanning Electron Microscope (SEM) Camera Imaging shows and helps analyze hidden organs in the human body. SEM image analysis provides in-depth and critical details of organ abnormalities. Similarly, the human tongue finds use in the detection of organ dysfunction with tongue reflexology.

Objective:

To detect diabetes at an early stage using a non-invasive method of diabetes detection through tongue images and to utilize the reasonable cost of modality (SEM camera) for capturing the tongue images instead of the existing and expensive imaging modalities like X-ray, Computed Tomography, Magnetic Resonance Imaging, Positron Emission Tomography, Single-Photon Emission Computed Tomography etc.

Methods:

The tongue image is captured via SEM camera, it is preprocessed to remove noise and resize the tongue such that it is suitable for segmentation. Greedy Snake Algorithm (GSA) is used to segment the tongue image. The texture features of the tongue are analyzed and finally it is classified as diabetic or normal.

Results:

Failure of organs stomach, intestine, liver and pancreas results in change of the color of the tongue, coating thickness and cracks on the tongue. Changes in pancreas proactive behavior also reflect on tongue coating. The tongue coating texture varies from white or vanilla to yellow also the tongue coating thickness also increases.

Conclusion:

In this paper, the author proposes to diagnose Diabetes Type2 (DT2) at an early stage from tongue digital image. The tongue image is acquired and processed with Greedy Snake Algorithm (GSA) to extract edge and texture features.

Management of Blood Pressure After Acute Ischemic Stroke

PURPOSE OF REVIEW

The present manuscript examines the significance of blood pressure elevation in patients with acute ischemic stroke, the physiologic principles worthy of consideration during its treatment, and the recent empirical evidence that should guide management protocols. It also provides a sound and practical approach to treatment along the time continuum, with particular relevance to reperfusion strategies.

RECENT FINDINGS

The existing evidence shows that both insufficient and excessive blood pressures are detrimental to the outcome of patients with acute ischemic stroke. This "U-shaped" relation, however, relates to measurements at the time of presentation, and clinical studies lack detail and specificity relative to differential measurements along the time continuum, particularly prior to and following reperfusion. Extrapolating from recent series, it is possible to construct treatment protocols balanced for effectiveness and safety. The management of blood pressure after acute ischemic stroke is an important, complex, and challenging aspect of care, requiring a thorough understanding of cerebrovascular physiology. Along the time continuum, the therapeutic priorities start with the preservation of penumbral tissue prior to reperfusion and then follow with the limitation of the damaging effects of excessive blood pressure readings after reperfusion, optimizing the chances of improved outcomes.

Non-Invasive Assessment of Intravascular Pressure Gradients: A Review of Current and Proposed Novel Methods

Invasive catheterization is associated with a low risk of serious complications. However, although it is the gold standard for measuring pressure gradients, it induces changes to blood flow and requires significant resources. Therefore, non-invasive alternatives are urgently needed. Pressure gradients are routinely estimated non-invasively in clinical settings using ultrasound and calculated with the simplified Bernoulli equation, a method with several limitations. A PubMed literature search on validation of non-invasive techniques was conducted, and studies were included if non-invasively estimated pressure gradients were compared with invasively measured pressure gradients in vivo. Pressure gradients were mainly estimated from velocities obtained with Doppler ultrasound or magnetic resonance imaging. Most studies used the simplified Bernoulli equation, but more recent studies have employed the expanded Bernoulli and Navier⁻Stokes equations. Overall, the studies reported good correlation between non-invasive estimation of pressure gradients and catheterization. Despite having strong correlations, several studies reported the non-invasive techniques to either overestimate or underestimate the invasive measurements, thus questioning the accuracy of the non-invasive methods. In conclusion, more advanced imaging techniques may be needed to overcome the shortcomings of current methods.

Blood homocysteine levels could predict major adverse cardiac events in patients with acute coronary syndrome: A STROBE-compliant observational study

The Global Registry of Acute Coronary Events (GRACE) risk score independently predicts major adverse cardiac events (MACEs) in patients with acute coronary syndrome (ACS). This study aims to evaluate whether the level of plasma homocysteine in addition to the GRACE score enhances the predictive value for MACEs in patients with acute coronary syndrome.A total of 361 patients with ACS evaluated at our hospital were included in the study and tested for blood homocysteine levels. We recorded 40 (11.1%) instances of MACE during a median follow-up of 43.3 months (quartile 40.6-44.4 months), including 29 cases (8.0%) of all-cause death and 11 cases (3.1%) of nonfatal myocardial infarction.The GRACE score was significantly associated with homocysteine levels, and multivariate Cox regression analysis showed that both the GRACE risk score and homocysteine content were independent predictors of MACEs (HR 2.63; 95% confidence interval (CI) 1.54 to 4.49; P < .001 and 2.27; 1.06 to 4.86; P = .035, respectively). Moreover, meta-analysis showed that as the homocysteine level increased, the incidence of MACEs also increased (log-rank 8.41; P = .015). GRACE scores adjusted by homocysteine level increased the area under the curve (AUC) from 0.78 to 0.83 (P = 0.006).Blood homocysteine levels are significantly associated with the GRACE risk score, and using both parameters can further improve risk stratification in patients with acute coronary syndrome.

The effect of inlet and outlet boundary conditions in image-based CFD modeling of aortic flow

Background

Computational modeling of cardiovascular flow is a growing and useful field, but such simulations usually require the researcher to guess the flow’s inlet and outlet conditions since they are difficult and expensive to measure. It is critical to determine the amount of uncertainty introduced by these assumptions in order to evaluate the degree to which cardiovascular flow simulations are accurate. Our work begins to address this question by examining the sensitivity of flow to several different assumed velocity inlet and outlet conditions in a patient-specific aorta model.

Methods

We examined the differences between plug flow, parabolic flow, linear shear flows, skewed cubic flow profiles, and Womersley flow at the inlet. Only the shape of the inlet velocity profile was varied—all other parameters were identical among these simulations. Secondary flow in the form of a counter-rotating pair of vortices was also added to parabolic axial flow to study its effect on the solution. In addition, we examined the differences between two-element Windkessel, three element Windkessel and the outflow boundary conditions. In these simulations, only the outlet boundary condition was varied.

Results

The results show axial and in-plane velocities are considerably different close to the inlet for the cases with different inlet velocity profile shapes. However, the solutions are qualitatively similar beyond 1.75D, where D is the inlet diameter. This trend is also observed in other quantities such as pressure and wall shear stress. Normalized root-mean-square deviation, a measure of axial velocity magnitude differences between the different cases, generally decreases along the streamwise coordinate. The linear shear inlet velocity boundary condition and plug velocity boundary condition solution exhibit the highest time-averaged wall shear stress, approximately \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8\%$$\end{document}8% higher than the parabolic inlet velocity boundary condition. Upstream of 1D from the inlet, adding secondary flow has a significant impact on temporal wall shear stress distributions. This is especially observable during diastole, when integrated wall shear stress magnitude varies about \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$26\%$$\end{document}26% between simulations with and without secondary flow. The results from the outlet boundary condition study show the Windkessel models differ from the outflow boundary condition by as much as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$18\%$$\end{document}18% in terms of time-averaged wall shear stress. Furthermore, normalized root-mean-square deviation of axial velocity magnitude, a measure of deviation between Windkessel and the outflow boundary condition, increases along the streamwise coordinate indicating larger variations near outlets.

Conclusion

It was found that the selection of inlet velocity conditions significantly affects only the flow region close to the inlet of the aorta. Beyond two diameters distal to the inlet, differences in flow solution are small. Although additional studies must be performed to verify this result, the data suggest that it is important to use patient-specific inlet conditions primarily if the researcher is concerned with the details of the flow very close to the inlet. Similarly, the selection of outlet conditions significantly affects the flow in the vicinity of the outlets. Upstream of five diameters proximal to the outlet, deviations between the outlet boundary conditions examined are insignificant. Although the inlet and outlet conditions only affect the flow significantly in their respective neighborhoods, our study indicates that outlet conditions influence a larger percentage of the solution domain.

Clinical validation and assessment of aortic hemodynamics using computational fluid dynamics simulations from computed tomography angiography

Background

Hemodynamic information including peak systolic pressure (PSP) and peak systolic velocity (PSV) carry an important role in evaluation and diagnosis of congenital heart disease (CHD). Since MDCTA cannot evaluate hemodynamic information directly, the aim of this study is to provide a noninvasive method based on a computational fluid dynamics (CFD) model, derived from multi-detector computed tomography angiography (MDCTA) raw data, to analyze the aortic hemodynamics in infants with CHD, and validate these results against echocardiography and cardiac catheter measurements.

Methods

This study included 25 patients (17 males, and 8 females; a median age of 2 years, range: 4 months–4 years) with CHD. All patients underwent both transthoracic echocardiography (TTE) and MDCTA within 2 weeks prior to cardiac catheterization. CFD models were created from MDCTA raw data. Boundary conditions were confirmed by lumped parameter model and transthoracic echocardiography (TTE). Peak systolic velocity derived from CFD models (PSV_CFD) was compared to TTE measurements (PSV_TTE), while the peak systolic pressure derived from CFD (PSP_CFD) was compared to catheterization (PSP_CC). Regions with low and high peak systolic wall shear stress (PSWSS) were also evaluated.

Results

PSV_CFD and PSP_CFD showed good agreements between PSV_TTE (r = 0.968, p < 0.001; mean bias = − 7.68 cm/s) and PSP_CC (r = 0.918, p < 0.001; mean bias = 1.405 mmHg). Regions with low and high PSWSS) can also be visualized. Skewing of velocity or helical blood flow was also observed at aortic arch in patients.

Conclusions

Our result demonstrated that CFD scheme based on MDCTA raw data is an accurate and convenient method in obtaining the velocity and pressure from aorta and displaying the distribution of PSWSS and flow pattern of aorta. The preliminary results from our study demonstrate the capability in combining clinical imaging data and novel CFD tools in infants with CHD and provide a noninvasive approach for diagnose of CHD such as coarctation of aorta in future.

Carotid DSA based CFD simulation in assessing the patient with asymptomatic carotid stenosis: a preliminary study

Background

Cerebrovascular events are frequently associated with hemodynamic disturbance caused by internal carotid artery (ICA) stenosis. It is challenging to determine the ischemia-related carotid stenosis during the intervention only using digital subtracted angiography (DSA). Inspired by the performance of well-established FFRct technique in hemodynamic assessment of significant coronary stenosis, we introduced a pressure-based carotid arterial functional assessment (CAFA) index generated from computational fluid dynamic (CFD) simulation in DSA data, and investigated its feasibility in the assessment of hemodynamic disturbance preliminarily using pressure-wired measurement and arterial spin labeling (ASL) MRI as references.

Methods

The cerebral multi-delay multi-parametric ASL-MRI and carotid DSA including trans-stenotic pressure-wired measurement were implemented on a 65-year-old man with asymptomatic unilateral (left) ICA stenosis. A CFD simulation using simplified boundary condition was performed in DSA data to calculate the CAFA index. The cerebral blood flow (CBF) and arterial transit time (ATT) of ICA territories were acquired.

Results

CFD simulation showed good correlation (r = 0.839, P = 0.001) with slight systematic overestimation (mean difference − 0.007, standard deviation 0.017) compared with pressure-wired measurement. No significant difference was observed between them (P = 0.09). Though the narrowing degree of in the involved ICA was about 70%, the simulated and measured CAFA (0.942/0.937) revealed a functionally nonsignificant stenosis which was also verified by a compensatory final CBF (fronto-temporal/fronto-parietal region: 51.58/45.62 ml/100 g/min) and slightly prolonged ATT (1.23/1.4 s) in the involved territories, together with a normal left–right percentage difference (2.1–8.85%).

Conclusions

The DSA based CFD simulation showed good consistence with invasive approach and could be used as a cost-saving and efficient way to study the relationship between hemodynamic disorder caused by ICA stenosis and subsequent perfusion variations in brain. Further research should focus on the role of noninvasive pressure-based CAFA in screening asymptomatic ischemia-causing carotid stenosis.

A patient-specific lumped-parameter model of coronary circulation

A new lumped-parameter model for coronary hemodynamics is developed. This model is developed for the whole coronary network based on CT scans of a patient-specific geometry including the right coronary tree, which is absent in many previous mathematical models. The model adopts the structured tree model boundary conditions similar to the work of Olufsen et al., thus avoiding the necessity of invasive perfusion measurements. In addition, we also incorporated the effects of the head loss at the two inlets of the large coronary arteries for the first time. The head loss could explain the phenomenon of a sudden increase of the resistance at the inlet of coronary vessel. The estimated blood pressure and flow rate results from the model agree well with the clinical measurements. The computed impedances also match the experimental perfusion measurement. The effects of coronary arterial stenosis are considered and the fractional flow reserve and relative flow in the coronary vessels for a stenotic vessel computed in this model show good agreement with published experimental data. It is believed that the approach could be readily translated to clinical practice to facilitate real time clinical diagnosis.