Christopher M. Putman

Efficient pipeline for image-based patient-specific analysis of cerebral aneurysm hemodynamics: technique and sensitivity

Hemodynamic factors are thought to be implicated in the progression and rupture of intracranial aneurysms. Current efforts aim to study the possible associations of hemodynamic characteristics such as complexity and stability of intra-aneurysmal flow patterns, size and location of the region of flow impingement with the clinical history of aneurysmal rupture. However, there are no reliable methods for measuring blood flow patterns in vivo. In this paper, an efficient methodology for patient-specific modeling and characterization of the hemodynamics in cerebral aneurysms from medical images is described. A sensitivity analysis of the hemodynamic characteristics with respect to variations of several variables over the expected physiologic range of conditions is also presented. This sensitivity analysis shows that although changes in the velocity fields can be observed, the characterization of the intra-aneurysmal flow patterns is not altered when the mean input flow, the flow division, the viscosity model, or mesh resolution are changed. It was also found that the variable that has the greater impact on the computed flow fields is the geometry of the vascular structures. We conclude that with the proposed modeling pipeline clinical studies involving large numbers cerebral aneurysms are feasible.

Association of Hemodynamic Characteristics and Cerebral Aneurysm Rupture

BACKGROUND AND PURPOSE: Hemodynamic factors are thought to play an important role in the initiation, growth, and rupture of cerebral aneurysms. This report describes a study of the associations between qualitative intra-aneurysmal hemodynamics and the rupture of cerebral aneurysms. MATERIALS AND METHODS: Two hundred ten consecutive aneurysms were analyzed by using patient-specific CFD simulations under pulsatile flow conditions. The aneurysms were classified into categories by 2 blinded observers, depending on the complexity and stability of the flow pattern, size of the impingement region, and inflow concentration. A statistical analysis was then performed with respect to the history of previous rupture. Interobserver variability analysis was performed. RESULTS: Ruptured aneurysms were more likely to have complex flow patterns (83%, P < .001), stable flow patterns (75%, P = .0018), concentrated inflow (66%, P = <.0001), and small impingement regions (76%, P = .0006) compared with unruptured aneurysms. Interobserver variability analyses indicated that all the classifications performed were in very good agreement—that is, well within the 95% CI. CONCLUSIONS: A qualitative hemodynamic analysis of cerebral aneurysms by using image-based patient-specific geometries has shown that concentrated inflow jets, small impingement regions, complex flow patterns, and unstable flow patterns are correlated with a clinical history of prior aneurysm rupture. These qualitative measures provide a starting point for more sophisticated quantitative analysis aimed at assigning aneurysm risk of future rupture. These analyses highlight the potential for CFD to play an important role in the clinical determination of aneurysm risks.

Aneurysm Rupture Following Treatment with Flow-Diverting Stents: Computational Hemodynamics Analysis of Treatment

BACKGROUND AND PURPOSE: Flow-diverting approaches to intracranial aneurysm treatment had many promising early results, but recent apparently successful treatments have been complicated by later aneurysm hemorrhage. We analyzed 7 cases of aneurysms treated with flow diversion to explore the possible rupture mechanisms. MATERIALS AND METHODS: CFD analysis of pre- and posttreatment conditions was performed on 3 giant aneurysms that ruptured after treatment and 4 successfully treated aneurysms. Pre- and posttreatment hemodynamics were compared including WSS, relative blood flows, vascular resistances, and pressures, to identify the effects of flow-diverter placements. RESULTS: Expected reductions in aneurysm velocity and WSS were obtained, indicating effective flow diversion from the sac into the parent artery, consistent with periprocedural observations. In each case with postaneurysm rupture, the result of flow diversion led to an increase in pressure within the aneurysm. This pressure increase is related to larger effective resistance in the parent artery from placement of the devices and, in 2 cases, the reduction of a preaneurysm stenosis. CONCLUSIONS: Flow-diversion devices can cause intra-aneurysmal pressure increases, which can potentially lead to rupture, especially for giant aneurysms. This relates both to changes in the parent artery configuration, such as reduction of a proximal stenosis, and to the flow diversion into higher resistance parent artery pathways combined with cerebral autoregulation, leading to higher pressure gradients. These may be important effects that should be considered when planning interventions. Potentially dangerous cases could be identified with angiography and/or patient-specific CFD models.

