Unraveling the relationship between arterial flow and intra-aneurysmal hemodynamics

A real-time patient-specific treatment strategy for enhanced external counterpulsation

Diastolic/systolic blood pressure ratio (D/S) ≥ 1.2 is the gold standard of enhanced external counterpulsation (EECP) treatment, but it does not show a clear clinical correspondence with the configuration of the EECP mode. As such, a single target results in different treatment effects in different individuals. The local haemodynamic effect (wall shear stress, WSS) of EECP on vascular endothelial cells is conducive to promote the growth of collateral circulation vessels and restore the blood supply distal to the stenosis lesion. Considering the haemodynamic effects of WSS on human arteries, this study developed a real-time patient-specific treatment strategy of EECP for patients with cardio-cerebrovascular diseases. Based on patient-specific haemodynamic data from 113 individuals, an optimization algorithm was developed to achieve the individualization of a 0D lumped-parameter model of the human circulatory system, thereby simulating the patient-specific global haemodynamic effects. 0D/3D coupled cardio-cerebrovascular models of two subjects were established to simulate the local WSS. We then established statistical models to evaluate clinically unmeasurable WSS based on measurable global haemodynamic indicators. With the aim of attaining appropriate area- and time-averaged WSS (ATAWSS, 4-7 Pa), as evaluated by global haemodynamic indicators, a closed-loop feedback tuning method was developed to provide patient-specific EECP treatment strategies. Results showed that for clinical data collected from 113 individuals, the individualized 0D model can accurately simulate patient-specific global haemodynamic effects (average error <5%). Based on two subjects, the statistical models can be used to evaluate local ATAWSS (error <6%) for coronary arteries and for cerebral arteries. An EECP mode planned by the patient-specific treatment strategy can promote an appropriate ATAWSS within a 16 s calculation time. The real-time patient-specific treatment strategy of EECP is expected to improve the long-term outcome for each patient and have potential clinical significance.

Sensitivity of hostile hemodynamics to aneurysm geometry via unsupervised shape interpolation

BACKGROUND AND OBJECTIVE

Vessel geometry and hemodynamics are intrinsically linked, whereby geometry determines hemodynamics, and hemodynamics influence vascular remodeling. Both have been used for testing clinical outcomes, but geometry/morphology generally has less uncertainty than hemodynamics derived from medical image-based computational fluid dynamics (CFD). To provide clinical utility, CFD-based hemodynamic parameters must be robust to modeling errors and/or uncertainties, but must also provide useful information not more-easily extracted from shape alone. The objective of this study was to methodically assess the response of hemodynamic parameters to gradual changes in shape created using an unsupervised 3D shape interpolation method.

METHODS

We trained the neural network NeuroMorph on 3 patient-derived intracranial aneurysm surfaces (labelled A, B, C), and then generated 3 distinct morph sequences (A→B, B→C, C→A) each containing 10 interpolated surfaces. From high-fidelity CFD simulation of these, we calculated a variety of common reduced hemodynamic parameters, including many previously associated with aneurysm rupture, and analyzed their responses to changes in shape, and their correlations.

RESULTS

The interpolated surfaces demonstrate complex, gradual changes in branch angles, vessel diameters, and aneurysm morphology. CFD simulation showed gradual changes in aneurysm jetting characteristics and wall-shear stress (WSS) patterns, but demonstrated a range of responses from the reduced hemodynamic parameters. Spatially and temporally averaged parameters including time-averaged WSS, time-averaged velocity, and low-shear area (LSA) showed low variation across all morph sequences, while parameters of flow complexity such as oscillatory shear, spectral broadening, and spectral bandedness indices showed high variation between slightly-altered neighboring surfaces. Correlation analysis revealed a great deal of mutual information with easier-to-measure shape-based parameters.

CONCLUSIONS

In the absence of large clinical datasets, unsupervised shape interpolation provides an ideal laboratory for exploring the delicate balance between robustness and sensitivity of nominal hemodynamic predictors of aneurysm rupture. Parameters like time-averaged WSS and LSA that are highly "robust" may, as a result, be effectively redundant to morphological predictors, whereas more sensitive parameters may be too uncertain for practical clinical use. Understanding these sensitivities may help identify parameters that are capable of providing added value to rupture risk assessment.

Subarachnoid Hemorrhage from Ruptured Internal Carotid Artery Aneurysm: Association with Arterial Tortuosity

OBJECTIVE

Many researchers have found a correlation between tortuous arteries and development of aneurysms in cerebral arteries. We decided to determine whether tortuosity of the internal carotid artery can be related to its aneurysm rupture.

METHODS

