How do coil configuration and packing density influence intra-aneurysmal hemodynamics?

CFD: computational fluid dynamics or confounding factor dissemination? The role of hemodynamics in intracranial aneurysm rupture risk assessment

Image-based computational fluid dynamics holds a prominent position in the evaluation of intracranial aneurysms, especially as a promising tool to stratify rupture risk. Current computational fluid dynamics findings correlating both high and low wall shear stress with intracranial aneurysm growth and rupture puzzle researchers and clinicians alike. These conflicting findings may stem from inconsistent parameter definitions, small datasets, and intrinsic complexities in intracranial aneurysm growth and rupture. In Part 1 of this 2-part review, we proposed a unifying hypothesis: both high and low wall shear stress drive intracranial aneurysm growth and rupture through mural cell-mediated and inflammatory cell-mediated destructive remodeling pathways, respectively. In the present report, Part 2, we delineate different wall shear stress parameter definitions and survey recent computational fluid dynamics studies, in light of this mechanistic heterogeneity. In the future, we expect that larger datasets, better analyses, and increased understanding of hemodynamic-biologic mechanisms will lead to more accurate predictive models for intracranial aneurysm risk assessment from computational fluid dynamics.

Numerical blood flow simulation in surgical corrections: what do we need for an accurate analysis?

BACKGROUND

Computational fluid dynamics has been increasingly used in congenital heart surgery to simulate pathophysiological blood flow, investigate surgical options, or design medical devices. Several commercial and research computational or numerical codes have been developed. They present different approaches to numerically solve the blood flow equations, raising the question whether these numerical codes are equally reliable to achieve accurate simulation results. Accordingly, we sought to examine the influence of numerical code selection in several complex congenital cardiac operations.

MATERIAL AND METHODS

The main steps of blood flow simulations are detailed (geometrical mesh, boundary conditions, and solver numerical methods) for congenital cardiac operations of increasing complexity. The first case tests different numerical solutions against an analytical, or exact, solution. In the second case, the three-dimensional domain is a patient-specific superior cavopulmonary anastomosis. As an analytical solution does not exist in such a complex geometry, different numerical solutions are compared. Finally, a realistic case of a systemic-to-pulmonary shunt is presented with both geometrically and physiologically challenging conditions. For all, solutions from a commercially available code and an open-source research code are compared.

RESULTS

In the first case, as the mesh or solver numerical method is refined, the simulation results for both codes converged to the analytical solution. In the second example, velocity differences between the two codes are greater when the resolution of the mesh were lower and less refined. The third case with realistic anatomy reveals that the pulsatile complex flow is very similar for both codes.

CONCLUSIONS

The precise setup of the numerical cases has more influence on the results than the choice of numerical codes. The need for detailed construction of the numerical model that requires high computational cost depends on the precision needed to answer the biomedical question at hand and should be assessed for each problem on a combination of clinically relevant patient-specific geometry and physiological conditions.

Understanding the role of hemodynamics in the initiation, progression, rupture, and treatment outcome of cerebral aneurysm from medical image-based computational studies

About a decade ago, the first image-based computational hemodynamic studies of cerebral aneurysms were presented. Their potential for clinical applications was the result of a right combination of medical image processing, vascular reconstruction, and grid generation techniques used to reconstruct personalized domains for computational fluid and solid dynamics solvers and data analysis and visualization techniques. A considerable number of studies have captivated the attention of clinicians, neurosurgeons, and neuroradiologists, who realized the ability of those tools to help in understanding the role played by hemodynamics in the natural history and management of intracranial aneurysms. This paper intends to summarize the most relevant results in the field reported during the last years.

Computational fluid dynamics of blood flow in coil-embolized aneurysms: effect of packing density on flow stagnation in an idealized geometry

Coil embolization is performed to induce flow stagnation in cerebral aneurysms and enhance blood clot formation, thus preventing rupture and further growth. We investigated hemodynamics in differently positioned aneurysms coiled at various packing densities to determine the effective packing density in terms of flow stagnation. As a first step, hemodynamic simulations were conducted for idealized geometries of both terminal- and sidewall-type aneurysms. Porous media modeling was employed to describe blood flow in coil-embolized aneurysms. The stagnant volume ratio (SVR) was analyzed to quantify the efficacy of coil embolization. Regardless of aneurysm type and angle, SVR increased with increasing packing density, but the increase in SVR varied depending on type. For sidewall-type aneurysms, the packing density required to achieve 60 % SVR was 20 %, roughly independent of aneurysm angle; flow stagnation was achieved at low packing density. In contrast, in terminal-type aneurysms, the packing density required to achieve 60 % SVR was highly dependent on aneurysm angle, accomplishing a 20 % packing density only for lower angles. Indications are that a relatively high packing density would be required, particularly when these aneurysms are angled against the parent artery. The packing density required for flow stagnation varies depending on aneurysm type and relative position.

