Large Neck and Strong Ostium Inflow as the Potential Causes for Delayed Occlusion of Unruptured Sidewall Intracranial Aneurysms Treated by Flow Diverter.

BACKGROUND AND PURPOSE: Flow diverter-induced hemodynamic change plays an important role in the mechanism of intracranial aneurysm occlusion. Our aim was to explore the relationship between aneurysm features and flow-diverter treatment of unruptured sidewall intracranial aneurysms. MATERIALS AND METHODS: MR imaging, 4D phase-contrast, was prospectively performed before flow diverter implantation in each patient with unruptured intracranial aneurysm. Two postprocedure follow-ups were scheduled at 6 and 12 months. Responses were grouped according to whether the aneurysms were occluded or remnant. Preprocedural aneurysm geometries and ostium hemodynamics in 38 patients were compared between the 2 groups at 6 and 12 months. Receiver operating characteristic curve analyses were performed for significant geometric and hemodynamic continuous parameters. RESULTS: After the 6-month assessment, 21 of 41 intracranial aneurysms were occluded, and 9 additional aneurysms were occluded at 12 months. Geometrically, the ostium maximum diameter was significantly larger in the remnant group at 6 and 12 months (both P < .001). Hemodynamically, the proximal inflow zone was more frequently observed in the remnant group at 6 months. Several preprocedural ostium hemodynamic parameters were significantly higher in the remnant group. As a prediction for occlusion, the areas under the curve of the ostium maximum diameter (for 6 and 12 months), systolic inflow rate ratio (for 6 months), and systolic inflow area (for 12 months) reached 0.843, 0.883, 0.855, and 0.860, respectively. CONCLUSIONS: Intracranial aneurysms with a large ostium and strong ostium inflow may need a longer time for occlusion. Preprocedural 4D flow MR imaging can well illustrate ostium hemodynamics and characterize aneurysm treatment responses.

How Flow Reduction Influences the Intracranial Aneurysm Occlusion: A Prospective 4D Phase-Contrast MRI Study.

BACKGROUND AND PURPOSE: Flow-diverter stents are widely used for the treatment of wide-neck intracranial aneurysms. Various parameters may influence intracranial aneurysm thrombosis, including the flow reduction induced by flow-diverter stent implantation, which is assumed to play a leading role. However, its actual impact remains unclear due to the lack of detailed intra-aneurysmal flow measurements. This study aimed to clarify this relationship by quantitatively measuring the intra-aneurysmal flow using 4D phase-contrast MR imaging. MATERIALS AND METHODS: We acquired prospective pre- and post-stent implantation 4D phase-contrast MR imaging data of a consecutive series of 23 patients treated with flow-diverter stents. Velocity field data were combined with the intraprocedural 3D angiogram vessel geometries for precise intracranial aneurysm extraction and partial volume correction. Intra-aneurysmal hemodynamic modifications were compared with occlusion outcomes at 6 and 12 months. RESULTS: The averaged velocities at systole were lower after flow-diverter stent implantation for all patients and ranged from 21.7 ± 7.1 cm/s before to 7.2 ± 2.9 cm/s after stent placement. The velocity reduction was more important for the group of patients with aneurysm thrombosis at 6 months (68.8%) and decreased gradually from 66.2% to 55% for 12-month thrombosis and no thrombosis, respectively (P = .08). CONCLUSIONS: We propose an innovative approach to measure intracranial flow changes after flow-diverter stent implantation. We identified a trend between flow reduction and thrombosis outcome that brings a new insight into current understanding of the flow-diversion treatment response.

Errors in power-law estimations of inflow rates for intracranial aneurysm CFD.