We retrospectively analyzed the internal carotid artery anatomy of 149 patients with internal carotid artery aneurysms. For each patient, we calculated relative length (RL), sum of angle metrics (SOAM), triangular index (TI), product of angle distance (PAD), and inflection count metrics (ICM).

RESULTS

A total of 33 patients (22.15%) had subarachnoid hemorrhage. These patients had significantly lower SOAM (0.31 ± 0.17 vs. 0.42 ± 0.21; P < 0.01), TI (0.27 ± 0.09 vs. 0.31 ± 0.11; P = 0.03) and ICM (0.25 ± 0.11 vs. 0.31 ± 0.17; P = 0.04). In multivariate logistic regression analysis, higher SOAM (odds ratio, 0.780; 95% confidence interval, 0.619-0.961; P = 0.025) remained independently associated with lower risk of internal carotid artery aneurysm rupture. In addition, we found significant positive correlation of aneurysm dome size with SOAM (R = 0.224; P = 0.013) and PAD (0.269; P < 0.01). Our study also showed that age (R = 0.252; P = 0.036), Glasgow Coma Scale score (R = -0.706; P < 0.01), and TI (R = -0.249; P = 0.042) were independently correlated with modified Rankin Scale score on discharge.

CONCLUSIONS

Lower tortuosity might be a protective factor against internal carotid artery aneurysm rupture and poor outcome after subarachnoid hemorrhage. Higher tortuosity is correlated with internal carotid artery aneurysm growth.

Evaluation of aneurysm rupture risk based upon flowrate-independent hemodynamic parameters: a multi-center pilot study

BACKGROUND

Specifying generic flow boundary conditions in aneurysm hemodynamic simulations yields a great degree of uncertainty for the evaluation of aneurysm rupture risk. Herein, we proposed the use of flowrate-independent parameters in discriminating unstable aneurysms and compared their prognostic performance against that of conventional absolute parameters.

METHODS

This retrospective study included 186 aneurysms collected from three international centers, with the stable aneurysms having a minimum follow-up period of 24 months. The flowrate-independent aneurysmal wall shear stress (WSS) and energy loss (EL) were defined as the coefficients of the second-order polynomials characterizing the relationships between the respective parameters and the parent-artery flows. Performance of the flowrate-independent parameters in discriminating unstable aneurysms with the logistic regression, Adaboost, and support-vector machine (SVM) methods was quantified and compared against that of the conventional parameters, in terms of sensitivity, specificity, and area under the curve (AUC).

RESULTS

In discriminating unstable aneurysms, the proposed flowrate-independent EL achieved the highest sensitivity (0.833, 95% CI 0.586 to 0.964) and specificity (0.833, 95% CI 0.672 to 0.936) on the SVM, with the AUC outperforming the conventional EL by 0.133 (95% CI 0.039 to 0.226, p=0.006). Likewise, the flowrate-independent WSS outperformed the conventional WSS in terms of the AUC (difference: 0.137, 95% CI 0.033 to 0.241, p=0.010).

CONCLUSION

The flowrate-independent hemodynamic parameters surpassed their conventional counterparts in predicting the stability of aneurysms, which may serve as a promising set of hemodynamic metrics to be used for the prediction of aneurysm rupture risk when physiologically real vascular boundary conditions are unavailable.

How patient-specific do internal carotid artery inflow rates need to be for computational fluid dynamics of cerebral aneurysms?

BACKGROUND

Computational fluid dynamics (CFD) has become a popular tool for studying 'patient-specific' blood flow dynamics in cerebral aneurysms; however, rarely are the inflow boundary conditions patient-specific. We aimed to test the impact of widespread reliance on generalized inflow rates.

METHODS

Internal carotid artery (ICA) flow rates were measured via 2D cine phase-contrast MRI for 24 patients scheduled for endovascular therapy of an ICA aneurysm. CFD models were constructed from 3D rotational angiography, and pulsatile inflow rates imposed as measured by MRI or estimated using an average older-adult ICA flow waveform shape scaled by a cycle-average flow rate (Q_avg) derived from the patient's ICA cross-sectional area via an assumed inlet velocity.

RESULTS

There was good overall qualitative agreement in the magnitudes and spatial distributions of time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and spectral power index (SPI) using generalized versus patient-specific inflows. Sac-averaged quantities showed moderate to good correlations: R²=0.54 (TAWSS), 0.80 (OSI), and 0.68 (SPI). Using patient-specific Q_avg to scale the generalized waveform shape resulted in near-perfect agreement for TAWSS, and reduced bias, but not scatter, for SPI. Patient-specific waveform had an impact only on OSI correlations, which improved to R²=0.93.

