Comapping Cellular Content and Extracellular Matrix with Hemodynamics in Intact Arterial Tissues Using Scanning Immunofluorescent Multiphoton Microscopy

Deviation of blood flow from an optimal range is known to be associated with the initiation and progression of vascular pathologies. Important open questions remain about how the abnormal flow drives specific wall changes in pathologies such as cerebral aneurysms where the flow is highly heterogeneous and complex. This knowledge gap precludes the clinical use of readily available flow data to predict outcomes and improve treatment of these diseases. As both flow and the pathological wall changes are spatially heterogeneous, a crucial requirement for progress in this area is a methodology for acquiring and comapping local vascular wall biology data with local hemodynamic data. Here, we developed an imaging pipeline to address this pressing need. A protocol that employs scanning multiphoton microscopy was developed to obtain three-dimensional (3D) datasets for smooth muscle actin, collagen, and elastin in intact vascular specimens. A cluster analysis was introduced to objectively categorize the smooth muscle cells (SMC) across the vascular specimen based on SMC actin density. Finally, direct quantitative comparison of local flow and wall biology in 3D intact specimens was achieved by comapping both heterogeneous SMC data and wall thickness to patient-specific hemodynamic results.

Evaluation of predictive models of aneurysm focal growth and bleb development using machine learning techniques

BACKGROUND

The presence of blebs increases the rupture risk of intracranial aneurysms (IAs).

OBJECTIVE

To evaluate whether cross-sectional bleb formation models can identify aneurysms with focalized enlargement in longitudinal series.

METHODS

Hemodynamic, geometric, and anatomical variables derived from computational fluid dynamics models of 2265 IAs from a cross-sectional dataset were used to train machine learning (ML) models for bleb development. ML algorithms, including logistic regression, random forest, bagging method, support vector machine, and K-nearest neighbors, were validated using an independent cross-sectional dataset of 266 IAs. The models' ability to identify aneurysms with focalized enlargement was evaluated using a separate longitudinal dataset of 174 IAs. Model performance was quantified by the area under the receiving operating characteristic curve (AUC), the sensitivity and specificity, positive predictive value, negative predictive value, F1 score, balanced accuracy, and misclassification error.

RESULTS

The final model, with three hemodynamic and four geometrical variables, along with aneurysm location and morphology, identified strong inflow jets, non-uniform wall shear stress with high peaks, larger sizes, and elongated shapes as indicators of a higher risk of focal growth over time. The logistic regression model demonstrated the best performance on the longitudinal series, achieving an AUC of 0.9, sensitivity of 85%, specificity of 75%, balanced accuracy of 80%, and a misclassification error of 21%.

CONCLUSIONS

Models trained with cross-sectional data can identify aneurysms prone to future focalized growth with good accuracy. These models could potentially be used as early indicators of future risk in clinical practice.

Phenotyping calcification in vascular tissues using artificial intelligence

Vascular calcification is implicated as an important factor in major adverse cardiovascular events (MACE), including heart attack and stroke. A controversy remains over how to integrate the diverse forms of vascular calcification into clinical risk assessment tools. Even the commonly used calcium score for coronary arteries, which assumes risk scales positively with total calcification, has important inconsistencies. Fundamental studies are needed to determine how risk is influenced by the diverse calcification phenotypes. However, studies of these kinds are hindered by the lack of high-throughput, objective, and non-destructive tools for classifying calcification in imaging data sets. Here, we introduce a new classification system for phenotyping calcification along with a semi-automated, non-destructive pipeline that can distinguish these phenotypes in even atherosclerotic tissues. The pipeline includes a deep-learning-based framework for segmenting lipid pools in noisy μ-CT images and an unsupervised clustering framework for categorizing calcification based on size, clustering, and topology. This approach is illustrated for five vascular specimens, providing phenotyping for thousands of calcification particles across as many as 3200 images in less than seven hours. Average Dice Similarity Coefficients of 0.96 and 0.87 could be achieved for tissue and lipid pool, respectively, with training and validation needed on only 13 images despite the high heterogeneity in these tissues. By introducing an efficient and comprehensive approach to phenotyping calcification, this work enables large-scale studies to identify a more reliable indicator of the risk of cardiovascular events, a leading cause of global mortality and morbidity.

Computational fluid dynamics-based virtual angiograms for the detection of flow stagnation in intracranial aneurysms

The goal of this study was to test if CFD-based virtual angiograms could be used to automatically discriminate between intracranial aneurysms (IAs) with and without flow stagnation. Time density curves (TDC) were extracted from patient digital subtraction angiography (DSA) image sequences by computing the average gray level intensity inside the aneurysm region and used to define injection profiles for each subject. Subject-specific 3D models were reconstructed from 3D rotational angiography (3DRA) and computational fluid dynamics (CFD) simulations were performed to simulate the blood flow inside IAs. Transport equations were solved numerically to simulate the dynamics of contrast injection into the parent arteries and IAs and then the contrast retention time (RET) was calculated. The importance of gravitational pooling of contrast agent within the aneurysm was evaluated by modeling contrast agent and blood as a mixture of two fluids with different densities and viscosities. Virtual angiograms can reproduce DSA sequences if the correct injection profile is used. RET can identify aneurysms with significant flow stagnation even when the injection profile is not known. Using a small sample of 14 IAs of which seven were previously classified as having flow stagnation, it was found that a threshold RET value of 0.46 s can successfully identify flow stagnation. CFD-based prediction of stagnation was in more than 90% agreement with independent visual DSA assessment of stagnation in a second sample of 34 IAs. While gravitational pooling prolonged contrast retention time it did not affect the predictive capabilities of RET. CFD-based virtual angiograms can detect flow stagnation in IAs and can be used to automatically identify aneurysms with flow stagnation even without including gravitational effects on contrast agents.

Understanding development of jugular bulb stenosis in vein of galen malformations: identifying metrics of complex flow dynamics in the cerebral venous vasculature of infants

Introduction: Computational fluid dynamics (CFD) assess biological systems based on specific boundary conditions. We propose modeling more advanced hemodynamic metrics, such as core line length (CL) and critical points which characterize complexity of flow in the context of cerebral vasculature, and specifically cerebral veins during the physiologically evolving early neonatal state of vein of Galen malformations (VOGM). CFD has not been applied to the study of arteriovenous shunting in Vein of Galen Malformations but could help illustrate the pathophysiology of this malformation. Methods: Three neonatal patients with VOGM at Boston Children's Hospital met inclusion criteria for this study. Structural MRI data was segmented to generate a mesh of the VOGM and venous outflow. Boundary condition flow velocity was derived from PC-MR sequences with arterial and venous dual velocity encoding. The mesh and boundary conditions were applied to model the cerebral venous flow. We computed flow variables including mean wall shear stress (WSSmean), mean OSI, CL, and the mean number of critical points (nCrPointsmean) for each patient specific model. A critical point is defined as the location where the shear stress vector field is zero (stationary point) and can be used to describe complexity of flow. Results: The division of flow into the left and right venous outflow was comparable between PC-MR and CFD modeling. A high complexity recirculating flow pattern observed on PC-MR was also identified on CFD modeling. Regions of similar WSSmean and OSImean (<1.3 fold) in the left and right venous outflow channels of a single patient have several-fold magnitude difference in higher order hemodynamic metrics (> 3.3 fold CL, > 1.7 fold nCrPointsmean). Specifically, the side which developed JBS in each model had greater nCrPointsmean compared to the jugular bulb with no stenosis (VOGM1: 4.49 vs. 2.53, VOGM2: 1.94 vs. 0, VOGM3: 1 vs. 0). Biologically, these regions had subsequently divergent development, with increased complexity of flow associating with venous stenosis. Discussion: Advanced metrics of flow complexity identified in computational models may reflect observed flow phenomena not fully characterized by primary or secondary hemodynamic parameters. These advanced metrics may indicate physiological states that impact development of jugular bulb stenosis in VOGM.

Differences Between Ruptured Aneurysms With and Without Blebs: Mechanistic Implications

PURPOSE

Blebs are known risk factors for intracranial aneurysm (IA) rupture. We analyzed differences between IAs that ruptured with blebs and those that ruptured without developing blebs to identify distinguishing characteristics among them and suggest possible mechanistic implications.

METHODS

Using image-based models, 25 hemodynamic and geometric parameters were compared between ruptured IAs with and without blebs (n = 673), stratified by location. Hemodynamic and geometric differences between bifurcation and sidewall aneurysms and for aneurysms at five locations were also analyzed.

RESULTS

Ruptured aneurysms harboring blebs were exposed to higher flow conditions than aneurysms that ruptured without developing blebs, and this was consistent across locations. Bifurcation aneurysms were exposed to higher flow conditions than sidewall aneurysms. They had larger maximum wall shear stress (WSS), more concentrated WSS distribution, and larger numbers of critical points than sidewall aneurysms. Additionally, bifurcation aneurysms were larger, more elongated, and had more distorted shapes than sidewall aneurysms. Aneurysm morphology was associated with aneurysm location (p < 0.01). Flow conditions were different between aneurysm locations.

CONCLUSION

Aneurysms at different locations are likely to develop into varying morphologies and thus be exposed to diverse flow conditions that may predispose them to follow distinct pathways towards rupture with or without bleb development. This could explain the diverse rupture rates and bleb presence in aneurysms at different locations.

Flow reversal in distal collaterals as a possible mechanism of delayed intraparenchymal hemorrhage after flow diversion treatment of cerebral aneurysms

Background and Purpose: Delayed intraparenchymal hemorrhages (DIPHs) are one of the most serious complications of cerebral aneurysm treatment with flow diverters (FD), yet their causes are largely unknown. This study analyzes distal hemodynamic alterations induced by the treatment of intracranial aneurysms with FDs.