A virtual coiling technique for image-based aneurysm models by dynamic path planning

Computational algorithms modeling the insertion of endovascular devices, such as coil or stents, have gained an increasing interest in recent years. This scientific enthusiasm is due to the potential impact that these techniques have to support clinicians by understanding the intravascular hemodynamics and predicting treatment outcomes. In this work, a virtual coiling technique for treating image-based aneurysm models is proposed. A dynamic path planning was used to mimic the structure and distribution of coils inside aneurysm cavities, and to reach high packing densities, which is desirable by clinicians when treating with coils. Several tests were done to evaluate the performance on idealized and image-based aneurysm models. The proposed technique was validated using clinical information of real coiled aneurysms. The virtual coiling technique reproduces the macroscopic behavior of inserted coils and properly captures the densities, shapes and coil distributions inside aneurysm cavities. A practical application was performed by assessing the local hemodynamic after coiling using computational fluid dynamics (CFD). Wall shear stress and intra-aneurysmal velocities were reduced after coiling. Additionally, CFD simulations show that coils decrease the amount of contrast entering the aneurysm and increase its residence time.

AngioLab--a software tool for morphological analysis and endovascular treatment planning of intracranial aneurysms

Determining whether and how an intracranial aneurysm should be treated is a tough decision that clinicians face everyday. Emerging computational tools could help clinicians analyze clinical data and make these decisions. AngioLab is a single graphical user interface, developed on top of the open source framework GIMIAS, that integrates some of the latest image analysis and computational modeling tools for intracranial aneurysms. Two workflows are available: Advanced Morphological Analysis (AMA) and Endovascular Treatment Planning (ETP). AngioLab has been evaluated by a total of 62 clinicians, who considered the information provided by AngioLab relevant and meaningful. They acknowledged the emerging need of these type of tools and the potential impact they might have on the clinical decision-making process.

Analysis and quantification of endovascular coil distribution inside saccular aneurysms using histological images

OBJECTIVE

Endovascular coiling is often performed by first placing coils along the aneurysm wall to create a frame and then by filling up the aneurysm core. However, little attention has been paid to quantifying this filling strategy and to see how it changes for different packing densities. The purpose of this work is to analyze and quantify endovascular coil distribution inside aneurysms based on serial histological images of experimental aneurysms.

METHOD

Seventeen histological images from 10 elastase-induced saccular aneurysms in rabbits treated with coils were studied. In-slice coil density, defined as the area taken up by coil winds, was calculated on each histological image. Images were analyzed by partitioning the aneurysm along its longitudinal and radial axes. Coil distribution was quantified by measuring and comparing the in-slice coil density of each partition.

RESULTS

Mean total in-slice coil density was 22.0 ± 6.2% (range 10.1-30.2%). The density was non-significantly different (p = 0.465) along the longitudinal axis. A significant difference (p < 0.001) between peripheral and core densities was found. Additionally, the peripheral-core density ratio was observed to be inversely proportional to the total in-slice coil density (R(2)=0.57, p <0.001). This ratio was near unity for high in-slice coil density (around 30%).

CONCLUSIONS

These findings demonstrate and confirm that coils tend to be located near the aneurysm periphery when few are inserted. However, when more coils are added, the radial distribution becomes more homogeneous. Coils are homogeneously distributed along the longitudinal axis.

Haemodynamic changes induced by intrasaccular packing on intracranial aneurysms: A computational fluid dynamic study

Endovascular treatment of intracranial aneurysms is a routine medical practice. The most widely used technique is the packing the aneurysm sac with an embolic material. To gain deeper understanding in the effects of specific treatment methods, the intra-aneurysmal haemodynamics are studied with the help of patient-specific computational models. Numerical simulations demonstrated that embolisation with liquid polymer results in an overall decrease of the wall shear stress and pressure in the aneurysm region. Within the range of clinically relevant packing density, simulation of coil embolisation showed homogenisation and decrease of the wall loads on the aneurysm sac. Increasing the packing density above 20% produces little or no further reduction of intra-aneurysmal flow. Sufficient packing of the aneurysm sac results in significant intra-aneurysmal flow decrease associated with reduced wall loads but locally increased pressure or wall shear stress zones may appear depending on the specific vessel geometry.