Generalized versus Patient-Specific Inflow Boundary Conditions in Computational Fluid Dynamics Simulations of Cerebral Aneurysmal Hemodynamics

I.G.H. Jansen; J.J. Schneiders; W.V. Potters; P. van Ooij; R. van den Berg; E. van Bavel; H.A. Marquering; C.B.L.M. Majoie

doi:10.3174/ajnr.A3901

Abstract

BACKGROUND AND PURPOSE: Attempts have been made to associate intracranial aneurysmal hemodynamics with aneurysm growth and rupture status. Hemodynamics in aneurysms is traditionally determined with computational fluid dynamics by using generalized inflow boundary conditions in a parent artery. Recently, patient-specific inflow boundary conditions are being implemented more frequently. Our purpose was to compare intracranial aneurysm hemodynamics based on generalized versus patient-specific inflow boundary conditions.

MATERIALS AND METHODS: For 36 patients, geometric models of aneurysms were determined by using 3D rotational angiography. 2D phase-contrast MR imaging velocity measurements of the parent artery were performed. Computational fluid dynamics simulations were performed twice: once by using patient-specific phase-contrast MR imaging velocity profiles and once by using generalized Womersley profiles as inflow boundary conditions. Resulting mean and maximum wall shear stress and oscillatory shear index values were analyzed, and hemodynamic characteristics were qualitatively compared.

RESULTS: Quantitative analysis showed statistically significant differences for mean and maximum wall shear stress values between both inflow boundary conditions (P < .001). Qualitative assessment of hemodynamic characteristics showed differences in 21 cases: high wall shear stress location (n = 8), deflection location (n = 3), lobulation wall shear stress (n = 12), and/or vortex and inflow jet stability (n = 9). The latter showed more instability for the generalized inflow boundary conditions in 7 of 9 patients.

CONCLUSIONS: Using generalized and patient-specific inflow boundary conditions for computational fluid dynamics results in different wall shear stress magnitudes and hemodynamic characteristics. Generalized inflow boundary conditions result in more vortices and inflow jet instabilities. This study emphasizes the necessity of patient-specific inflow boundary conditions for calculation of hemodynamics in cerebral aneurysms by using computational fluid dynamics techniques.

ABBREVIATIONS:

CFD: computational fluid dynamics
PC-MR imaging: 2D phase-contrast MR imaging
WSS: wall shear stress

It has been estimated that the prevalence of intracranial aneurysms in the adult population is between 1% and 5%.¹ Although most aneurysms go undetected, acute rupture resulting in subarachnoid hemorrhage is associated with high morbidity and fatality rates.^2,3 Ruptured aneurysms are treated by coiling or clipping to prevent rebleed. The indication for preventive treatment of unruptured aneurysms is, however, not straightforward.^4,5 The risk of treatment has to be carefully balanced against the risk of rupture. At present, rupture-risk assessment of unruptured intracranial aneurysms and the decision to treat or wait and scan are mainly based on size, location, and growth of the aneurysm.⁶ It is, however, clear that the predictive value of these characteristics is limited.^1,6⇓–8 It is therefore crucial to search for additional and more predictive parameters for aneurysm rupture risk assessment.

Aneurysmal hemodynamics, in particular wall shear stress (WSS) and vortex instability, have been proposed as additional risk factors for aneurysm growth and rupture.^9,10 It has been shown that the combination of vortex instability and high or low WSS within the aneurysm is more prevalent in ruptured cases.^11⇓–13

In many studies, computational fluid dynamics (CFD) is used to simulate aneurysmal hemodynamics. CFD is traditionally performed by using generalized inflow boundary conditions based on typical flow rates in a healthy adult.^{14⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓–27} Recently, several studies have replaced these generalized inflow boundary conditions by patient-specific velocity measurements in the vessels proximal to the aneurysm.^{10,12,28⇓⇓⇓–32} In these studies, either 2D phase-contrast MR imaging (PC-MR imaging) or transcranial Doppler sonography was used to measure the flow. So far, only 3 studies have compared patient-specific with generalized inflow boundary conditions in a total of 14 aneurysms.^28⇓–30 Evidently, the necessity of using patient-specific inflow boundary conditions has not been elucidated to the full extent. In this study, we assessed the effects of patient-specific inflow boundary conditions in a group of 36 patients.

Materials and Methods

Patient Selection

Image data of 36 aneurysms in 36 patients who presented at the Academic Medical Center, Amsterdam, The Netherlands with ruptured or unruptured aneurysms from January 2009 to October 2011 were retrospectively selected from a cohort of 164 patients in an ongoing study of aneurysm hemodynamics. This was done on the basis of the high signal quality of inflow velocity measurements. Of the selected aneurysms, 9 were located in the medial cerebral artery; 7 in the carotid artery; 6 in the anterior communicating artery; 6 in the posterior communicating artery; 3 in the basilar artery; 2 in the pericallosal artery; 1 in the anterior cerebral artery; 1 in the vertebral artery; and 1 in the ophthalmic artery. Six aneurysms were ruptured. Aneurysm size ranged from 3.2 to 12.4 mm. Dome-to-neck ratio ranged from 0.72 to 2.32. A Glasgow Outcome Score of ≥4 for patients with ruptured aneurysms was mandatory.³³ Exclusion criteria were contraindications for MR imaging, including treatment of the aneurysm by surgical clipping. The study was approved by the local ethics committee. Written informed consent was obtained from all patients.

Imaging

