Multi-Dimensional Modeling of Cerebral Hemodynamics: A Systematic Review.

The Circle of Willis (CoW) describes the arterial system in the human brain enabling the neurovascular blood supply. Neurovascular diseases like intracranial aneurysms (IAs) can occur within the CoW and carry the risk of rupture, which can lead to subarachnoid hemorrhage. The assessment of hemodynamic information in these pathologies is crucial for their understanding regarding detection, diagnosis and treatment. Multi-dimensional in silico approaches exist to evaluate these hemodynamics based on patient-specific input data. The approaches comprise low-scale (zero-dimensional, one-dimensional) and high-scale (three-dimensional) models as well as multi-scale coupled models. The input data can be derived from medical imaging, numerical models, literature-based assumptions or from measurements within healthy subjects. Thus, the most realistic description of neurovascular hemodynamics is still controversial. Within this systematic review, first, the models of the three scales (0D, 1D, 3D) and second, the multi-scale models, which are coupled versions of the three scales, were discussed. Current best practices in describing neurovascular hemodynamics most realistically and their clinical applicablility were elucidated. The performance of 3D simulation entails high computational expenses, which could be reduced by analyzing solely the region of interest in detail. Medical imaging to establish patient-specific boundary conditions is usually rare, and thus, lower dimensional models provide a realistic mimicking of the surrounding hemodynamics. Multi-scale coupling, however, is computationally expensive as well, especially when taking all dimensions into account. In conclusion, the 0D-1D-3D multi-scale approach provides the most realistic outcome; nevertheless, it is least applicable. A 1D-3D multi-scale model can be considered regarding a beneficial trade-off between realistic results and applicable performance.

Assessment of intracranial aneurysm neck deformation after contour deployment.

PURPOSE: The contour neurovascular system (CNS) is a novel device to treat intracranial wide-necked bifurcation aneurysms, with few studies assessing its long-term effects. Particularly its impact on aneurysm morphology has not been explored yet. We present a preliminary study to explore this impact for the first time, focusing on the neck curve and ostium of the aneurysm. METHODS: We investigated seven aneurysms treated with the CNS to assess ostium deformation after CNS deployment by comparing models extracted from in vivo medical pre-treatment and follow-up scans via morphological analysis. Time between pre- and follow-up scans was ten months on average. Size and shape indices like area, neck diameter, ellipticity index, undulation index, and more were assessed. RESULTS: Ostium size was reduced after treatment. On average, ostium area was reduced at a rate of - 0.58 (± 4.88) mm2 per year, from 15.52 (± 3.51) mm2 to 13.30 (± 2.27) mm2, and ostium width from 5.01 (± 0.54) mm to 4.49 (± 0.45) mm, with an average reduction of - 0.59 (± 0.87) mm. This shrinking positively correlated with time passing. Shape deformation was low, though notably mean ellipticity index was reduced by 0.06 (± 0.15) on average, indicating ostia were less elongated after treatment. CONCLUSION: We interpret the shrinking of the ostium as part of the healing process. Shape changes were found to be small enough to conclude no shape deformation of the ostium from CNS deployment, but the analysis of more cases with more parameters and information is necessary.

Image-based hemodynamic simulations for intracranial aneurysms: the impact of complex vasculature.

PURPOSE: Hemodynamics play an important role in the assessment of intracranial aneurysm (IA) development and rupture risk. The purpose of this study was to examine the impact of complex vasculatures onto the intra-vessel and intra-aneurysmal blood flow. METHODS: Complex segmentation of a subject-specific, 60-outlet and 3-inlet circle of Willis model captured with 7T magnetic resonance imaging was performed. This model was trimmed to a 10-outlet model version. Two patient-specific IAs were added onto both models yielding two pathological versions, and image-based blood flow simulations of the four resulting cases were carried out. To capture the differences between complex and trimmed model, time-averaged and centerline velocities were compared. The assessment of intra-saccular blood flow within the IAs involved the evaluation of wall shear stresses (WSS) at the IA wall and neck inflow rates (NIR). RESULTS: Lower flow values are observed in the majority of the complex model. However, at specific locations (left middle cerebral artery 0.5 m/s, left posterior cerebral artery 0.25 m/s), higher flow rates were visible when compared to the trimmed counterpart. Furthermore, at the centerlines the total velocity values reveal differences up to 0.15 m/s. In the IAs, the reduction in the neck inflow rate and WSS in the complex model was observed for the first IA (IA-A δNIRmean = - 0.07ml/s, PCA.l δWSSmean = - 0.05 Pa). The second IA featured an increase in the neck inflow rate and WSS (IA-B δNIRmean = 0.04 ml/s, PCA.l δWSSmean = 0.07 Pa). CONCLUSION: Both the magnitude and shape of the flow distribution vary depending on the model's complexity. The magnitude is primarily influenced by the global vessel model, while the shape is determined by the local structure. Furthermore, intra-aneurysmal flow strongly depends on the location in the vessel tree, emphasizing the need for complex model geometries for realistic hemodynamic assessment and rupture risk analysis.

