(Stroke. 2004;35:2500.)
© 2004 American Heart Association, Inc.
Original Contributions |
From the Department of Neurosurgery (M.S., A.M., T.K.), Faculty of Medicine, and the Institute of Industrial Science (M.S., M.O., R.T.), University of Tokyo, Japan; the Department of Neurosurgery (K.T.), Faculty of Medicine, Teikyo University, Japan; and the Departments of Neurosurgery (M.H.) and Radiology (K.K.), Faculty of Medicine, Fujita Health University, Japan.
Correspondence to Dr Masaaki Shojima, Department of Neurosurgery, Faculty of Medicine, University of Tokyo, Japan. 7-3-1, Hongo, Bunkyo-ku, Tokyo, Japan 113-0033. E-mail mshoji-tky{at}umin.ac.jp
| Abstract |
|---|
|
|
|---|
Methods Twenty mathematical models of MCA vessels with aneurysms were created by 3-dimensional computed tomographic angiography. CFD calculations were performed by using our original finite-element solver with the assumption of Newtonian fluid property for blood and the rigid wall property for the vessel and the aneurysm.
Results The maximum WSS in the calculated region tended to occur near the neck of the aneurysm, not in its tip or bleb. The magnitude of the maximum WSS was 14.39±6.21 N/m2, which was 4-times higher than the average WSS in the vessel region (3.64±1.25 N/m2). The average WSS of the aneurysm region (1.64±1.16 N/m2) was significantly lower than that of the vessel region (P<0.05). The WSSs at the tip of ruptured aneurysms were markedly low.
Conclusions These results suggest that in contrast to the pathogenic effect of a high WSS in the initiating phase, a low WSS may facilitate the growing phase and may trigger the rupture of a cerebral aneurysm by causing degenerative changes in the aneurysm wall. The WSS of the aneurysm region may be of some help for the prediction of rupture.
Key Words: aneurysm biomechanics hemodynamics shear strength
| Introduction |
|---|
|
|
|---|
Hemodynamic stresses are considered to have profound effects on the development of cerebral aneurysms.7,8 One of these, the wall shear stress (WSS), acts directly on the vascular endothelium as a biological stimulator that modulates the cellular function of the endothelium.9,10 Thus, the focus of the study presented here was aneurysm WSS. The close relationship between high WSS and the initiation of cerebral aneurysm formation has already been demonstrated in animal experiments.11 The WSS may also play an important role in the growth and rupture of cerebral aneurysms.
The measurement of WSS in vivo is becoming feasible; however, it remains very difficult, especially in small and tortuous intracranial arteries.12,13 With recent advances in computer technology, the magnitude and distribution of WSSs in complex arterial models have been observed with the aid of the computational fluid dynamics (CFD) technique, which is also a useful clinical tool for planning endovascular treatment.14,15 The use of this technique has been limited to just a few cases; therefore, a statistical analysis of the results has not been possible. Here, we present the results of a statistical analysis of the magnitude of WSSs in and around saccular cerebral aneurysms.
| Materials and Methods |
|---|
|
|
|---|
From January 2001 to December 2002, there were 42 middle cerebral artery (MCA) aneurysms that were diagnosed with 3-dimensional (3-D) computed tomographic (CT) angiography in 40 patients. Among them, 20 aneurysms of 19 patients had an adequate image quality for CFD calculation and were analyzed in this study. Patient population consisted of 7 men and 12 women, with a mean age of 61.5 years (range, 51 to 75 years). Three aneurysms were ruptured, and one of them was accompanied with an unruptured aneurysm. The other 16 aneurysms were detected before rupture by the screening examinations. All the aneurysm studied here were saccular aneurysms originating at the first major bifurcation of MCA.
Modeling of the Aneurysms
Three-dimensional CT angiography data with a voxel size of 0.21 mmx0.21 mmx0.50 mm were obtained with the aid of a multislice CT system,16 the Aquilion multi 16 (Toshiba). The 16 central rows of 0.5-mm detector elements that were used had the following parameters: 0.75-second rotation, a scanning pitch of 0.69, 135 kV, and 260 mA. A total of 80 mL nonionic contrast medium (300 mg/mL) was injected at a rate of 3 mL/s via the median cubital vein. Digital images were transferred to a Unix workstation and vessel surfaces were constructed with the Fly-through mode of Alatoview (Toshiba). The fine irregularities of the original models resulting from partial volume effects or slice gaps were refined without changing the comprehensive geometry using our original smoothing software based on the algorithm of Garland mesh simplification17 and Taubin mesh smoothing.18 The aneurysm and the 20 mm of vessel surrounding it were trimmed out for the analysis. When possible, an intersecting plane dividing the aneurysm volume from that of the parent artery was made to allow a comparison of the WSS of the aneurysm and the vessel (Figure 1). Computational meshes were generated for these models with
60 000 hexahedral elements.
|
Numerical Simulation
CFD simulations were performed using our original finite-element solver under the governing equations of mass conservation and NavieStokes.1921 The boundary conditions were applied as follows. Blood was assumed to be an incompressible isothermal Newtonian fluid22 with a specific gravity of 1000 kg/m3 and a viscosity of 4.0x103 N/m2 per second. The viscoelastic properties of the vessel wall were neglected and a rigid wall with no-slip condition was applied.15 For the inlet condition, a pulsatile flow with a Womersley velocity profile was simulated,23 with the typical MCA velocity obtained by transcranial Doppler scanning (mean velocity, 0.6 m/s; maximum velocity, 0.81 m/s; heart rate, 80 bpm). This inlet velocity condition was applied to all aneurysms (mean Reynolds number, 413; mean Womersley number, 3.99). A traction-free boundary condition24 was applied to the outlets. The width of the time step for calculation was set at 0.0001 seconds.
To confirm numerical stability, the calculation was performed for 5 cardiac cycles and the result at the fifth cardiac cycle was used for the analysis. The calculation time for one aneurysm was
48 hours. WSS distributions were calculated from the 3-D velocity field data.21
Statistical Analysis
The maximum WSS region and its value were recorded for all aneurysms. When it was possible to divide the aneurysm from the parent artery with an intersecting plane, we calculated the spatially averaged WSS for each region at the peak systole. The WSS of the aneurysm was then compared with that of the vessel (paired Student t test). A nonpaired Student t test was applied to a comparison of the WSS of the ruptured aneurysms with that of the unruptured aneurysms.
For each aneurysm, the diameter of the inlet of the aneurysm and the maximum height of the aneurysmal sac from the inlet plane were measured from the 3-D model, and the aspect ratio was determined by dividing the latter by the former. Pearson correlation coefficients were calculated among the aspect ratio, the WSS of the aneurysm region, and the volumetric flow into the aneurysm. Statistical significance was taken as P<0.05.
| Results |
|---|
|
|
|---|
|
Intersecting planes dividing the aneurysm from the parent artery were created successfully for all but 3 small unruptured aneurysms. For these 17 cases, the WSS of the vessel region (spatially averaged at the peak systole) was 3.64±1.25 N/m2. In contrast, the WSS of the aneurysm region (spatially averaged at the peak systole) was 1.64±1.16 N/m2, which is significantly lower than that of the vessel (P<0.05; Figure 3).
|
Flow Structure in the Aneurysms
The shape of the aneurysm had a profound impact on the flow structure within it. In case 1 (Figure 4, left), a small aneurysm with a smooth contour and a low aspect ratio of 0.56, the intra-aneurysmal flow velocity did not decrease particularly. The averaged WSS of the aneurysm region at the peak systole was 3.35 N/m2, one of the highest among our cases. However, in case 2, which had a large aneurysm with an aspect ratio of 1.58 (Figure 4, right), the intra-aneurysmal flow velocity was markedly low and the flow field showed a conspicuous recirculating zone. In this case, the spatially averaged WSS of the aneurysm region was markedly low (0.14 N/m2) even at the peak systole because of this large area of the flow stasis.
|
The aspect ratio of the aneurysm had a significant negative correlation (r=0.67, P<0.05; Figure 5) with the WSS of the aneurysm region (spatially averaged at the peak systole). A mild positive correlation was observed between the spatially averaged WSS of the aneurysm region and the volumetric flow into the aneurysm (r=0.50, P=0.06). The correlation between the volumetric flow into the aneurysm and the aspect ratio of the aneurysm was weak (r=0.36, P=0.06).
|
Ruptured Versus Unruptured Aneurysms
The mean size (diameter) and the aspect ratio of the aneurysms were 3.36 mm and 0.73, respectively, for the ruptured cases (n=3) and 4.31 mm and 0.92, respectively, for the unruptured cases (n=17). The difference between the respective figure for the ruptured and unruptured cases was not statistically significant.
When the spatially averaged WSS of aneurysm region at the peak systole was compared between ruptured and unruptured cases, it was found to be significantly higher for ruptured cases (2.92 N/m2 versus 1.48 N/m2, P<0.05).
In all ruptured cases, high and low WSS were mixed in the small aneurysm area (Figure 6). The blood of parent artery flowed into the aneurysm more directly in the ruptured cases and high WSSs appeared at the body or the neck of aneurysm. This resulted in a higher averaged WSS of the aneurysm region in ruptured cases. However, at the tip, the stasis of the blood with recirculating zones was observed. This caused markedly low WSS at the tip (<0.5 N/m2). These findings make fine contrast to the flow structure of the unruptured aneurysm, shown in the left side of Figure 4, which also had a high averaged WSS of the aneurysm region. The intra-aneurysmal velocity was not delayed, particularly near the wall. The WSS at the tip of this unruptured aneurysm was not decreased severely (1.7 N/m2).
|
| Discussion |
|---|
|
|
|---|
WSS is a flow-induced stress that can be described as the frictional force of viscous blood.10 The 3-D geometry and the 3-D velocity field of vessels are indispensable to the establishment of spatial distribution of WSS and the flow structure. In vivo measurements of the 3-D velocity field and WSS have become possible with the development of phase contrast magnetic resonance velocimetry for large and simple arteries like the aorta25 or the carotid bifurcation.13 However, in small and tortuous vessels like the intracranial arteries, the currently available techniques, like phase contrast magnetic resonance velocimetry, cannot be applied in calculating the spatial distribution of the WSS. To investigate the flow dynamics in cerebral arteries, simulation in vitro or by computer is necessary. The major difference between in vitro fluid experiments and CFD simulations may be the fact that the quality of computational mesh generation has some effect on the results of computer simulations. However, the reliability of computer simulations with proper mesh generation has been established.24 Comparisons between our results and those of in vitro experiments are ongoing in our laboratory.20
Recent studies have indicated the involvement of WSS in the formation of saccular cerebral aneurysms.11 A prolonged high WSS fragments the internal elastic lamina of vessels26 and gives rise to the initial change involved in the formation of a cerebral aneurysm. Our results have established that the magnitude of the WSS of well-developed aneurysms is very low, in accordance with the previous hypothesis that the strength of the WSS of the aneurysm region is not sufficient to mechanically tear the wall of the aneurysm.7 The WSS is converted to biological signals via mechanoreceptors on endothelial cells, and it modulates gene expressions and the cellular functions of the vessel wall.9,10 It is assumed that a WSS of
2.0 N/m2 is suitable for maintaining the structure of arterial vessels and a WSS lower than 1.5 N/m2 will degenerate endothelial cells via the apoptotic cell cycle.10 The WSS of the aneurysm region was barely 1.64 N/m2 even in the peak systole and seems to be too low to maintain the regular cellar functions of endothelial cells. This excessively low WSS may be one of the main factors underlying the degeneration, indicating the structural fragility of the aneurysmal wall. Although a high WSS plays an essential role in the initiation of cerebral aneurysms,11 a low WSS might be a major factor for its growth.
Our 3 ruptured aneurysms had higher averaged WSS of aneurysm region than unruptured aneurysms and they had markedly low WSS in their tip or bleb with high WSS in the body or fundus of aneurysm. We speculate that this low WSS at the tip or the bleb might be responsible for the fragile change of the aneurysm and led to the rupture. Endothelial cells react differently to the high and low WSS.10 The proximity of high and low WSS in a small aneurysm region might enhance the degenerative change of the aneurysm wall.
Until now, there has been no study to our knowledge that demonstrated the changes of the size and the shape of the aneurysm immediately before and after the rupture. These changes might have affected our results. The current and previous4,5 results of ruptured aneurysms might not characterize the aneurysm with high risk of rupture, but may only document the feature of the aneurysm after rupture. The application of these results to the clinical materials will prove the validity.
Conclusions
The results of this study suggest that the CFD technique has the potential to be a useful clinical tool for the prediction of the initiation, growth, and rupture of cerebral aneurysms.
Received May 4, 2004; revision received July 2, 2004; accepted August 10, 2004.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. S. E. Bor, B. K. Velthuis, C. B. Majoie, and G. J.E. Rinkel Configuration of intracranial arteries and development of aneurysms: A follow-up study Neurology, February 26, 2008; 70(9): 700 - 705. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Boussel, M. Wintermark, A. Martin, B. Dispensa, R. VanTijen, J. Leach, V. Rayz, G. Acevedo-Bolton, M. Lawton, R. Higashida, et al. Monitoring Serial Change in the Lumen and Outer Wall of Vertebrobasilar Aneurysms AJNR Am. J. Neuroradiol., February 1, 2008; 29(2): 259 - 264. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Moftakhar, B. Aagaard-Kienitz, K. Johnson, P.A. Turski, A.S. Turk, D.B. Niemann, D. Consigny, J. Grinde, O. Wieben, and C.A. Mistretta Noninvasive Measurement of Intra-Aneurysmal Pressure and Flow Pattern Using Phase Contrast with Vastly Undersampled Isotropic Projection Imaging AJNR Am. J. Neuroradiol., October 1, 2007; 28(9): 1710 - 1714. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S. Alnaes, J. Isaksen, K.-A. Mardal, B. Romner, M. K. Morgan, and T. Ingebrigtsen Computation of Hemodynamics in the Circle of Willis Stroke, September 1, 2007; 38(9): 2500 - 2505. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Tateshima, K. Tanishita, H. Omura, J.P. Villablanca, and F. Vinuela Intra-Aneurysmal Hemodynamics during the Growth of an Unruptured Aneurysm: In Vitro Study Using Longitudinal CT Angiogram Database AJNR Am. J. Neuroradiol., April 1, 2007; 28(4): 622 - 627. [Abstract] [Full Text] [PDF] |
||||
![]() |
M.A. Castro, C.M. Putman, and J.R. Cebral Patient-Specific Computational Fluid Dynamics Modeling of Anterior Communicating Artery Aneurysms: A Study of the Sensitivity of Intra-Aneurysmal Flow Patterns to Flow Conditions in the Carotid Arteries AJNR Am. J. Neuroradiol., November 1, 2006; 27(10): 2061 - 2068. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Mantha, C. Karmonik, G. Benndorf, C. Strother, and R. Metcalfe Hemodynamics in a Cerebral Artery before and after the Formation of an Aneurysm. AJNR Am. J. Neuroradiol., May 1, 2006; 27(5): 1113 - 1118. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Isoda, M. Hirano, H. Takeda, T. Kosugi, M.T. Alley, M. Markl, N.J. Pelc, and H. Sakahara Visualization of Hemodynamics in a Silicon Aneurysm Model Using Time-Resolved, 3D, Phase-Contrast MRI. AJNR Am. J. Neuroradiol., May 1, 2006; 27(5): 1119 - 1122. [Abstract] [Full Text] [PDF] |
||||
![]() |
L.-D. Jou, G. Wong, B. Dispensa, M. T. Lawton, R. T. Higashida, W. L. Young, and D. Saloner Correlation between Lumenal Geometry Changes and Hemodynamics in Fusiform Intracranial Aneurysms AJNR Am. J. Neuroradiol., October 1, 2005; 26(9): 2357 - 2363. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Satoh, M. Omi, C. Ohsako, A. Katsumata, Y. Yoshimoto, S. Tsuchimoto, K. Onoda, K. Tokunaga, K. Sugiu, and I. Date Influence of Perianeurysmal Environment on the Deformation and Bleb Formation of the Unruptured Cerebral Aneurysm: Assessment with Fusion Imaging of 3D MR Cisternography and 3D MR Angiography AJNR Am. J. Neuroradiol., September 1, 2005; 26(8): 2010 - 2018. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Shojima, M. Oshima, K. Takagi, R. Torii, K. Nagata, I. Shirouzu, A. Morita, and T. Kirino Role of the Bloodstream Impacting Force and the Local Pressure Elevation in the Rupture of Cerebral Aneurysms Stroke, September 1, 2005; 36(9): 1933 - 1938. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Stroke Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2004 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |