Patient-specific computational modeling of cerebral aneurysms with multiple avenues of flow from 3D rotational angiography images

Hormonal influences on cerebral aneurysms: unraveling the complex connections

INTRODUCTION

Intracranial aneurysms (IAs) occur in 3-5% of the general population and are characterized by localized structural deterioration of the arterial wall with loss of internal elastic lamina and disruption of the media. The risk of incidence and rupture of aneurysms depends on age, sex, ethnicity, and other different factors, indicating the influence of genetic and environmental factors. When an aneurysm ruptures, there is an estimated 20% mortality rate, along with an added 30-40% morbidity in survivors. The alterations in hormonal levels can influence IAs, while the rupture of an aneurysm can have various impacts on endocrine pathways and affect their outcome.

AREA COVERED

This review explores the reciprocal relationship between endocrinological changes (estrogen, growth hormone, and thyroid hormones) and IAs, as well as the effects of aneurysm ruptures on endocrine fluctuations.

EXPERT OPINION

Based on the data presented in this paper, we recommend further exploration into the influence of hormones on aneurysm formation and rupture. Additionally, we propose conducting endocrine assessments for patients who have experienced a rupture of IAs. Monitoring hormonal changes in patients with IAs could serve as a potential risk factor for rupture, leading to interventions in the approach to managing IAs.

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.

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.

Virtual endovascular treatment of intracranial aneurysms: models and uncertainty

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

Computational Fluid Dynamics modeling of contrast transport in basilar aneurysms following flow-altering surgeries

In vivo measurement of blood velocity fields and flow descriptors remains challenging due to image artifacts and limited resolution of current imaging methods; however, in vivo imaging data can be used to inform and validate patient-specific computational fluid dynamics (CFD) models. Image-based CFD can be particularly useful for planning surgical interventions in complicated cases such as fusiform aneurysms of the basilar artery, where it is crucial to alter pathological hemodynamics while preserving flow to the distal vasculature. In this study, patient-specific CFD modeling was conducted for two basilar aneurysm patients considered for surgical treatment. In addition to velocity fields, transport of contrast agent was simulated for the preoperative and postoperative conditions using two approaches. The transport of a virtual contrast passively following the flow streamlines was simulated to predict post-surgical flow regions prone to thrombus deposition. In addition, the transport of a mixture of blood with an iodine-based contrast agent was modeled to compare and verify the CFD results with X-ray angiograms. The CFD-predicted patterns of contrast flow were qualitatively compared to in vivo X-ray angiograms acquired before and after the intervention. The results suggest that the mixture modeling approach, accounting for the flow rates and properties of the contrast injection, is in better agreement with the X-ray angiography data. The virtual contrast modeling assessed the residence time based on flow patterns unaffected by the injection procedure, which makes the virtual contrast modeling approach better suited for prediction of thrombus deposition, which is not limited to the peri-procedural state.

Understanding the role of hemodynamics in the initiation, progression, rupture, and treatment outcome of cerebral aneurysm from medical image-based computational studies

About a decade ago, the first image-based computational hemodynamic studies of cerebral aneurysms were presented. Their potential for clinical applications was the result of a right combination of medical image processing, vascular reconstruction, and grid generation techniques used to reconstruct personalized domains for computational fluid and solid dynamics solvers and data analysis and visualization techniques. A considerable number of studies have captivated the attention of clinicians, neurosurgeons, and neuroradiologists, who realized the ability of those tools to help in understanding the role played by hemodynamics in the natural history and management of intracranial aneurysms. This paper intends to summarize the most relevant results in the field reported during the last years.

Physical factors effecting cerebral aneurysm pathophysiology

Many factors that are either blood-, wall-, or hemodynamics-borne have been associated with the initiation, growth, and rupture of intracranial aneurysms. The distribution of cerebral aneurysms around the bifurcations of the circle of Willis has provided the impetus for numerous studies trying to link hemodynamic factors (flow impingement, pressure, and/or wall shear stress) to aneurysm pathophysiology. The focus of this review is to provide a broad overview of such hemodynamic associations as well as the subsumed aspects of vascular anatomy and wall structure. Hemodynamic factors seem to be correlated to the distribution of aneurysms on the intracranial arterial tree and complex, slow flow patterns seem to be associated with aneurysm growth and rupture. However, both the prevalence of aneurysms in the general population and the incidence of ruptures in the aneurysm population are extremely low. This suggests that hemodynamic factors and purely mechanical explanations by themselves may serve as necessary, but never as necessary and sufficient conditions of this disease's causation. The ultimate cause is not yet known, but it is likely an additive or multiplicative effect of a handful of biochemical and biomechanical factors.