Methods: A realistic model of the brain arterial network was constructed from MRA images and extended with a constrained constructive optimization technique down to vessel diameters of approximately 50μm. Different variants of the circle of Willis were created by alternatively occluding communicating arteries. Collateral vessels connecting different arterial trees were then added to the model, and a distributed lumped parameter approach was used to model the pulsatile blood flow in the arterial network. The treatment of an ICA aneurysm was modeled by changing the local resistance, flow inertia, and compliance of the aneurysmal segment.

Results: The maximum relative change in distal pressure induced by the aneurysm treatment was below 1%. However, for certain combinations of the circle of Willis and distal collateralization, important flow reversals (with a wall shear stress larger than approximately 1.0 dyne/cm2) were observed in collateral vessels, both ipsilaterally and contralaterally to the treated aneurysm.

Conclusion: This study suggests the hypothesis that flow diverters treatment of intracranial aneurysms could cause important flow reversal in distal collaterals. Flow reversal has previously been shown to be pro-inflammatory and pro-atherogenic and could therefore have a detrimental effect on these collateral vessels, and thus could be a suitable explanation of DIPHs, while the small distal pressure increase is not.

Identification of Small, Regularly Shaped Cerebral Aneurysms Prone to Rupture

BACKGROUND AND PURPOSE

Many small, regularly shaped cerebral aneurysms rupture; however, they usually receive a low score based on current risk-assessment methods. Our goal was to identify patient and aneurysm characteristics associated with rupture of small, regularly shaped aneurysms and to develop and validate predictive models of rupture in this aneurysm subpopulation.

MATERIALS AND METHODS

Cross-sectional data from 1079 aneurysms smaller than 7 mm with regular shapes (without blebs) were used to train predictive models for aneurysm rupture using machine learning methods. These models were based on the patient population, aneurysm location, and hemodynamic and geometric characteristics derived from image-based computational fluid dynamics models. An independent data set with 102 small, regularly shaped aneurysms was used for validation.

RESULTS

Adverse hemodynamic environments characterized by strong, concentrated inflow jets, high speed, complex and unstable flow patterns, and concentrated, oscillatory, and heterogeneous wall shear stress patterns were associated with rupture in small, regularly shaped aneurysms. Additionally, ruptured aneurysms were larger and more elongated than unruptured aneurysms in this subset. A total of 5 hemodynamic and 6 geometric parameters along with aneurysm location, multiplicity, and morphology, were used as predictive variables. The best machine learning rupture prediction-model achieved a good performance with an area under the curve of 0.84 on the external validation data set.

CONCLUSIONS

This study demonstrated the potential of using predictive machine learning models based on aneurysm-specific hemodynamic, geometric, and anatomic characteristics for identifying small, regularly shaped aneurysms prone to rupture.

Evaluation of Outcome Prediction of Flow Diversion for Intracranial Aneurysms

BACKGROUND AND PURPOSE

Identifying and predicting which aneurysms are likely to quickly occlude and which ones are likely to remain open following treatment with flow-diverting devices is important to develop optimal patient management strategies. The purpose of this study was to evaluate predictions based on computational fluid dynamics models using the elastase rabbit aneurysm model.

MATERIALS AND METHODS

A series of 13 aneurysms created in rabbits were treated with flow diverters, and outcomes were angiographically assessed at 8 weeks' follow-up. Computational fluid dynamics models were constructed from pretreatment 3D rotational angiograms and Doppler ultrasound flow velocity measurements. Postimplantation mean aneurysm inflow rate and flow velocity were used to prospectively predict aneurysm occlusion blinded to the actual outcomes. Specifically, if both variables were below their corresponding thresholds, fast occlusion was predicted, while if one of them was above the threshold, slow or incomplete occlusion was predicted.

RESULTS

Of the 13 aneurysms included, 8 were incompletely occluded 8 weeks after treatment, and 5 were completely occluded. A total of 10 computational fluid dynamics-based predictions agreed with the angiographic outcome, reaching 77% accuracy, 80% sensitivity, and 75% specificity. Posttreatment mean velocity alone was able to achieve the same predictive power as the combination of inflow rate and velocity.

CONCLUSIONS

Subject-specific computational fluid dynamics models of the hemodynamic conditions created immediately after implantation of flow-diverting devices in experimental aneurysms created in rabbits are capable of prospectively predicting, with a reasonable accuracy, which aneurysms will completely occlude and which ones will remain incompletely occluded.

Prediction of bleb formation in intracranial aneurysms using machine learning models based on aneurysm hemodynamics, geometry, location, and patient population

BACKGROUND

Bleb presence in intracranial aneurysms (IAs) is a known indication of instability and vulnerability.

OBJECTIVE

To develop and evaluate predictive models of bleb development in IAs based on hemodynamics, geometry, anatomical location, and patient population.

METHODS

Cross-sectional data (one time point) of 2395 IAs were used for training bleb formation models using machine learning (random forest, support vector machine, logistic regression, k-nearest neighbor, and bagging). Aneurysm hemodynamics and geometry were characterized using image-based computational fluid dynamics. A separate dataset with 266 aneurysms was used for model evaluation. Model performance was quantified by the area under the receiving operating characteristic curve (AUC), true positive rate (TPR), false positive rate (FPR), precision, and balanced accuracy.

RESULTS

The final model retained 18 variables, including hemodynamic, geometrical, location, multiplicity, and morphology parameters, and patient population. Generally, strong and concentrated inflow jets, high speed, complex and unstable flow patterns, and concentrated, oscillatory, and heterogeneous wall shear stress patterns together with larger, more elongated, and more distorted shapes were associated with bleb formation. The best performance on the validation set was achieved by the random forest model (AUC=0.82, TPR=91%, FPR=36%, misclassification error=27%).

CONCLUSIONS

Based on the premise that aneurysm characteristics prior to bleb formation resemble those derived from vascular reconstructions with their blebs virtually removed, machine learning models can identify aneurysms prone to bleb development with good accuracy. Pending further validation with longitudinal data, these models may prove valuable for assessing the propensity of IAs to progress to vulnerable states and potentially rupturing.

Quantification of the Rupture Potential of Patient-Specific Intracranial Aneurysms under Contact Constraints

Intracranial aneurysms (IAs) are localized enlargements of cerebral blood vessels that cause substantial rates of mortality and morbidity in humans. The rupture possibility of these aneurysms is a critical medical challenge for physicians during treatment planning. This treatment planning while assessing the rupture potential of aneurysms becomes more complicated when they are constrained by an adjacent structure such as optic nerve tissues or bones, which is not widely studied yet. In this work, we considered and studied a constitutive model to investigate the bio-mechanical response of image-based patient-specific IA data using cardiovascular structural mechanics equations. We performed biomechanical modeling and simulations of four different patient-specific aneurysms’ data (three middle cerebral arteries and one internal carotid artery) to assess the rupture potential of those aneurysms under a plane contact constraint. Our results suggest that aneurysms with plane contact constraints produce less or almost similar maximum wall effective stress compared to aneurysms with no contact constraints. In our research findings, we observed that a plane contact constraint on top of an internal carotid artery might work as a protective wall due to the 16.6% reduction in maximum wall effective stress than that for the case where there is no contact on top of the aneurysm.

Analysis of hemodynamic changes from aneurysm inception to large sizes

While previous studies have identified many risk factors for the progression and rupture of cerebral aneurysms, the changes in aneurysm flow characteristics during its evolution are not fully understood. This work analyzes the changes in the aneurysm hemodynamic environment from its initial development to later stages when the aneurysm has substantially enlarged. A total of 88 aneurysms at four locations were studied with image based computational fluid dynamics (CFD). Two synthetic sequences representing the aneurysm geometry at three earlier stages were generated by shrinking the aneurysm sac while keeping the neck fixed or shrinking the neck simultaneously. The flow conditions were then quantitatively compared between these two modes of evolution. As aneurysms enlarged, the inflow rate increased in growing neck sequences, but decreased in fixed neck sequences. The inflow jet became more concentrated in both sequences. The mean aneurysm flow velocity and wall shear stress decreased in both sequences, but they decreased faster in enlarging aneurysms if the neck was fixed. Additionally, the intra-aneurysmal flows became more complex and more unstable, wall shear stress distribution became more oscillatory, and the area under low wall shear stress increased for both sequences. The evolution of flow characteristics of aneurysms with fixed and growing necks are different. The observed trends suggest that fixed neck aneurysms may evolve towards a flow environment characteristic of stable aneurysms faster than aneurysms with growing necks, which could also evolve towards a more disfavorable environment.

Regional Aneurysm Wall Enhancement is Affected by Local Hemodynamics: A 7T MRI Study

BACKGROUND AND PURPOSE

Aneurysm wall enhancement has been proposed as a biomarker for inflammation and instability. However, the mechanisms of aneurysm wall enhancement remain unclear. We used 7T MR imaging to determine the effect of flow in different regions of the wall.

MATERIALS AND METHODS

Twenty-three intracranial aneurysms imaged with 7T MR imaging and 3D angiography were studied with computational fluid dynamics. Local flow conditions were compared between aneurysm wall enhancement and nonenhanced regions. Aneurysm wall enhancement regions were subdivided according to their location on the aneurysm and relative to the inflow and were further compared.

RESULTS