Quantitative Characterization of the Hemodynamic Environment in Ruptured and Unruptured Brain Aneurysms

BACKGROUND AND PURPOSE: Hemodynamics are thought to play an important role in the mechanisms of aneurysm pathogenesis, progression, and rupture. The purpose of this study was to define quantitative measures related to qualitative flow characteristics previously analyzed and to investigate their relationship to aneurysm rupture. MATERIALS AND METHODS: The hemodynamic environments in 210 cerebral aneurysms were analyzed by using image-based CFD under different flow conditions. Quantitative hemodynamic variables were defined and extracted from the simulation results. A statistical analysis of the relationship to the previous history of aneurysm rupture was performed, and the variability with flow conditions was assessed. RESULTS: Ruptured aneurysms were more likely to have larger inflow concentrations, larger MWSS, larger shear concentrations, and lower viscous dissipation ratios than unruptured aneurysms. Areas under low WSS and measures of abnormally low shear force distributions of ruptured and unruptured aneurysms were not statistically different. Although the values of hemodynamic quantities changed with different flow conditions, the statistical differences or ratios between their mean values over the ruptured and unruptured groups were maintained, for both pulsatile and steady flows. CONCLUSIONS: Concentrated inflow streams and WSS distributions with elevated levels of MWSS and low aneurysmal viscous dissipation are statistically associated with a clinical history of prior aneurysm rupture. In contrast, the area and total viscous shear force applied in the aneurysm region subjected to abnormally low WSS levels are not. This study highlights the potential for image-based CFD for investigating aneurysm-evolution mechanisms and for clinical assessment of aneurysm risks.

CT angiography in the rapid triage of patients with hyperacute stroke to intraarterial thrombolysis: accuracy in the detection of large vessel thrombus.

Purpose The purpose of this work was to evaluate the accuracy of CT angiography (CTA) for the detection of large vessel intracranial thrombus in clinically suspected hyperacute (<6 h) stroke patients. Method Forty-four consecutive intraarterial thrombolysis candidates underwent noncontrast CT followed immediately by CTA. Axial source and two-dimensional collapsed maximum intensity projection reformatted CTA images were rated for the presence or absence of large vessel occlusion. Five hundred seventy-two circle-of-Willis vessels were reviewed; arteriographic correlation was available for 224 of these. Results Sensitivity and specificity for the detection of large vessel occlusion were 98.4 and 98.1%; accuracy, calculated using receiver operating characteristic analysis, was 99%. Mean time for acquisition, reconstruction, and analysis of CTA images was approximately 15 min. Conclusion CTA is highly accurate for the detection and exclusion of large vessel intracranial occlusion and may therefore be valuable in the rapid triage of hyperacute stroke patients to intraarterial thrombolytic treatment.

Results after Surgical and Endovascular Treatment of Paraclinoid Aneurysms by a Combined Neurovascular Team