Suggested connections between risk factors of intracranial aneurysms: a review

The purpose of this article is to review studies of aneurysm risk factors and the suggested hypotheses that connect the different risk factors and the underlying mechanisms governing the aneurysm natural history. The result of this work suggests that at the center of aneurysm evolution there is a cycle of wall degeneration and weakening in response to changing hemodynamic loading and biomechanic stress. This progressive wall degradation drives the geometrical evolution of the aneurysm until it stabilizes or ruptures. Risk factors such as location, genetics, smoking, co-morbidities, and hypertension seem to affect different components of this cycle. However, details of these interactions or their relative importance are still not clearly understood.

Magnetic resonance fluid dynamics for intracranial aneurysms--comparison with computed fluid dynamics

BACKGROUND

Hemodynamics in intracranial aneurysms is thought to play an important role in their growth and rupture. Usual computed fluid dynamics (CFD) based on three-dimensional (3D) computed tomographic (CT) angiography requires a time-consuming process for analysis. Magnetic resonance fluid dynamics (MRFD) based on MR images is a new tool for analyzing flow dynamics and a promising method for obtaining such information more easily. We compared the data from MRFD and CFD and studied the clinical feasibility of MRFD.

METHODS

A total of 15 aneurysms, including two ruptured ones, in 15 patients were investigated with MR imaging and 3D-CT angiography. The flow data of MRFD and CFD, 3D stream lines, flow velocity profile and wall shear stress (WSS) were extracted from the image reconstruction and were compared each other.

RESULTS

Both flow dynamics images showed quite similar 3D flow pattern and WSS map. However, the calculated value of maximum WSS was quite different and there was no significant correlation. Further, in one ruptured case, CFD showed less visualization to evaluate the intra-aneurysmal flow. Interestingly, one delayed rupture case showed a particular flow pattern with abnormal secondary flow in the bottom of the aneurysm before rupture, which might suggest the specific finding of rupture risk.

CONCLUSION

MRFD is a valuable and less invasive tool to evaluate aneurysmal fluid dynamics. It can be obtained from the usual MRI examination without contrast medium and exposure to radiation. Although there is a problem of consistency of the absolute value of WSS between MRFD and conventional CFD, it may be useful to predict the risk of enlargement or rupture of aneurysms based on the information of the similar distribution of WSS and flow patterns. The quantifiable analysis and establishment of a meaningful threshold for high risk should be further studied.

Sensitivity to outflow boundary conditions and level of geometry description for a cerebral aneurysm

Mathematical models, namely the flow boundary conditions, as well as the detail of the bounding geometry, can highly influence the computed flow field. In this work, an anatomically realistic portion of cerebral vasculature with a saccular aneurysm, and its geometric idealisation, are considered. The importance of the geometric description, namely including the side branches or modelling them as holes in the main vessel, is studied. Several approaches to prescribe the outflow boundary conditions at the side branches are analysed, including the traction-free condition, zero velocity (hence neglecting the side-branch), and the coupling with simple zero-dimensional and one-dimensional models. Results of the effects of outflow boundary modelling choice on computed haemodynamic parameters are used to identify appropriateness of the models based on the physical interpretation. Estimated range of error-bars associated to outflow boundary model choice and the level of geometric details are presented for patient-specific computational haemodynamics, and can serve as invitation for future studies. The zero-dimensional and one-dimensional models are shown to provide good representations of the side branches in the case of the clipped geometry.

Clinical Application of Image-Based CFD for Cerebral Aneurysms