On average, wall shear stress was lower in enhanced than in nonenhanced regions (P = .05). Aneurysm wall enhancement regions at the neck had higher wall shear stress gradients (P = .05) with lower oscillations (P = .05) than nonenhanced regions. In contrast, aneurysm wall enhancement regions at the aneurysm body had lower wall shear stress (P = .01) and wall shear stress gradients (P = .008) than nonenhanced regions. Aneurysm wall enhancement regions far from the inflow had lower wall shear stress (P = .006) than nonenhanced regions, while aneurysm wall enhancement regions close to the inflow tended to have higher wall shear stress than the nonenhanced regions, but this association was not significant.

CONCLUSIONS

Aneurysm wall enhancement regions tend to have lower wall shear stress than nonenhanced regions of the same aneurysm. Moreover, the association between flow conditions and aneurysm wall enhancement seems to depend on the location of the region on the aneurysm sac. Regions at the neck and close to the inflow tend to be exposed to higher wall shear stress and wall shear stress gradients. Regions at the body, dome, or far from the inflow tend to be exposed to uniformly low wall shear stress and have more aneurysm wall enhancement.

Hemodynamics in aneurysm blebs with different wall characteristics

BACKGROUND

Blebs are important secondary structures of intracranial aneurysms associated with increased rupture risk and can affect local wall stress and hemodynamics. Mechanisms of bleb development and evolution are not clearly understood. We investigate the relationship between blebs with different wall characteristics and local hemodynamics and rupture sites.

METHODS

Blebs with different wall appearances in intra-operative videos were analyzed with image-based computational fluid dynamics. Thin red blebs were compared against thick atherosclerotic/hyperplastic white/yellow blebs. Rupture points were identified in videos of ruptured aneurysms harboring blebs.

RESULTS

Thin blebs tended to be closer to the inflow than atherosclerotic blebs of the same aneurysm (P=0.0234). Blebs near the inflow had higher velocity (P=0.0213), vorticity (P=0.0057), shear strain rate (P=0.0084), wall shear stress (WSS) (P=0.0085), and WSS gradient (P=0.0151) than blebs far from the inflow. In a subset of 12 ruptured aneurysms harboring blebs, rupture points were associated with thin blebs in 42% of aneurysms, atherosclerotic blebs in 25%, and were away from blebs in the remaining 33%.

CONCLUSIONS

Not all blebs are equal; some have thin translucent walls while others have thick atherosclerotic walls. Thin blebs tend to be located closer to the inflow than atherosclerotic blebs. Blebs near the inflow are exposed to stronger flows with higher and spatially variable WSS than blebs far from the inflow which tend to have uniformly lower WSS. Aneurysms can rupture at thin blebs, atherosclerotic blebs, and even away from blebs. Further study of wall failure in aneurysms with different bleb types is needed.

Extending statistical learning for aneurysm rupture assessment to Finnish and Japanese populations using morphology, hemodynamics, and patient characteristics

OBJECTIVE

Incidental aneurysms pose a challenge for physicians, who need to weigh the rupture risk against the risks associated with treatment and its complications. A statistical model could potentially support such treatment decisions. A recently developed aneurysm rupture probability model performed well in the US data used for model training and in data from two European cohorts for external validation. Because Japanese and Finnish patients are known to have a higher aneurysm rupture risk, the authors' goals in the present study were to evaluate this model using data from Japanese and Finnish patients and to compare it with new models trained with Finnish and Japanese data.

METHODS

Patient and image data on 2129 aneurysms in 1472 patients were used. Of these aneurysm cases, 1631 had been collected mainly from US hospitals, 249 from European (other than Finnish) hospitals, 147 from Japanese hospitals, and 102 from Finnish hospitals. Computational fluid dynamics simulations and shape analyses were conducted to quantitatively characterize each aneurysm's shape and hemodynamics. Next, the previously developed model's discrimination was evaluated using the Finnish and Japanese data in terms of the area under the receiver operating characteristic curve (AUC). Models with and without interaction terms between patient population and aneurysm characteristics were trained and evaluated including data from all four cohorts obtained by repeatedly randomly splitting the data into training and test data.

RESULTS

The US model's AUC was reduced to 0.70 and 0.72, respectively, in the Finnish and Japanese data compared to 0.82 and 0.86 in the European and US data. When training the model with Japanese and Finnish data, the average AUC increased only slightly for the Finnish sample (to 0.76 ± 0.16) and Finnish and Japanese cases combined (from 0.74 to 0.75 ± 0.14) and decreased for the Japanese data (to 0.66 ± 0.33). In models including interaction terms, the AUC in the Finnish and Japanese data combined increased significantly to 0.83 ± 0.10.

CONCLUSIONS

Developing an aneurysm rupture prediction model that applies to Japanese and Finnish aneurysms requires including data from these two cohorts for model training, as well as interaction terms between patient population and the other variables in the model. When including this information, the performance of such a model with Japanese and Finnish data is close to its performance with US or European data. These results suggest that population-specific differences determine how hemodynamics and shape associate with rupture risk in intracranial aneurysms.

Hemodynamic conditions that favor bleb formation in cerebral aneurysms

BACKGROUND

Although it is generally believed that blebs represent weaker spots in the walls of intracranial aneurysms (IAs), it is largely unknown which aneurysm characteristics favor their development.

OBJECTIVE

To investigate possible associations between aneurysm hemodynamic and geometric characteristics and the development of blebs in intracranial aneurysms.

METHODS

A total of 270 IAs in 199 patients selected for surgical clipping were studied. Blebs were visually identified and interactively marked on patient-specific vascular models constructed from presurgical images. Blebs were then deleted from the vascular reconstruction to approximate the aneurysm before bleb formation. Computational fluid dynamics studies were performed in these models and in cases without blebs. Hemodynamic and geometric characteristics of aneurysms with and without blebs were compared.

RESULTS

A total of 173 aneurysms had no blebs, while 97 aneurysms had a total of 122 blebs. Aneurysms favoring bleb formation had stronger (p<0.0001) and more concentrated inflow jets (p<0.0001), higher flow velocity (p=0.0061), more complex (p<0.0001) and unstable (p=0.0157) flow patterns, larger maximum wall shear stress (WSS; p<0.0001), more concentrated (p=0.0005) and oscillatory (p=0.0004) WSS distribution, and a more heterogeneous WSS field (p<0.0001), than aneurysms without blebs. They were also larger (p<0.0001), more elongated (p<0.0001), had wider necks (p=0.0002), and more distorted and irregular shapes (p<0.0001).

CONCLUSIONS

Strong and concentrated inflow jets, high-speed, complex, and unstable flow patterns, and concentrated, oscillatory, and heterogeneous WSS patterns favor the formation of blebs in IAs. Blebs are more likely to form in large, elongated, and irregularly shaped aneurysms. These adverse characteristics could be considered signs of aneurysm instability when evaluating aneurysms for conservative observation or treatment.

Blebs in intracranial aneurysms: prevalence and general characteristics

BACKGROUND

Blebs are rupture risk factors in intracranial aneurysms (IAs), but their prevalence, distribution, and associations with clinical factors as well as their causes and effects on aneurysm vulnerability remain unclear.

METHODS

A total of 122 blebs in 270 IAs selected for surgery were studied using patient-specific vascular reconstructions from 3D angiographic images. Bleb geometry, location on the aneurysm, and frequency of occurrence in aneurysms at different locations were analyzed. Associations between gender, age, smoking, hypertension, hormone therapy, dental infection, and presence of blebs were investigated.

RESULTS

Of all aneurysms with blebs, 77% had a single bleb and 23% had multiple blebs. Only 6% of blebs were at the neck, while 46% were in the body and 48% in the dome. Aneurysms with blebs were larger (p<0.0001), more elongated (p=0.0002), and with wider necks than aneurysms without blebs. Bleb presence was associated with dental infection (p=0.0426) and negatively associated with hormone therapy (p=0.0426) in women. Anterior and posterior communicating arteries had larger percentages of aneurysms with blebs than internal carotid arteries. Patients with a history of hypertension tended to have a larger percentage of aneurysms with blebs. However, these trends did not reach significance in this sample.

CONCLUSIONS

Blebs are common in IAs, and most aneurysms harboring blebs have a single bleb. Blebs in the aneurysm neck are rare, but they are equally common in the body and dome. The presence of blebs in IAs was associated with dental infection, and negatively associated with hormone replacement therapy.

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.

Comparing Morphology and Hemodynamics of Stable-versus-Growing and Grown Intracranial Aneurysms

BACKGROUND AND PURPOSE

Aneurysm growth has been related to higher rupture risk. A better understanding of the characteristics related to growth may assist in the treatment decisions of unruptured intracranial aneurysms. This study aimed to identify morphologic and hemodynamic characteristics associated with aneurysm growth and to determine whether these characteristics deviate further from those of stable aneurysms after growth.

MATERIALS AND METHODS

We included 81 stable and 56 growing aneurysms. 3D vascular models were segmented on CTA, MRA, or 3D rotational angiographic images. With these models, we performed computational fluid dynamics simulations. Morphologic (size, size ratios, and shape) and hemodynamic (inflow, vorticity, shear stress, oscillatory shear index, flow instability) characteristics were automatically calculated. We compared the characteristics between aneurysms that were stable and those that had grown at baseline and final imaging. The significance level after Bonferroni correction was P < .002.

RESULTS

At baseline, no significant differences between aneurysms that were stable and those that had grown were detected (P > .002). Significant differences between aneurysms that were stable and those that had grown were seen at the final imaging for shear rate, aneurysm velocity, vorticity, and mean wall shear stress (P < .002). The latter was 11.5 (interquartile range, 5.4-18.8 dyne/cm²) compared with 17.5 (interquartile range, 11.2-29.9 dyne/cm²) in stable aneurysms (P = .001). Additionally, a trend toward lower area weighted average Gaussian curvature in aneurysms that had grown was observed with a median of 6.0 (interquartile range, 3.2-10.7 cm^-2) compared with 10.4 (interquartile range, 5.0-21.2 cm^-2) in stable aneurysms (P = .004).