Fusiform versus Saccular Intracranial Aneurysms-Hemodynamic Evaluation of the Pre-Aneurysmal, Pathological, and Post-Interventional State.

Minimally-invasive therapies are well-established treatment methods for saccular intracranial aneurysms (SIAs). Knowledge concerning fusiform IAs (FIAs) is low, due to their wide and alternating lumen and their infrequent occurrence. However, FIAs carry risks like ischemia and thus require further in-depth investigation. Six patient-specific IAs, comprising three position-identical FIAs and SIAs, with the FIAs showing a non-typical FIA shape, were compared, respectively. For each model, a healthy counterpart and a treated version with a flow diverting stent were created. Eighteen time-dependent simulations were performed to analyze morphological and hemodynamic parameters focusing on the treatment effect (TE). The stent expansion is higher for FIAs than SIAs. For FIAs, the reduction in vorticity is higher (Δ35-75% case 2/3) and the reduction in the oscillatory velocity index is lower (Δ15-68% case 2/3). Velocity is reduced equally for FIAs and SIAs with a TE of 37-60% in FIAs and of 41-72% in SIAs. Time-averaged wall shear stress (TAWSS) is less reduced within FIAs than SIAs (Δ30-105%). Within this study, the positive TE of FDS deployed in FIAs is shown and a similarity in parameters found due to the non-typical FIA shape. Despite the higher stent expansion, velocity and vorticity are equally reduced compared to identically located SIAs.

Is Accurate Lumen Segmentation More Important than Outlet Boundary Condition in Image-Based Blood Flow Simulations for Intracranial Aneurysms?

PURPOSE: Image-based blood flow simulations are increasingly used to investigate the hemodynamics in intracranial aneurysms (IAs). However, a strong variability in segmentation approaches as well as the absence of individualized boundary conditions (BCs) influence the quality of these simulation results leading to imprecision and decreased reliability. This study aims to analyze these influences on relevant hemodynamic parameters within IAs. METHODS: As a follow-up study of an international multiple aneurysms challenge, the segmentation results of five IAs differing in size and location were investigated. Specifically, five possible outlet BCs were considered in each of the IAs. These are comprised of the zero-pressure condition (BC1), a flow distribution based on Murray's law with the exponents n = 2 (BC2) and n = 3 (BC3) as well as two advanced flow-splitting models considering the real vessels by including circular cross sections (BC4) or anatomical cross sections (BC5), respectively. In total, 120 time-dependent blood flow simulations were analyzed qualitatively and quantitatively, focusing on five representative intra-aneurysmal flow and five shear parameters such as vorticity and wall shear stress. RESULTS: The outlet BC variation revealed substantial differences. Higher shear stresses (up to Δ9.69 Pa), intrasaccular velocities (up to Δ0.15 m/s) and vorticities (up to Δ629.22 1/s) were detected when advanced flow-splitting was applied compared to the widely used zero-pressure BC. The tendency of outlets BCs to over- or underestimate hemodynamic parameters is consistent across different segmentations of a single aneurysm model. Segmentation-induced variability reaches Δ19.58 Pa, Δ0.42 m/s and Δ957.27 1/s, respectively. Excluding low fidelity segmentations, however, (a) reduces the deviation drastically (>43%) and (b) leads to a lower impact of the outlet BC on hemodynamic predictions. CONCLUSION: With a more realistic lumen segmentation, the influence of the BC on the resulting hemodynamics is decreased. A realistic lumen segmentation can be ensured, e.g., by using high-resolved 2D images. Furthermore, the selection of an advanced outflow-splitting model is advised and the use of a zero-pressure BC and BC based on Murray's law with exponent n = 3 should be avoided.

In vitro and in silico assessment of flow modulation after deploying the Contour Neurovascular System in intracranial aneurysm models.