CONCLUSIONS

Aneurysm CFD demonstrates the ability to stratify cases by nominal hemodynamic 'risk' factors when employing an age- and vascular-territory-specific recipe for generalized inflow rates. Q_avg has a greater influence than waveform shape, suggesting some improvement could be achieved by including measurement of patient-specific Q_avg into aneurysm imaging protocols.

Incorporating variability of patient inflow conditions into statistical models for aneurysm rupture assessment

BACKGROUND

Hemodynamic patterns have been associated with cerebral aneurysm instability. For patient-specific computational fluid dynamics (CFD) simulations, the inflow rates of a patient are typically not known. The aim of this study was to analyze the influence of inter- and intra-patient variations of cerebral blood flow on the computed hemodynamics through CFD simulations and to incorporate these variations into statistical models for aneurysm rupture prediction.

METHODS

Image data of 1820 aneurysms were used for patient-specific steady CFD simulations with nine different inflow rates per case, capturing inter- and intra-patient flow variations. Based on the computed flow fields, 17 hemodynamic parameters were calculated and compared for the different flow conditions. Next, statistical models for aneurysm rupture were trained in 1571 of the aneurysms including hemodynamic parameters capturing the flow variations either by defining hemodynamic "response variables" (model A) or repeatedly randomly selecting flow conditions by patients (model B) as well as morphological and patient-specific variables. Both models were evaluated in the remaining 249 cases.

RESULTS

All hemodynamic parameters were significantly different for the varying flow conditions (p < 0.001). Both the flow-independent "response" model A and the flow-dependent model B performed well with areas under the receiver operating characteristic curve of 0.8182 and 0.8174 ± 0.0045, respectively.

CONCLUSIONS

The influence of inter- and intra-patient flow variations on computed hemodynamics can be taken into account in multivariate aneurysm rupture prediction models achieving a good predictive performance. Such models can be applied to CFD data independent of the specific inflow boundary conditions.

Role of Salivary Duct Morphology in the Etiology of Chronic Obstructive Parotitis: Statistical Analysis of Sialographic Features and Computational Fluid Dynamics Analysis of Salivary Flow

PURPOSE

The causes of some cases of chronic obstructive parotitis (COP) without obstructive factors are still unclear. The authors hypothesized that some morphologic features of salivary ducts might contribute to the development of COP. This study investigated the role of salivary duct morphology in the etiology of COP.

MATERIALS AND METHODS

The authors designed and implemented a case-and-control study. Cases were defined as patients with COP, diagnosed from September 2014 to May 2017 at the Affiliated Hospital of Stomatology of the Sun Yat-sen University (Guangzhou, China), and controls were healthy participants. The primary predictor variables were the occurrence of an accessory duct (AD), the number of branches uniting to form the Stensen duct (SD), the angle between the AD and the SD, and the angle between branches identified on sialographic computed tomograms. Data from the 2 groups were compared to investigate the association between these variables and COP. The χ² test, Student t test, and logistic regression were computed, with significance set at a P value less than .05. Fluid dynamics analysis was used to analyze salivary flow field in models of salivary ducts with different morphologic features reconstructed from sialographic computed tomograms.

RESULTS

The sample was composed of 39 patients with COP and 18 controls without COP. The 2 groups were not similar for incidences of an AD (71.8 vs 38.9%) and the angle between branches (96.5 ± 26.0° vs 71.5 ± 21.2°). There was no relevant difference between groups in the number of branches and the angle between the AD and the SD. The area of low velocity was larger in the model with the wider angle between branches.

CONCLUSIONS

The results suggest that the presence of an AD and a wider angle between duct branches are associated with COP.

Analysis of Anterior Cerebral Artery Tortuosity: Association with Anterior Communicating Artery Aneurysm Rupture

BACKGROUND

Many researchers have found a correlation between tortuous arteries and development of aneurysms in cerebral arteries. However, there are no studies analyzing the impact of tortuosity on risk of subarachnoid hemorrhage (SAH) occurrence. Therefore, we decided to determine whether tortuosity of the anterior cerebral artery can be related to the rupture of anterior communicating artery aneurysm and to severity and treatment outcome of SAH.

METHODS

We retrospectively analyzed anterior cerebral artery anatomy of 121 patients with anterior communicating artery aneurysms. From patients' medical records, we obtained their history including previous and current diseases and medications. For each patient we calculated relative length, sum of angle metrics, triangular index, product of angle distance, and inflection count metrics.

RESULTS