Patient-specific inflow rates are rarely available for computational fluid dynamics (CFD) studies of intracranial aneurysms. Instead, inflow rates are often estimated from parent artery diameters via power laws, i.e. Q ∝ Dn, reflecting adaptation of conduit arteries to demanded flow. The present study aimed to validate the accuracy of these power laws. Internal carotid artery (ICA) flow rates were measured from 25 ICA aneurysm patients via 2D phase contrast MRI. ICA diameters, derived from 3D segmentation of rotational angiograms, were used to estimate inflow rates via power laws from the aneurysm CFD literature assuming the same inlet wall shear stress (WSS) (n = 3), velocity (n = 2) or flow rate (n = 0) for all cases. To illustrate the potential impact of errors in flow rate estimates, pulsatile CFD was carried out for four cases having large errors for at least one power law. Flow rates estimated by n = 3 and n = 0 power laws had significant (p < 0.01) mean biases of -22% to +32%, respectively, but with individual errors ranging from -78% to +120%. The n = 2 power law had no significant bias, but had non-negligible individual errors of -58% to +71%. CFD showed similarly large errors for time-averaged sac WSS; however, these were reduced after normalizing by parent artery WSS. High frequency WSS fluctuations, evident in 2/4 aneurysms, were also sensitive to inflow rate errors. Care should therefore be exercised in the interpretation of aneurysm CFD studies that rely on power law estimates of inflow rates, especially if absolute (vs. normalized) WSS, or WSS instabilities, are of interest.

Better Than Nothing: A Rational Approach for Minimizing the Impact of Outflow Strategy on Cerebrovascular Simulations.

BACKGROUND AND PURPOSE: Computational fluid dynamics simulations of neurovascular diseases are impacted by various modeling assumptions and uncertainties, including outlet boundary conditions. Many studies of intracranial aneurysms, for example, assume zero pressure at all outlets, often the default ("do-nothing") strategy, with no physiological basis. Others divide outflow according to the outlet diameters cubed, nominally based on the more physiological Murray's law but still susceptible to subjective choices about the segmented model extent. Here we demonstrate the limitations and impact of these outflow strategies, against a novel "splitting" method introduced here. MATERIALS AND METHODS: With our method, the segmented lumen is split into its constituent bifurcations, where flow divisions are estimated locally using a power law. Together these provide the global outflow rate boundary conditions. The impact of outflow strategy on flow rates was tested for 70 cases of MCA aneurysm with 0D simulations. The impact on hemodynamic indices used for rupture status assessment was tested for 10 cases with 3D simulations. RESULTS: Differences in flow rates among the various strategies were up to 70%, with a non-negligible impact on average and oscillatory wall shear stresses in some cases. Murray-law and splitting methods gave flow rates closest to physiological values reported in the literature; however, only the splitting method was insensitive to arbitrary truncation of the model extent. CONCLUSIONS: Cerebrovascular simulations can depend strongly on the outflow strategy. The default zero-pressure method should be avoided in favor of Murray-law or splitting methods, the latter being released as an open-source tool to encourage the standardization of outflow strategies.

Vessel calibre and flow splitting relationships at the internal carotid artery terminal bifurcation.

OBJECTIVE: Vessel lumen calibres and flow rates are thought to be related by mathematical power laws, reflecting the optimization of cardiac versus metabolic work. While these laws have been confirmed indirectly via measurement of branch calibres, there is little data confirming power law relationships of flow distribution to branch calibres at individual bifurcations. APPROACH: Flow rates and diameters of parent and daughter vessels of the internal carotid artery terminal bifurcation were determined, via robust and automated methods, from 4D phase-contrast magnetic resonance imaging and 3D rotational angiography of 31 patients. MAIN RESULTS: Junction exponents were 2.06 ± 0.44 for relating parent to daughter branch diameters (geometrical exponent), and 2.45 ± 0.75 for relating daughter branch diameters to their flow division (flow split exponent). These exponents were not significantly different, but showed large inter- and intra-individual variations, and with confidence intervals excluding the theoretical optimum of 3. Power law fits of flow split versus diameter ratio and pooled flow rates versus diameters showed exponents of 2.17 and 1.96, respectively. A significant negative correlation was found between age and the geometrical exponent (r = -0.55, p = 0.003) but not the flow split exponent. We also found a dependence of our results on how lumen diameter is measured, possibly explaining some of the variability in the literature. SIGNIFICANCE: Our study confirms that, on average, division of flow to the middle and anterior cerebral arteries is related to these vessels' relative calibres via a power law, but it is closer to a square law than a cube law as commonly assumed.

