Skip to main content
Advertisement

Main menu

  • Home
  • Content
    • Current Issue
    • Accepted Manuscripts
    • Article Preview
    • Past Issue Archive
    • Video Articles
    • AJNR Case Collection
    • Case of the Week Archive
    • Case of the Month Archive
    • Classic Case Archive
  • Special Collections
    • AJNR Awards
    • ASNR Foundation Special Collection
    • Most Impactful AJNR Articles
    • Photon-Counting CT
    • Spinal CSF Leak Articles (Jan 2020-June 2024)
  • Multimedia
    • AJNR Podcasts
    • AJNR SCANtastic
    • Trainee Corner
    • MRI Safety Corner
    • Imaging Protocols
  • For Authors
    • Submit a Manuscript
    • Submit a Video Article
    • Submit an eLetter to the Editor/Response
    • Manuscript Submission Guidelines
    • Statistical Tips
    • Fast Publishing of Accepted Manuscripts
    • Graphical Abstract Preparation
    • Imaging Protocol Submission
    • Author Policies
  • About Us
    • About AJNR
    • Editorial Board
    • Editorial Board Alumni
  • More
    • Become a Reviewer/Academy of Reviewers
    • Subscribers
    • Permissions
    • Alerts
    • Feedback
    • Advertisers
    • ASNR Home

User menu

  • Alerts
  • Log in

Search

  • Advanced search
American Journal of Neuroradiology
American Journal of Neuroradiology

American Journal of Neuroradiology

ASHNR American Society of Functional Neuroradiology ASHNR American Society of Pediatric Neuroradiology ASSR
  • Alerts
  • Log in

Advanced Search

  • Home
  • Content
    • Current Issue
    • Accepted Manuscripts
    • Article Preview
    • Past Issue Archive
    • Video Articles
    • AJNR Case Collection
    • Case of the Week Archive
    • Case of the Month Archive
    • Classic Case Archive
  • Special Collections
    • AJNR Awards
    • ASNR Foundation Special Collection
    • Most Impactful AJNR Articles
    • Photon-Counting CT
    • Spinal CSF Leak Articles (Jan 2020-June 2024)
  • Multimedia
    • AJNR Podcasts
    • AJNR SCANtastic
    • Trainee Corner
    • MRI Safety Corner
    • Imaging Protocols
  • For Authors
    • Submit a Manuscript
    • Submit a Video Article
    • Submit an eLetter to the Editor/Response
    • Manuscript Submission Guidelines
    • Statistical Tips
    • Fast Publishing of Accepted Manuscripts
    • Graphical Abstract Preparation
    • Imaging Protocol Submission
    • Author Policies
  • About Us
    • About AJNR
    • Editorial Board
    • Editorial Board Alumni
  • More
    • Become a Reviewer/Academy of Reviewers
    • Subscribers
    • Permissions
    • Alerts
    • Feedback
    • Advertisers
    • ASNR Home
  • Follow AJNR on Twitter
  • Visit AJNR on Facebook
  • Follow AJNR on Instagram
  • Join AJNR on LinkedIn
  • RSS Feeds

AJNR is seeking candidates for the AJNR Podcast Editor. Read the position description.

Research ArticleInterventional

How Do Coil Configuration and Packing Density Influence Intra-Aneurysmal Hemodynamics?

H.G. Morales, M. Kim, E.E. Vivas, M.-C. Villa-Uriol, I. Larrabide, T. Sola, L. Guimaraens and A.F. Frangi
American Journal of Neuroradiology November 2011, 32 (10) 1935-1941; DOI: https://doi.org/10.3174/ajnr.A2635
H.G. Morales
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
M. Kim
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
E.E. Vivas
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
M.-C. Villa-Uriol
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
I. Larrabide
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
T. Sola
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
L. Guimaraens
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
A.F. Frangi
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • Responses
  • References
  • PDF
Loading

Abstract

BACKGROUND AND PURPOSE: Endovascular coiling is a well-established therapy for treating intracranial aneurysms. Nonetheless, postoperative hemodynamic changes induced by this therapy remain not fully understood. The purpose of this work is to assess the influence of coil configuration and packing density on intra-aneurysmal hemodynamics.

MATERIALS AND METHODS: Three 3D rotational angiography images of 3 intracranial aneurysms before and after endovascular coiling were used. For each aneurysm, a 3D representation of the vasculature was obtained after the segmentation of the images. Afterward, a virtual coiling technique was used to treat the aneurysm geometries with coil models. The aneurysms were coiled with 5 packing densities, and each was generated by using 3 coil configurations. Computational fluid dynamics analyses were carried out in both untreated and treated aneurysm geometries. Statistical tests were performed to evaluate the relative effect of coil configuration on local hemodynamics.

RESULTS: The intra-aneurysmal blood flow velocity and wall shear stress were diminished as packing density increased. Aneurysmal flow velocity was reduced >50% due to the first inserted coils (packing density <12%) but with a high dependency on coil configuration. Nonsignificant differences (P > .01) were found in the hemodynamics due to coil configuration for high packing densities (near 30%). A damping effect was observed on the intra-aneurysmal blood flow waveform after coiling.

CONCLUSIONS: Intra-aneurysmal hemodynamics are altered by coils. Coil configuration might reduce its influence on intra-aneurysmal hemodynamics as the packing density increases until an insignificant influence could be achieved for high packing densities.

Abbreviations

3DRA
3D rotational angiography
ANOVA
analysis of variance
CFD
computational fluid dynamics
DSA
digital subtraction angiography
WSS
wall shear stress

Intracranial aneurysms are balloon-like dilations of arteries often occurring at branching points in the circle of Willis. In most cases, these stay asymptomatic and are only diagnosed through medical imaging or after spontaneous rupture. Aneurysm rupture typically leads to subarachnoid hemorrhage producing rates of mortality of approximately 50%.1 The existing therapeutic options seek isolating the aneurysm from the bloodstream and preventing its rupture. The most important therapies are surgical clipping and endovascular interventions,2 such as coiling and stent placement.