CONCLUSIONS

Morphologic and hemodynamic characteristics at baseline were not associated with aneurysm growth in our population. After growth, almost all indices increase toward values associated with higher rupture risks. Therefore, we stress the importance of longitudinal imaging and repeat risk assessment in unruptured aneurysms.

Comparison of statistical learning approaches for cerebral aneurysm rupture assessment

PURPOSE

Incidental aneurysms pose a challenge to physicians who need to decide whether or not to treat them. A statistical model could potentially support such treatment decisions. The aim of this study was to compare a previously developed aneurysm rupture logistic regression probability model (LRM) to other machine learning (ML) classifiers for discrimination of aneurysm rupture status.

METHODS

Hemodynamic, morphological, and patient-related information of 1631 cerebral aneurysms characterized by computational fluid dynamics simulations were used to train support vector machines (SVMs) with linear and RBF kernel (RBF-SVM), k-nearest neighbors (kNN), decision tree, random forest, and multilayer perceptron (MLP) neural network classifiers for predicting the aneurysm rupture status. The classifiers' accuracy, sensitivity, specificity, and area under the receiver operating characteristic curve (AUC) were evaluated and compared to the LRM using 249 test cases obtained from two external cohorts. Additionally, important variables were determined based on the random forest and weights of the linear SVM.

RESULTS

The AUCs of the MLP, LRM, linear SVM, RBF-SVM, kNN, decision tree, and random forest were 0.83, 0.82, 0.80, 0.81, 0.76, 0.70, and 0.79, respectively. The accuracy ranged between 0.76 (decision tree,) and 0.79 (linear SVM, RBF-SVM, and MLP). Important variables for predicting the aneurysm rupture status included aneurysm location, the mean surface curvature, and maximum flow velocity.

CONCLUSION

The performance of the LRM was overall comparable to that of the other ML classifiers, confirming its potential for aneurysm rupture assessment. To further improve the predictions, additional information, e.g., related to the aneurysm wall, might be needed.

Calcification in Human Intracranial Aneurysms Is Highly Prevalent and Displays Both Atherosclerotic and Nonatherosclerotic Types

OBJECTIVE

Although the clinical and biological importance of calcification is well recognized for the extracerebral vasculature, its role in cerebral vascular disease, particularly, intracranial aneurysms (IAs), remains poorly understood. Extracerebrally, 2 distinct mechanisms drive calcification, a nonatherosclerotic, rapid mineralization in the media and a slower, inflammation driven, atherosclerotic mechanism in the intima. This study aims to determine the prevalence, distribution, and type (atherosclerotic, nonatherosclerotic) of calcification in IAs and assess differences in occurrence between ruptured and unruptured IAs. Approach and Results: Sixty-five 65 IA specimens (48 unruptured, 17 ruptured) were resected perioperatively. Calcification and lipid pools were analyzed nondestructively in intact samples using high resolution (0.35 μm) microcomputed tomography. Calcification is highly prevalent (78%) appearing as micro (<500 µm), meso (500 µm-1 mm), and macro (>1 mm) calcifications. Calcification manifests in IAs as both nonatherosclerotic (calcification distinct from lipid pools) and atherosclerotic (calcification in the presence of lipid pools) with 3 wall types: Type I-only calcification, no lipid pools (20/51, 39%), Type II-calcification and lipid pools, not colocalized (19/51, 37%), Type III-calcification colocalized with lipid pools (12/51, 24%). Ruptured IAs either had no calcifications or had nonatherosclerotic micro- or meso-calcifications (Type I or II), without macro-calcifications.

CONCLUSIONS

Calcification in IAs is substantially more prevalent than previously reported and presents as both nonatherosclerotic and atherosclerotic types. Notably, ruptured aneurysms had only nonatherosclerotic calcification, had significantly lower calcification fraction, and did not contain macrocalcifications. Improved understanding of the role of calcification in IA pathology should lead to new therapeutic targets.

Analysis of Flow Dynamics and Outcomes of Cerebral Aneurysms Treated with Intrasaccular Flow-Diverting Devices

BACKGROUND AND PURPOSE

Intrasaccular flow diversion offers a promising treatment option for complex bifurcation aneurysms. The purpose of this study was to compare the flow conditions between successfully occluded and incompletely occluded aneurysms treated with intrasaccular devices.

MATERIALS AND METHODS

The hemodynamics in 18 completely occluded aneurysms after treatment with intrasaccular devices was compared against 18 that were incompletely occluded at follow-up. Hemodynamic and geometric parameters were obtained from computational fluid dynamics models constructed from 3D angiographies. Models of the intrasaccular devices were created and interactively deployed within the vascular models using posttreatment angiography images for guidance. Hemodynamic and geometric variables were compared using the Mann-Whitney test and univariate logistic regression analysis.

RESULTS

Incomplete occlusion was associated with large posttreatment mean aneurysm inflows (P = .02) and small reductions in the mean inflow rate (P = .01) and inflow concentration index (P = .03). Incompletely occluded aneurysms were larger (P = .002) and had wider necks (P = .004) than completely occluded aneurysms and tended to have more complex flow patterns, though this trend was not significant after adjusting for multiple testing.

CONCLUSIONS

The outcome of cerebral aneurysm treatment with intrasaccular flow diverters is associated with flow conditions created immediately after device implantation. Flow conditions unfavorable for immediate and complete occlusion seem to be created by improper positioning or orientation of the device. Complete occlusion is more difficult to achieve in larger aneurysms, aneurysms with wider necks, and aneurysms with stronger and more complex flows.

Image-based modeling of blood flow in cerebral aneurysms treated with intrasaccular flow diverting devices

Modeling the flow dynamics in cerebral aneurysms after the implantation of intrasaccular devices is important for understanding the relationship between flow conditions created immediately posttreatment and the subsequent outcomes. This information, ideally available a priori based on computational modeling prior to implantation, is valuable to identify which aneurysms will occlude immediately and which aneurysms will likely remain patent and would benefit from a different procedure or device. In this report, a methodology for modeling the hemodynamics in intracranial aneurysms treated with intrasaccular flow diverting devices is described. This approach combines an image-guided, virtual device deployment within patient-specific vascular models with an immersed boundary method on adaptive unstructured grids. A partial mesh refinement strategy that reduces the number of mesh elements near the aneurysm dome where the flow conditions are largely stagnant was compared with the full refinement strategy that refines the mesh everywhere around the device wires. The results indicate that using the partial mesh refinement approach is adequate for analyzing the posttreatment hemodynamics, at a reduced computational cost. The results obtained on a series of four cerebral aneurysms treated with different intrasaccular devices were in good qualitative agreement with angiographic observations. Promising results were obtained relating posttreatment flow conditions and outcomes of treatments with intrasaccular devices, which need to be confirmed on larger series.

Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH)-phase II: rupture risk assessment

PURPOSE

Assessing the rupture probability of intracranial aneurysms (IAs) remains challenging. Therefore, hemodynamic simulations are increasingly applied toward supporting physicians during treatment planning. However, due to several assumptions, the clinical acceptance of these methods remains limited.

METHODS

To provide an overview of state-of-the-art blood flow simulation capabilities, the Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH) was conducted. Seventeen research groups from all over the world performed segmentations and hemodynamic simulations to identify the ruptured aneurysm in a patient harboring five IAs. Although simulation setups revealed good similarity, clear differences exist with respect to the analysis of aneurysm shape and blood flow results. Most groups (12/71%) included morphological and hemodynamic parameters in their analysis, with aspect ratio and wall shear stress as the most popular candidates, respectively.

RESULTS

The majority of groups (7/41%) selected the largest aneurysm as being the ruptured one. Four (24%) of the participating groups were able to correctly select the ruptured aneurysm, while three groups (18%) ranked the ruptured aneurysm as the second most probable. Successful selections were based on the integration of clinically relevant information such as the aneurysm site, as well as advanced rupture probability models considering multiple parameters. Additionally, flow characteristics such as the quantification of inflow jets and the identification of multiple vortices led to correct predictions.

CONCLUSIONS

MATCH compares state-of-the-art image-based blood flow simulation approaches to assess the rupture risk of IAs. Furthermore, this challenge highlights the importance of multivariate analyses by combining clinically relevant metadata with advanced morphological and hemodynamic quantification.

Local Hemodynamic Conditions Associated with Focal Changes in the Intracranial Aneurysm Wall

BACKGROUND AND PURPOSE

Aneurysm hemodynamics has been associated with wall histology and inflammation. We investigated associations between local hemodynamics and focal wall changes visible intraoperatively.

MATERIALS AND METHODS

Computational fluid dynamics models were constructed from 3D images of 65 aneurysms treated surgically. Aneurysm regions with different visual appearances were identified in intraoperative videos: 1) "atherosclerotic" (yellow), 2) "hyperplastic" (white), 3) "thin" (red), 4) rupture site, and 5) "normal" (similar to parent artery), They were marked on 3D reconstructions. Regional hemodynamics was characterized by the following: wall shear stress, oscillatory shear index, relative residence time, wall shear stress gradient and divergence, gradient oscillatory number, and dynamic pressure; these were compared using the Mann-Whitney test.

RESULTS

