Improving accuracy for finite element modeling of endovascular coiling of intracranial aneurysm

Background: Computer modeling of endovascular coiling intervention for intracranial aneurysm could enable a priori patient-specific treatment evaluation. To that end, we previously developed a finite element method (FEM) coiling technique, which incorporated simplified assumptions. To improve accuracy in capturing real-life coiling, we aimed to enhance the modeling strategies and experimentally test whether improvements lead to more accurate coiling simulations. Methods: We previously modeled coils using a pre-shape based on mathematical curves and mechanical properties based on those of platinum wires. In the improved version, to better represent the physical properties of coils, we model coil pre-shapes based on how they are manufactured, and their mechanical properties based on their spring-like geometric structures. To enhance the deployment mechanics, we include coil advancement to the aneurysm in FEM simulations. To test if these new strategies produce more accurate coil deployments, we fabricated silicone phantoms of 2 patient-specific aneurysms in duplicate, deployed coils in each, and quantified coil distributions from intra-aneurysmal cross-sections using coil density (CD) and lacunarity (L). These deployments were simulated 9 times each using the original and improved techniques, and CD and L were calculated for cross-sections matching those in the experiments. To compare the 2 simulation techniques, Euclidean distances (dMin, dMax, and dAvg) between experimental and simulation points in standardized CD-L space were evaluated. Univariate tests were performed to determine if these distances were significantly different between the 2 simulations. Results: Coil deployments using the improved technique agreed better with experiments than the original technique. All dMin, dMax, and dAvg values were smaller for the improved technique, and the average values across all simulations for the improved technique were significantly smaller than those from the original technique (dMin: p = 0.014, dMax: p = 0.013, dAvg: p = 0.045). Conclusion: Incorporating coil-specific physical properties and mechanics improves accuracy of FEM simulations of endovascular intracranial aneurysm coiling.