Coiling has been the most popular endovascular option for the last 15 years.1–3 This minimally invasive treatment consists of the placement of biocompatible metal alloy coils to block the blood flow into the aneurysm and to prevent rupture. Although coiling is associated to lower mortality and morbidity rates compared with surgical clipping,2 its outcome is not always predictable.1–5 Packing density is thought to have a strong influence in postcoiling outcome, and usually the aneurysms are packed as dense as possible.6–8 Packing density is defined as the ratio between the inserted coils and aneurysm volume.

Moreover, the success of coiling is intimately related to the biologic responses to hemodynamics inside the aneurysm.9 Certainly, CFD is an effective and well-accepted technique to investigate intravascular hemodynamics in silico.10 However, the use of this technique in the presence of endovascular coils is still challenging, due to the geometric complexity of the devices, their small size (diameters of 0.010–0.015 inch),11 and their unpredictable distribution inside the aneurysm.

To simulate the hemodynamics in coiled aneurysm geometries, researchers have modeled explicitly the coils, including straight12 or helical cylinders,13 big hollow spheres,14 or small overlapping spheres along the axis of curved filaments.15 Other approaches assume that coils act as a porous medium inside the aneurysm.16–19 Most of these studies have taken into account the influence of packing density on the hemodynamics. However, little is known about how coil configuration, ie, the arrangement of coils inside the aneurysm, could affect the blood flow. Recently, an in vitro experiment has studied qualitatively the influence of both coil configuration and packing density, but only on near-wall flow.20 In addition, a CFD-based study on an idealized sidewall aneurysm geometry found that the orientation of helical coils has influence on intra-aneurysmal hemodynamics.13

The purpose of this work is to qualitatively and quantitatively assess the influence of coil configuration and packing density on intra-aneurysmal hemodynamics. Three patient-specific aneurysm geometries were virtually treated with coils, and CFD simulations were performed in both untreated and treated geometries. The main hypothesis behind our CFD experiments was that intra-aneurysmal hemodynamics strongly depend on the mechanical resistance produced by the coils as was reported in the literature.21

Materials and Methods

Vascular Models

Three patient-specific aneurysms at different locations in the human brain vasculature were chosen. The locations for cases 1, 2, and 3 were respectively left internal carotid sinus, right middle cerebral artery bifurcation, and basilar tip. All 3 cases were selected for treatment with endovascular coils. 3DRA images of the aneurysms were obtained by using an Integris System (Philips Healthcare, Best, the Netherlands). The size of the images in voxels was 256 × 256 × 256, with an FOV of 134.95 mm. A 3D representation of the luminal surfaces of the vessels and aneurysms was obtained after segmentation of the medical images by using an iso-intensity surface extraction approach. Afterward, the extracted 3D geometries were cleaned, cropped, and smoothed.22

Virtual Coiling

The 3 cases were treated by using a virtual coiling technique.23,24 This technique was designed to place coil models inside aneurysm geometries. The coils are characterized by their geometry, namely, wire diameter and length.

Each coil is progressively placed by advancing its tip, which starts from an initial position inside the aneurysm. The coil tip motion is constrained by the aneurysm dome and its ostium. If the tip of the coil is close to either the surface of the aneurysm or another coil previously inserted, the tip smoothly changes its direction to remain inside the aneurysm and to avoid collision with the other coils. When the coil tip cannot move further because it is blocked by other coils previously placed, the virtual coiling technique partially retracts and redirects the tip toward another region. For a more detailed explanation of the method, see the Appendix.

The technique sequentially places the coils inside the aneurysms until the desired packing density is achieved. As a result, a plausible representation of the distribution of the coils is produced as it could be obtained after a real intervention (Fig 1).

Fig 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig 1.

Results of virtual coiling for case 2. Untreated model (A) and treated models corresponding to coil packing densities of 5.6% (B), 11.2% (C), 16.7% (D), 22.2% (E), and 30.0% (F).

This virtual coiling technique was used in our 3 cases. Five coil packing densities were produced (Table 1). Packing densities were calculated from the volumes of the 3D geometries. Also, for each packing density, 3 coil configurations were generated and labeled as I, II, and III. In this way, the influence of geometric dissimilarities due to the coil configuration on intra-aneurysmal hemodynamics could be evaluated.

View this table:
  • View inline
  • View popup
Table 1:

Design parameters of virtual coiling experiments

Hemodynamic Modeling

To solve the governing equations of fluid flow, body-fitted unstructured tetrahedral meshes were generated with the vascular geometries and coils by using the commercial software ANSYS ICEM CFD version 11 (Ansys, Canonsburg, Pennsylvania). In the treated models, the volumetric mesh near the coils was refined, and a convergence test was performed to obtain mesh-independent CFD simulations. The test criterion was set to be below 5% of differences of the averaged WSS at the aneurysm between consecutive mesh refinements. Average WSS was found to fulfill this criterion when the coil cross-section perimeter was enclosed with approximately 10 mesh elements.

These volumetric meshes were exported in ANSYS CFX version 11 (Ansys) for the CFD simulations. Blood was modeled as a Newtonian fluid with viscosity of 0.00375 Pa × s. Flow was assumed to be laminar and incompressible (blood density of 1060 kg/m3). Vessel wall and coil surfaces were imposed to be rigid with no-slip boundary condition on the flow. At the inlet of the models, a parabolic velocity profile (spatial condition) was imposed as well as a physiologic blood flow waveform (temporal condition) from a normal individual. These flow waveforms were scaled to ensure a time-averaged Reynolds number of 300 at the inlet of each model. At the outlets, a traction-free boundary condition was applied.