All patients underwent high-resolution 3D rotational angiography as part of the standard clinical work-up with de novo aneurysms. This was done in either the awake state (in the case of unruptured aneurysms) or during endovascular treatment with the patient under general anesthesia (in the case of ruptured aneurysms). A single-plane angiographic unit was used (Integris Allura Neuro; Philips Healthcare, Best, the Netherlands). Twenty-one milliliters of contrast agent was administered at 3 mL/s (iodixanol, Visipaque; GE Healthcare, Cork, Ireland). This resulted in a 256³ isotropic image volume. Following this, 3D velocity measurements proximal to the aneurysm were obtained with PC-MR imaging. A single-section PC-MR imaging was performed on a 3T scanner (Intera; Philips Healthcare). Scan resolution was 0.64 × 0.65 × 3 mm. Further imaging parameters were the following: TE/TR/flip angle, 5.7 ms/8.5 ms/10°; receiver bandwidth, 172 kHz; imaging volume, 200 × 200 × 3 mm in 1 section; parallel imaging factor, 2. The velocity-encoding was 100 cm/s in all directions. The number of measured cardiac phases (ie, temporal resolution) depended on the heart rate and ranged between 23 and 36 cardiac phases, to keep the scanning time close to 3 minutes 30 seconds. The view-sharing factor for the retrospective sorting of acquired k-lines was set to 1.8.

Because patients with ruptured aneurysms were treated within 24 hours after onset, MR imaging velocity measurements to assess inflow boundary conditions could not be performed before the coiling procedure. Therefore, for this patient group, postprocedural PC-MR images were obtained at the standard follow-up 6 months after coiling. Patients with unruptured aneurysms were requested to undergo an additional preprocedural PC-MR imaging study at admission.

Geometric Vascular Models