Hyperplastic regions had lower average wall shear stress (P = .005) and pressure (P = .009) than normal regions. Flow conditions in atherosclerotic and hyperplastic regions were similar but had higher average relative residence time (P = .03) and oscillatory shear index (P = .04) than thin regions. Hyperplastic regions also had a higher average gradient oscillatory number (P = .002) than thin regions. Thin regions had lower average relative residence time (P < .001), oscillatory shear index (P = .006), and gradient oscillatory number (P < .001) than normal regions, and higher average wall shear stress (P = .006) and pressure (P = .009) than hyperplastic regions. Thin regions tended to be aligned with the flow stream, while atherosclerotic and hyperplastic regions tended to be aligned with recirculation zones.

CONCLUSIONS

Local hemodynamics is associated with visible focal wall changes. Slow swirling flow with low and oscillatory wall shear stress was associated with atherosclerotic and hyperplastic changes. High flow conditions prevalent in regions near the flow impingement site characterized by higher and less oscillatory wall shear stress were associated with local "thinning" of the wall.

External validation of cerebral aneurysm rupture probability model with data from two patient cohorts

BACKGROUND

For a treatment decision of unruptured cerebral aneurysms, physicians and patients need to weigh the risk of treatment against the risk of hemorrhagic stroke caused by aneurysm rupture. The aim of this study was to externally evaluate a recently developed statistical aneurysm rupture probability model, which could potentially support such treatment decisions.

METHODS

Segmented image data and patient information obtained from two patient cohorts including 203 patients with 249 aneurysms were used for patient-specific computational fluid dynamics simulations and subsequent evaluation of the statistical model in terms of accuracy, discrimination, and goodness of fit. The model's performance was further compared to a similarity-based approach for rupture assessment by identifying aneurysms in the training cohort that were similar in terms of hemodynamics and shape compared to a given aneurysm from the external cohorts.

RESULTS

When applied to the external data, the model achieved a good discrimination and goodness of fit (area under the receiver operating characteristic curve AUC = 0.82), which was only slightly reduced compared to the optimism-corrected AUC in the training population (AUC = 0.84). The accuracy metrics indicated a small decrease in accuracy compared to the training data (misclassification error of 0.24 vs. 0.21). The model's prediction accuracy was improved when combined with the similarity approach (misclassification error of 0.14).

CONCLUSIONS

The model's performance measures indicated a good generalizability for data acquired at different clinical institutions. Combining the model-based and similarity-based approach could further improve the assessment and interpretation of new cases, demonstrating its potential use for clinical risk assessment.

Associations of hemodynamics, morphology, and patient characteristics with aneurysm rupture stratified by aneurysm location

PURPOSE

The mechanisms of cerebral aneurysm rupture are not fully understood. We analyzed the associations of hemodynamics, morphology, and patient age and gender with aneurysm rupture stratifying by location.

METHODS

Using image-based models, 20 hemodynamic and 17 morphological parameters were compared in 1931 ruptured and unruptured aneurysms with univariate logistic regression. Rupture rates were compared between males and females as well as younger and older patients and bifurcation versus sidewall aneurysms for different aneurysm locations. Subsequently, associations between hemodynamics and morphology and patient as well as aneurysm characteristics were analyzed for aneurysms at five locations.

RESULTS

Compared to unruptured aneurysms, ruptured aneurysms were characterized by a more irregular shape and were exposed to a more adverse hemodynamic environment described by faster flow, higher wall shear stress, more oscillatory shear, and more unstable and complex flows. These associations with rupture status were consistent for different aneurysm locations. Rupture rates were significantly higher in males at the internal carotid artery (ICA) bifurcation, ophthalmic ICA, and the middle cerebral artery (MCA) bifurcation. At the anterior communicating artery (ACOM) and MCA bifurcation, they were significantly higher for younger patients. Bifurcation aneurysms had significantly larger rupture rates at the MCA and posterior communicating artery (PCOM). In these groups with higher rupture rates, aneurysms were characterized by adverse hemodynamics and more complex shapes.

CONCLUSION

Hemodynamic and morphological differences between ruptured and unruptured aneurysms are consistent across locations. Adverse morphology and hemodynamics are related to rupture as well as younger age, male gender, and bifurcation aneurysms.

Downstream vascular changes after flow-diverting device deployment in a rabbit model

BACKGROUND

Flow diverters (FDs) are increasingly used in the treatment of intracranial aneurysms, and carry the risk of thromboembolic complications, even in patients treated with dual antiplatelet therapy. The effect of FDs on the downstream vascular is unknown. The aim of the study was to investigate vascular wall pulse wave velocity (PWV) and contractility changes following FD treatment in a rabbit model.

METHODS

FDs (Pipeline Embolic Device, Medtronic Inc., Irvine, California, USA) were implanted in the aorta of normal rabbits and sham-operated aorta were used as controls (n=6 per group). Pulse wave imaging with ultra-fast ultrasound at 1600 frames per second (Vantage, Verasonics, Inc., Kirkland, WA) was performed in the vessel wall distal to FD prior to device implantation and at 8- week follow-up to measure the PWV. Force contraction vascular reactivity studies were conducted in the aortic rings using an organ bath.

RESULTS

The difference in mean PWV in the follow-up compared with pre-implantation was significantly higher in the distal vessels compared with sham controls (1.18 m/s [SD=0.54] vs. 0.37 m/s [SD=1.09], P=0.03). Conversely, the aortic segments distal to the FD exhibited a 55% increase in vascular contractility compared with proximal segments (P=0.002). We observed a significant positive correlation between mean PWV and mean vascular contractility.

CONCLUSION

Implantation of FD was associated with increased PWV and vascular contractility, suggesting that FD implantation causes changes to the vascular wall. Further studies are needed to understand the clinical implication of changes in vascular PWV and contractility.

Combining data from multiple sources to study mechanisms of aneurysm disease: Tools and techniques

INTRODUCTION

Connecting local hemodynamics, biomechanics, and tissue properties in cerebral aneurysms is important for understanding the processes of wall degeneration and subsequent aneurysm progression and rupture. This challenging problem requires integration of data from multiple sources.

METHODS

This paper describes the tools and techniques developed to integrate data from multiple sources, including clinical information, 3D imaging, intraoperative videos, ex vivo micro-computed tomography (CT), and multiphoton microscopy. Central to this approach is a 3D tissue model constructed from micro-CT images of aneurysm samples resected during neurosurgery. This model is aligned to vascular models constructed from 3D clinical images and is used to map and compare flow, biomechanics, and tissue data.

RESULTS

The approach is illustrated with data of three human intracranial aneurysms. These case studies demonstrated the ability of this approach to study relationships between different factors affecting the aneurysm wall and produced provocative observations that will be further studied with larger series. For instance, "atherosclerotic" and "hyperplastic" looking parts of the aneurysm corresponded to thicker walls and occurred in regions of recirculating flow and low wall shear stress (WSS); thin regions were associated with inflow jets, flow impingement, and high WSS; blebs had walls of varying structures, including calcified, thin, or hyperplastic walls.

CONCLUSIONS

The current approach enables the study of interactions of multiple factors thought to be responsible for the progressive degradation and weakening of the aneurysm wall during its evolution.

Identification of Hostile Hemodynamics and Geometries of Cerebral Aneurysms: A Case-Control Study

BACKGROUND AND PURPOSE

Hostile hemodynamic conditions and geometries are thought to predispose aneurysms for instability and rupture. This study compares stable, unstable, and ruptured aneurysms while controlling for location and patient characteristics.

MATERIALS AND METHODS

The hemodynamics and geometries of 165 stable, 65 unstable, and 554 ruptured aneurysms were compared. Hemodynamics was modeled using image-based computational fluid dynamics. Case-control pairs were selected matching aneurysm location, patient age, and sex. Paired Wilcoxon tests were used to compare hemodynamic and geometric variables among different aneurysm groups. The pairing was repeated 100 times, and the combined P values were calculated and adjusted for multiple testing.

RESULTS

Ruptured aneurysms had lower minimum wall shear stress (P = .03), higher maximum wall shear stress (P = .03), more concentrated (P = .03) and mean oscillatory shear stress (P = .03), higher maximum velocity (P = .03), and more complex flows (vortex core-line length, P = .03) than stable aneurysms. Similarly, unstable aneurysms had more concentrated shear stress (P = .04) and more complex flows (vortex core-line length, P = .04) than stable aneurysms. Compared with stable aneurysms, ruptured aneurysms were larger (size ratio, aneurysm size/vessel size, P = .03), more elongated (aspect ratio, P = .03), and irregular (nonsphericity index, P = .03). Similarly, unstable aneurysms were larger (size ratio, P = .04), more elongated (aspect ratio, P = .04), and irregular (bulge location, P = .04; area-weighted Gaussian curvature; P = .04) than stable aneurysms. No significant differences were found between unstable and ruptured aneurysms.

CONCLUSIONS

Unstable and ruptured aneurysms have more complex flows with concentrated wall shear stress and are larger, more elongated, and irregular than stable aneurysms, independent of aneurysm location and patient sex and age.

Development and internal validation of an aneurysm rupture probability model based on patient characteristics and aneurysm location, morphology, and hemodynamics

PURPOSE

Unruptured cerebral aneurysms pose a dilemma for physicians who need to weigh the risk of a devastating subarachnoid hemorrhage against the risk of surgery or endovascular treatment and their complications when deciding on a treatment strategy. A prediction model could potentially support such treatment decisions. The aim of this study was to develop and internally validate a model for aneurysm rupture based on hemodynamic and geometric parameters, aneurysm location, and patient gender and age.

METHODS

Cross-sectional data from 1061 patients were used for image-based computational fluid dynamics and shape characterization of 1631 aneurysms for training an aneurysm rupture probability model using logistic group Lasso regression. The model's discrimination and calibration were internally validated based on the area under the curve (AUC) of the receiver operating characteristic and calibration plots.