Two entire cardiac cycles were simulated, and the results from the second cycle were analyzed as time-averaged, and at peak systole and end diastole. The blood flow velocity, the average WSS, and the areas of low (<0.4 Pa) and high WSS (>1.5 Pa) were evaluated in the aneurysm. These WSS thresholds were selected because of their association with pathologic states of the arterial wall.25

Statistical Evaluation

ANOVA tests were performed to evaluate the relative effect of the coil configuration on the intra-aneurysmal hemodynamics at different packing densities. The results for an entire cardiac cycle were taken into account as samples for the tests. Before applying the tests, the hypotheses of normal distribution, homogeneity of the variance, and independence were confirmed for each sample.

Results

Qualitative Findings

Figure 2 depicts velocity magnitudes at peak systole on a fixed aneurysm cross-section for case 2. Packing densities are shown in Fig 2A, and Fig 2B compares different coil configurations. The untreated aneurysms showed high-velocity regions with strong flow impingement on the distal wall and a vortical flow structure inside the aneurysm. For packing densities <11.2%, inflow was disturbed and aneurysmal flow velocity was decreased, especially at the aneurysm fundus. As the number of coils increased, a progressive reduction in the flow velocity inside the aneurysm was observed. Nonetheless, small additional velocity reductions were associated with coil packing densities >22.2%. Furthermore, Fig 2B shows how 2 coil configurations for case 2 were able to visually induce substantial differences on the flow fields at a low packing density, such as 11.2%.

Fig 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig 2.

Blood flow velocity magnitudes in gray-scale for an aneurysm cross-section at peak systole for case 2. A, Changes in the velocity are presented with a fixed coil configuration and increasing packing density. B, Two coil configurations are compared for packing densities of 11.2% and 16.7%. The coils are presented in yellow and the percentages represent the packing densities achieved for each model.

Visualization of WSS distributions at peak systole and end-diastole are presented in Fig 3A and B, respectively. Here, simulations for case 2 are shown with and without coils. A reduction of WSS was already observed for a low packing density of 5.6%, whereas above 16.7%, a stable low WSS throughout the cardiac cycle was found at the aneurysm dome. A remarkably high WSS at peak systole near the aneurysms' neck were also observed. Similar WSS distributions as well as flow structures and velocity reduction trends were visualized for cases 1 and 3.

Fig 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig 3.

WSS distributions for the untreated and several treated models for case 2, calculated at peak systole (A) and end diastole (B).

Quantitative Findings

Time-averaged velocity reduction inside the aneurysm based on the untreated model was computed for the different coil configurations and packing densities (Fig 4A). The velocity reduction was almost independent of the coil configuration for packing densities >25%. In general, velocity reduction increased as packing density rose. For high packing densities (near 30%), the coils induced a maximum velocity reduction of approximately 70% for case 1, 65% for case 2, and almost 90% for case 3.

Fig 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig 4.

A, Time-averaged intra-aneurysmal mean flow velocity reduction rate. B, Time-averaged WSS in the aneurysm dome.

Figure 4B shows time-averaged aneurysmal WSS. For all cases, a reduction in WSS was found after coiling, even for those with lower packing densities. A progressive decrease of the average WSS was observed as packing density increased. For low packing densities (<10%), a large variability in the average WSS was induced by the coil configuration, especially for case 2, with 5.6% of packing density (mean WSS of 1.65 Pa and SD of 0.4 Pa).

To analyze the evolution of the WSS in a cardiac cycle, the areas of low (<0.4 Pa) and high (>1.5 Pa) WSS were computed for each aneurysm. In Fig 5, these areas were presented for 3 CFD models of case 2, including the untreated model (Fig 5A). A clear difference in these areas is visualized between the simulations with and without coils. After coiling, low WSS areas increased and reached values of approximately 75% for case 1, 60% for case 2%, and 95% for case 3. The high WSS areas decreased and were approximately 20% for cases 1 and 2, and <10% for case 3. Also, there were noticeable differences in the WSS areas due to the coil configuration for a packing density <22.2% (Fig 5B), especially near peak systole. However, these differences were less evident as packing density increased (Fig 5C).

Fig 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig 5.

Variation of the aneurysm areas with low and high WSS expressed as a percentage of the whole aneurysm dome during one cardiac cycle for case 2. A, Untreated model. B and C, treated models with packing densities of 22.2% and 30.0%, respectively.

Analysis of Statistical Data

The ANOVA test results (Table 2) showed that there were nonsignificant differences (P > .01) on the average WSS and high WSS areas due to the coil configuration when packing density was >20%. However, there were significant differences (P < .01) in the velocity reduction and low WSS areas when the packing density was <22.2% in all 3 cases (case 2 had the highest packing density). The only exception to this result was for the low WSS areas in case 1, in which the differences induced by the coil configuration were negligible. Table 2 also shows that, for the studied cases, the hemodynamic variables derived from different coil configurations were nonsignificant for packing densities of approximately 30%.

View this table:
  • View inline
  • View popup
Table 2:

P values of ANOVA tests performed for each of the cases at different packing densities

Discussion

The present study aims at elucidating how packing density and coil configuration might influence the intra-aneurysmal hemodynamics. Four main findings were obtained from our numeric modeling and data analyses. First, the intra-aneurysmal velocity and WSS decreased as packing density increased. Second, the influence of coil configuration on intra-aneurysmal hemodynamics proved to be negligible for elevated packing densities (near 30%). However for lower packing densities, this influence was larger and statistically significant. Third, aneurysmal flow velocity was reduced by >50% after the insertion of the first coils (packing density <12%). Fourth, evidence was found that the blood flow waveform is damped near the aneurysm wall as packing density increases, which was reflected in the flattening of the areas with low and high WSS.

These findings were obtained by qualitative observations, analyses of quantitative data, and statistical evaluations of our CFD simulations. The first hemodynamic variable analyzed was blood flow velocity inside the aneurysm. In Fig 2, a progressive reduction of the velocity inside the aneurysm was observed as packing density increased. This finding is consistent with previous CFD-based studies and in vitro experiments.13,15,19,20,26 Particular flow patterns were visualized due to the coil configuration (Fig 2B), which are associated with the special location of the coils at the aneurysm ostium. Specific flow characteristics also were reported by an in vivo study, in which short and disturbed streamlines were found near the aneurysm wall.20