During the last decade, the convergence of medical imaging and computational modeling technologies has enabled tremendous progress in the development and application of image-based computational fluid dynamics modeling of patient-specific blood flows. These techniques have been used for studying the basic mechanisms involved in the initiation and progression of vascular diseases, for studying possible ways to improve the diagnosis and evaluation of patients by incorporating hemodynamics information to the anatomical data typically available, and for the development of computational tools that can be used to improve surgical and endovascular treatment planning. However, before these technologies can have a significant impact on the routine clinical practice, it is still necessary to demonstrate the connection between the extra information provided by the models and the natural progression of vascular diseases and the outcome of interventions. This paper summarizes some of our contributions in this direction, focusing in particular on cerebral aneurysms.

Automated segmentation of cerebral vasculature with aneurysms in 3DRA and TOF-MRA using geodesic active regions: an evaluation study

PURPOSE

To evaluate the suitability of an improved version of an automatic segmentation method based on geodesic active regions (GAR) for segmenting cerebral vasculature with aneurysms from 3D x-ray reconstruction angiography (3DRA) and time of flight magnetic resonance angiography (TOF-MRA) images available in the clinical routine.

METHODS

Three aspects of the GAR method have been improved: execution time, robustness to variability in imaging protocols, and robustness to variability in image spatial resolutions. The improved GAR was retrospectively evaluated on images from patients containing intracranial aneurysms in the area of the Circle of Willis and imaged with two modalities: 3DRA and TOF-MRA. Images were obtained from two clinical centers, each using different imaging equipment. Evaluation included qualitative and quantitative analyses of the segmentation results on 20 images from 10 patients. The gold standard was built from 660 cross-sections (33 per image) of vessels and aneurysms, manually measured by interventional neuroradiologists. GAR has also been compared to an interactive segmentation method: isointensity surface extraction (ISE). In addition, since patients had been imaged with the two modalities, we performed an intermodality agreement analysis with respect to both the manual measurements and each of the two segmentation methods.

RESULTS

Both GAR and ISE differed from the gold standard within acceptable limits compared to the imaging resolution. GAR (ISE) had an average accuracy of 0.20 (0.24) mm for 3DRA and 0.27 (0.30) mm for TOF-MRA, and had a repeatability of 0.05 (0.20) mm. Compared to ISE, GAR had a lower qualitative error in the vessel region and a lower quantitative error in the aneurysm region. The repeatability of GAR was superior to manual measurements and ISE. The intermodality agreement was similar between GAR and the manual measurements.

CONCLUSIONS

The improved GAR method outperformed ISE qualitatively as well as quantitatively and is suitable for segmenting 3DRA and TOF-MRA images from clinical routine.

Cerebrovascular biomodeling for aneurysm surgery: simulation-based training by means of rapid prototyping technologies

OBJECTIVE

Opportunities for developing procedural skills are progressively rare. Therefore, sophisticated educational tools are highly warranted.

METHODS

This study compared stereolithography and 3-dimensional printing for simulating cerebral aneurysm surgery. The latter jets multiple materials simultaneously and thus has the ability to print assemblies of multiple materials with different features. The authors created the solid skull and the cerebral vessels in different materials to simulate the real aneurysm when clipped.

RESULTS

Precise plastic replicas of complex anatomical data provide intuitive tactile views that can be scrutinized from any perspective. Hollowed out vessel sections allow serial clipping efforts, evaluation of different clips, and clip positions. The models can be used for accurate prediction of vascular anatomy, for optimization of teaching surgical skills, for advanced procedural competency training, and for patient counseling.

CONCLUSION

Simultaneous 3-dimensional printing is the most promising rapid prototyping technique to produce biomodels that meet the high demands of neurovascular surgery.

Association of hemodynamic characteristics and cerebral aneurysm rupture

BACKGROUND AND PURPOSE

Hemodynamic factors are thought to play an important role in the initiation, growth, and rupture of cerebral aneurysms. This report describes a study of the associations between qualitative intra-aneurysmal hemodynamics and the rupture of cerebral aneurysms.

MATERIALS AND METHODS

Two hundred ten consecutive aneurysms were analyzed by using patient-specific CFD simulations under pulsatile flow conditions. The aneurysms were classified into categories by 2 blinded observers, depending on the complexity and stability of the flow pattern, size of the impingement region, and inflow concentration. A statistical analysis was then performed with respect to the history of previous rupture. Interobserver variability analysis was performed.

RESULTS