RESULTS

The final model retained 11 hemodynamic and 12 morphological variables, aneurysm location, as well as patient age and gender. An adverse hemodynamic environment characterized by a higher maximum oscillatory shear index, higher kinetic energy and smaller low shear area as well as a more complex aneurysm shape, male gender and younger age were associated with an increased rupture risk. The corresponding AUC of the model was 0.86 (95% CI [0.85, 0.86], after correction for optimism 0.84).

CONCLUSION

The model combining variables from various domains was able to discriminate between ruptured and unruptured aneurysms with an AUC of 86%. Internal validation indicated potential for the application of this model in clinical practice after evaluation with longitudinal data.

Hemodynamic characteristics of stable and unstable vertebrobasilar dolichoectatic and fusiform aneurysms

Background and purpose

Vertebrobasilar dolichoectatic and fusiform aneurysms (VBDAs) are known to have a poor natural history, with high rates of growth, rupture, and stroke. The purpose of this study was to identify hemodynamic characteristics that differ between VBDAs associated with growth, rupture, and stroke.

Materials and methods

VBDAs with CT angiography or MR angiography followed longitudinally without treatment were studied. Unstable aneurysms were defined as those that grew or ruptured during follow-up. Aneurysms associated with stroke were defined as those associated with posterior circulation infarct at follow-up. Baseline data, including demographics, comorbidities, and aneurysm morphology and size were collected. Image based computational fluid dynamics models were created and run under pulsatile flow conditions. Relevant hemodynamic and geometric variables were calculated and compared between groups (stable vs unstable and no stroke vs stroke) using the Wilcoxon test.

Results

A total of 37 VBDAs were included (24 stable, 13 unstable; 30 no stroke, 7 stroke). Unstable aneurysms had lower shear rates (P=0.05), blood flow velocity (P=0.03), and lower vorticity (P=0.049) than stable aneurysms. In addition, unstable aneurysms had higher mean oscillatory shear indices (P=0.001). There were no differences in the hemodynamic characteristics of aneurysms in the stroke group compared with the non-stroke group.

Conclusion

This small study suggests there may be hemodynamic differences between unstable and stable VBDAs. Unstable VBDAs appear to be under lower flow conditions with lower velocity, vorticity, and shear rates, and have more oscillatory flow. There was no difference in the hemodynamic characteristics of aneurysms in the stroke and no stroke group.

A data-driven approach for addressing the lack of flow waveform data in studies of cerebral arterial flow in older adults

OBJECTIVE

Blood flow waveforms-essential data for hemodynamic modeling-are often in practice unavailable to researchers. The objectives of this work were to assess the variability among the waveforms for a clinically relevant older population, and develop data-based methods for addressing the missing waveform data for hemodynamic studies.

APPROACH

We analyzed 272 flow waveforms from the internal carotid arteries of older patients (73  ±  13 yr) with moderate cardiovascular disease, and used these data to develop methods to guide new approaches for hemodynamic studies.

MAIN RESULTS

Profound variations in waveform parameters were found within the aged population that were not seen in published data for young subjects. Common features in the aged population relative to the young included a larger systole-to-diastole flow rate ratio, increased flow during late systole, and absence of a dicrotic notch. Eight waveforms were identified that collectively represent the range of waveforms in the older population. A relationship between waveform shape and flow rate was obtained that, in conjunction with equations relating flow rate to diameter, can be used to provide individualized waveforms for patient-specific geometries. The dependence of flow rate on diameter was statistically different between male and female patients.

SIGNIFICANCE

It was shown that a single archetypal waveform cannot well-represent the diverse waveforms found within an aged population, although this approach is frequently used in studies of flow in the cerebral vasculature. Motivated by these results, we provided a set of eight waveforms that can be used to assess the hemodynamic uncertainty associated with the lack of patient-specific waveform data. We also provided a methodology for generating individualized waveforms when patient gender, age, and cardiovascular disease state are known. These data-driven approaches can be used to devise more relevant in vitro or in silico intra-cranial hemodynamic studies for older patients.

Hemodynamic Characteristics of Ruptured and Unruptured Multiple Aneurysms at Mirror and Ipsilateral Locations

BACKGROUND AND PURPOSE

Different hemodynamic patterns have been associated with aneurysm rupture. The objective was to test whether hemodynamic characteristics of the ruptured aneurysm in patients with multiple aneurysms were different from those in unruptured aneurysms in the same patient.

MATERIALS AND METHODS

Twenty-four mirror and 58 ipsilateral multiple aneurysms with 1 ruptured and the others unruptured were studied. Computational fluid dynamics models were created from 3D angiographies. Case-control studies of mirror and ipsilateral aneurysms were performed with paired Wilcoxon tests.

RESULTS

In mirror pairs, the ruptured aneurysm had more oscillatory wall shear stress (P = .007) than the unruptured one and tended to be more elongated (higher aspect ratio), though this trend achieved only marginal significance (P = .03, 1-sided test). In ipsilateral aneurysms, ruptured aneurysms had larger maximum wall shear (P = .05), more concentrated (P < .001) and oscillatory wall shear stress (P < .001), stronger (P < .001) and more concentrated inflow jets (P < .001), larger maximum velocity (P < .001), and more complex flow patterns (P < .001) compared with unruptured aneurysms. Additionally, ruptured aneurysms were larger (P < .001) and more elongated (P < .001) and had wider necks (P < .001) and lower minimum wall shear stress (P < .001) than unruptured aneurysms.

CONCLUSIONS

High wall shear stress oscillations and larger aspect ratios are associated with rupture in mirror aneurysms. Adverse flow conditions characterized by high and concentrated inflow jets; high, concentrated, and oscillatory wall shear stress; and strong, complex and unstable flow patterns are associated with rupture in ipsilateral multiple aneurysms. In multiple ipsilateral aneurysms, these unfavorable flow conditions are more likely to develop in larger, more elongated, more wide-necked, and more distal aneurysms.

Angioarchitectures and Hemodynamic Characteristics of Posterior Communicating Artery Aneurysms and Their Association with Rupture Status

BACKGROUND AND PURPOSE

Intracranial aneurysms originating at the posterior communicating artery are known to have high rupture risk compared with other locations. We tested the hypothesis that different angioarchitectures (ie, branch point configuration) of posterior communicating artery aneurysms are associated with aneurysm hemodynamics, which in turn predisposes aneurysms to rupture.

MATERIALS AND METHODS

A total of 313 posterior communicating artery aneurysms (145 ruptured, 168 unruptured) were studied with image-based computational fluid dynamics. Aneurysms were classified into different angioarchitecture types depending on the location of the aneurysm with respect to parent artery bifurcation. Hemodynamic characteristics were compared between ruptured and unruptured aneurysms, as well as among aneurysms with different angioarchitectures.

RESULTS

Angioarchitecture was associated with rupture (P = .003). Ruptured aneurysms had higher, more concentrated, and more oscillatory wall shear stress distributions (maximum wall shear stress, P < .001; shear concentration index, P < .001; mean oscillatory shear index, P < .001), stronger and more concentrated inflow jets (represented as Q, P = .01; inflow concentration index, P < .001), and more complex and unstable flow patterns (vortex core length, P < .001; proper orthogonal decomposition entropy, P < .001) compared with unruptured aneurysms. These adverse conditions were more common in aneurysms with bifurcation-type angioarchitectures compared with those with lateral or sidewall angioarchitectures. Interestingly, ruptured aneurysms also had lower normalized mean wall shear stress (P = .02) and minimum wall shear stress (P = .002) than unruptured aneurysms.

CONCLUSIONS

High-flow intrasaccular hemodynamic characteristics, commonly found in bifurcation-type angioarchitectures, are associated with the posterior communicating artery aneurysm rupture status. These characteristics include strong and concentrated inflow jets, concentrated regions of elevated wall shear stress, oscillatory wall shear stress, lower normalized wall shear stress, and complex and unstable flow patterns.

Understanding Angiography-Based Aneurysm Flow Fields through Comparison with Computational Fluid Dynamics

BACKGROUND AND PURPOSE

Hemodynamics is thought to be an important factor for aneurysm progression and rupture. Our aim was to evaluate whether flow fields reconstructed from dynamic angiography data can be used to realistically represent the main flow structures in intracranial aneurysms.

MATERIALS AND METHODS

DSA-based flow reconstructions, obtained during interventional treatment, were compared qualitatively with flow fields obtained from patient-specific computational fluid dynamics models and quantitatively with projections of the computational fluid dynamics fields (by computing a directional similarity of the vector fields) in 15 cerebral aneurysms.

RESULTS

The average similarity between the DSA and the projected computational fluid dynamics flow fields was 78% in the parent artery, while it was only 30% in the aneurysm region. Qualitatively, both the DSA and projected computational fluid dynamics flow fields captured the location of the inflow jet, the main vortex structure, the intrasaccular flow split, and the main rotation direction in approximately 60% of the cases.

CONCLUSIONS

Several factors affect the reconstruction of 2D flow fields from dynamic angiography sequences. The most important factors are the 3-dimensionality of the intrasaccular flow patterns and inflow jets, the alignment of the main vortex structure with the line of sight, the overlapping of surrounding vessels, and possibly frame rate undersampling. Flow visualization with DSA from >1 projection is required for understanding of the 3D intrasaccular flow patterns. Although these DSA-based flow quantification techniques do not capture swirling or secondary flows in the parent artery, they still provide a good representation of the mean axial flow and the corresponding flow rate.

Differences in Hemodynamics and Rupture Rate of Aneurysms at the Bifurcation of the Basilar and Internal Carotid Arteries

BACKGROUND AND PURPOSE

Cerebral aneurysms in the posterior circulation are known to have a higher rupture risk than those in the anterior circulation. We sought to test the hypothesis that differences in hemodynamics can explain the difference in rupture rates.

MATERIALS AND METHODS

A total of 117 aneurysms, 63 at the tip of the basilar artery (27 ruptured, 36 unruptured, rupture rate = 43%) and 54 at the bifurcation of the internal carotid artery (11 ruptured, 43 unruptured, rupture rate = 20%) were analyzed with image-based computational fluid dynamics. Several hemodynamic variables were compared among aneurysms at each location and between ruptured and unruptured aneurysms at each location.

RESULTS

On average, aneurysms at the basilar tip had more concentrated inflow (P < .001), a larger inflow rate (P < .001), a larger maximum oscillatory shear index (P = .003), more complex flows (P = .033), and smaller areas under low wall shear stress (P < .001) than aneurysms at the bifurcation of the internal carotid artery. In general, ruptured aneurysms had larger inflow concentration (P = .02), larger shear concentration (P = .02), more complex flows (P < .001), and smaller minimum wall shear stress (P = .003) than unruptured aneurysms.

CONCLUSIONS

High flow conditions, characterized by large and concentrated inflow jets, complex and oscillatory flow patterns, and wall shear stress distributions with focalized regions of high shear and large regions of low shear, are associated with aneurysm rupture, especially for basilar tip aneurysms. The higher flow conditions in basilar tip aneurysms could explain their increased rupture risk compared with internal carotid bifurcation aneurysms.

Hemodynamic differences between unstable and stable unruptured aneurysms independent of size and location: a pilot study

BACKGROUND

While clinical and angiographic risk factors for intracranial aneurysm instability are well established, it is reasonable to postulate that intra-aneurysmal hemodynamics also have a role in aneurysm instability.

OBJECTIVE

To identify hemodynamic characteristics that differ between radiologically unstable and stable unruptured intracranial aneurysms.

MATERIALS AND METHODS

12 pairs of unruptured intracranial aneurysms with a 3D rotational angiographic set of images and followed up longitudinally without treatment were studied. Each pair consisted of one stable aneurysm (no change on serial imaging) and one unstable aneurysm (demonstrated growth of at least 1 mm diameter or ruptured during follow-up) of matching size (within 10%) and locations. Patient-specific computational fluid dynamics models were created and run under pulsatile flow conditions. Relevant hemodynamic and geometric variables were calculated and compared between groups using the paired Wilcoxon test.

RESULTS

The area of the aneurysm under low wall shear stress (low shear stress area (LSA)) was 2.26 times larger in unstable aneurysms than in stable aneurysms (p=0.0499). The mean aneurysm vorticity was smaller by a factor of 0.57 in unstable aneurysms compared with stable aneurysms (p=0.0499). No statistically significant differences in geometric variables or shape indices were found.

CONCLUSIONS

This pilot study suggests there may be hemodynamic differences between unstable and stable unruptured cerebral aneurysms. In particular, the area under low wall shear stress was larger in unstable aneurysms. These findings should be considered tentative until confirmed by future larger studies.

Cluster Analysis of Vortical Flow in Simulations of Cerebral Aneurysm Hemodynamics

Computational fluid dynamic (CFD) simulations of blood flow provide new insights into the hemodynamics of vascular pathologies such as cerebral aneurysms. Understanding the relations between hemodynamics and aneurysm initiation, progression, and risk of rupture is crucial in diagnosis and treatment. Recent studies link the existence of vortices in the blood flow pattern to aneurysm rupture and report observations of embedded vortices -a larger vortex encloses a smaller one flowing in the opposite direction -whose implications are unclear. We present a clustering-based approach for the visual analysis of vortical flow in simulated cerebral aneurysm hemodynamics. We show how embedded vortices develop at saddle-node bifurcations on vortex core lines and convey the participating flow at full manifestation of the vortex by a fast and smart grouping of streamlines and the visualization of group representatives. The grouping result may be refined based on spectral clustering generating a more detailed visualization of the flow pattern, especially further off the core lines. We aim at supporting CFD engineers researching the biological implications of embedded vortices.

The Computational Fluid Dynamics Rupture Challenge 2013—Phase II: Variability of Hemodynamic Simulations in Two Intracranial Aneurysms

With the increased availability of computational resources, the past decade has seen a rise in the use of computational fluid dynamics (CFD) for medical applications. There has been an increase in the application of CFD to attempt to predict the rupture of intracranial aneurysms, however, while many hemodynamic parameters can be obtained from these computations, to date, no consistent methodology for the prediction of the rupture has been identified. One particular challenge to CFD is that many factors contribute to its accuracy; the mesh resolution and spatial/temporal discretization can alone contribute to a variation in accuracy. This failure to identify the importance of these factors and identify a methodology for the prediction of ruptures has limited the acceptance of CFD among physicians for rupture prediction. The International CFD Rupture Challenge 2013 seeks to comment on the sensitivity of these various CFD assumptions to predict the rupture by undertaking a comparison of the rupture and blood-flow predictions from a wide range of independent participants utilizing a range of CFD approaches. Twenty-six groups from 15 countries took part in the challenge. Participants were provided with surface models of two intracranial aneurysms and asked to carry out the corresponding hemodynamics simulations, free to choose their own mesh, solver, and temporal discretization. They were requested to submit velocity and pressure predictions along the centerline and on specified planes. The first phase of the challenge, described in a separate paper, was aimed at predicting which of the two aneurysms had previously ruptured and where the rupture site was located. The second phase, described in this paper, aims to assess the variability of the solutions and the sensitivity to the modeling assumptions. Participants were free to choose boundary conditions in the first phase, whereas they were prescribed in the second phase but all other CFD modeling parameters were not prescribed. In order to compare the computational results of one representative group with experimental results, steady-flow measurements using particle image velocimetry (PIV) were carried out in a silicone model of one of the provided aneurysms. Approximately 80% of the participating groups generated similar results. Both velocity and pressure computations were in good agreement with each other for cycle-averaged and peak-systolic predictions. Most apparent “outliers” (results that stand out of the collective) were observed to have underestimated velocity levels compared to the majority of solutions, but nevertheless identified comparable flow structures. In only two cases, the results deviate by over 35% from the mean solution of all the participants. Results of steady CFD simulations of the representative group and PIV experiments were in good agreement. The study demonstrated that while a range of numerical schemes, mesh resolution, and solvers was used, similar flow predictions were observed in the majority of cases. To further validate the computational results, it is suggested that time-dependent measurements should be conducted in the future. However, it is recognized that this study does not include the biological aspects of the aneurysm, which needs to be considered to be able to more precisely identify the specific rupture risk of an intracranial aneurysm.

Wall Mechanical Properties and Hemodynamics of Unruptured Intracranial Aneurysms

BACKGROUND AND PURPOSE

Aneurysm progression and rupture is thought to be governed by progressive degradation and weakening of the wall in response to abnormal hemodynamics. Our goal was to investigate the relationship between the intra-aneurysmal hemodynamic conditions and wall mechanical properties in human aneurysms.

MATERIALS AND METHODS

A total of 8 unruptured aneurysms were analyzed. Computational fluid dynamics models were constructed from preoperative 3D rotational angiography images. The aneurysms were clipped, and the domes were resected and mechanically tested to failure with a uniaxial testing system under multiphoton microscopy. Linear regression analysis was performed to explore possible correlations between hemodynamic quantities and the failure characteristics and stiffness of the wall.

RESULTS

The ultimate strain was correlated negatively to aneurysm inflow rate (P = .021), mean velocity (P = .025), and mean wall shear stress (P = .039). It was also correlated negatively to inflow concentration, oscillatory shear index, and measures of the complexity and instability of the flow; however, these trends did not reach statistical significance. The wall stiffness at high strains was correlated positively to inflow rate (P = .014), mean velocity (P = .008), inflow concentration (P = .04), flow instability (P = .006), flow complexity (P = .019), wall shear stress (P = .002), and oscillatory shear index (P = .004).

CONCLUSIONS

In a study of 8 unruptured intracranial aneurysms, ultimate strain was correlated negatively with aneurysm inflow rate, mean velocity, and mean wall shear stress. Wall stiffness was correlated positively with aneurysm inflow rate, mean velocity, wall shear stress, flow complexity and stability, and oscillatory shear index. These trends and the impact of hemodynamics on wall structure and mechanical properties should be investigated further in larger studies.

Diversity in the Strength and Structure of Unruptured Cerebral Aneurysms

Intracranial aneurysms are pathological enlargements of brain arteries that are believed to arise from progressive wall degeneration and remodeling. Earlier work using classical histological approaches identified variability in cerebral aneurysm mural content, ranging from layered walls with intact endothelium and aligned smooth muscle cells, to thin, hypocellular walls. Here, we take advantage of recent advances in multiphoton microscopy, to provide novel results for collagen fiber architecture in 15 human aneurysm domes without staining or fixation as well as in 12 control cerebral arteries. For all aneurysm samples, the elastic lamina was absent and the abluminal collagen fibers had similar diameters to control arteries. In contrast, the collagen fibers on the luminal side showed great variability in both diameter and architecture ranging from dense fiber layers to sparse fiber constructs suggestive of ineffective remodeling efforts. The mechanical integrity of eight aneurysm samples was assessed using uniaxial experiments, revealing two sub-classes (i) vulnerable unruptured aneurysms (low failure stress and failure pressure), and (ii) strong unruptured aneurysms (high failure stress and failure pressure). These results suggest a need to refine the end-point of risk assessment studies that currently do not distinguish risk levels among unruptured aneurysms. We propose that a measure of wall integrity that identifies this vulnerable wall subpopulation will be useful for interpreting future biological and structural data.