OBJECTIVE Advances in surgical and endovascular techniques have improved treatment for paraclinoid aneurysms. A combined surgical and endovascular team can formulate individualized treatment strategies for patients with paraclinoid aneurysms. Patients who are considered to be at high surgical risk can be treated endovascularly to minimize morbidity. We reviewed the clinical and radiographic outcomes of 238 paraclinoid aneurysms treated by our combined surgical and endovascular unit. METHODS From 1991 to 1999, the neurovascular team treated 238 paraclinoid aneurysms in 216 patients at the Massachusetts General Hospital. The modality of treatment for each aneurysm was chosen based on anatomic and clinical risk factors, with endovascular treatment offered to patients considered to have higher surgical risks. One hundred eighty aneurysms were treated by direct surgery, 57 were treated by endovascular occlusion, and one was treated by surgical extracranial-intracranial bypass and endovascular internal carotid artery balloon occlusion. Locations were transitional, 12 (5%); carotid cave, 11 (5%); ophthalmic, 131 (55%); posterior carotid wall, 38 (16%); and superior hypophyseal 46 (19%). Lesions contained completely within the cavernous sinus were excluded from this analysis. RESULTS Using the Glasgow Outcome Scale (GOS), overall clinical outcomes were excellent or good (GOS 5 or 4), 86%; fair (GOS 3), 7%; poor (GOS 2), 4%; and death (GOS 1), 3%. Among the surgically treated patients, 90% experienced excellent or good outcomes (GOS 5 or 4), 6% had fair outcomes (GOS 3), 2% had poor outcomes (GOS 2), and 3% died (GOS 1). Among the endovascularly treated patients, 74% had excellent or good outcomes (GOS 5 or 4), 12% had fair outcomes (GOS 3), 10% had poor outcomes (GOS 2), and 4% died (GOS 1). The overall major and minor complication rate from surgery was 29%, with a 6% surgery-related permanent morbidity rate and a mortality rate of 0%. The overall major and minor complication rate from endovascular treatment was 21%, with a 3% endovascular-related permanent morbidity rate and a 2% mortality rate. Visual outcomes for patients who presented with visual symptoms were as follows: improved, 69%; no change, 25%; worsened, 6%; and new visual deficits, 3%. In general, angiographic efficacy was lower in the endovascular treatment group. CONCLUSION A combined team approach of direct surgery and endovascular coiling can lead to good outcomes in the treatment for paraclinoid aneurysms, including high-risk lesions that might not have been treated in previous surgical series.

Hemodynamic Patterns of Anterior Communicating Artery Aneurysms: A Possible Association with Rupture

BACKGROUND AND PURPOSE: The anterior communicating artery (AcomA) is a predilect location of aneurysms which typically carry higher rupture risks than other locations in the anterior circulation. The purpose of this study was to characterize the different flow types present in AcomA aneurysms and to investigate possible associations with rupture. MATERIALS AND METHODS: Patient-specific computational models of 26 AcomA aneurysms were constructed from 3D rotational angiography images. Bilateral images were acquired in 15 patients who had both A1 segments of the anterior cerebral arteries, and models of the whole anterior circulation were created by fusing the reconstructed left and right arterial trees. Computational fluid dynamics simulations were performed under pulsatile flow conditions measured on a healthy subject. Visualizations of flow velocity, instantaneous streamlines, and wall shear stress (WSS) were performed. These were analyzed for flow patterns, size of the impaction zone, and peak WSS and then correlations were made with prior history of rupture. RESULTS: Aneurysms with small impaction zones were more likely to have ruptured than those with large impaction zones (83% versus 63%). Maximum intra-aneurysmal WSS (MWSS) for the unruptured aneurysms ranged from 10 to 230 dyne/cm2 (mean, 114 dyne/cm2) compared with ruptured aneurysms, which ranged from 35 to 1500 dyne/cm2 (mean, 271 dyne/cm2). This difference in MWSS was statistically significant at 90% confidence levels (P = .10). CONCLUSIONS: Aneurysms with small impaction zones, higher flow rates entering the aneurysm, and elevated MWSS are associated with a clinical history of previous rupture.

Complications of Cerebral Arteriovenous Malformation Embolization: Multivariate Analysis of Predictive Factors