The second evaluated hemodynamic variable was average WSS. In Fig 3, a stable (systole and diastole) low WSS distribution (blue areas) was found as packing densities were >16.7%. This was mainly observed in the aneurysm fundus. Furthermore, high WSS (red areas) was localized at the aneurysm neck, which is associated with the blockage and reflection of the flow stream due to the coils. These WSS distributions were previously observed in a CFD-based study modeling the coils as a porous medium.18 Also, the reduction of average WSS as packing density increased was noticeable in Fig 3 by the increment of the blue area.

The third hemodynamic variable was the intra-aneurysmal velocity reduction. In Fig 4A, the previous qualitative observations were confirmed as for the influence of packing density and coil configuration on the local hemodynamics. The efficiency of the coils reducing the intra-aneurysmal flow velocity depends on coil configuration and packing density as is visualized among cases. Furthermore, the dependency of the velocity reduction on coil configuration was lower for high packing densities (>25%).

For our 3 cases, the insertion of the first coils (packing densities <12%) produced the most important velocity reduction (>50%) inside the aneurysm. However, more variability of the velocity reduction was induced by low packing densities (see cases 2 and 3 in Fig 4A). In general, an asymptotic trend of the velocity reduction was observed as packing density increased. In other words, the last coils have a smaller role in the velocity reduction inside the aneurysm. Additionally, the maximum velocity reduction observed in Fig 4A for our cases (70%–90%), corresponds well to previous CFD-based studies18 and in vivo animal experiments.27

The last hemodynamic variables studied were the aneurysmal areas of low (<0.4 Pa) and high (>1.5 Pa) WSS, presented in Fig 5. Here, the high/low WSS areas were damped by the coils, suggesting that the hemodynamics near the wall tend to be steady as packing density increases. However, the damping was strongly related to packing density. Furthermore, most of the variations in WSS areas happened near the aneurysm neck (Fig 3), thereby WSS at the aneurysm dome was almost time-invariant. This steady hemodynamic condition at the aneurysm wall combined with the previously mentioned low-shear stresses and low-velocity flow might be favorable hemodynamic conditions to isolate the aneurysm from the bloodstream.9,28 In addition, it was observed how coil configuration reduces its influence on hemodynamics when packing density increases (Fig 5).

The statistical tests were performed to study the variability induced by the coil configuration as packing density increased on the studied hemodynamic variables. We found that coil configuration has an influence on intra-aneurysmal hemodynamics when packing density is low. However, the hemodynamic dependence on coil configuration was reduced as packing density increased. This dependence was statistically nonsignificant (P > .01) for packing densities of approximately 30%.

From our results, we observe that intra-aneurysmal hemodynamics tend to stabilize as packing density increases and becomes independent of coil configuration. However, the reported packing densities do not ensure the success of coiling, mainly because the required hemodynamic conditions to warranty aneurysm occlusion by coiling are unknown. Also, the presented results are based on the flow waveforms imposed as boundary conditions, which were physiologic and measured under rest condition but neither patient-specific nor generic. As it has been found that local hemodynamic forces are associated with heart rate and aneurysm geometry,29 further investigation is required to generalize or extrapolate our findings to other physiologic flow waveforms (different heart rate, flow rate, or waveform shape) and different vascular morphologies.

Our results have proved to be comparable to previous CFD-based studies with coils, but our virtual coiling technique presents some advantages over these. The first advantage is that our technique has the potential of becoming relevant in clinical practice. It is able to reproduce the wide range of packing densities as it is obtained in real interventions6–8 and in vitro experiments30 with coils. Furthermore, it considers the geometric specifications of coils, such as wire diameter and length, and can be used in patient-specific geometries. The technique uses medical images before treatment to extract the 3D representation of the vasculature, and the manual definition of the ostium to constrain the space occupied by the coil models provides a tool for preoperatory hemodynamic assessment. In addition, the use of images after treatment allows us to extract similar regions to insert the coils as the regions visualized with medical images, which provides postoperative hemodynamic assessment.

The other advantage is the explicit modeling of the coils, which allowed us to investigate the influence of coil configuration and packing density on the intra-aneurysmal hemodynamics. Certainly, packing density can be modeled with a porous medium approach by setting the porosity,16,18,19 but the simulation of different coil configurations is not trivial due to the assumptions of homogeneity and isotropy of the medium. In contrast, a study by Cebral and Löhner15 proposed an explicit modeling of the coils for further CFD analysis. However, it is not clear the packing densities that can be achieved with the generation of arbitrary objects and more specifically, the coil distribution near the aneurysm walls. Finally, the helical coils presented by Schirmer and Malek13 in an idealized aneurysm geometry are fully user-dependent and did not achieve high packing densities.

Despite these advantages and our interesting findings, this work has some limitations. One limitation is that the coils were assumed as rigid bodies and their motion was not allowed in our CFD simulations, which means that the coils were not affected by the flow. Although this assumption considerably simplifies the numeric solution, coil deformation as well as a more precise interaction with the flow should be taken into account to investigate other relevant issues related to this therapy, such as coil compaction. Also, another limitation is that the vessel and aneurysm walls were considered as rigid bodies due to the lack of information regarding wall thickness and mechanical properties of the tissue. Nevertheless, such choice is supported by several CFD-based studies that assumed the rigid wall condition and claimed that this effect is minor compared with others, such as lumen morphology.15,16,18,19,31