Hemodynamics in growing and stable cerebral aneurysms

Objective

The detailed mechanisms of cerebral aneurysm evolution are poorly understood but are important for objective aneurysm evaluation and improved patient management. The purpose of this study was to identify hemodynamic conditions that may predispose aneurysms to growth.

Methods

A total of 33 intracranial unruptured aneurysms longitudinally followed with three-dimensional imaging were studied. Patient-specific computational fluid dynamics models were constructed and used to quantitatively characterize the hemodynamic environments of these aneurysms. Hemodynamic characteristics of growing (n=16) and stable (n=17) aneurysms were compared. Logistic regression statistical models were constructed to test the predictability of aneurysm growth by hemodynamic features.

Results

Growing aneurysms had significantly smaller shear rate ratios (p=0.01), higher concentration of wall shear stress (p=0.03), smaller vorticity ratios (p=0.01), and smaller viscous dissipation ratios (p=0.01) than stable aneurysms. They also tended to have larger areas under low wall shear stress (p=0.06) and larger aspect ratios (p=0.18), but these trends were not significant. Mean wall shear stress was not significantly different between growing and stable aneurysms. Logistic regression models based on hemodynamic variables were able to discriminate between growing and stable aneurysms with a high degree of accuracy (94–100%).

Conclusions

Growing aneurysms tend to have complex intrasaccular flow patterns that induce non-uniform wall shear stress distributions with areas of concentrated high wall shear stress and large areas of low wall shear stress. Statistical models based on hemodynamic features seem capable of discriminating between growing and stable aneurysms.

Hemodynamic analysis of fast and slow aneurysm occlusions by flow diversion in rabbits

PURPOSE

To assess hemodynamic differences between aneurysms that occlude rapidly and those occluding in delayed fashion after flow diversion in rabbits.

METHODS

Thirty-six elastase-induced aneurysms in rabbits were treated with flow diverting devices. Aneurysm occlusion was assessed angiographically immediately before they were sacrificed at 1 (n=6), 2 (n=4), 4 (n=8) or 8 weeks (n=18) after treatment. The aneurysms were classified into a fast occlusion group if they were completely or near completely occluded at 4 weeks or earlier and a slow occlusion group if they remained incompletely occluded at 8 weeks. The immediate post-treatment flow conditions in aneurysms of each group were quantified using subject-specific computational fluid dynamics and statistically compared.

RESULTS

Nine aneurysms were classified into the fast occlusion group and six into the slow occlusion group. Aneurysms in the fast occlusion group were on average significantly smaller (fast=0.9 cm, slow=1.393 cm, p=0.024) and had smaller ostia (fast=0.144 cm2, slow=0.365 cm2, p=0.015) than aneurysms in the slow occlusion group. They also had a lower mean post-treatment inflow rate (fast=0.047 mL/s, slow=0.155 mL/s, p=0.0239), kinetic energy (fast=0.519 erg, slow=1.283 erg, p=0.0468), and velocity (fast=0.221 cm/s, slow=0.506 cm/s, p=0.0582). However, the differences in the latter two variables were only marginally significant.

CONCLUSIONS

Hemodynamic conditions after flow diversion treatment of cerebral aneurysms in rabbits are associated with the subsequent aneurysm occlusion time. Specifically, smaller inflow rate, kinetic energy, and velocity seem to promote faster occlusions, especially in smaller and small-necked aneurysms. These results are consistent with previous studies based on clinical series.

Strategy for analysis of flow diverting devices based on multi-modality image-based modeling

Quantification and characterization of the hemodynamic environment created after flow diversion treatment of cerebral aneurysms is important to understand the effects of flow diverters and their interactions with the biology of the aneurysm wall and the thrombosis process that takes place subsequently. This paper describes the construction of multi-modality image-based subject-specific CFD models of experimentally created aneurysms in rabbits and subsequently treated with flow diverters. Briefly, anatomical models were constructed from 3D rotational angiography images, flow conditions were derived from Doppler ultrasound measurements, stent models were created and virtually deployed, and the results were compared with in vivo digital subtraction angiography and Doppler ultrasound images. The models were capable of reproducing in vivo observations, including velocity waveforms measured in the parent artery, peak velocity values measured in the aneurysm, and flow structures observed with digital subtraction angiography before and after deployment of flow diverters. The results indicate that regions of aneurysm occlusion after flow diversion coincide with slow and smooth flow patterns, whereas regions still permeable at the time of animal sacrifice were observed in parts of the aneurysm exposed to larger flow activity, that is, higher velocities, more swirling, and more complex flow structures.

Analysis of flow changes in side branches jailed by flow diverters in rabbit models

Understanding the flow alteration in side branches during flow diversion treatment of cerebral aneurysms is important to prevent ischemic complications and improve device designs. Flow diverters were placed in the aorta of four rabbits crossing the origin of side arteries. Subject-specific computational models were constructed from 3D angiographies and Doppler ultrasounds (DUSs). Flow simulations were run before and after virtually deploying the flow diverters, assuming distal resistances remained unchanged after treatment. All jailed arteries remained patent angiographically 8 weeks after treatment. The computational models estimated decreases compared to pretreatment in the mean flow rates between 2% and 20% and in peak flow rates between 5% and 36%. The major changes were observed during systole. Flow patterns did not exhibit recirculation zones before treatment. Implantation of the flow diverters altered the flow structure only locally near the device wires. No major recirculation regions were created or destroyed. Flow diverters seem safe with respect to perforator or side branch occlusion. Relatively small changes in flow rates through jailed arteries are expected, even for moderate to large degrees of coverage of their origins. These results seem consistent with previous clinical experiences where no or very few complications related to perforator occlusion have been reported.

Strategy for modeling flow diverters in cerebral aneurysms as a porous medium

Simulations using the patient-specific geometry of the aneurysm may help in a better planning of the treatment and in a consequent reduction of the associated risks. We propose, evaluate, and implement a methodology for the simulation of flow diverter (FD) devices in intracranial aneurysms by using a porous medium method (PMM), which greatly reduces the computational cost of these simulations compared with immersed method (IMM) approaches used to model complex FDs. The method relies on parameters from an empirical correlation derived from experimental observations in wire screens, consistent with CFD simulations. The verification of our PMM strategy was carried out by comparing the results of simulations in different (patient-specific) geometries and FDs, to those obtained under identical conditions by the IMM. Overall, both quantitative and qualitative results are consistent between IMM and PMM in cases where the local porosity remains roughly uniform throughout the neck, with differences in the reduction of the observables lower than 10%. This PMM strategy is up to 10 times faster than the IMM, which allows for a runtime of hours instead of days, bringing it closer for its application in the clinic.

Unsteady wall shear stress analysis from image-based computational fluid dynamic aneurysm models under Newtonian and Casson rheological models

The aim of this work was to determine whether or not Newtonian rheology assumption in image-based patient-specific computational fluid dynamics (CFD) cerebrovascular models harboring cerebral aneurysms may affect the hemodynamics characteristics, which have been previously associated with aneurysm progression and rupture. Ten patients with cerebral aneurysms with lobulations were considered. CFD models were reconstructed from 3DRA and 4DCTA images by means of region growing, deformable models, and an advancing front technique. Patient-specific FEM blood flow simulations were performed under Newtonian and Casson rheological models. Wall shear stress (WSS) maps were created and distributions were compared at the end diastole. Regions of lower WSS (lobulation) and higher WSS (neck) were identified. WSS changes in time were analyzed. Maximum, minimum and time-averaged values were calculated and statistically compared. WSS characterization remained unchanged. At high WSS regions, Casson rheology systematically produced higher WSS minimum, maximum and time-averaged values. However, those differences were not statistically significant. At low WSS regions, when averaging over all cases, the Casson model produced higher stresses, although in some cases the Newtonian model did. However, those differences were not significant either. There is no evidence that Newtonian model overestimates WSS. Differences are not statistically significant.

Analysis of hemodynamics and aneurysm occlusion after flow-diverting treatment in rabbit models

BACKGROUND AND PURPOSE

Predicting the outcome of flow diversion treatment of cerebral aneurysms remains challenging. Our aim was to investigate the relationship between hemodynamic conditions created immediately after flow diversion and subsequent occlusion of experimental aneurysms in rabbits.

MATERIALS AND METHODS

The hemodynamic environment before and after flow-diversion treatment of elastase-induced aneurysms in 20 rabbits was modeled by using image-based computational fluid dynamics. Local aneurysm occlusion was quantified by using a voxelization technique on 3D images acquired 8 weeks after treatment. Global and local voxel-by-voxel hemodynamic variables were used to statistically compare aneurysm regions that later thrombosed to regions that remained patent.

RESULTS

Six aneurysms remained patent at 8 weeks, while 14 were completely or nearly completely occluded. Patent aneurysms had statistically larger neck sizes (P = .0015) and smaller mean transit times (P = .02). The velocity, vorticity, and shear rate were approximately 2.8 times (P < .0001) larger in patent regions-that is, they had larger "flow activity" than regions that progressed to occlusion. Statistical models based on local hemodynamic variables were capable of predicting local occlusion with good precision (84% accuracy), especially away from the neck (92%-94%). Predictions near the neck were poorer (73% accuracy).

CONCLUSIONS

These results suggests that the dominant healing mechanism of occlusion within the aneurysm dome is related to slow-flow-induced thrombosis, while near the neck, other processes could be at play simultaneously.