BACKGROUND: The novel Contour Neurovascular System (Contour) has been reported to be efficient and safe for the treatment of intracranial, wide-necked bifurcation aneurysms. Flow in the aneurysm and posterior cerebral arteries (PCAs) after Contour deployment has not been analyzed in detail yet. However, this information is crucial for predicting aneurysm treatment outcomes. METHODS: Time-resolved three-dimensional velocity maps in 14 combinations of patient-based basilar tip aneurysm models with and without Contour devices (sizes between 5 and 14 mm) were analyzed using four-dimensionsal (4D) flow MRI and numerical/image-based flow simulations. A complex virtual processing pipeline was developed to mimic the experimental shape and position of the Contour together with the simulations. RESULTS: On average, the Contour significantly reduced intra-aneurysmal flow velocity by 67% (mean w/ = 0.03m/s; mean w/o = 0.12m/s; p-value=0.002), and the time-averaged wall shear stress by more than 87% (mean w/ = 0.17Pa; mean w/o = 1.35Pa; p-value=0.002), as observed by numerical simulations. Furthermore, a significant reduction in flow (P<0.01) was confirmed by the neck inflow rate, kinetic energy, and inflow concentration index after Contour deployment. Notably, device size has a stronger effect on reducing flow than device positioning. However, positioning affected flow in the PCAs, while being robust in effectively reducing flow. CONCLUSIONS: This study showed the high efficacy of the Contour device in reducing flow within aneurysms regardless of the exact position. However, we observed an effect on the flow in PCAs, which needs to be investigated further.

Intraventricular Flow Simulations in Singular Right Ventricles Reveal Deteriorated Washout and Low Vortex Formation.

PURPOSE: Patients with a functionally univentricular heart represent one of the most common severe cardiac lesions with a prevalence of 3 per 10,000 live births. Hemodynamics of the singular ventricle is a major research topic in cardiology and there exists a relationship between fluid dynamical features and cardiac behavior in health and disease. The aim of the present work was to compare intraventricular flow in single right ventricle (SRV) patients and subjects with healthy left hearts (LV) through patient-specific CFD simulations. METHODS: Three-dimensional real-time echocardiographic images were obtained for five SRV patients and two healthy subjects and CFD simulations with a moving mesh methodology were performed. Intraventricular vortex formation and vortex formation time (VFT) as well as the turbulent kinetic energy (TKE) and ventricular washout were evaluated. RESULTS: The results show significantly lower values for the VFT and the TKE in SRV patients compared with healthy LV subjects. Furthermore, vortex formation does not progress to the apex in SRV patients. These findings were confirmed by a significantly lower washout in SRV patients. CONCLUSIONS: The study pinpoints the intriguing role of intraventricular flows to characterize performance of SRVs that goes beyond standard clinical metrics such as ejection fraction.

Can Endovascular Treatment of Fusiform Intracranial Aneurysms Restore the Healthy Hemodynamic Environment?-A Virtual Pilot Study.

Numerous studies assess intracranial aneurysm rupture risk based on morphological and hemodynamic parameter analysis in addition to clinical information such as aneurysm localization, age, and sex. However, intracranial aneurysms mostly occur with a saccular shape located either lateral to the parent artery or at a bifurcation. In contrast, fusiform intracranial aneurysms (FIAs), i.e., aneurysms with a non-saccular, dilated form, occur in approximately 3-13% of all cases and therefore have not yet been as thoroughly studied. To improve the understanding of FIA hemodynamics, this pilot study contains morphological analyses and image-based blood flow simulations in three patient-specific cases. For a precise and realistic comparison to the pre-pathological state, each dilation was manually removed and the time-dependent blood flow simulations were repeated. Additionally, a validated fast virtual stenting approach was applied to evaluate the effect of virtual endovascular flow-diverter deployment focusing on relevant hemodynamic quantities. For two of the three patients, post-interventional information was available and included in the analysis. The results of this numerical pilot study indicate that complex flow structures, i.e., helical flow phenomena and the presence of high oscillating flow features, predominantly occur in FIAs with morphologically differing appearances. Due to the investigation of the individual healthy states, the original flow environment could be restored which serves as a reference for the virtual treatment target. It was shown that the realistic deployment led to a considerable stabilization of the individual hemodynamics in all cases. Furthermore, a quantification of the stent-induced therapy effect became feasible for the treating physician. The results of the morphological and hemodynamic analyses in this pilot study show that virtual stenting can be used in FIAs to quantify the effect of the planned endovascular treatment.