A CFD simulation able to assess the success or failure of coiling an aneurysm should consider the thrombus formation inside the aneurysm. In this work, thrombus formation was not tackled because efforts were mainly focused on evaluating the hemodynamic resistance of the coils as these have been found to be relevant experimentally.21 In addition, thrombus formation has been related to hemodynamics and blood viscosity.9 Our study was limited to the use of a Newtonian model for blood viscosity, and further investigation is required to establish the influence of the non-Newtonian behavior of blood on the hemodynamics of coiled aneurysms. The assumption of Newtonian flow is justified because it has been reported to be as good as non-Newtonian models in similar studies in arteries with aneurysms compared with the uncertainties introduced by geometric models.31 This aspect also has been recently used to assess the effect of flow diverters on intra-aneurysmal hemodynamics.32 However, the consideration of a non-Newtonian model opens questions such as which non-Newtonian model (eg, Casson, Herschel-Bulkley, or Power law) represents more accurately the behavior of blood and why, what are the differences and similarities compared with a Newtonian model, and whether these differences are meaningful.

Conclusions

We have studied the influence of coiling treatment on intra-aneurysmal hemodynamics in 3 patient-specific aneurysm geometries by using a virtual coiling technique and CFD modeling. Four main findings were extracted from our results. First, coils attenuate blood flow velocity and WSS with increasing packing density. Second, the influence of coil configuration on intra-aneurysmal hemodynamics is reduced with increasing packing density, such that no significant influence (P > .01) is found for packing densities of 30%. Third, the first inserted coils (packing density <12%) cause >50% of the velocity reduction inside the aneurysm. Finally, a damping effect is observed on WSS areas due to the coils, which creates a stable hemodynamic condition inside the aneurysm.

Appendix

Virtual Coiling Technique

The virtual coiling technique uses a dynamic path planning approach to insert one-by-one the coils inside a closed geometry (in this case, aneurysm geometry was defined by its wall and ostium). Inserting a coil means consecutively selecting the next position of the coil tip. This position is selected among a set of candidate positions around the tip by evaluating a potential field φ for each of them. The candidate position with the lowest φ is chosen as the next tip position. In this way, the algorithm inserts a coil (creates a path) by advancing the coil tip.

The potential field φ has 4 rules (or subpotential fields) that mimic the observed behavior of coils during their placement. Embedded Image The first rule φA constrains the region where the coil tip is allowed to move, and it is limited by the aneurysm wall and the ostium. It ensures that the coil tip remains inside the aneurysm, and it is defined as 0 inside the aneurysm and 1 elsewhere.

The second rule φC avoids the collision between coils and keeps the coil tip away from highly dense coiled areas. The rule φC takes values between 0 and 1, being 1 when the evaluated position is occupied by a previously inserted coil.

The third rule φD considers the deformation of the coils. Basically, the larger the angle of the evaluated candidate is, compared with the previous direction, the larger φD will be. This rule presents values between 0 and 0.5, being 0 for an angle equal to 0 and 0.5 when the angle of the candidate position is π.

Finally, the rule φPB is meant to solve the situation where the tip is trapped by the aneurysm wall, other coil(s), or both, which means that all candidate locations have φ ≥ 1 (dead-end). In this situation, the algorithm pulls back and advances the coil tip in a different direction.

To initialize a coil, a set of candidate initial positions is created and the minimizing φC + φD is selected. For the first inserted coil, only φD is considered, which is compared with a direction provided by the user. This guarantees that by changing this initial parameter, the coil configuration is changed.

Virtual Coiling Technique Assessment

To assess the effect of the virtual coiling technique on hemodynamics in comparison with the real coils,24 a CFD model was generated from medical images available for a coiled aneurysm. The aneurysm was located at the right middle cerebral artery, and treated with a 43% packing density.

To construct the CFD model, a 3DRA image of the patient was segmented, cropped, cleaned, and smoothed to obtain the vascular model of the aneurysm and parent artery. Then, the virtual coiling technique was used to insert in the vascular model the same number of coils with the same geometric characteristics as in the real treatment. To set the boundary conditions of the CFD model, pseudopatient-specific conditions were imposed. A time-intensity curve was extracted from the DSA image sequence at the location of the model inlet. This curve was calculated as the average pixel intensity in a region near the inlet for each time instant. This information was used as boundary condition to solve an advection-diffusion equation to calculate the transport of a contrast material, modeled as a mass-less scalar.

After the CFD simulation, a set of virtual angiographic images were created. These virtual images were qualitatively compared with the DSA image sequence of the patient and a good agreement between these was found. Time-intensity curves also were measured at several control regions along the vasculature on the DSA and virtual images (CFD-based) for a quantitative comparison, including locations downstream of the coiled aneurysm. The comparison of both time-intensity curves showed an excellent correlation in all the evaluated control regions.

Footnotes

  • Disclosures: Hernán Morales, Research Support (including provision of equipment or materials): Universitat Pompeu Fabra/CISTIB, Details: This work was partially supported within the CENIT program, as part of CDTEAM and cvREMOD projects funded by the Spanish CDTI and partly within the framework of the @neurIST Project (IST-2005-027703), which is cofinanced by the European Commission within the IST Program of the Sixth Framework Program; Minsuok Kim, Research Support (including provision of equipment or materials): Spanish CDTI, European Commission, Details: This work was partially supported within the CENIT program, as part of CDTEAM and cvREMOD projects funded by the Spanish CDTI and partly within the framework of the @neurIST Project (IST-2005-027703), which is cofinanced by the European Commission within the IST Program of the Sixth Framework Program; Ignacio Larrabide, Research Support (including provision of equipment or materials): CIBER-BBN, Details: Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina. This researcher is funded by this institution; Alejandro Frangi, Research Support (including provision of equipment or materials): ANSYS, Details: The company provides us with free software licenses for numeric simulation of flow dynamics. This researcher does not anticipate any conflict of interest.

  • All authors have substantially contributed to the preparation of this manuscript. Preliminary results for the virtual coiling technique and further CFD modeling previously presented at World Congress on Medical Physics and Biomedical Engineering, September 7–12, 2009, Munich, Germany.

References

  1. 1.↵
    1. van Gijn J,
    2. Rinkel GJ
    . Subarachnoid haemorrhage: diagnosis, causes and management. Brain 2001;124:249–78
    Abstract/FREE Full Text
  2. 2.↵
    1. Molyneux AJ,
    2. Kerr RS,
    3. Yu LM,
    4. et al
    . International subarachnoid aneurysm trial (ISAT) of neurosurgical clipping versus endovascular coiling in 2143 patients with ruptured intracranial aneurysms: a randomised comparison of effects on survival, dependency, seizures, rebleeding, subgroups, and aneurysm occlusion. Lancet 2005;366:809–17
    CrossRefPubMed
  3. 3.↵
    1. Henkes H,
    2. Fischer S,
    3. Weber W,
    4. et al
    . Endovascular coil occlusion of 1811 intracranial aneurysms: early angiographic and clinical results. Neurosurgery 2004;54:268–80; discussion 80–85
    CrossRefPubMed
  4. 4.↵
    1. Park SH,
    2. Lee CY,
    3. Yim MB
    . The merits of endovascular coil surgery for patients with unruptured intracranial aneurysms. J Korean Neurosurg Soc 2008;43:270–74
    CrossRefPubMed
  5. 5.↵
    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
  6. 6.↵
    1. Kawanabe Y,
    2. Sadato A,
    3. Taki W,
    4. et al
    . Endovascular occlusion of intracranial aneurysms with Guglielmi detachable coils: correlation between coil packing density and coil compaction. Acta Neurochir (Wien) 2001;143:451–55
    CrossRefPubMed
  7. 7.↵
    1. Sluzewski M,
    2. van Rooij WJ,
    3. Slob MJ,
    4. et al
    . Relation between aneurysm volume, packing, and compaction in 145 cerebral aneurysms treated with coils. Radiology 2004;231:653–58
    CrossRefPubMed
  8. 8.↵
    1. Tamatani S,
    2. Ito Y,
    3. Abe H,
    4. et al
    . Evaluation of the stability of aneurysms after embolization using detachable coils: correlation between stability of aneurysms and embolized volume of aneurysms. AJNR Am J Neuroradiol 2002;23:762–67
    Abstract/FREE Full Text
  9. 9.↵
    1. Wootton DM,
    2. Ku DN
    . Fluid mechanics of vascular systems, diseases, and thrombosis. Annu Rev Biomed Eng 2003;1:299–329
  10. 10.↵
    1. Alnaes MS,
    2. Isaksen J,
    3. Mardal KA,
    4. et al
    . Computation of hemodynamics in the circle of Willis. Stroke 2007;38:2500–05
    Abstract/FREE Full Text
  11. 11.↵
    1. White JB,
    2. Ken CG,
    3. Cloft HJ,
    4. et al
    . Coils in a nutshell: a review of coil physical properties. AJNR Am J Neuroradiol 2008;29:1242–46
    Abstract/FREE Full Text
  12. 12.↵
    1. Narracott A,
    2. Smith S,
    3. Lawford P,
    4. et al
    . Development and validation of models for the investigation of blood clotting in idealized stenoses and cerebral aneurysms. J Artif Organs 2005;8:56–62
    CrossRefPubMed
  13. 13.↵
    1. Schirmer CM,
    2. Malek AM
    . Critical influence of framing coil orientation on intra-aneurysmal and neck region hemodynamics in a sidewall aneurysm model. Neurosurgery 2010;67:1692–702
    CrossRefPubMed
  14. 14.↵
    1. Byun HS,
    2. Rhee K
    . CFD modeling of blood flow following coil embolization of aneurysms. Med Eng Phys 2004;26:755–61
    CrossRefPubMed
  15. 15.↵
    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
    CrossRefPubMed
  16. 16.↵
    1. Cha KS,
    2. Balaras E,
    3. Lieber BB,
    4. et al
    . Modeling the interaction of coils with the local blood flow after coil embolization of intracranial aneurysms. J Biomech Eng 2007;129:873–79
    CrossRefPubMed
  17. 17.↵
    1. Groden C,
    2. Eckert B,
    3. Ries T,
    4. et al
    . Angiographic follow-up of vertebrobasilar artery aneurysms treated with detachable coils. Neuroradiology 2003;45:435–40
    CrossRefPubMed
  18. 18.↵
    1. Jou L-D,
    2. Saloner D,
    3. Higashida R
    . Determining intra-aneurysmal flow for coiled cerebral aneurysm with digital fluoroscopy. Biomed Eng Appl Basis Comm 2004;16:6
  19. 19.↵
    1. Kakalis NM,
    2. Mitsos AP,
    3. Byrne JV,
    4. et al
    . The haemodynamics of endovascular aneurysm treatment: a computational modelling approach for estimating the influence of multiple coil deployment. IEEE Trans Med Imaging 2008;27:814–24
    CrossRefPubMed
  20. 20.↵
    1. Goubergrits L,
    2. Thamsen B,
    3. Berthe A,
    4. et al
    . In vitro study of near-wall flow in a cerebral aneurysm model with and without coils. AJNR Am J Neuroradiol 2010;31:1521–28
    Abstract/FREE Full Text
  21. 21.↵
    1. Sorteberg A,
    2. Sorteberg W,
    3. Aagaard BD,
    4. et al
    . Hemodynamic versus hydrodynamic effects of Guglielmi detachable coils on intra-aneurysmal pressure and flow at varying pulse rate and systemic pressure. AJNR Am J Neuroradiol 2004;25:1049–57
    Abstract/FREE Full Text
  22. 22.↵
    1. Villa-Uriol MC,
    2. Larrabide I,
    3. Geers AJ,
    4. et al
    . AngioLab: integrated technology for patient-specific management of intracranial aneurysms. Conf Proc IEEE Eng Med Biol Soc 2010;1:6801–04
  23. 23.↵
    1. Morales HG,
    2. Kim M,
    3. Villa-Uriol M-C,
    4. et al
    . Influence of coil packing rate and configuration on intracranial aneurysm hemodynamics. World Congress on Medical Physics and Biomedical Engineering, September 7–12, 2009, Munich, Germany 2009;25/4
  24. 24.↵
    1. Morales HG,
    2. Larrabide I,
    3. Villa-Uriol M-C,
    4. et al
    . Towards the validation of a virtual coiling technique using a real versus a simulated bolus injection. Second International Conference on Computational & Mathematical Biomedical Engineering, March 30-April 1, 2011, Washington, DC; 512–15
  25. 25.↵
    1. Malek AM,
    2. Alper SL,
    3. Izumo S
    . Hemodynamic shear stress and its role in atherosclerosis. JAMA 1999;282:2035–42
    CrossRefPubMed
  26. 26.↵
    1. Babiker MH,
    2. Gonzalez LF,
    3. Albuquerque F,
    4. et al
    . Quantitative effects of coil packing density on cerebral aneurysm fluid dynamics: an in vitro steady flow study. Ann Biomed Eng 2010;38:2293–301
    CrossRefPubMed
  27. 27.↵
    1. Sorteberg A,
    2. Sorteberg W,
    3. Rappe A,
    4. et al
    . Effect of Guglielmi detachable coils on intraaneurysmal flow: experimental study in canines. AJNR Am J Neuroradiol 2002;23:288–94
    Abstract/FREE Full Text
  28. 28.↵
    1. Mori K,
    2. Nakao Y,
    3. Horinaka N,
    4. et al
    . Cerebral aneurysm regrowth and coil unraveling after incomplete Guglielmi detachable coil embolization: serial angiographical and histological findings. Neurol Med Chir (Tokyo) 2003;43:293–97
    CrossRefPubMed
  29. 29.↵
    1. Jiang J,
    2. Strother C
    . Computational fluid dynamics simulations of intracranial aneurysms at varying heart rates: a “patient-specific” study. J Biomech Eng 2009;131:091001
  30. 30.↵
    1. Piotin M,
    2. Mandai S,
    3. Murphy KJ,
    4. et al
    . Dense packing of cerebral aneurysms: an in vitro study with detachable platinum coils. AJNR Am J Neuroradiol 2000;21:757–60
    Abstract/FREE Full Text
  31. 31.↵
    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
    CrossRefPubMed
  32. 32.↵
    1. Cebral JR,
    2. Mut F,
    3. Raschi M,
    4. et al
    . Aneurysm rupture following treatment with flow-diverting stents: computational hemodynamics analysis of treatment. AJNR Am J Neuroradiol 2011;32:27–33
    Abstract/FREE Full Text
  • Received December 14, 2010.
  • Accepted after revision March 5, 2011.
  • © 2011 by American Journal of Neuroradiology
