Intracranial aneurysmal pulsatility as a new individual criterion for rupture risk evaluation: biomechanical and numeric approach (IRRAs Project)

Modeling and evaluation of biomechanics and hemodynamic based on patient-specific small intracranial aneurysm using fluid-structure interaction

BACKGROUND AND OBJECTIVE

Rupture of small intracranial aneurysm (IA) often leads to the development of highly fatal clinical syndromes such as subarachnoid hemorrhage. Due to the patient specificity of small IA, there are many difficulties in evaluating the rupture risk of small IA such as multiple influencing factors, high clinical experience requirements and poor reusability.

METHODS

In this study, clinical methods such as transcranial doppler (TCD) and magnetic resonance imaging (MRI) are used to obtain patient-specific parameters, and the fluid-structure interaction method (FSI) is used to model and evaluate the biomechanics and hemodynamics of patient-specific small IA.

RESULTS

The results show that a spiral vortex stably exists in the patient-specific small IA. Due to the small size of the patient-specific small IA, the blood flow velocity still maintains a high value with maximum reaching 3 m/s. The inertial impact of blood flow and vortex convection have certain influence on hemodynamic and biomechanics parameters. They cause three high value areas of WSSM on the patient-specific small IA with maximum of 180 Pa, 130 Pa and 110 Pa, respectively. They also cause two types of WSS concentration points, positive normal stress peak value areas and negative normal stress peak value areas to appear.

CONCLUSION

This paper found that the factors affecting hemodynamic parameters and biomechanical parameters are different. Unlike hemodynamic parameters, biomechanical parameters are also affected by blood pressure in addition to blood flow velocity. This study reveals the relationship between the flow field distribution and changes of patient-specific small IA, biomechanics and hemodynamics.

Assessing a heterogeneous model for accounting for endovascular devices in hemodynamic simulations of cerebral aneurysms

The heterogeneous model developed by Berod et al [Int J Numer Method Biomed Eng 38, 2021] for representing the hemodynamic effects of endovascular prostheses is applied to a series of 10 patient specific cerebral aneurysms, 6 being treated by flow diverters, 4 being equipped with WEBs. Two markers correlated with the medical outcome of the treatment are used to assess the potential of the model, namely the saccular mean velocity and the inflow rate at the neck of the aneurysm. The comparison with the corresponding wire-resolved simulations is very favorable in both cases, and the model-based simulations also retrieve the jetting-type flows generated downstream of the struts. Noteworthy, the very same model was used for representing the flow diverters and the WEBs, showing the versatility and robustness of the heterogeneous modeling of the devices.

Endothelial Responses to Curvature-Induced Flow Patterns in Engineered Cerebral Aneurysms

Hemodynamic factors have long been associated with clinical outcomes in the treatment of cerebral aneurysms. Computational studies of cerebral aneurysm hemodynamics have provided valuable estimates of the mechanical environment experienced by the endothelium in both the parent vessel and aneurysmal dome walls and have correlated them with disease state. These computational-clinical studies have recently been correlated with the response of endothelial cells (EC) using either idealized or patient-specific models. Here, we present a robust workflow for generating anatomic-scale aneurysm models, establishing luminal cultures of ECs at physiological relevant flow profiles, and comparing EC responses to curvature mediated flow. We show that flow patterns induced by parent vessel curvature produce changes in wall shear stress (WSS) and wall shear stress gradients (WSSG) that are correlated with differences in cell morphology and cellular protein localization. Cells in higher WSS regions align better with the flow and display strong Notch1-extracellular domain (ECD) polarization, while, under low WSS, differences in WSSG due to curvature change were associated with less alignment and attenuation of Notch1-ECD polarization in ECs of the corresponding regions. These proof-of-concept results highlight the use of engineered cellularized aneurysm models for connecting computational fluid dynamics to the underlying endothelial biology that mediates disease.

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.

The quantitative comparison between high wall shear stress and high strain in the formation of paraclinoid aneurysms