Improved reduced-order modelling of cerebrovascular flow distribution by accounting for arterial bifurcation pressure drops.

Reduced-order modelling offers the possibility to study global flow features in cardiovascular networks. In order to validate these models, previous studies have been conducted in which they compared 3D computational fluid dynamics simulations with reduced-order simulations. Discrepancies have been reported between the two methods. The loss of energy at the bifurcations is usually neglected and has been pointed out as a possible explanation for these discrepancies. We present distributed lumped models of cerebrovasculatures created automatically from 70 cerebrovascular networks segmented from 3D angiograms. The outflow rate repartitions predicted with and without modelling the energy loss at the bifurcations are compared against 3D simulations. When neglecting the energy loss at the bifurcations, the flow rates though the anterior cerebral arteries are overestimated by 4.7±6.8% (error relative to the inlet flow rate, mean ± standard deviation), impacting the remaining volume of flow going to the other vessels. When the energy loss is modelled, this error is dropping to 0.1±3.2%. Overall, over the total of 337 outlet vessels, when the energy losses at the bifurcations are not modelled the 95% of agreement is in the range of ±13.5% and is down to ±6.5% when the energy losses are considered. With minimal input and computational resources, the presented method can estimate the outflow rates reliably. This study constitutes the largest validation of a reduced-order flow model against 3D simulations. The impact of the energy loss at the bifurcations is here demonstrated for cerebrovasculatures but can be applied to other physiological networks.

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

Virtual-versus-Real Implantation of Flow Diverters: Clinical Potential and Influence of Vascular Geometry.

BACKGROUND AND PURPOSE: Intracranial stents have become extremely important in the endovascular management of complex intracranial aneurysms. Sizing and landing zone predictions are still very challenging steps in the procedure. Virtual stent deployment may help therapeutic planning, device choice, and hemodynamic simulations. We aimed to assess the predictability of our recently developed virtual deployment model by comparing in vivo and virtual stents implanted in a consecutive series of patients presenting with intracranial aneurysms. MATERIALS AND METHODS: Virtual stents were implanted in patient-specific geometries of intracranial aneurysms treated with the Pipeline Embolization Device. The length and cross-section of virtual and real stents measured with conebeam CT were compared. The influence of vessel geometry modifications occurring during the intervention was analyzed. RESULTS: The virtual deployment based on pre- and poststent implantation 3D rotational angiography overestimated (underestimated) the device length by 13% ± 11% (-9% ± 5%). These differences were highly correlated (R2 = 0.67) with the virtual-versus-real stent radius differences of -6% ± 7% (5% ± 4%) for predictions based on pre- and poststent implantation 3D rotational angiography. These mismatches were due principally to implantation concerns and vessel-shape modifications. CONCLUSIONS: The recently proposed geometric model was shown to predict accurately the deployment of Pipeline Embolization Devices when the stent radius was well-assessed. However, unpredictable delivery manipulations and variations of vessel geometry occurring during the intervention might impact the stent implantation.

Improving workflows of neuro-interventional procedures with autostereoscopic 3D visualization of multi-modality imaging in hybrid interventional suites.