OBJECTIVE:Embolization is an important therapeutic modality in the multidisciplinary management of arteriovenous malformations (AVM); however, prior series have reported a wide variability in overall complication rates caused by embolization (10–50% neurological deficit, 1–4% mortality). In this study, we reviewed our experience with AVM embolization and analyzed factors that might predict complications and clinical outcomes after AVM embolization. METHODS:We analyzed our combined neurovascular unit’s results with AVM embolization from 1993 to 2004 for the following outcomes measures: 1) clinically significant complications, 2) technical complications without clinical sequelae, 3) discharge Glasgow Outcome Scale score, and 4) death. To determine embolization efficacy, we analyzed perioperative blood transfusion and rate of AVM obliteration. Univariate and multivariate analyses were performed for patient age, sex, history of rupture, history of seizure, associated aneurysms, AVM size, deep venous drainage, eloquent location, Spetzler-Martin grade, number of embolization stages, number of pedicles embolized, and primary treatment modality. RESULTS:Over an 11 year period, 295 embolization procedures (761 pedicles embolized) were performed in 168 patients with embolization as the primary treatment modality (n = 16) or as an adjunct to surgery (n = 124) or radiosurgery (n = 28). There were a total of 27 complications in this series, of which 11 were clinically significant (6.5% of patients, 3.7% per procedure), and 16 were technical complications (9.5% of patients, 5.4% per procedure). Excellent or good outcomes (Glasgow Outcome Scale ≥ 4) were observed in 152 (90.5%) patients. Unfavorable outcomes (Glasgow Outcome Scale 1–3) as a direct result of embolization were both 3.0% at discharge and at follow-up, with a 1.2% embolization-related mortality. In the 124 surgical patients, 96.8% had complete AVM obliteration after initial resection, and 31% received perioperative transfusion (mean 1.4 units packed red blood cells per surgical patient). Predictors of unfavorable outcome caused by embolization by univariate analysis were deep venous drainage (P < 0.05), Spetzler-Martin Grade III to V (P < 0.05), and periprocedural hemorrhage (P < 0.0001) and by multivariate analysis were Spetzler-Martin III to V (odds ratio 10.6, P = 0.03) and periprocedural hemorrhage (odds ratio 17, P = 0.004). CONCLUSION:In a single-center, retrospective, nonrandomized study, 90.5% of patients had excellent or good outcomes after AVM embolization, with a complication rate lower than previously reported. Spetzler-Martin grade III to V and periprocedural hemorrhage were the most important predictive factors in determining outcome after embolization.

Hemodynamics and Bleb Formation in Intracranial Aneurysms

BACKGROUND AND PURPOSE: Intracranial aneurysms with irregular shapes and blebs or secondary outpouchings have been correlated with increased rupture risk. The purpose of this study was to investigate possible associations between the local hemodynamics and the formation of blebs in cerebral aneurysms. MATERIALS AND METHODS: Computational models of 20 cerebral aneurysms harboring 30 well-defined blebs were constructed from 3D rotational angiographies. Models representing the aneurysm before bleb formation were constructed by virtually removing the blebs from the anatomic models. Computational fluid dynamics simulations of the aneurysm before and after bleb formation were performed under pulsatile flows. Flow and WSS visualizations were used to analyze the local hemodynamics in the region of the aneurysm that developed the bleb. RESULTS: Most blebs (80%) occurred at or adjacent to the aneurysm region with the highest WSS before bleb formation, and near the flow impaction zone. Most blebs (83%) were found in regions of the aneurysm previously subjected to high or moderate WSS and progressed to low WSS states after the blebs were formed. Most blebs (77%) were aligned or adjacent to the inflow jet, whereas 17% were aligned with the outflow jet, and only 6% were not aligned with the flow direction. In addition, 90% of the aneurysms had maximal WSS higher than or similar to the WSS in the parent artery. CONCLUSIONS: Blebs form at or adjacent to regions of high WSS and are aligned with major intra-aneurysmal flow structures. Formation of blebs results in a lower WSS state with formation of a counter current vortex. These findings imply that locally elevated WSS could contribute to the focalized wall damage that formed these structures.

CFD analysis incorporating the influence of wall motion: application to intracranial aneurysms

Haemodynamics, and in particular wall shear stress, is thought to play a critical role in the progression and rupture of intracranial aneurysms. A novel method is presented that combines image-based wall motion estimation obtained through non-rigid registration with computational fluid dynamics (CFD) simulations in order to provide realistic intra-aneurysmal flow patterns and understand the effects of deforming walls on the haemodynamic patterns. In contrast to previous approaches, which assume rigid walls or ad hoc elastic parameters to perform the CFD simulations, wall compliance has been included in this study through the imposition of measured wall motions. This circumvents the difficulties in estimating personalized elasticity properties. Although variations in the aneurysmal haemodynamics were observed when incorporating the wall motion, the overall characteristics of the wall shear stress distribution do not seem to change considerably. Further experiments with more cases will be required to establish the clinical significance of the observed variations.