Patients with SAH had significantly higher relative length (0.70 ± 0.19 vs. 0.63 ± 0.22; P = 0.03) and significantly lower inflection count metrics (0.10 ± 0.08 vs. 0.16 ± 0.19; P < 0.01), respectively. In multivariate logistic regression analysis, after adjustment of all possible confounders, diabetes mellitus (odds ratio [OR], 0.154; 95% confidence interval [CI], 0.032-0.553; P < 0.01) and higher inflection count metrics (OR, 0.604; 95% CI, 0.357-0.909; P = 0.042) remained independently associated with lower risk of SAH. We also found an independent correlation between aneurysm dome size (R = -0.289; P = 0.02) and triangular index (R = 0.273; P = 0.03) and Glasgow Coma Scale score on admission.

CONCLUSIONS

Higher anterior cerebral artery tortuosity might be a protective factor against anterior communicating artery aneurysm rupture.

Real-World Variability in the Prediction of Intracranial Aneurysm Wall Shear Stress: The 2015 International Aneurysm CFD Challenge

Purpose

Image-based computational fluid dynamics (CFD) is widely used to predict intracranial aneurysm wall shear stress (WSS), particularly with the goal of improving rupture risk assessment. Nevertheless, concern has been expressed over the variability of predicted WSS and inconsistent associations with rupture. Previous challenges, and studies from individual groups, have focused on individual aspects of the image-based CFD pipeline. The aim of this Challenge was to quantify the total variability of the whole pipeline.

Methods

3D rotational angiography image volumes of five middle cerebral artery aneurysms were provided to participants, who were free to choose their segmentation methods, boundary conditions, and CFD solver and settings. Participants were asked to fill out a questionnaire about their solution strategies and experience with aneurysm CFD, and provide surface distributions of WSS magnitude, from which we objectively derived a variety of hemodynamic parameters.

Results

A total of 28 datasets were submitted, from 26 teams with varying levels of self-assessed experience. Wide variability of segmentations, CFD model extents, and inflow rates resulted in interquartile ranges of sac average WSS up to 56%, which reduced to < 30% after normalizing by parent artery WSS. Sac-maximum WSS and low shear area were more variable, while rank-ordering of cases by low or high shear showed only modest consensus among teams. Experience was not a significant predictor of variability.

Conclusions

Wide variability exists in the prediction of intracranial aneurysm WSS. While segmentation and CFD solver techniques may be difficult to standardize across groups, our findings suggest that some of the variability in image-based CFD could be reduced by establishing guidelines for model extents, inflow rates, and blood properties, and by encouraging the reporting of normalized hemodynamic parameters.

Synergistic Integration of Laboratory and Numerical Approaches in Studies of the Biomechanics of Diseased Red Blood Cells

In red blood cell (RBC) disorders, such as sickle cell disease, hereditary spherocytosis, and diabetes, alterations to the size and shape of RBCs due to either mutations of RBC proteins or changes to the extracellular environment, lead to compromised cell deformability, impaired cell stability, and increased propensity to aggregate. Numerous laboratory approaches have been implemented to elucidate the pathogenesis of RBC disorders. Concurrently, computational RBC models have been developed to simulate the dynamics of RBCs under physiological and pathological conditions. In this work, we review recent laboratory and computational studies of disordered RBCs. Distinguished from previous reviews, we emphasize how experimental techniques and computational modeling can be synergically integrated to improve the understanding of the pathophysiology of hematological disorders.

Non-Newtonian versus numerical rheology: Practical impact of shear-thinning on the prediction of stable and unstable flows in intracranial aneurysms

Computational fluid dynamics (CFD) shows promise for informing treatment planning and rupture risk assessment for intracranial aneurysms. Much attention has been paid to the impact on predicted hemodynamics of various modelling assumptions and uncertainties, including the need for modelling the non-Newtonian, shear-thinning rheology of blood, with equivocal results. Our study clarifies this issue by contextualizing the impact of rheology model against the recently demonstrated impact of CFD solution strategy on the prediction of aneurysm flow instabilities. Three aneurysm cases were considered, spanning a range of stable to unstable flows. Simulations were performed using a high-resolution/accuracy solution strategy with Newtonian and modified-Cross rheology models and compared against results from a so-called normal-resolution strategy. Time-averaged and instantaneous wall shear stress (WSS) distributions, as well as frequency content of flow instabilities and dome-averaged WSS metrics, were minimally affected by the rheology model, whereas numerical solution strategy had a demonstrably more marked impact when the rheology model was fixed. We show that point-wise normalization of non-Newtonian by Newtonian WSS values tended to artificially amplify small differences in WSS of questionable physiological relevance in already-low WSS regions, which might help to explain the disparity of opinions in the aneurysm CFD literature regarding the impact of non-Newtonian rheology. Toward the goal of more patient-specific aneurysm CFD, we conclude that attention seems better spent on solution strategy and other likely "first-order" effects (eg, lumen segmentation and choice of flow rates), as opposed to "second-order" effects such as rheology.

