A comprehensive numerical approach to coil placement in cerebral aneurysms: mathematical modeling and in silico occlusion classification

Endovascular coil embolization is one of the primary treatment techniques for cerebral aneurysms. Although it is a well-established and minimally invasive method, it bears the risk of suboptimal coil placement which can lead to incomplete occlusion of the aneurysm possibly causing recurrence. One of the key features of coils is that they have an imprinted natural shape supporting the fixation within the aneurysm. For the spatial discretization, our mathematical coil model is based on the discrete elastic rod model which results in a dimension-reduced 1D system of differential equations. We include bending and twisting responses to account for the coils natural curvature and allow for the placement of several coils having different material parameters. Collisions between coil segments and the aneurysm wall are handled by an efficient contact algorithm that relies on an octree based collision detection. In time, we use a standard symplectic semi-implicit Euler time stepping method. Our model can be easily incorporated into blood flow simulations of embolized aneurysms. In order to differentiate optimal from suboptimal placements, we employ a suitable in silico Raymond-Roy-type occlusion classification and measure the local packing density in the aneurysm at its neck, wall region and core. We investigate the impact of uncertainties in the coil parameters and embolization procedure. To this end, we vary the position and the angle of insertion of the micro-catheter, and approximate the local packing density distributions by evaluating sample statistics.