In the hemodynamic study, computational fluid dynamics (CFD) analysis has shown that high wall shear stress (WSS) is an important parameter in cerebral aneurysm formation. However, CFD analysis is not more realistic than fluid–structure interaction (FSI) analysis given its lack of considering the involvement of vascular structures. To investigate the relationship between the hemodynamic parameters and the aneurysm formation, the locations of high WSS and high strain were extracted from the CFD and FSI analyses, respectively. Then the distances between the aneurysm formation site and the locations of high WSS or high strain were calculated. A total of 37 intracranial paraclinoid aneurysms were enrolled for quantitative comparison. Additionally, the dura mater was modeled to facilitate realistic results in FSI analysis. The average distance from the location of the aneurysm formation site to the high strain (1.74 mm \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm $$\end{document}± 1.04 mm) was smaller than the average distance to the high WSS (3.33 mm \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm $$\end{document}± 1.18 mm). The presence of dura mater also influenced the findings in the aneurysm formation site. High strain extracted by FSI analysis is an important hemodynamic factor related to the formation of cerebral aneurysms. Strain parameter could help to predict the formation of aneurysms and elucidate the appropriate treatment.

Flow-splitting-based computation of outlet boundary conditions for improved cerebrovascular simulation in multiple intracranial aneurysms

PURPOSE

Image-based hemodynamic simulations have great potential for precise blood flow predictions in intracranial aneurysms. Due to model assumptions and simplifications with respect to boundary conditions, clinical acceptance remains limited.

METHODS

Within this study, we analyzed the influence of outflow-splitting approaches on multiple aneurysm studies and present a new outflow-splitting approach that takes the precise morphological vessel cross sections into account. We provide a detailed comparison of five outflow strategies considering eight intracranial aneurysms: zero-pressure configuration (1), a flow splitting inspired by Murray's law with a square (2) and a cubic (3) vessel diameter, a flow splitting incorporating vessel bifurcations based on circular vessel cross sections (4) and our novel flow splitting including vessel bifurcations and anatomical vessel cross sections (5). Other boundary conditions remain constant. For each simulation and each aneurysm, we conducted an evaluation based on common hemodynamic parameters, e.g., normalized wall shear stress and inflow concentration index.

RESULTS

The comparison of five outflow strategies for image-based simulations shows a large variability regarding the parameters of interest. Qualitatively, our strategy based on anatomical cross sections yields a more uniform flow rate distribution with increased aneurysm inflow rates. The commonly used zero-pressure approach shows the largest variations, especially for more distal aneurysms. A rank ordering of multiple aneurysms in one patient might still be possible, since the ordering appeared to be independent of the outflow strategy.

CONCLUSIONS

The results reveal that outlet boundary conditions have a crucial impact on image-based blood flow simulations, especially for multiple aneurysm studies. We could confirm the advantages of the more complex outflow-splitting model (4) including an incremental improvement (5) compared to strategies (1), (2) and (3) for this application scenario. Furthermore, we discourage from using zero-pressure configurations that lack a physiological basis.

Backward sensitivity analysis and reduced-order covariance estimation in noninvasive parameter identification for cerebral arteries

Using a previously developed inversion platform for functional cerebral medical imaging with ensemble Kalman filters, this work analyzes the sensitivity of the results with respect to different parameters entering the physical model and inversion procedure, such as the inlet flow rate from the heart, the choice of the boundary conditions, and the nonsymmetry in the network terminations. It also proposes an alternative low complexity construction for the covariance matrix of the hemodynamic parameters of a network of arteries including the circle of Willis. The platform takes as input patient-specific blood flow rates extracted from magnetic resonance angiography and magnetic resonance imaging (dicom files) and is applied to several available patients data. The paper presents full analysis of the results for one of these patients, including a sensitivity study with respect to the proximal and distal boundary conditions. The results notably show that the uncertainties on the inlet flow rate led to uncertainties of the same order of magnitude in the estimated parameters (blood pressure and elastic parameters) and that three-lumped parameters boundary conditions are necessary for a correct retrieval of the target signals.

Hemodynamics in a giant intracranial aneurysm characterized by in vitro 4D flow MRI

Experimental and computational data suggest that hemodynamics play a critical role in the development, growth, and rupture of cerebral aneurysms. The flow structure, especially in aneurysms with a large sac, is highly complex and three-dimensional. Therefore, volumetric and time-resolved measurements of the flow properties are crucial to fully characterize the hemodynamics. In this study, phase-contrast Magnetic Resonance Imaging is used to assess the fluid dynamics inside a 3D-printed replica of a giant intracranial aneurysm, whose hemodynamics was previously simulated by multiple research groups. The physiological inflow waveform is imposed in a flow circuit with realistic cardiovascular impedance. Measurements are acquired with sub-millimeter spatial resolution for 16 time steps over a cardiac cycle, allowing for the detailed reconstruction of the flow evolution. Moreover, the three-dimensional and time-resolved pressure distribution is calculated from the velocity field by integrating the fluid dynamics equations, and is validated against differential pressure measurements using precision transducers. The flow structure is characterized by vortical motions that persist within the aneurysm sac for most of the cardiac cycle. All the main flow statistics including velocity, vorticity, pressure, and wall shear stress suggest that the flow pattern is dictated by the aneurysm morphology and is largely independent of the pulsatility of the inflow, at least for the flow regimes investigated here. Comparisons are carried out with previous computational simulations that used the same geometry and inflow conditions, both in terms of cycle-averaged and systolic quantities.