View Abstract
PreviousNext
Back to top

In this issue

American Journal of Neuroradiology: 32 (10)
American Journal of Neuroradiology
Vol. 32, Issue 10
1 Nov 2011
  • Table of Contents
  • Index by author
Advertisement
Print
Download PDF
Email Article

Thank you for your interest in spreading the word on American Journal of Neuroradiology.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
How Do Coil Configuration and Packing Density Influence Intra-Aneurysmal Hemodynamics?
(Your Name) has sent you a message from American Journal of Neuroradiology
(Your Name) thought you would like to see the American Journal of Neuroradiology web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Cite this article
H.G. Morales, M. Kim, E.E. Vivas, M.-C. Villa-Uriol, I. Larrabide, T. Sola, L. Guimaraens, A.F. Frangi
How Do Coil Configuration and Packing Density Influence Intra-Aneurysmal Hemodynamics?
American Journal of Neuroradiology Nov 2011, 32 (10) 1935-1941; DOI: 10.3174/ajnr.A2635

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
0 Responses
Respond to this article
Share
Bookmark this article
How Do Coil Configuration and Packing Density Influence Intra-Aneurysmal Hemodynamics?
H.G. Morales, M. Kim, E.E. Vivas, M.-C. Villa-Uriol, I. Larrabide, T. Sola, L. Guimaraens, A.F. Frangi
American Journal of Neuroradiology Nov 2011, 32 (10) 1935-1941; DOI: 10.3174/ajnr.A2635
del.icio.us logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One
Purchase

Jump to section

  • Article
    • Abstract
    • Abbreviations
    • Materials and Methods
    • Results
    • Discussion
    • Conclusions
    • Appendix
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • Responses
  • References
  • PDF

Related Articles

  • No related articles found.
  • PubMed
  • Google Scholar

Cited By...

  • Improving Computational Fluid Dynamics Simulations of Coiled Aneurysms Using Finite Element Modeling
  • How to perform intra-aneurysmal coil embolization after Pipeline deployment: a study from a hemodynamic viewpoint
  • Differences in hemodynamic characteristics under high packing density between the porous media model and finite element analysis in computational fluid dynamics of intracranial aneurysm virtual treatment
  • Optimal first coil selection to avoid aneurysmal recanalization in endovascular intracranial aneurysmal coiling
  • Risk Factor Analysis of Recanalization Timing in Coiled Aneurysms: Early versus Late Recanalization
  • Computational fluid dynamics of cerebral aneurysm coiling using high-resolution and high-energy synchrotron X-ray microtomography: comparison with the homogeneous porous medium approach
  • Effect of hemodynamics on outcome of subtotally occluded paraclinoid aneurysms after stent-assisted coil embolization
  • Volume versus standard coils in the treatment of intracranial aneurysms
  • Fate of Coiled Aneurysms with Minor Recanalization at 6 Months: Rate of Progression to Further Recanalization and Related Risk Factors
  • Aneurysm permeability following coil embolization: packing density and coil distribution
  • HydroCoils Are Associated with Lower Angiographic Recurrence Rates Than Are Bare Platinum Coils in Treatment of "Difficult-to-Treat" Aneurysms: A Post Hoc Subgroup Analysis of the HELPS Trial
  • CFD: Computational Fluid Dynamics or Confounding Factor Dissemination? The Role of Hemodynamics in Intracranial Aneurysm Rupture Risk Assessment
  • Initial experience with Penumbra Coil 400 versus standard coils in embolization of cerebral aneurysms: a retrospective review
  • Analysis and quantification of endovascular coil distribution inside saccular aneurysms using histological images
  • Crossref (75)
  • Google Scholar