Ruptured aneurysms were more likely to have complex flow patterns (83%, P < .001), stable flow patterns (75%, P = .0018), concentrated inflow (66%, P = <.0001), and small impingement regions (76%, P = .0006) compared with unruptured aneurysms. Interobserver variability analyses indicated that all the classifications performed were in very good agreement-that is, well within the 95% CI.

CONCLUSIONS

A qualitative hemodynamic analysis of cerebral aneurysms by using image-based patient-specific geometries has shown that concentrated inflow jets, small impingement regions, complex flow patterns, and unstable flow patterns are correlated with a clinical history of prior aneurysm rupture. These qualitative measures provide a starting point for more sophisticated quantitative analysis aimed at assigning aneurysm risk of future rupture. These analyses highlight the potential for CFD to play an important role in the clinical determination of aneurysm risks.

Temporal variations of wall shear stress parameters in intracranial aneurysms--importance of patient-specific inflow waveforms for CFD calculations

PURPOSE

To assess reliability of wall shear stress (WSS)calculations using computational fluid dynamics (CFD) dependent on inflow in internal carotid artery aneurysms (ICA).

MATERIALS AND METHODS

Six unruptured ICA aneurysms were studied. 3D computational meshes were created from 3D digital subtraction angiographic images (Axiom Artis dBA, Siemens Medical Solutions). Transient CFD simulations(Fluent, ANSYS Inc.) were performed for two inflow conditions: (1) idealized averaged waveform from normal subjects (ID) and (2) patient-specific waveform (PS)measured with 2D phase contrast magnetic resonance imaging. Stability of calculation was assessed by comparing mean WSS (<WSS>), temporal wall shear stress magnitude variation (Delta WSS), and oscillatory shear index(OSI, a measure of variation in the WSS direction) on the aneurysmal wall for both conditions.

RESULTS

For all cases, mean relative difference (PS-ID) of WSS (<WSS>) was -15% (range -32% to 11%). Mean Delta WSS difference was -29.3% ( -100% to 67%). Mean OSI difference was 7.5% (-12% to 40%). Large variations in histograms of these parameters were noted.

CONCLUSION

For accurate calculations of WSS parameters,patient-specific information on physiological flow may be necessary. Results obtained with averaged or idealized flow waveforms may have to be interpreted with caution.

Toward integrated management of cerebral aneurysms

In the last few years, some of the visionary concepts behind the virtual physiological human began to be demonstrated on various clinical domains, showing great promise for improving healthcare management. In the current work, we provide an overview of image- and biomechanics-based techniques that, when put together, provide a patient-specific pipeline for the management of intracranial aneurysms. The derivation and subsequent integration of morphological, morphodynamic, haemodynamic and structural analyses allow us to extract patient-specific models and information from which diagnostic and prognostic descriptors can be obtained. Linking such new indices with relevant clinical events should bring new insights into the processes behind aneurysm genesis, growth and rupture. The development of techniques for modelling endovascular devices such as stents and coils allows the evaluation of alternative treatment scenarios before the intervention takes place and could also contribute to the understanding and improved design of more effective devices. A key element to facilitate the clinical take-up of all these developments is their comprehensive validation. Although a number of previously published results have shown the accuracy and robustness of individual components, further efforts should be directed to demonstrate the diagnostic and prognostic efficacy of these advanced tools through large-scale clinical trials.

Computer-assisted extraction of intracranial aneurysms on 3D rotational angiograms for computational fluid dynamics modeling

PURPOSE

Three-dimensional rotational angiography (3DRA) is an evolving imaging procedure from traditional digital subtraction angiography and is gaining much interest for detecting intracranial aneurysms. Computational fluid dynamics (CFD) modeling plays an important role in understanding the biomechanical properties and in facilitating the prediction of aneurysm rupture. A successful computational study relies on an accurate description of the vascular geometry that is obtained from volumetric images.

METHODS

The authors propose a new aneurysm segmentation algorithm to facilitate the study of CFD. This software combines a region-growing segmentation method with the 3D extension of a deformable contour based on a charged fluid model. A charged fluid model essentially consists of a set of charged elements that are governed by the nature of electrostatics. The approach requires no prior knowledge of anatomic structures and automatically segments the vasculature after the end-user selects a vessel section in a plane image.

RESULTS