PURPOSE: Recent developments in interventional neuroradiology techniques, medical imaging modalities, endovascular stenting and embolization materials lead to an increasing number of patients with cerebral aneurysms and arteriovenous malformations that are eligible for endovascular treatment and have opened new perspectives for novel ways for patient treatment in general. In this paper, we describe a software tool for 3D image fusion of multi-modal acquisitions to assist endovascular treatment of cerebral malformations. The software and an autostereoscopic 3D display were implemented and tested in clinical applications in a hybrid interventional suite that is used for radio-interventional as well as neurosurgical procedures. Our hypothesis is that fusion of image data acquired prior to intervention procedures with images acquired during those procedures should allow better visualizing and navigating through complex cerebral vasculature. This should also improve workflows of neuro-interventional procedures. METHODS: Preoperative and intra-operative acquisitions of vascular images of the brain were performed and transferred to a dedicated imaging workstation to be processed with our image fusion and visualization software tool. The tool was developed as a plugin extension to the open-source DICOM viewer OsiriX and is based on a modular and scalable architecture. Several processing modules were implemented to allow spatial co-registration and fusion of preoperative and intra-operative modalities. A special extension was also implemented for interactive autostereosopic, glass-free 3D visualization of fused results. RESULTS: The software platform was validated and evaluated in nine in vivo procedures by expert users. All patient cases were related to interventional treatment of neuro-vascular diseases. The emphasis was laid on the added value of spatial co-registration and fusion of preoperative and intra-operative modalities, as well as the overall impact on workflow during the intervention. The co-registered and fused images were visualized on an autostereoscopic 3D monitor installed in hybrid interventional suite. All experiments were evaluated and scored by interventional physicians and technicians. CONCLUSIONS: Displaying 3D-4D representations of brain vascular anomalies based on multi-modal acquisitions on a 3D autostereoscopic display is beneficial for the workflow and efficiency of interventional radiologists. The implemented software tool fulfills the premise of applicability of an open-source platform for more advanced, multi-modal visualization and processing of brain vascular structures for image-guided therapeutic interventions.

Intra-aneurysmal flow patterns: illustrative comparison among digital subtraction angiography, optical flow, and computational fluid dynamics.

BACKGROUND AND PURPOSE: Digital subtraction angiography is the gold standard vascular imaging and it is used for all endovascular treatment of intracranial anerysms. Optical flow imaging has been described as a potential method to evaluate cerebral hemodynamics through DSA. In this study, we aimed to compare the flow patterns measured during angiography, by using an optical flow method, with those measured by using computational fluid dynamics in intracranial aneurysms. MATERIALS AND METHODS: A consecutive series of 21 patients harboring unruptured saccular intracranial aneurysms who underwent diagnostic angiography before treatment was considered. High-frame-rate digital subtraction angiography was performed to obtain an intra-aneurysmal velocity field by following the cardiac-modulated contrast wave through the vascular structures by using optical flow principles. Additionally, computational fluid dynamics modeling was performed for every case by using patient-specific inlet-boundary conditions measured with the optical flow method from both DSA and 3D rotational angiography datasets. Three independent observers compared qualitatively both the inflow direction and the apparent recirculation in regular DSA, optical flow images, and computational fluid dynamics flow patterns for each patient; κ statistics were estimated. RESULTS: We included 21 patients. In 14 of these 21, the flow patterns were conclusive and matching between the optical flow images and computational fluid dynamics within the same projection view (κ = .91). However, in only 8 of these 14 patients the optical flow images were conclusive and matching regular DSA images (observer κ = 0.87). In 7 of the 21 patients, the flow patterns in the optical flow images were inconclusive, possibly due to improper projection angles. CONCLUSIONS: The DSA-based optical flow technique was considered qualitatively consistent with computational fluid dynamics outcomes in evaluating intra-aneurysmal inflow direction and apparent recirculation. Moreover, the optical flow technique may provide the premises for new solutions for improving the visibility of flow patterns when contrast motion in DSA is not apparent. This technique is a diagnostic method to evaluate intra-aneurysmal flow patterns and could be used in the future for validation and patient evaluation.

Quantification of internal carotid artery flow with digital subtraction angiography: validation of an optical flow approach with Doppler ultrasound.

