On the major role played by the lumen curvature of intracranial aneurysms walls in determining their mechanical response, local hemodynamics, and rupture likelihood

The properties of intracranial aneurysms (IAs) walls are known to be driven by the underlying hemodynamics adjacent to the IA sac. Different pathways exist explaining the connections between hemodynamics and local tissue properties. The emergence of such theories is essential if one wishes to compute the mechanical response of a patient-specific IA wall and predict its rupture. Apart from the hemodynamics and tissue properties, one could assume that the mechanical response also depends on the local morphology, more specifically, the curvature of the luminal surface, with larger values at highly-curved wall portions. Nonetheless, this contradicts observations of IA rupture sites more often found at the dome, where the curvature is lower. This seeming contradiction indicates a complex interaction between the hemodynamics adjacent to the aneurysm wall, its morphology, and mechanical response, which warrants further investigation. This was the main goal of this work. We accomplished this by analyzing the stress and stretch fields in different regions of the wall for a sample of IAs, which have been classified based on particular hemodynamics conditions and lumen curvature. Pulsatile numerical simulations were performed using the one-way fluid-solid interaction strategy implemented in OpenFOAM (solids4foam toolbox). We found that the variable best correlated with regions of high stress and stretch was the lumen curvature. Additionally, our data suggest a connection between the local curvature and particular hemodynamics conditions adjacent to the wall, indicating that the lumen curvature is a property that could be used to assess both mechanical response and hemodynamic conditions, and, moreover, suggest new rupture indicators based on the curvature.

A numerical investigation of the mechanics of intracranial aneurysms walls: Assessing the influence of tissue hyperelastic laws and heterogeneous properties on the stress and stretch fields

Numerical simulations have been extensively used in the past two decades for the study of intracranial aneurysms (IAs), a dangerous disease that occurs in the arteries that reach the brain and affect overall 3.2% of a population without comorbidity with up to 60% mortality rate, in case of rupture. The majority of those studies, though, assumed a rigid-wall model to simulate the blood flow. However, to also study the mechanics of IAs walls, it is important to assume a fluid-solid interaction (FSI) modeling. Progress towards more reliable FSI simulations is limited because FSI techniques pose severe numerical difficulties, but also due to scarce data on the mechanical behavior and material constants of IA tissue. Additionally, works that have investigated the impact of different wall modeling choices for patient-specific IAs geometries are a few and often with limited conclusions. Thus our present study investigated the effect of different modeling approaches to simulate the motion of an IA. We used three hyperelastic laws - the Yeoh law, the three-parameter Mooney-Rivlin law, and a Fung-like law with a single parameter - and two different ways of modeling the wall thickness and tissue mechanical properties - one assumed that both were uniform while the other accounted for the heterogeneity of the wall by using a "hemodynamics-driven" approach in which both thickness and material constants varied spatially with the cardiac-cycle-averaged hemodynamics. Pulsatile numerical simulations, with patient-specific vascular geometries harboring IAs, were carried out using the one-way fluid-solid interaction solution strategy implemented in solids4foam, an extension of OpenFOAM®, in which the blood flow is solved and applied as the driving force of the wall motion. We found that different wall morphology models yield smaller absolute differences in the mechanical response than different hyperelastic laws. Furthermore, the stretch levels of IAs walls were more sensitive to the hyperelastic and material constants than the stress. These findings could be used to guide modeling decisions on IA simulations, since the computational behavior of each law was different, for example, with the Yeoh law being the fastest to converge.