Experimental results on 15 cases indicate that aneurysm structures were effectively segmented and in good agreement with manual delineation outcomes. In comparison with the existing methods, the algorithm provided a much higher overlap index with respect to the ground truth. Furthermore, the outcomes of the proposed approach achieved a clean representation of vascular structures that is advantageous for hemodynamics analyses.

CONCLUSIONS

A new aneurysm segmentation framework in an attempt to automatically segment vascular structures in 3DRA image volumes has been developed. The proposed algorithm demonstrated promising performance and unique characteristics to adequately segment aneurysms in 3DRA image volumes for further study in computational fluid dynamics.

Hemodynamics and bleb formation in intracranial aneurysms

BACKGROUND AND PURPOSE

Intracranial aneurysms with irregular shapes and blebs or secondary outpouchings have been correlated with increased rupture risk. The purpose of this study was to investigate possible associations between the local hemodynamics and the formation of blebs in cerebral aneurysms.

MATERIALS AND METHODS

Computational models of 20 cerebral aneurysms harboring 30 well-defined blebs were constructed from 3D rotational angiographies. Models representing the aneurysm before bleb formation were constructed by virtually removing the blebs from the anatomic models. Computational fluid dynamics simulations of the aneurysm before and after bleb formation were performed under pulsatile flows. Flow and WSS visualizations were used to analyze the local hemodynamics in the region of the aneurysm that developed the bleb.

RESULTS

Most blebs (80%) occurred at or adjacent to the aneurysm region with the highest WSS before bleb formation, and near the flow impaction zone. Most blebs (83%) were found in regions of the aneurysm previously subjected to high or moderate WSS and progressed to low WSS states after the blebs were formed. Most blebs (77%) were aligned or adjacent to the inflow jet, whereas 17% were aligned with the outflow jet, and only 6% were not aligned with the flow direction. In addition, 90% of the aneurysms had maximal WSS higher than or similar to the WSS in the parent artery.

CONCLUSIONS

Blebs form at or adjacent to regions of high WSS and are aligned with major intra-aneurysmal flow structures. Formation of blebs results in a lower WSS state with formation of a counter current vortex. These findings imply that locally elevated WSS could contribute to the focalized wall damage that formed these structures.

Vascular graph model to simulate the cerebral blood flow in realistic vascular networks

At its most fundamental level, cerebral blood flow (CBF) may be modeled as fluid flow driven through a network of resistors by pressure gradients. The composition of the blood as well as the cross-sectional area and length of a vessel are the major determinants of its resistance to flow. Here, we introduce a vascular graph modeling framework based on these principles that can compute blood pressure, flow and scalar transport in realistic vascular networks. By embedding the network in a computational grid representative of brain tissue, the interaction between the two compartments can be captured in a truly three-dimensional manner and may be applied, among others, to simulate oxygen extraction from the vessels. Moreover, we have devised an upscaling algorithm that significantly reduces the computational expense and eliminates the need for detailed knowledge on the topology of the capillary bed. The vascular graph framework has been applied to investigate the effect of local vascular dilation and occlusion on the flow in the surrounding network.

Quantitative hemodynamic analysis of brain aneurysms at different locations

BACKGROUND AND PURPOSE

Studies have shown that the occurrence of brain aneurysms and risk of rupture vary between locations. However, the reason that aneurysms at different branches of the cerebral arteries have different clinical presentations is not clear. Because research has indicated that aneurysm hemodynamics may be one of the important factors related to aneurysm growth and rupture, our aim was to analyze and compare the flow parameters in aneurysms at different locations.

MATERIALS AND METHODS

A total of 24 patient-specific aneurysm models were constructed by using 3D rotational angiographic data for the hemodynamic simulation. Previously developed computational fluid dynamics software was applied to each aneurysm to simulate the blood-flow properties. Hemodynamic data at peak pulsatile flow were recorded, and wall shear stress (WSS) and flow rate in the aneurysms and parent arteries were quantitatively compared. To validate our method, a comparison with a previously reported technique was also performed.

RESULTS

WSS and flow rate in the aneurysms at the peak of the cardiac cycle were found to differ in magnitude between different locations. Multiple comparisons among locations showed higher WSS and flow rate in middle cerebral artery aneurysms and lower WSS and flow rate in basilar artery and anterior communicating artery aneurysms.

CONCLUSIONS

We observed changes in hemodynamic values that may be related to aneurysm location. Further study of aneurysm locations with a large number of cases is needed to test this hypothesis.

Hemodynamic patterns of anterior communicating artery aneurysms: a possible association with rupture

BACKGROUND AND PURPOSE

The anterior communicating artery (AcomA) is a predilect location of aneurysms which typically carry higher rupture risks than other locations in the anterior circulation. The purpose of this study was to characterize the different flow types present in AcomA aneurysms and to investigate possible associations with rupture.

MATERIALS AND METHODS

Patient-specific computational models of 26 AcomA aneurysms were constructed from 3D rotational angiography images. Bilateral images were acquired in 15 patients who had both A1 segments of the anterior cerebral arteries, and models of the whole anterior circulation were created by fusing the reconstructed left and right arterial trees. Computational fluid dynamics simulations were performed under pulsatile flow conditions measured on a healthy subject. Visualizations of flow velocity, instantaneous streamlines, and wall shear stress (WSS) were performed. These were analyzed for flow patterns, size of the impaction zone, and peak WSS and then correlations were made with prior history of rupture.

RESULTS

Aneurysms with small impaction zones were more likely to have ruptured than those with large impaction zones (83% versus 63%). Maximum intra-aneurysmal WSS (MWSS) for the unruptured aneurysms ranged from 10 to 230 dyne/cm(2) (mean, 114 dyne/cm(2)) compared with ruptured aneurysms, which ranged from 35 to 1500 dyne/cm(2) (mean, 271 dyne/cm(2)). This difference in MWSS was statistically significant at 90% confidence levels (P = .10).

CONCLUSIONS

Aneurysms with small impaction zones, higher flow rates entering the aneurysm, and elevated MWSS are associated with a clinical history of previous rupture.

Patient-specific modeling of cardiovascular mechanics

Advances in numerical methods and three-dimensional imaging techniques have enabled the quantification of cardiovascular mechanics in subject-specific anatomic and physiologic models. Patient-specific models are being used to guide cell culture and animal experiments and test hypotheses related to the role of biomechanical factors in vascular diseases. Furthermore, biomechanical models based on noninvasive medical imaging could provide invaluable data on the in vivo service environment where cardiovascular devices are employed and on the effect of the devices on physiologic function. Finally, patient-specific modeling has enabled an entirely new application of cardiovascular mechanics, namely predicting outcomes of alternate therapeutic interventions for individual patients. We review methods to create anatomic and physiologic models, obtain properties, assign boundary conditions, and solve the equations governing blood flow and vessel wall dynamics. Applications of patient-specific models of cardiovascular mechanics are presented, followed by a discussion of the challenges and opportunities that lie ahead.

The haemodynamics of endovascular aneurysm treatment: a computational modelling approach for estimating the influence of multiple coil deployment

This paper proposes a novel computational methodology for modelling the haemodynamic effects of endovascular coil embolization for cerebral aneurysms. We employ high-resolution 3-D angiographic data to reconstruct the intracranial geometry and we model the coiled part of the aneurysm as a porous medium, with porosity decreasing as coils are inserted. The actual dimensions of the coils employed are used to determine the characteristics of the porous medium. Simulation results for saccular aneurysms from the anterior communicating and middle cerebral arteries show that insertion of coils rapidly changes intraaneurysmal blood flow and causes reduction in mural pressure and blood velocity up to stagnation, providing favorable conditions for thrombus formation and obliteration of the aneurysm.

BOTTLENECK FACTOR AND HEIGHT-WIDTH RATIO

OBJECTIVE

Determining factors predictive of the natural risk of rupture of cerebral aneurysms is difficult because of the need to control for confounding variables. We studied factors associated with rupture in a study model of patients with multiple cerebral aneurysms, one aneurysm that had ruptured and one or more that had not, in which each patient served as their own internal control.

METHODS