BACKGROUND AND PURPOSE: Digital subtraction angiography is the reference standard technique to evaluate intracranial vascular anatomy and used on the endovascular treatment of vascular diseases. A dedicated optical flow-based algorithm was applied to DSA to measure arterial flow. The first quantification results of internal carotid artery flow validated with Doppler sonography are reported. MATERIALS AND METHODS: We included 22 consecutive patients who underwent endovascular procedures. To assess the sensitivity of the algorithm to contrast agent-blood mixing dynamics, we acquired high-frame DSA series (60 images/s) with different injection rates: 1.5 mL/s (n = 19), 2.0 mL/s (n = 18), and 3.0 mL/s (n = 13). 3D rotational angiography was used to extract the centerline of the vessel and the arterial section necessary for volume flow calculation. Optical flow was used to measure flow velocities in straight parts of the ICAs; these data were further compared with Doppler sonography data. DSA mean flow rates were linearly regressed on Doppler sonography measurements, and regression slope coefficient bias from value 1 was analyzed within the 95% confidence interval. RESULTS: DSA mean flow rates measured with the optical flow approach significantly matched Doppler sonography measurements (slope regression coefficient, b = 0.83 ± 0.19, P = .05) for injection rate = 2.0 mL/s and circulating volumetric blood flow <6 mL/s. For injection rate = 1.5 mL/s, volumetric blood flow <3 mL/s correlated well with Doppler sonography (b = 0.67 ± 0.33, P = .05). Injection rate = 3.0 mL/s failed to provide DSA-optical flow measurements correlating with Doppler sonography because of the lack of measurable pulsatility. CONCLUSIONS: A new model-free optical flow technique was tested reliably on the ICA. DSA-based blood flow velocity measurements were essentially validated with Doppler sonography whenever the conditions of measurable pulsatility were achieved (injection rates = 1.5 and 2.0 mL/s).

Evaluation of the influence of inlet boundary conditions on computational fluid dynamics for intracranial aneurysms: a virtual experiment.

Inlet boundary conditions (BCs) are important inputs of computational fluid dynamics (CFD) in intracranial aneurysms (IAs). We performed sensibility analysis of CFD to different inlet BCs applied to illustrative patient-specific aneurysm-vessel geometry. BCs corresponding to generic and patient-specific pulsatile flow curves were applied to three vascular geometry models of carotid ophthalmic aneurysm-vessel geometry, in which the inlet lengths were different. CFD outcomes were compared to high frame rate Digital Subtraction Angiography (DSA) sequences. The streamlines were found to match contrast agent (CA) motion pattern in the case where the non-truncated inlet vessel model was coupled to generic Womersley BC solution. Even though dynamic pressure loss (55%) was equal for all models and different BCs, the minimum distance to wall of the fastest velocity fields for the non-truncated model was significantly larger (p=0.002) and mean vorticity sign was different. Significant difference in spatial distributions of wall shear stress (WSS) and oscillating shear stress index (OSI) was found in aneurysm between Womersley and Plugflow BC conditions, only. Reliable CFD for carotid ophthalmic aneurysm would require avoiding truncation of the inlet vessel to be independent of the solution applied to generate CFD.

A DSA-based method using contrast-motion estimation for the assessment of the intra-aneurysmal flow changes induced by flow-diverter stents.

BACKGROUND AND PURPOSE: Flow-diverter stents modify intra-aneurysmal blood flow and induce the progressive thrombosis of intracranial aneurysms followed by stable vascular reconstruction. The aim of this study was to report a new method for the appraisal of intracranial blood flow from DSA performed during endovascular treatment procedures. MATERIALS AND METHODS: A cohort of 24 patients with unruptured IAs who underwent FDS implantation was prospectively recruited. Pre- and post-DSA sequences in combination with 3D rotational angiography were acquired. The quantification of arterial and intra-aneurysmal flow was accomplished by using an optical flow approach. Flow reduction was assessed by using a new metric termed the mean aneurysm flow amplitude ratio. The correlation between the MAFA ratio and the incidence of aneurysm thrombosis was assessed by using receiver operating characteristic analysis and the Fisher exact test when the optimum Youden index was found. RESULTS: The quantification of flow was successfully achieved in 21 of 24 patients (87.5%). On the imaging follow-up, 18 aneurysms developed complete thrombosis (87.5%) and 3 displayed residual circulation (12.5%). The threshold analysis of the MAFA ratio significantly predicted thrombosis at 12 months below a threshold of 1.03 (P=.035). There was no significant correlation between the time for complete occlusion of the aneurysm and contrast stagnation inside the aneurysm after treatment (P>.05). CONCLUSIONS: The MAFA ratio based on DSA flow quantification appears to be a reliable predictor for the assessment of stent treatment outcomes in this small study. These results open the door for perioperative flow quantification and provide indices that may help clinicians make appropriate intraprocedural decisions.