This article has been cited by the following articles in journals that are participating in Crossref Cited-by Linking.

  • CFD: Computational Fluid Dynamics or Confounding Factor Dissemination? The Role of Hemodynamics in Intracranial Aneurysm Rupture Risk Assessment
    J. Xiang, V.M. Tutino, K.V. Snyder, H. Meng
    American Journal of Neuroradiology 2014 35 10
  • Newtonian and non-Newtonian blood flow in coiled cerebral aneurysms
    Hernán G. Morales, Ignacio Larrabide, Arjan J. Geers, Martha L. Aguilar, Alejandro F. Frangi
    Journal of Biomechanics 2013 46 13
  • Finite element modeling of endovascular coiling and flow diversion enables hemodynamic prediction of complex treatment strategies for intracranial aneurysm
    Robert J. Damiano, Ding Ma, Jianping Xiang, Adnan H. Siddiqui, Kenneth V. Snyder, Hui Meng
    Journal of Biomechanics 2015 48 12
  • Finite element modeling of embolic coil deployment: Multifactor characterization of treatment effects on cerebral aneurysm hemodynamics
    M. Haithem Babiker, Brian Chong, L. Fernando Gonzalez, Sachmanik Cheema, David H. Frakes
    Journal of Biomechanics 2013 46 16
  • Overlapping stents for blood blister-like aneurysms of the internal carotid artery
    Yi-Bin Fang, Qiang Li, Yi-Na Wu, Qi Zhang, Peng-Fei Yang, Wen-Yuan Zhao, Qing-Hai Huang, Bo Hong, Yi Xu, Jian-Min Liu
    Clinical Neurology and Neurosurgery 2014 123
  • Intrasaccular flow-diversion for treatment of intracranial aneurysms: the Woven EndoBridge
    Benjamin Mine, Laurent Pierot, Boris Lubicz
    Expert Review of Medical Devices 2014 11 3
  • Hemodynamic Effect of Flow Diverter and Coils in Treatment of Large and Giant Intracranial Aneurysms
    Linkai Jing, Jingru Zhong, Jian Liu, Xinjian Yang, Nikhil Paliwal, Hui Meng, Shengzhang Wang, Ying Zhang
    World Neurosurgery 2016 89
  • A Virtual Coiling Technique for Image-Based Aneurysm Models by Dynamic Path Planning
    H. G. Morales, I. Larrabide, A. J. Geers, Luis San Roman, J. Blasco, J. M. Macho, A. F. Frangi
    IEEE Transactions on Medical Imaging 2013 32 1
  • Computational fluid dynamics of blood flow in coil-embolized aneurysms: effect of packing density on flow stagnation in an idealized geometry
    Tomohiro Otani, Masanori Nakamura, Toshiyuki Fujinaka, Masayuki Hirata, Junko Kuroda, Katsuhiko Shibano, Shigeo Wada
    Medical & Biological Engineering & Computing 2013 51 8
  • Numerical simulation of patient-specific endovascular stenting and coiling for intracranial aneurysm surgical planning
    Xiaochang Leng, Yang Wang, Jing Xu, Yeqing Jiang, Xiaolong Zhang, Jianping Xiang
    Journal of Translational Medicine 2018 16 1

More in this TOC Section

  • SAVE vs. Solumbra Techniques for Thrombectomy
  • Contrast-Induced Encephalopathy after NeuroIR
  • CT Perfusion&Reperfusion in Acute Ischemic Stroke
Show more Interventional

Similar Articles

Advertisement

Indexed Content

  • Current Issue
  • Accepted Manuscripts
  • Article Preview
  • Past Issues
  • Editorials
  • Editor's Choice
  • Fellows' Journal Club
  • Letters to the Editor
  • Video Articles

Cases

  • Case Collection
  • Archive - Case of the Week
  • Archive - Case of the Month
  • Archive - Classic Case

Special Collections

  • AJNR Awards
  • ASNR Foundation Special Collection
  • Most Impactful AJNR Articles
  • Photon-Counting CT
  • Spinal CSF Leak Articles (Jan 2020-June 2024)

More from AJNR

  • Trainee Corner
  • Imaging Protocols
  • MRI Safety Corner

Multimedia

  • AJNR Podcasts
  • AJNR Scantastics

Resources

  • Turnaround Time
  • Submit a Manuscript
  • Submit a Video Article
  • Submit an eLetter to the Editor/Response
  • Manuscript Submission Guidelines
  • Statistical Tips
  • Fast Publishing of Accepted Manuscripts
  • Graphical Abstract Preparation
  • Imaging Protocol Submission
  • Evidence-Based Medicine Level Guide
  • Publishing Checklists
  • Author Policies
  • Become a Reviewer/Academy of Reviewers
  • News and Updates

About Us

  • About AJNR
  • Editorial Board
  • Editorial Board Alumni
  • Alerts
  • Permissions
  • Not an AJNR Subscriber? Join Now
  • Advertise with Us
  • Librarian Resources
  • Feedback
  • Terms and Conditions
  • AJNR Editorial Board Alumni

American Society of Neuroradiology

  • Not an ASNR Member? Join Now

© 2025 by the American Society of Neuroradiology All rights, including for text and data mining, AI training, and similar technologies, are reserved.
Print ISSN: 0195-6108 Online ISSN: 1936-959X

Powered by HighWire