We collected aneurysm location, one-dimensional measurements, and two-dimensional indices from the computed tomographic angiograms of patients in the proposed study model and compared ruptured versus unruptured aneurysms. Bivariate statistics were supplemented with multivariable logistic regression analysis to model ruptured status. A total of 40 candidate models were evaluated for predictive power and fit with Wald scoring, Cox and Snell R2, Hosmer and Lemeshow tests, case classification counting, and residual analysis to determine which of the computed tomographic angiographic measurements or indices were jointly associated with and predictive of aneurysm rupture.

RESULTS

Thirty patients with 67 aneurysms (30 ruptured, 37 unruptured) were studied. Maximum diameter, height, maximum width, bulge height, parent artery diameter, aspect ratio, bottleneck factor, and aneurysm/parent artery ratio were significantly (P &lt; 0.05) associated with ruptured aneurysms on bivariate analysis. When best subsets and stepwise multivariable logistic regression was performed, bottleneck factor (odds ratio = 1.25, confidence interval = 1.11–1.41 for every 0.1 increase) and height-width ratio (odds ratio = 1.23, confidence interval = 1.03–1.47 for every 0.1 increase) were the only measures that were significantly predictive of rupture.

CONCLUSION

In a case-control study of patients with multiple cerebral aneurysms, increased bottleneck factor and height-width ratio were consistently associated with rupture.

Noninvasive measurement of intra-aneurysmal pressure and flow pattern using phase contrast with vastly undersampled isotropic projection imaging

BACKGROUND AND PURPOSE

Currently, more reliable parameters to predict the risk of aneurysmal rupture are needed. Intra-aneurysmal pressure gradients and flow maps could provide additional information regarding the risk of rupture. Our hypothesis was that phase contrast with vastly undersampled isotropic projection reconstruction (PC-VIPR), a novel 3D MR imaging sequence, could accurately assess intra-aneurysmal pressure gradients in a canine aneurysmal model when compared with invasive measurements.

MATERIALS AND METHODS

A total of 13 surgically created aneurysms in 8 canines were included in this study. Pressure measurements were performed in the parent vessel, aneurysm neck, and 5 regions within the aneurysmal sac with a microcatheter. PC-VIPR sequence was used to obtain cardiac-gated velocity measurements in a region covering the entire aneurysm. The velocity and pressure gradient maps derived from the PC-VIPR data were then coregistered with the anatomic DSA images and compared with catheter measurements.

RESULTS

In 7 of the bifurcation aneurysms, the velocity flow maps demonstrated a recirculation flow pattern with a small neck-to-dome pressure gradient (mean, +0.5 mm Hg). In 1 bifurcation aneurysm, a flow jet extending from the neck to the dome with significantly greater pressure gradient (+50.2 mm Hg) was observed. All sidewall aneurysms had low flow in the sac with intermediate pressure gradients (mean, +8.3 mm Hg). High statistical correlation existed between PC-VIPR aneurysmal pressures and microcatheter pressure measurements (R = 0.82, P < .01).

CONCLUSION

PC-VIPR can provide anatomic as well as noninvasive quantitative and qualitative hemodynamic information in the canine aneurysm model. The PC-VIPR intra-aneurysmal pressure measurements correlated well with catheter measurements.

Computational fluid dynamics modeling of intracranial aneurysms: qualitative comparison with cerebral angiography

RATIONALE AND OBJECTIVE

The purpose of this study is to determine whether computational fluid dynamics modeling can correctly predict the location of the major intra-aneurysmal flow structures that can be identified by conventional angiography.

MATERIALS AND METHODS

Patient-specific models of three cerebral aneurysms were constructed from three-dimensional rotational angiography images and computational fluid dynamic simulations performed. Using these velocity fields, contrast transport was simulated and visualizations constructed to provide a "virtual" angiogram. These models were then compared to images from high frame rate conventional angiography to compare flow structures.

RESULTS

Computational fluid dynamics simulations showed three distinct flow types ranging from simple to complex. Virtual angiographic images showed good agreement with images from conventional angiography for all three aneurysms with analogous size and orientation of the inflow jet, regions of impaction, and flow type. Large intra-aneurysmal vortices and regions of outflow also corresponded between the images.

CONCLUSIONS

Patient-specific image-based computational models of cerebral aneurysms can realistically reproduce the major intra-aneurysmal flow structures observed with conventional angiography. The agreement between computational models and angiographic structures is less for slower zones of recirculation later in the cardiac cycle.