Virtual endovascular treatment of intracranial aneurysms: models and uncertainty

Virtual endovascular treatment models (VETMs) have been developed with the view to aid interventional neuroradiologists and neurosurgeons to pre-operatively analyze the comparative efficacy and safety of endovascular treatments for intracranial aneurysms. Based on the current state of VETMs in aneurysm rupture risk stratification and in patient-specific prediction of treatment outcomes, we argue there is a need to go beyond personalized biomechanical flow modeling assuming deterministic parameters and error-free measurements. The mechanobiological effects associated with blood clot formation are important factors in therapeutic decision making and models of post-treatment intra-aneurysmal biology and biochemistry should be linked to the purely hemodynamic models to improve the predictive power of current VETMs. The influence of model and parameter uncertainties associated to each component of a VETM is, where feasible, quantified via a random-effects meta-analysis of the literature. This allows estimating the pooled effect size of these uncertainties on aneurysmal wall shear stress. From such meta-analyses, two main sources of uncertainty emerge where research efforts have so far been limited: (1) vascular wall distensibility, and (2) intra/intersubject systemic flow variations. In the future, we suggest that current deterministic computational simulations need to be extended with strategies for uncertainty mitigation, uncertainty exploration, and sensitivity reduction techniques. WIREs Syst Biol Med 2017, 9:e1385. doi: 10.1002/wsbm.1385 For further resources related to this article, please visit the WIREs website.

Uncertainty quantification of wall shear stress in intracranial aneurysms using a data-driven statistical model of systemic blood flow variability

Adverse wall shear stress (WSS) patterns are known to play a key role in the localisation, formation, and progression of intracranial aneurysms (IAs). Complex region-specific and time-varying aneurysmal WSS patterns depend both on vascular morphology as well as on variable systemic flow conditions. Computational fluid dynamics (CFD) has been proposed for characterising WSS patterns in IAs; however, CFD simulations often rely on deterministic boundary conditions that are not representative of the actual variations in blood flow. We develop a data-driven statistical model of internal carotid artery (ICA) flow, which is used to generate a virtual population of waveforms used as inlet boundary conditions in CFD simulations. This allows the statistics of the resulting aneurysmal WSS distributions to be computed. It is observed that ICA waveform variations have limited influence on the time-averaged WSS (TAWSS) on the IA surface. In contrast, in regions where the flow is locally highly multidirectional, WSS directionality and harmonic content are strongly affected by the ICA flow waveform. As a consequence, we argue that the effect of blood flow variability should be explicitly considered in CFD-based IA rupture assessment to prevent confounding the conclusions.

Does Arterial Flow Rate Affect the Assessment of Flow-Diverter Stent Performance?

BACKGROUND AND PURPOSE

Our aim was to assess the performance of flow-diverter stents. The pre- and end-of-treatment angiographies are commonly compared. However, the arterial flow rate may change between acquisitions; therefore, a better understanding of its influence on the local intra-aneurysmal hemodynamics before and after flow-diverter stent use is required.

MATERIALS AND METHODS

Twenty-five image-based aneurysm models extracted from 3D rotational angiograms were conditioned for computational fluid dynamics simulations. Pulsatile simulations were performed at different arterial flow rates, covering a wide possible range of physiologic flows among 1-5 mL/s. The effect of flow-diverter stents on intra-aneurysmal hemodynamics was numerically simulated with a porous medium model. Spatiotemporal-averaged intra-aneurysmal flow velocity and flow rate were calculated for each case to quantify the hemodynamics after treatment. The short-term flow-diverter stent performance was characterized by the relative velocity reduction inside the aneurysm.

RESULTS

Spatiotemporal-averaged intra-aneurysmal flow velocity before and after flow-diverter stent use is linearly proportional to the mean arterial flow rate (minimum R² > 0.983 of the linear regression models for untreated and stented models). Relative velocity reduction asymptotically decreases with increasing mean arterial flow rate. When the most probable range of arterial flow rate was considered (3-5 mL/s), instead of the wide possible flow range, the mean SD of relative velocity reduction was reduced from 3.6% to 0.48%.

CONCLUSIONS

Both intra-aneurysmal aneurysm velocity and flow-diverter stent performance depend on the arterial flow rate. The performance could be considered independent of the arterial flow rates within the most probable range of physiologic flows.