Combined Visualization of Vessel Deformation and Hemodynamics in Cerebral Aneurysms

We present the first visualization tool that combines patient-specific hemodynamics with information about the vessel wall deformation and wall thickness in cerebral aneurysms. Such aneurysms bear the risk of rupture, whereas their treatment also carries considerable risks for the patient. For the patient-specific rupture risk evaluation and treatment analysis, both morphological and hemodynamic data have to be investigated. Medical researchers emphasize the importance of analyzing correlations between wall properties such as the wall deformation and thickness, and hemodynamic attributes like the Wall Shear Stress and near-wall flow. Our method uses a linked 2.5D and 3D depiction of the aneurysm together with blood flow information that enables the simultaneous exploration of wall characteristics and hemodynamic attributes during the cardiac cycle. We thus offer medical researchers an effective visual exploration tool for aneurysm treatment risk assessment. The 2.5D view serves as an overview that comprises a projection of the vessel surface to a 2D map, providing an occlusion-free surface visualization combined with a glyph-based depiction of the local wall thickness. The 3D view represents the focus upon which the data exploration takes place. To support the time-dependent parameter exploration and expert collaboration, a camera path is calculated automatically, where the user can place landmarks for further exploration of the properties. We developed a GPU-based implementation of our visualizations with a flexible interactive data exploration mechanism. We designed our techniques in collaboration with domain experts, and provide details about the evaluation.

Fluid-Structure Simulations of a Ruptured Intracranial Aneurysm: Constant versus Patient-Specific Wall Thickness

Computational Fluid Dynamics is intensively used to deepen the understanding of aneurysm growth and rupture in order to support physicians during therapy planning. However, numerous studies considering only the hemodynamics within the vessel lumen found no satisfactory criteria for rupture risk assessment. To improve available simulation models, the rigid vessel wall assumption has been discarded in this work and patient-specific wall thickness is considered within the simulation. For this purpose, a ruptured intracranial aneurysm was prepared ex vivo, followed by the acquisition of local wall thickness using μCT. The segmented inner and outer vessel surfaces served as solid domain for the fluid-structure interaction (FSI) simulation. To compare wall stress distributions within the aneurysm wall and at the rupture site, FSI computations are repeated in a virtual model using a constant wall thickness approach. Although the wall stresses obtained by the two approaches—when averaged over the complete aneurysm sac—are in very good agreement, strong differences occur in their distribution. Accounting for the real wall thickness distribution, the rupture site exhibits much higher stress values compared to the configuration with constant wall thickness. The study reveals the importance of geometry reconstruction and accurate description of wall thickness in FSI simulations.

A NUMERICAL STUDY ON THE HEMODYNAMIC AND SHEAR STRESS OF DOUBLE ANEURYSM THROUGH S-SHAPED VESSEL

Abdominal aortic aneurysm (AAA) is a degenerative disease defined as the abnormal ballooning of the abdominal aorta (AA) wall which is usually caused by atherosclerosis. The aneurysm grows larger and eventually ruptures if it is not diagnosed and treated. Aneurysms occur mostly in the aorta, the main artery of the chest and abdomen. The aorta carries blood flow from the heart to all parts of the body, including the vital organs, the legs, and feet. The objective of the present study is to investigate the combined effects of aneurysm and curvature on flow characteristics in S-shaped bends with sweep angle of 90° at Reynolds number of 900. The fluid mechanics of blood flow in a curved artery with abnormal aortic is studied through a mathematical analysis and employing Cosmos flow simulation. Blood is modeled as an incompressible non-Newtonian fluid and the flow is assumed to be steady and laminar. Hemodynamic characteristics are analyzed. Grid independence is tested on three successively refined meshes. It is observed that the abrupt expansion induced by AAA results in an immensely disturbed regime. The results may have implications not only for understanding the mechanical behavior of the blood flow inside an aneurysm artery but also for investigating the mechanical behavior of the blood flow in different arterial diseases, such as atherosclerosis.