To generate vascular models that were usable for CFD, we segmented the aneurysm and its connected arteries in 3D rotational angiography images by using a level-set algorithm by using the Vascular Modeling Toolkit, VMTK Version 0.9.0 (http://www.vmtk.org). Subsequently, the segmented volumes were converted to volumetric meshes consisting of approximately 1,000,000 tetrahedral elements.

Computational Fluid Dynamics

For each aneurysm, we performed 2 CFD simulations: one applying spatiotemporal patient-specific inflow boundary conditions, acquired by PC-MR imaging velocity measurements, and the other applying generalized inflow boundary conditions.³⁴ The generalized inflow velocity profile was defined by predetermined Womersley profiles for fully developed pulsatile flow.^35,36 The flow was scaled so that the total generalized inflow equaled the measured inflow rate as determined by PC-MR imaging. Zero pressure boundary conditions were prescribed at all outlets. A no-slip boundary was set, and rigid walls were assumed. Transient Navier-Stokes equations were solved by using a pressure-based, 3D double-precision solver following the Semi-Implicit Method for Pressure Linked Equations.³⁷ Blood was modeled with an attenuation of 1040 kg/m³ and a dynamic viscosity of 0.004 Pa/s. CFD simulations were performed with Fluent software (ANSYS, Canonsburg, Pennsylvania). Three cardiac cycles were simulated to account for the transient character. Only the third complete cycle was used for analysis.³⁶

Quantitative Assessment of Hemodynamic Features

The mean WSS, maximum WSS, mean oscillatory shear index, and maximum oscillatory shear index values were calculated. Paired differences between the generalized and patient-specific inflow boundary conditions were analyzed by using paired t test statistics. Results with P values of <.05 were considered statistically significant.

Qualitative Assessment of Hemodynamic Features

Velocity-based streamlines and WSS patterns within the aneurysm during 1 cardiac cycle were visualized as movie clips with ParaView software (Kitware; Los Alamos National Laboratory, Los Alamos, New Mexico). On the basis of these movie clips, hemodynamic characteristics were scored in consensus by 2 neuroradiologists with >10 years of experience. The cases were randomly presented to the observers, who were blinded to the method of inflow boundary conditions used.

We assessed the following hemodynamic characteristics: inflow jet concentration, inflow jet stability, number of vortices, vortex stability, location of the highest WSS, and location of the deflection zone. If the aneurysm contained a lobulation, WSS on the lobulation sac and on the ostium of the lobulation were also assessed, as well as the direction of the main inflow jet toward the lobulation.

The inflow jet was defined as an isosurface of 25% of the maximum velocity magnitude within the aneurysm. It was considered “concentrated” when interpreted as smaller than half the size of the aneurysm neck, and “diffuse” if larger. In the event of the flow pattern deteriorating or changing considerably during the cardiac cycle, it was labeled as “unstable.” The number of vortices was defined as the number of flow structures within the aneurysm and, if present, in the lobulation. The location of the highest WSS and deflection zone were classified according to their position as dome, body, or neck.³⁶ The “deflection zone” was defined as the area of divergence of the inflow jet on the aneurysm wall. For each hemodynamic parameter, differences between the 2 inflow boundary conditions were assessed and the amount of difference was rated between 1 and 5, with 5 representing a large difference. Differences between the scores of the observers were discussed in an additional meeting to reach a consensus.

Results

Figure 1 presents differences in hemodynamic characteristics between generalized or patient-specific inflow boundary conditions in the same aneurysm. For each scored hemodynamic characteristic, 1 example is given. An example of the flow rate curve for both methods is displayed in Fig 2.

Fig 1.

Examples of visualized differences in hemodynamic characteristics. The left column shows patient-specific inflow boundary conditions; the right column shows generalized boundary conditions. A–H, Differences in distribution for all assessed WSS characteristics. A, No difference in any characteristic. B, High WSS area on the ostium of lobulation. C, High WSS location on the primary aneurysm. D, High WSS on the lobulation sac. E, Deflection zone. F and G, Isosurface projections of inflow jet characteristics. F, Inflow jet concentration. G, Inflow jet aim in the lobulation. H, Visualization of intra-aneurysmal flow structures represented by velocity magnitude streamlines. On this figure, a difference in the number of vortices can be appreciated, with an additional vortex for the patient-specific inflow boundary conditions (blue arrow).

Fig 2.

Examples of the applied flow-rate pattern during 1 heart cycle for both methods. PS indicates patient-specific flow rate curves of 3 separate cases (blue, purple, and green); G, generalized flow-rate curve (red).

Quantitative Assessment of Hemodynamic Features

Table 1 displays the mean and maximum WSS and oscillatory shear index values per method, including their relative differences. Mean WSS was larger for the generalized inflow boundary conditions, with an average of 3.5 Pa (P < .001). The maximum WSS was also larger for the generalized approach with an average of 65 Pa (P = .0013). The differences for mean or maximum oscillatory shear index values were not statistically significant (P = .42 and P = .65, respectively). In On-line Table 1, all measured values per case used in the quantitative assessment, including their relative differences, are displayed.

View this table:

Table 1:

Quantitative values for mean and maximum WSS and OSI^{^a}

Qualitative Assessment of Hemodynamic Features

Table 2 displays the number of differences in hemodynamic characteristics as scored during the qualitative assessment. In 21 of 36 aneurysms, at least 1 of the characteristics was scored differently for the 2 inflow boundary conditions. Seven aneurysms with differences in WSS characteristics showed no differences in vortex and inflow jet characteristics (eg, patient 16). Six patients with no differences in WSS characteristics showed distinct differences in vortex and inflow jet characteristics (eg, patient 27). Furthermore, aneurysms showing differences in vortex and inflow jet characteristics were more often assessed as unstable for the generalized inflow boundary conditions. This was seen in 7 of 9 aneurysms with differences in vortex characteristics and all aneurysms with differences in inflow jet characteristics. In On-line Table 2, a detailed overview of the qualitative assessment is given, in which each case can be appreciated for which hemodynamic characteristic differences between the 2 methods were scored.

View this table:

Table 2:

Qualitative assessment: number of cases for each scored hemodynamic characteristic showing a difference between the patient-specific inflow boundary conditions and the generalized inflow boundary conditions

Discussion

We showed that the choice of general or patient-specific inflow boundary conditions results in large differences in WSS magnitude and distribution. Differences in vortex and inflow jet characteristics occurred less frequently. The approach by using generalized inflow boundary conditions led to considerably more unstable vortices and inflow jets, suggesting that inflow jet and vortex instability are sensitive to inflow boundary conditions. These data indicate that previous findings on vortex instability by using generalized inflow boundary conditions should be considered with caution.³⁸

Both patient-specific and generalized inflow boundary conditions are traditionally used for CFD analysis of cerebral aneurysm hemodynamics. However, only 2 authors have so far compared both approaches.^28⇓–30 The large differences in WSS magnitude and WSS distribution found in our study are in line with findings of a smaller study in 6 patients by Karmonik et al,²⁸ who pointed out differences in both WSS values and in WSS distribution between the 2 inflow boundary conditions. In addition, Venugopal et al³⁹ showed that WSS distributions are sensitive to changes in flow-rate distribution in the proximal artery. Marzo et al³⁰ also reported differences in WSS magnitude but did not report changes in WSS distribution outside the order of physiologic variations.

Twenty-one of 36 aneurysms showed at least 1 difference in a hemodynamic characteristic for the 2 inflow boundary conditions. Most of these were WSS-related, such as high WSS location and level of WSS on the lobulation sac. Fewer differences were found in vortex- or inflow jet–related characteristics. This discrepancy could be explained by findings of Cebral et al,³⁶ who have shown that variations of up to 25% in blood flow rate do not affect the flow patterns inside the aneurysm.

There are several limitations related to the design of this study. It did not address the role of outflow boundary conditions and flow divisions distal to the aneurysm, which could influence the accuracy of the results. In addition, for patients with a ruptured aneurysm, postprocedural PC-MR imaging velocity measurements were used for the inflow boundary conditions to simulate preprocedural hemodynamics. The embolization itself may alter the local hemodynamics, resulting in an inadequate representation of pre-embolization hemodynamics. However, we have compared pre- and postprocedural PC-MR imaging measurements of the parent artery in a small number of patients (with unruptured aneurysms) and did not detect notable differences. Finally, vascular compliance and fluid-wall interaction were not incorporated in this model. Disregarding these effects may have influenced the resulting hemodynamic characteristics and may have led to an overestimation of mean and maximum WSS values.¹⁰

The findings of this study emphasize that running simulations with generalized boundary conditions may result in variations in WSS magnitude and distribution and may overestimate vortex instability. These variations could greatly influence the association of hemodynamics with the rupture of cerebral aneurysms. Therefore, interpretation of WSS profiles should be applied with great caution when generalized inflow boundary conditions are used for this purpose. More research is also needed to further investigate the influence of other boundary conditions on CFD in aneurysms, ideally incorporating the effect of wall and outflow boundary conditions.

Conclusions

Patient-specific and generalized inflow boundary conditions in CFD-based simulations of aneurysmal hemodynamics resulted in large differences in WSS magnitudes. In addition, 21 of 36 aneurysms showed differences in hemodynamics characteristics. Aneurysms showing differences in vortex and inflow jet characteristics were more often assessed as unstable when generalized inflow boundary conditions were applied. This study emphasizes the necessity of the use of patient-specific inflow boundary conditions for the calculation of hemodynamics in cerebral aneurysms by using CFD techniques.

Acknowledgments

This study was supported by a grant from the Nuts Ohra Foundation, Amsterdam, The Netherlands, and the Dutch technology foundation STW under project number 11629.

Footnotes

H.A. Marquering and C.B.L.M. Majoie shared senior authorship of this article.
Disclosures: Wouter V. Potters—RELATED: Other: Dutch Technology Foundation STW (Carisma 111629),* Comments: Salary was paid by a grant from a government granting agency (the Dutch Technology Foundation STW) for this submitted work. Charles B.L. Majoie—RELATED: Grant: Nuts Ohra Foundation,* UNRELATED: Grants/Grants Pending: Dutch Heart Foundation.* *Money paid to the institution.
This study was supported by a grant from the Nuts Ohra Foundation, Amsterdam, the Netherlands.

Indicates open access to non-subscribers at www.ajnr.org

REFERENCES

1.↵
1. Brisman JL,
2. Song JK,
3. Newell DW
. Cerebral aneurysms. N Engl J Med 2006;355:928–39
CrossRef PubMed
2.↵
1. Taylor TN
. The medical economics of stroke. Drugs 1997;54(suppl 3):51–57
CrossRef PubMed
3.↵
1. Taylor TN,
2. Davis PH,
3. Torner JC,
4. et al
. Lifetime cost of stroke in the United States. Stroke 1996;27:1459–66
Abstract/FREE Full Text
4.↵
1. Raymond J,
2. Guilbert F,
3. Weill A,
4. et al
. Long-term angiographic recurrences after selective endovascular treatment of aneurysms with detachable coils. Stroke 2003;34:1398–403
Abstract/FREE Full Text
5.↵
1. Murayama Y,
2. Nien YL,
3. Duckwiler G,
4. et al
. Guglielmi detachable coil embolization of cerebral aneurysms: 11 years' experience. J Neurosurg 2003;98:959–66
CrossRef PubMed
6.↵
The International Study of Unruptured Intracranial Aneurysms Investigators. Unruptured intracranial aneurysms: risk of rupture and risks of surgical interventions. N Engl J Med 1998;339:1725–33
CrossRef PubMed
7.↵
1. Tsutsumi K,
2. Ueki K,
3. Morita A,
4. et al
. Risk of rupture from incidental cerebral aneurysms. J Neurosurg 2000;93:550–53
CrossRef PubMed
8.↵
1. Juvela S,
2. Porras M,
3. Poussa K
. Natural history of unruptured intracranial aneurysms: probability of and risk factors for aneurysm rupture. J Neurosurg 2000;93:379–87
CrossRef PubMed
9.↵
1. Metaxa E,
2. Tremmel M,
3. Natarajan SK,
4. et al
. Characterization of critical hemodynamics contributing to aneurysmal remodeling at the basilar terminus in a rabbit model. Stroke 2010;41:1774–82
Abstract/FREE Full Text
10.↵
1. Jou LD,
2. Quick CM,
3. Young WL,
4. et al
. Computational approach to quantifying hemodynamic forces in giant cerebral aneurysms. AJNR Am J Neuroradiol 2003;24:1804–10
Abstract/FREE Full Text
11.↵
1. Meng H,
2. Tutino VM,
3. Xiang J,
4. et al
. High WSS or low WSS? Complex interactions of hemodynamics with intracranial aneurysm initiation, growth, and rupture: toward a unifying hypothesis. AJNR Am J Neuroradiol 2014;35:1254–62
Abstract/FREE Full Text
12.↵
1. Boussel L,
2. Rayz V,
3. Mcculloch C,
4. et al
. Aneurysm growth occurs at region of low wall shear stress: patient-specific correlation of hemodynamics and growth in a longitudinal study. Stroke 2008;39:2997–3002
Abstract/FREE Full Text
13.↵
1. Xiang J,
2. Natarajan SK,
3. Tremmel M,
4. et al
. Hemodynamic-morphologic discriminants for intracranial aneurysm rupture. Stroke 2011;42:144–52
Abstract/FREE Full Text
14.↵
1. Ford MD,
2. Alperin N
. Characterization of volumetric flow rate waveforms in the normal internal carotid and vertebral arteries. Physiol Meas 2005;26:477–88
CrossRef PubMed
15.↵
1. Castro MA,
2. Putman CM,
3. Cebral JR
. Computational fluid dynamics modeling of intracranial aneurysms: effects of parent artery segmentation on intra-aneurysmal hemodynamics. AJNR Am J Neuroradiol 2006;27:1703–09
Abstract/FREE Full Text
16.↵
1. Castro MA,
2. Putman CM,
3. Cebral JR
. Patient-specific computational fluid dynamics modelling of anterior communicating artery aneurysms: a study of sensitivity of intra-aneurysmal flow patterns to flow conditions in the carotid arteries. AJNR Am J Neuroradiol 2006;27:2061–68
Abstract/FREE Full Text
17.↵
1. Cebral J,
2. Sheridan M,
3. Putman CM
. Hemodynamics and bleb formation in intracranial aneurysms. AJNR Am J Neuroradiol 2010;31:304–10
Abstract/FREE Full Text
18.↵
1. Cebral JR,
2. Castro MA,
3. Burgess JE,
4. et al
. Characterization of cerebral aneurysms for assessing risk of rupture by using patient-specific computational hemodynamics models. AJNR Am J Neuroradiol 2005;26:2550–59
Abstract/FREE Full Text
19.↵
1. Cebral JR,
2. Löhner R
. Efficient simulation of blood flow past complex endovascular devices using an adaptive embedding technique. IEEE Trans Med Imaging 2005;24:468–76
CrossRef PubMed
20.↵
1. Cebral JR,
2. Pergolizzi RS,
3. Putman CM
. Computational fluid dynamics modeling of intracranial aneurysms: qualitative comparison with cerebral angiography. Acad Radiol 2007;14:804–13
CrossRef PubMed
21.↵
1. Dempere-Marco L,
2. Oubel E,
3. Castro M,
4. et al
. CFD analysis incorporating the influence of wall motion: application to intracranial aneurysms. Med Image Comput Comput Assist Interv 2006;9(pt 2):438–45
PubMed
22.↵
1. Ford MD,
2. Stuhne GR,
3. Nikolov HN,
4. et al
. Virtual angiography for visualization and validation of computational models of aneurysm hemodynamics. IEEE Trans Med Imaging 2005;24:1586–92
CrossRef PubMed
23.↵
1. Hoi Y,
2. Woodward SH,
3. Kim M,
4. et al
. Validation of CFD simulations of cerebral aneurysms with implication of geometric variations. J Biomech Eng 2006;128:844–51
CrossRef PubMed
24.↵
1. Imai Y,
2. Sato K,
3. Ishikawa T,
4. et al
. Inflow into saccular cerebral aneurysms at arterial bends. Ann Biomed Eng 2008;36:1489–95
CrossRef PubMed
25.↵
1. Mantha A,
2. Karmonik C,
3. Benndorf G,
4. et al
. Hemodynamics in a cerebral artery before and after the formation of an aneurysm. AJNR Am J Neuroradiol 2006;27:1113–18
Abstract/FREE Full Text
26.↵
1. Shojima M,
2. Nemoto S,
3. Morita A,
4. et al
. Role of shear stress in the blister formation of cerebral aneurysms. Neurosurgery 2010;67:1268–74
CrossRef PubMed
27.↵
1. Shojima M,
2. Oshima M,
3. Takagi K,
4. et al
. Magnitude and role of wall shear stress on cerebral aneurysm: computational fluid dynamic study of 20 middle cerebral artery aneurysms. Stroke 2004;35:2500–05
Abstract/FREE Full Text
28.↵
1. Karmonik C,
2. Yen C,
3. Diaz O,
4. et al
. Temporal variations of wall shear stress parameters in intracranial aneurysms: importance of patient-specific inflow waveforms for CFD calculations. Acta Neurochir (Wien) 2010;152:1391–98, discussion 1398
CrossRef PubMed
29.↵
1. Karmonik C,
2. Yen C,
3. Grossman RG,
4. et al
. Intra-aneurysmal flow patterns and wall shear stresses calculated with computational flow dynamics in an anterior communicating artery aneurysm depend on knowledge of patient-specific inflow rates. Acta Neurochir (Wien) 2009;151:479–85, discussion 485
CrossRef PubMed
30.↵
1. Marzo A,
2. Singh P,
3. Larrabide I,
4. et al
. Computational hemodynamics in cerebral aneurysms: the effects of modeled versus measured boundary conditions. Ann Biomed Eng 2011;39:884–96
CrossRef PubMed
31.↵
1. Omodaka S,
2. Sugiyama S,
3. Inoue T,
4. et al
. Local hemodynamics at the rupture point of cerebral aneurysms determined by computational fluid dynamics analysis. Cerebrovasc Dis 2012;34:121–29
CrossRef PubMed
32.↵
1. Hassan T,
2. Ezura M,
3. Timofeev EV,
4. et al
. Computational simulation of therapeutic parent artery occlusion to treat giant vertebrobasilar aneurysm. AJNR Am J Neuroradiol 2004;25:63–68
Abstract/FREE Full Text
33.↵
1. Mitchell P,
2. Kerr R,
3. Mendelow AD,
4. et al
. Could late rebleeding overturn the superiority of cranial aneurysm coil embolization over clip ligation seen in the International Subarachnoid Aneurysm Trial? J Neurosurg 2008;108:437–42
CrossRef PubMed
34.↵
1. van Ooij P,
2. Schneiders JJ,
3. Marquering HA,
4. et al
. 3D cine phase-contrast MRI at 3T in intracranial aneurysms compared with patient-specific computational fluid dynamics. AJNR Am J Neuroradiol 2013;34:1785–91
Abstract/FREE Full Text
35.↵
1. Steinman DA,
2. Hoi Y,
3. Fahy P,
4. et al
. Variability of computational fluid dynamics solutions for pressure and flow in a giant aneurysm: the ASME 2012 Summer Bioengineering Conference CFD Challenge. J Biomech Eng 2013;135:021016
36.↵
1. Cebral JR,
2. Castro MA,
3. Appanaboyina S,
4. et al
. Efficient pipeline for image-based patient-specific analysis of cerebral aneurysm hemodynamics: technique and sensitivity. IEEE Trans Med Imaging 2005;24:457–67
CrossRef PubMed
37.↵
1. Patankar SV,
2. Spalding DB
. A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows. Int J Heat Mass Transf 1972;15:1787–806
CrossRef
38.↵
1. Byrne G,
2. Mut F,
3. Cebral JR
. Quantifying the large-scale hemodynamics of intracranial aneurysms. AJNR Am J Neuroradiol 2014;35:333–38
Abstract/FREE Full Text
39.↵
1. Venugopal P,
2. Valentino D,
3. Schmitt H,
4. et al
. Sensitivity of patient-specific numerical simulation of cerebral aneurysm hemodynamics to inflow boundary conditions. J Neurosurg 2007;106:1051–60
CrossRef PubMed
40.
1. González-Alonso J,
2. Dalsgaard MK,
3. Osada T,
4. et al
. Brain and central haemodynamics and oxygenation during maximal exercise in humans. J Physiol 2004;557(pt 1):331–42
CrossRef PubMed

Received November 15, 2013.
Accepted after revision December 30, 2013.

View Abstract

[1] 1.↵
Brisman JL,
Song JK,
Newell DW
. Cerebral aneurysms. N Engl J Med 2006;355:928–39
CrossRef PubMed

[2] Brisman JL,

[3] Song JK,

[4] Newell DW

[5] 2.↵
Taylor TN
. The medical economics of stroke. Drugs 1997;54(suppl 3):51–57
CrossRef PubMed

[6] Taylor TN

[7] 3.↵
Taylor TN,
Davis PH,
Torner JC,
et al
. Lifetime cost of stroke in the United States. Stroke 1996;27:1459–66
Abstract/FREE Full Text

[8] Taylor TN,

[9] Davis PH,

[10] Torner JC,

[11] et al

[12] 4.↵
Raymond J,
Guilbert F,
Weill A,
et al
. Long-term angiographic recurrences after selective endovascular treatment of aneurysms with detachable coils. Stroke 2003;34:1398–403
Abstract/FREE Full Text

[13] Raymond J,

[14] Guilbert F,

[15] Weill A,

[16] et al

[17] 5.↵
Murayama Y,
Nien YL,
Duckwiler G,
et al
. Guglielmi detachable coil embolization of cerebral aneurysms: 11 years' experience. J Neurosurg 2003;98:959–66
CrossRef PubMed

[18] Murayama Y,

[19] Nien YL,

[20] Duckwiler G,

[21] et al

[22] 6.↵
The International Study of Unruptured Intracranial Aneurysms Investigators. Unruptured intracranial aneurysms: risk of rupture and risks of surgical interventions. N Engl J Med 1998;339:1725–33
CrossRef PubMed

[23] 7.↵
Tsutsumi K,
Ueki K,
Morita A,
et al
. Risk of rupture from incidental cerebral aneurysms. J Neurosurg 2000;93:550–53
CrossRef PubMed

[24] Tsutsumi K,

[25] Ueki K,

[26] Morita A,

[27] et al

[28] 8.↵
Juvela S,
Porras M,
Poussa K
. Natural history of unruptured intracranial aneurysms: probability of and risk factors for aneurysm rupture. J Neurosurg 2000;93:379–87
CrossRef PubMed

[29] Juvela S,

[30] Porras M,

[31] Poussa K

[32] 9.↵
Metaxa E,
Tremmel M,
Natarajan SK,
et al
. Characterization of critical hemodynamics contributing to aneurysmal remodeling at the basilar terminus in a rabbit model. Stroke 2010;41:1774–82
Abstract/FREE Full Text

[33] Metaxa E,

[34] Tremmel M,

[35] Natarajan SK,

[36] et al

[37] 10.↵
Jou LD,
Quick CM,
Young WL,
et al
. Computational approach to quantifying hemodynamic forces in giant cerebral aneurysms. AJNR Am J Neuroradiol 2003;24:1804–10
Abstract/FREE Full Text

[38] Jou LD,

[39] Quick CM,

[40] Young WL,

[41] et al

[42] 11.↵
Meng H,
Tutino VM,
Xiang J,
et al
. High WSS or low WSS? Complex interactions of hemodynamics with intracranial aneurysm initiation, growth, and rupture: toward a unifying hypothesis. AJNR Am J Neuroradiol 2014;35:1254–62
Abstract/FREE Full Text

[43] Meng H,

[44] Tutino VM,

[45] Xiang J,

[46] et al

[47] 12.↵
Boussel L,
Rayz V,
Mcculloch C,
et al
. Aneurysm growth occurs at region of low wall shear stress: patient-specific correlation of hemodynamics and growth in a longitudinal study. Stroke 2008;39:2997–3002
Abstract/FREE Full Text

[48] Boussel L,

[49] Rayz V,

[50] Mcculloch C,

[51] et al

[52] 13.↵
Xiang J,
Natarajan SK,
Tremmel M,
et al
. Hemodynamic-morphologic discriminants for intracranial aneurysm rupture. Stroke 2011;42:144–52
Abstract/FREE Full Text

[53] Xiang J,

[54] Natarajan SK,

[55] Tremmel M,

[56] et al

[57] 14.↵
Ford MD,
Alperin N
. Characterization of volumetric flow rate waveforms in the normal internal carotid and vertebral arteries. Physiol Meas 2005;26:477–88
CrossRef PubMed

[58] Ford MD,

[59] Alperin N

[60] 15.↵
Castro MA,
Putman CM,
Cebral JR
. Computational fluid dynamics modeling of intracranial aneurysms: effects of parent artery segmentation on intra-aneurysmal hemodynamics. AJNR Am J Neuroradiol 2006;27:1703–09
Abstract/FREE Full Text

[61] Castro MA,

[62] Putman CM,

[63] Cebral JR

[64] 16.↵
Castro MA,
Putman CM,
Cebral JR
. Patient-specific computational fluid dynamics modelling of anterior communicating artery aneurysms: a study of sensitivity of intra-aneurysmal flow patterns to flow conditions in the carotid arteries. AJNR Am J Neuroradiol 2006;27:2061–68
Abstract/FREE Full Text

[65] Castro MA,

[66] Putman CM,

[67] Cebral JR

[68] 17.↵
Cebral J,
Sheridan M,
Putman CM
. Hemodynamics and bleb formation in intracranial aneurysms. AJNR Am J Neuroradiol 2010;31:304–10
Abstract/FREE Full Text

[69] Cebral J,

[70] Sheridan M,

[71] Putman CM

[72] 18.↵
Cebral JR,
Castro MA,
Burgess JE,
et al
. Characterization of cerebral aneurysms for assessing risk of rupture by using patient-specific computational hemodynamics models. AJNR Am J Neuroradiol 2005;26:2550–59
Abstract/FREE Full Text

[73] Cebral JR,

[74] Castro MA,

[75] Burgess JE,

[76] et al

[77] 19.↵
Cebral JR,
Löhner R
. Efficient simulation of blood flow past complex endovascular devices using an adaptive embedding technique. IEEE Trans Med Imaging 2005;24:468–76
CrossRef PubMed

[78] Cebral JR,

[79] Löhner R

[80] 20.↵
Cebral JR,
Pergolizzi RS,
Putman CM
. Computational fluid dynamics modeling of intracranial aneurysms: qualitative comparison with cerebral angiography. Acad Radiol 2007;14:804–13
CrossRef PubMed

[81] Cebral JR,

[82] Pergolizzi RS,

[83] Putman CM

[84] 21.↵
Dempere-Marco L,
Oubel E,
Castro M,
et al
. CFD analysis incorporating the influence of wall motion: application to intracranial aneurysms. Med Image Comput Comput Assist Interv 2006;9(pt 2):438–45
PubMed

[85] Dempere-Marco L,

[86] Oubel E,

[87] Castro M,

[88] et al

[89] 22.↵
Ford MD,
Stuhne GR,
Nikolov HN,
et al
. Virtual angiography for visualization and validation of computational models of aneurysm hemodynamics. IEEE Trans Med Imaging 2005;24:1586–92
CrossRef PubMed

[90] Ford MD,

[91] Stuhne GR,

[92] Nikolov HN,

[93] et al

[94] 23.↵
Hoi Y,
Woodward SH,
Kim M,
et al
. Validation of CFD simulations of cerebral aneurysms with implication of geometric variations. J Biomech Eng 2006;128:844–51
CrossRef PubMed

[95] Hoi Y,

[96] Woodward SH,

[97] Kim M,

[98] et al

[99] 24.↵
Imai Y,
Sato K,
Ishikawa T,
et al
. Inflow into saccular cerebral aneurysms at arterial bends. Ann Biomed Eng 2008;36:1489–95
CrossRef PubMed

[100] Imai Y,

[101] Sato K,

[102] Ishikawa T,

[103] et al

[104] 25.↵
Mantha A,
Karmonik C,
Benndorf G,
et al
. Hemodynamics in a cerebral artery before and after the formation of an aneurysm. AJNR Am J Neuroradiol 2006;27:1113–18
Abstract/FREE Full Text

[105] Mantha A,

[106] Karmonik C,

[107] Benndorf G,

[108] et al

[109] 26.↵
Shojima M,
Nemoto S,
Morita A,
et al
. Role of shear stress in the blister formation of cerebral aneurysms. Neurosurgery 2010;67:1268–74
CrossRef PubMed

[110] Shojima M,

[111] Nemoto S,

[112] Morita A,

[113] et al

[114] 27.↵
Shojima M,
Oshima M,
Takagi K,
et al
. Magnitude and role of wall shear stress on cerebral aneurysm: computational fluid dynamic study of 20 middle cerebral artery aneurysms. Stroke 2004;35:2500–05
Abstract/FREE Full Text

[115] Shojima M,

[116] Oshima M,

[117] Takagi K,

[118] et al

[119] 28.↵
Karmonik C,
Yen C,
Diaz O,
et al
. Temporal variations of wall shear stress parameters in intracranial aneurysms: importance of patient-specific inflow waveforms for CFD calculations. Acta Neurochir (Wien) 2010;152:1391–98, discussion 1398
CrossRef PubMed

[120] Karmonik C,

[121] Yen C,

[122] Diaz O,

[123] et al

[124] 29.↵
Karmonik C,
Yen C,
Grossman RG,
et al
. Intra-aneurysmal flow patterns and wall shear stresses calculated with computational flow dynamics in an anterior communicating artery aneurysm depend on knowledge of patient-specific inflow rates. Acta Neurochir (Wien) 2009;151:479–85, discussion 485
CrossRef PubMed

[125] Karmonik C,

[126] Yen C,

[127] Grossman RG,

[128] et al

[129] 30.↵
Marzo A,
Singh P,
Larrabide I,
et al
. Computational hemodynamics in cerebral aneurysms: the effects of modeled versus measured boundary conditions. Ann Biomed Eng 2011;39:884–96
CrossRef PubMed

[130] Marzo A,

[131] Singh P,

[132] Larrabide I,

[133] et al

[134] 31.↵
Omodaka S,
Sugiyama S,
Inoue T,
et al
. Local hemodynamics at the rupture point of cerebral aneurysms determined by computational fluid dynamics analysis. Cerebrovasc Dis 2012;34:121–29
CrossRef PubMed

[135] Omodaka S,

[136] Sugiyama S,

[137] Inoue T,

[138] et al

[139] 32.↵
Hassan T,
Ezura M,
Timofeev EV,
et al
. Computational simulation of therapeutic parent artery occlusion to treat giant vertebrobasilar aneurysm. AJNR Am J Neuroradiol 2004;25:63–68
Abstract/FREE Full Text

[140] Hassan T,

[141] Ezura M,

[142] Timofeev EV,

[143] et al

[144] 33.↵
Mitchell P,
Kerr R,
Mendelow AD,
et al
. Could late rebleeding overturn the superiority of cranial aneurysm coil embolization over clip ligation seen in the International Subarachnoid Aneurysm Trial? J Neurosurg 2008;108:437–42
CrossRef PubMed

[145] Mitchell P,

[146] Kerr R,

[147] Mendelow AD,

[148] et al

[149] 34.↵
van Ooij P,
Schneiders JJ,
Marquering HA,
et al
. 3D cine phase-contrast MRI at 3T in intracranial aneurysms compared with patient-specific computational fluid dynamics. AJNR Am J Neuroradiol 2013;34:1785–91
Abstract/FREE Full Text

[150] van Ooij P,

[151] Schneiders JJ,

[152] Marquering HA,

[153] et al

[154] 35.↵
Steinman DA,
Hoi Y,
Fahy P,
et al
. Variability of computational fluid dynamics solutions for pressure and flow in a giant aneurysm: the ASME 2012 Summer Bioengineering Conference CFD Challenge. J Biomech Eng 2013;135:021016

[155] Steinman DA,

[156] Hoi Y,

[157] Fahy P,

[158] et al

[159] 36.↵
Cebral JR,
Castro MA,
Appanaboyina S,
et al
. Efficient pipeline for image-based patient-specific analysis of cerebral aneurysm hemodynamics: technique and sensitivity. IEEE Trans Med Imaging 2005;24:457–67
CrossRef PubMed

[160] Cebral JR,

[161] Castro MA,

[162] Appanaboyina S,

[163] et al

[164] 37.↵
Patankar SV,
Spalding DB
. A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows. Int J Heat Mass Transf 1972;15:1787–806
CrossRef

[165] Patankar SV,

[166] Spalding DB

[167] 38.↵
Byrne G,
Mut F,
Cebral JR
. Quantifying the large-scale hemodynamics of intracranial aneurysms. AJNR Am J Neuroradiol 2014;35:333–38
Abstract/FREE Full Text

[168] Byrne G,

[169] Mut F,

[170] Cebral JR

[171] 39.↵
Venugopal P,
Valentino D,
Schmitt H,
et al
. Sensitivity of patient-specific numerical simulation of cerebral aneurysm hemodynamics to inflow boundary conditions. J Neurosurg 2007;106:1051–60
CrossRef PubMed

[172] Venugopal P,

[173] Valentino D,

[174] Schmitt H,

[175] et al

[176] 40.
González-Alonso J,
Dalsgaard MK,
Osada T,
et al
. Brain and central haemodynamics and oxygenation during maximal exercise in humans. J Physiol 2004;557(pt 1):331–42
CrossRef PubMed

[177] González-Alonso J,

[178] Dalsgaard MK,

[179] Osada T,

[180] et al

Main menu

User menu

Search

American Journal of Neuroradiology

Generalized versus Patient-Specific Inflow Boundary Conditions in Computational Fluid Dynamics Simulations of Cerebral Aneurysmal Hemodynamics

Abstract

ABBREVIATIONS:

Materials and Methods

Patient Selection

Imaging

Geometric Vascular Models

Computational Fluid Dynamics

Quantitative Assessment of Hemodynamic Features

Qualitative Assessment of Hemodynamic Features

Results

Quantitative Assessment of Hemodynamic Features

Qualitative Assessment of Hemodynamic Features

Discussion

Conclusions

Acknowledgments

Footnotes

REFERENCES

In this issue

Citation Manager Formats

Related Articles

Cited By...

More in this TOC Section

Similar Articles

News and Updates

Resources

Opportunities

American Society of Neuroradiology

Main menu

User menu

Search

Generalized versus Patient-Specific Inflow Boundary Conditions in Computational Fluid Dynamics Simulations of Cerebral Aneurysmal Hemodynamics

Abstract

ABBREVIATIONS:

Materials and Methods

Patient Selection

Imaging

Geometric Vascular Models

Computational Fluid Dynamics

Quantitative Assessment of Hemodynamic Features

Qualitative Assessment of Hemodynamic Features

Results

Quantitative Assessment of Hemodynamic Features

Qualitative Assessment of Hemodynamic Features

Discussion

Conclusions

Acknowledgments

Footnotes

REFERENCES

In this issue

Citation Manager Formats

Jump to section

Related Articles

Cited By...

More in this TOC Section

Similar Articles

News and Updates

Resources

Opportunities

American Society of Neuroradiology