Aneurysm growth occurs at region of low wall shear stress: patient-specific correlation of hemodynamics and growth in a longitudinal study

Computational analysis to assess hemodynamic forces in descending thoracic aortic aneurysms

Descending thoracic aortic aneurysm (DTAA) is a life-threatening disorder, defined as a localized enlargement of the descending portion of the thoracic aorta. In this context, we develop a fluid-structure interaction (FSI) computational framework, with the inclusion of a turbulence model and different material properties for the healthy and the aneurysmatic portions of the vessel, to study the hemodynamics and its relationship with DTAA. We first provide an analysis on nine ideal scenarios, accounting for different aortic arch types and DTAA ubications, to study changes in blood pressure, flow patterns, turbulence, wall shear stress, drag forces and internal wall stresses. Our findings demonstrate that the hemodynamics in DTAA is profoundly disturbed, with the presence of flow re-circulation, formation of vortices and transition to turbulence. In particular, configurations with a steeper aortic arch exhibit a more chaotic hemodynamics. We notice also an increase in pressure values for configurations with less steep aortic arch and in drag forces for configurations with distal DTAA. Second, we replicate our analysis for three patient-specific cases (one for type of arch) obtaining comforting results in terms of accordance with the ideal scenarios. Finally, in a very preliminary way, we try to relate our findings to possible stent-graft migrations after TEVAR procedure to provide predictions on the postoperative state. KEY POINTS: This study employs computational methods to assess hemodynamic forces in descending thoracic aortic aneurysms. We consider ideal cases by varying aortic arch type and aneurysm location. Our results show: chaotic hemodynamics for steep aortic arches; increase in pressure values for less steep aortic arches; high risk of plaque deposition in the aneurysmal sac for proximal aneurysms and near the neck for distal aneurysms. We also analyse three patient-specific cases, confirming the major outcomes found for the ideal cases. We try to suggest how our preoperative findings may correlate to assess the risk of stent-graft migration of a possible TEVAR procedure.

Cerebrovascular morphology: Insights into normal variations, aging effects, and disease implications

Cerebrovascular morphology plays a critical role in brain health, influencing cerebral blood flow (CBF) and contributing to the pathogenesis of various neurological diseases. This review examines the anatomical structure of the cerebrovascular network and its variations in healthy and diseased populations and highlights age-related changes and their implications in various neurological conditions. Normal variations, including the completeness and anatomical anomalies of the Circle of Willis and collateral circulation, are discussed in relation to their impact on CBF and susceptibility to ischemic events. Age-related changes in the cerebrovascular system, such as alterations in vessel geometry and density, are explored for their contributions to age-related neurological disorders, including Alzheimer's disease and vascular dementia. Advances in medical imaging and computational methods have enabled automatic quantitative assessment of cerebrovascular structures, facilitating the identification of pathological changes in both acute and chronic cerebrovascular disorders. Emerging technologies, including machine learning and computational fluid dynamics, offer new tools for predicting disease risk and patient outcomes based on vascular morphology. This review underscores the importance of understanding cerebrovascular remodeling for early diagnosis and the development of novel therapeutic approaches in brain diseases.

Mechanisms of aortic dissection: From pathological changes to experimental and in silico models

Aortic dissection continues to be responsible for significant morbidity and mortality, although recent advances in medical data assimilation and in experimental and in silico models have improved our understanding of the initiation and progression of the accumulation of blood within the aortic wall. Hence, there remains a pressing necessity for innovative and enhanced models to more accurately characterize the associated pathological changes. Early on, experimental models were employed to uncover mechanisms in aortic dissection, such as hemodynamic changes and alterations in wall microstructure, and to assess the efficacy of medical implants. While experimental models were once the only option available, more recently they are also being used to validate in silico models. Based on an improved understanding of the deteriorated microstructure of the aortic wall, numerous multiscale material models have been proposed in recent decades to study the state of stress in dissected aortas, including the changes associated with damage and failure. Furthermore, when integrated with accessible patient-derived medical data, in silico models prove to be an invaluable tool for identifying correlations between hemodynamics, wall stresses, or thrombus formation in the deteriorated aortic wall. They are also advantageous for model-guided design of medical implants with the aim of evaluating the deployment and migration of implants in patients. Nonetheless, the utility of in silico models depends largely on patient-derived medical data, such as chosen boundary conditions or tissue properties. In this review article, our objective is to provide a thorough summary of medical data elucidating the pathological alterations associated with this disease. Concurrently, we aim to assess experimental models, as well as multiscale material and patient data-informed in silico models, that investigate various aspects of aortic dissection. In conclusion, we present a discourse on future perspectives, encompassing aspects of disease modeling, numerical challenges, and clinical applications, with a particular focus on aortic dissection. The aspiration is to inspire future studies, deepen our comprehension of the disease, and ultimately shape clinical care and treatment decisions.

The order of precedence in treatment of multiple intracranial aneurysms: insights from a fluid-structure interaction study

Treatment strategy for multiple intracranial aneurysms is challenging, as in many cases the choice of the order in which to treat aneurysms is not based on high-quality evidence. We aimed to digitally simulate clinical scenarios of two different orders in which multiple intracranial aneurysms were treated and analyze changes in hemodynamics after first stage of treatment. We prospectively included patients with two intracranial aneurysms, with order of treatment difficult to determine based on clinical data alone. For each patient we prepared three models of arteries harboring aneurysms: with both aneurysms present and with one of them removed. Computational modeling of blood flow using fluid-structure interaction methodology was performed for each model. Hemodynamic parameters of aneurysm domes were compared between models with both aneurysm present, and models in with aneurysms were removed in changing order. In 25 included patients, the calculated hemodynamic parameters such as Time-Averaged Wall Shear Stress (0.46 ± 0.40 vs. 0.54 ± 0.44 Pa; p < 0.01) and surface vortex fraction (12.73% ± 7.92% vs. 14.26% ± 7.46%; p = 0.02) decreased after first stage of treatment, while Time-Averaged Wall Shear Stress Gradient (1.44 ± 0.41 vs. 1.34 ± 0.46 Pa; p = 0.04) and percentage of wall shear stress < 0.5 Pa (50.13% ± 33.01% vs. 44.08% ± 34.16%; p < 0.01) increased. Changes of wall shear stress in remaining aneurysm dome were independently correlated with dome-to-neck ratio of both removed and remaining aneurysms. Hemodynamics of untreated aneurysm worsens after first stage of treatment. Dome-to-neck ratio of both treated and untreated aneurysm was the strongest and independent predictor of that change.

Aneurysm Wall Enhancement Can Predict Rupture Point in Intracranial Aneurysms With Multiple Blebs

BACKGROUND AND OBJECTIVES

Ruptured aneurysms visualized by vessel wall MRI (VW-MRI) exhibit characteristic aneurysm wall enhancement (AWE). A secondary bulge of the aneurysmal wall, called a bleb, is often the site of rupture in ruptured aneurysms. We hypothesized that a higher degree of AWE would identify the rupture point in aneurysms with multiple blebs.

METHODS

AWE was quantitatively analyzed in consecutive ruptured intracranial aneurysms with multiple blebs (31 aneurysms with a total of 72 blebs) using VW-MRI. A 3-dimensional T1-weighted fast spin-echo sequence was obtained after contrast media injection, and the contrast ratio of the aneurysm wall against the pituitary stalk (CR stalk ) was calculated as the AWE indicator. Bleb characteristics, including CR stalk and wall shear stress (WSS), were compared between ruptured and unruptured blebs. Odds ratios with 95% confidence intervals for ruptures were calculated by conditional logistic regression analysis.

RESULTS

Ruptured blebs had a higher CR stalk and lower WSS compared with unruptured blebs. CR stalk remained significantly associated with the bleb rupture status in the conditional logistic regression (adjusted odds ratio 3.9, 95% CIs 1.6-9.7).

CONCLUSION

AWE is associated with the bleb rupture status independent of WSS. Contrast-enhanced VW-MRI may be a useful noninvasive tool for identifying the rupture point and guiding the treatment strategy.

Hemodynamic Analysis in Aortic Dilatation after Arterial Switch Operation for Patients with Transposition of Great Arteries Using Computational Fluid Dynamics

After an arterial switch operation for complete transposition of the great arteries, neo-aortic root dilatation occurs, with unclear hemodynamic effects. This study analyzes three groups (severe dilation, mild dilation, and normal) using computational fluid dynamics (CFD) on cardiac CT scans. Aortic arch angles in severe (median 72.3, range: 68.5-77.2) and mild dilation (76.6, 71.1-85.2) groups are significantly smaller than the normal group (97.3, 87.4-99.0). In the normal and mild dilatation groups, Wall Shear Stress (WSS) exhibits a consistent pattern: it is lowest at the aortic root, gradually increases until just before the bend in the aortic arch, peaks, and then subsequently decreases. However, severe dilation shows disrupted WSS patterns, notably lower in the distal ascending aorta, attributed to local recirculation. This unique WSS pattern observed in severely dilated patients, especially in the transverse aorta. CFD plays an essential role in comprehensively studying the pathophysiology underlying aortic dilation in this population.

Hemodynamic Characteristics in Ruptured and Unruptured Intracranial Aneurysms: A Prospective Cohort Study Utilizing the AneurysmFlow Tool

BACKGROUND AND PURPOSE

Hemodynamic factors significantly influence the onset, progression, and rupture of intracranial aneurysms (IAs). Current rupture risk prediction scores focus primarily on the clinical, anatomic, and morphologic aspects. This study aimed to investigate the hemodynamic characteristics differences between ruptured and unruptured IAs.

MATERIALS AND METHODS

Conducted from July 2021 to July 2022, this prospective cohort study involved patients with ruptured and unruptured IAs undergoing DSA. Hemodynamic characteristics were assessed by using the AneurysmFlow tool. Hemodynamic, clinical, anatomic, and morphologic parameters were compared between ruptured and unruptured IA groups.

RESULTS

The study included 127 patients with 135 aneurysms (67 ruptured, 68 unruptured). Complex flow patterns (type 3 and 4) were observed more frequently in ruptured aneurysms compared with unruptured aneurysms (OR, 5.57; 95% CI, 2.49-12.45; P < .001) in univariate analysis, and were also more common in unruptured aneurysms associated with daughter sac features (P = .015). The mean aneurysm flow amplitude (MAFA) was lower in ruptured aneurysms, and associated with lower flow velocity in the parent artery related to vasospasm. MAFA in the aneurysmal dome or any additional daughter sacs was lowest compared with other regions inside the aneurysms. The technical failure rate of AneurysmFlow measurements was 8.5% (12 of 139 patients). Additionally, hypertension (OR, 0.42; 95% CI, 0.30-0.54; P < .001), bifurcation location (anterior communicating artery/anterior cerebral artery/MCA/posterior communicating artery/posterior circulation) (OR, 0.17; 95% CI, 0.05-0.29; P = .005), and irregular shape (OR, 0.19; 95% CI, 0.05-0.35; P = .012) were identified as independently associated with rupture.

CONCLUSIONS

Complex flow patterns identified on the AneurysmFlow tool are significantly more common in ruptured and unruptured aneurysms associated with daughter sac features. The lowest MAFA in the aneurysmal dome and daughter sacs likely indicates specific pathophysiologic changes within the aneurysm wall associated with rupture incidence. Hypertension, bifurcation location, and an irregular shape are independently associated with the risk of rupture. Further multicenter studies with larger sample sizes are needed to validate these findings.

Impact of A1 segment asymmetry on hemodynamic conditions around the circle of Willis and anterior communicating artery aneurysm formation

Background

This study aims to investigate how A1 segment asymmetry-also known as A1 dominancy-influences the development of the anterior communicating artery aneurysm (AcomA) as it affects hemodynamic conditions within the circle of Willis (COW). Using time-of-flight magnetic resonance angiography (TOF-MRA), the research introduces a novel approach to assessing shear stress in A1 segments to uncover the hemodynamic factors contributing to AcomA formation.

Method

An observational study was conducted over 6 years at a tertiary university hospital's outpatient clinic. Recruited patients who underwent TOF-MRA imaging were divided into AcomA and non-AcomA groups. MRA images were analyzed using semi-automatic software (VINT, Mediimg, Inc.) to calculate the signal intensity gradient (SIG), which reflects wall shear stress. The comparison metrics included general demographics, anatomical characteristics, and hemodynamic attributes of the COW, mainly focusing on A1 segment asymmetry.

Results

Among the 700 subjects, 106 were categorized into the AcomA group, while 594 were placed in the non-AcomA group. The AcomA group showed a more significant difference in the bilateral A1 diameter (49.0% vs. 20.8%, p < 0.001) and a greater prevalence of unilateral A1 aplasia (32.1% vs. 6.7%, p < 0.001) compared to the non-AcomA group. Increased bilateral A1 asymmetry in the AcomA group corresponded with notable variations in A1 SIG, indicating increased wall shear stress. The occurrence of AcomA is associated with both anatomical factors of the circle of Willis, represented by the bilateral A1 diameter ratio, and hemodynamic factors, represented by the bilateral A1 SIG ratio, suggesting that both factors are almost equally significant.

Conclusion

Our findings suggest that A1 segment asymmetry influences hemodynamic changes within the COW, contributing to AcomA formation. Hemodynamic factors provide an intuitive understanding of how anatomical characteristics within the COW can lead to aneurysm development.

Computational fluid dynamics and shape analysis enhance aneurysm rupture risk stratification

PURPOSE

Accurately quantifying the rupture risk of unruptured intracranial aneurysms (UIAs) is crucial for guiding treatment decisions and remains an unmet clinical challenge. Computational Flow Dynamics and morphological measurements have been shown to differ between ruptured and unruptured aneurysms. It is not clear if these provide any additional information above routinely available clinical observations or not. Therefore, this study investigates whether incorporating image-derived features into the established PHASES score can improve the classification of aneurysm rupture status.

METHODS

A cross-sectional dataset of 170 patients (78 with ruptured aneurysm) was used. Computational fluid dynamics (CFD) and shape analysis were performed on patients' images to extract additional features. These derived features were combined with PHASES variables to develop five ridge constrained logistic regression models for classifying the aneurysm rupture status. Correlation analysis and principal component analysis were employed for image-derived feature reduction. The dataset was split into training and validation subsets, and a ten-fold cross validation strategy with grid search optimisation and bootstrap resampling was adopted for determining the models' coefficients. Models' performances were evaluated using the area under the receiver operating characteristic curve (AUC).

RESULTS

The logistic regression model based solely on PHASES achieved AUC of 0.63. All models incorporating derived features from CFD and shape analysis demonstrated improved performance, reaching an AUC of 0.71. Non-sphericity index (shape variable) and maximum oscillatory shear index (CFD variable) were the strongest predictors of a ruptured status.

CONCLUSION

This study demonstrates the benefits of integrating image-based fluid dynamics and shape analysis with clinical data for improving the classification accuracy of aneurysm rupture status. Further evaluation using longitudinal data is needed to assess the potential for clinical integration.

Verifying the Accuracy of Hemodynamic Analysis Using High Spatial Resolution 3D Phase-contrast MR Imaging on a 7T MR System: Comparison with a 3T System

PURPOSE

Hemodynamics is important in the initiation, growth, and rupture of intracranial aneurysms. Since intracranial aneurysms are small, a high-field MR system with high spatial resolution and high SNR is desirable for this hemodynamic analysis. The purpose of this study was to investigate whether the accuracy of MR fluid dynamic (MRFD) results based on 3D phase-contrast MR (3D PC MR, non-electrocardiogram[ECG]-gated 4D Flow MRI) data from a human cerebrovascular phantom and human healthy subjects obtained by a 7T MR system was superior to those by a 3T MR system.

METHODS

3D PC MR and 3D time of flight MR angiography (3D TOF MRA) imaging were performed on a 3T MR system and a 7T MR system for a human cerebrovascular phantom and 10 healthy human subjects, and MRFD analysis was performed using these data. The MRFD results from each MR system were then compared with the following items based on the computational fluid dynamics (CFD) results: 3D velocity vector field; correlation coefficient (R), angular similarity index (ASI), and magnitude similarity index (MSI) of blood flow velocity vectors.

RESULTS

In the MRFD results of 3D velocity vectors of the cerebrovascular phantom, noise-like vectors were observed near the vascular wall on the 3T MR system, but no noise was observed on the 7T MR system, showing results similar to those of CFD. In the MRFD results of the cerebrovascular phantom and healthy subjects, the correlation coefficients R, ASI, and MSI of the 7T MR system were higher than those of the 3T MR system, and ASI and MSI of healthy human subjects were significantly different between the two systems.

CONCLUSIONS

The accuracy of high spatial resolution MRFD using the 7T MR system exceeded that of the 3T MR system.

Sequential Hemodynamic Analysis of Ruptured Posterior Communicating Artery Aneurysms Treated With Coil Embolization and Delayed Flow Diversion

BACKGROUND AND OBJECTIVES

Computational fluid dynamics has advanced our knowledge of the pathogenesis of intracranial aneurysms and the dynamic changes observed after treatment. Herein, we analyze hemodynamic changes throughout the intervention stages for ruptured posterior communicating artery (PComA) aneurysms, treated with acute coiling and delayed flow diversion (FD).

METHODS

We performed a retrospective analysis of ruptured PComA aneurysms treated with the acute coiling and delayed FD strategy between June 2013 to November 2022, using 3-dimensional reconstructions of digital subtraction angiographies. Hemodynamic simulations using ANSYS® calculated aneurysmal and adjacent arteries' wall shear stress (WSS), aneurysmal low shear areas (LSA), and mean velocities in the distal arteries.

RESULTS

Six of the 22 patients were selected for computational fluid dynamics analysis, including 4 females and 2 males with a median age of 60 years. The average aneurysm volume was 984.12 mm3, with an average surface area of 386.11 mm2; LSA was 22.90%, and the average WSS was 3.39 Pa. The 2 largest aneurysms also had the highest LSA values. After coiling, there was a reduction in the aneurysmal volume (-78.42%) and the average surface area (-55.28%), and aneurysmal WSS increased to 6.10 Pa (+79.90%). WSS values for the middle cerebral artery (MCA) increased to 10.76 Pa, while anterior cerebral artery (ACA) increased to 7.51 Pa. Complete occlusion was achieved with delayed FD at a median follow-up of 19.7 months. After FD, average WSS increased to 14.94 Pa for the MCA (+70.64%) and to 10.82 Pa for the ACA (+30.10%). The mean MCA velocity increased to 43.04 cm/s (+36.85%), and 3 cases showed an increase in ACA velocities.

CONCLUSION

LSA may have triggered rupture for the PComA aneurysms analyzed. After coiling, average WSS increased in the aneurysm wall and downstream vessels in the majority of cases analyzed. Delayed FD caused hemodynamic disturbances distal to deployment, reflected in the sequential increase in the WSS and velocities in both the ACA and MCA.

Cerebral Aneurysm: Filling the Gap Between Pathophysiology and Nanocarriers

Intracranial aneurysms, characterized by abnormal dilations of cerebral arteries, pose significant health risks due to their potential to rupture, leading to subarachnoid hemorrhage with high mortality and morbidity rates. This paper aim is to explore the innovative application of nanoparticles in treating intracranial aneurysms, offering a promising avenue for enhancing current therapeutic strategies. We took into consideration the pathophysiology of cerebral aneurysms, focusing on the role of hemodynamic stress, endothelial dysfunction, and inflammation in their development and progression. By comparing cerebral aneurysms with other types, such as aortic aneurysms, we identify pathophysiological similarities and differences that could guide the adaptation of treatment approaches. The review highlights the potential of nanoparticles to improve drug delivery, targeting, and efficacy while minimizing side effects. We discuss various nanocarriers, including liposomes and polymeric nanoparticles, and their roles in overcoming biological barriers and enhancing therapeutic outcomes. Additionally, we discuss the potential of specific compounds, such as Edaravone and Tanshinone IIA, when used in conjunction with nanocarriers, to provide neuroprotective and anti-inflammatory benefits. By extrapolating insights from studies on aortic aneurysms, new research directions and therapeutic strategies for cerebral aneurysms are proposed. This interdisciplinary approach underscores the potential of nanoparticles to positively influence the management of intracranial aneurysms, paving the way for personalized treatment options that could significantly improve patient outcomes.

Geometrical Factors Affect Wall Shear Stress in Saccular Aneurysms of the Infrarenal Abdominal Aorta

BACKGROUND

Low wall shear stress (WSS) is predictive of aortic aneurysm growth and rupture. Yet, estimating WSS in a clinical setting is impractical, whereas measuring aneurysm geometry is feasible. This study investigates the association between saccular aneurysm geometry of the infrarenal aorta and WSS.

METHODS

Starting with a nonaneurysmal, patient-specific, computational fluid dynamics model of the aorta, saccular aneurysms of varying geometry were created by incrementally increasing the neck width and sac depth from 1 cm to 4 cm. The aspect ratio (the ratio between sac depth and neck width) varied between 0.25 and 4. The peak WSS, time-averaged WSS (TAWSS), and oscillatory shear index (OSI) were measured within the aneurysm sac.

RESULTS

Decreasing the neck width from 4 cm to 1 cm decreased the peak WSS by 69% and the TAWSS by 83%. Increasing the sac depth from 1 cm to 4 cm decreased the peak WSS by 55% and the OSI by 37%. The aspect ratio was negatively correlated to peak WSS (Rs -0.85; P < 0.001).

CONCLUSIONS

In saccular aneurysms of the infrarenal aorta, a smaller neck width, deeper aneurysm sac, and larger aspect ratio are associated with lower peak WSS.

Hemodynamic analysis of a thrombosed bleb in an unruptured cerebral aneurysm

Background

The intricate hemodynamic mechanisms of thrombosis in the bleb and/or dome of cerebral aneurysms remain unresolved. We encountered a unique case where the bleb of an unruptured internal carotid-posterior communicating artery (IC-PC) aneurysm underwent thrombosis over 7 years. Complete spontaneous thrombosis of a bleb in an unruptured cerebral aneurysm has not been previously reported. Therefore, a hemodynamic evaluation using computational fluid dynamics (CFD) analysis was conducted to examine the thrombotic development within the bleb of this aneurysm.

Case Description

We observed a case in which thrombosis led to the disappearance of a bleb in the dome of an unruptured IC-PC aneurysm over 7 years. CFD analysis was employed to investigate the hemodynamics of bleb thrombosis and the thrombosed bleb-neck regions of the dome in this IC-PC aneurysm. The reduction and disappearance of the bleb were associated with a decreased flow rate within the bleb, reduced magnitude of wall shear stress (WSSm), a lower WSSm ratio between the bleb and dome, increased vector direction of wall shear stress (WSSv), and discrete streamlines entering the bleb-neck region, resulting in stasis and subsequent thrombosis within the bleb. Seven years later, the dome region corresponding to the thrombosed bleb-neck exhibited localized areas with low WSSm and high WSSv along the dome wall.

Conclusion

Hemodynamically, spontaneously thrombosed bleb and thrombosed post-bleb-neck dome walls were characterized by low WSSm and high WSSv. These findings underscore the importance of CFD analysis in predicting thrombotic events in cerebral aneurysms, which can inform better clinical management strategies.

Computational fluid dynamics; a new diagnostic tool in giant intracerebral aneurysm treatment

Giant intracerebral aneurysms (GIA) comprise up to 5 % of all intracranial aneurysms. The indirect surgical strategy, which leaves the GIA untouched but reverses the blood flow by performing a bypass in combination with proximal parent artery occlusion is a useful method to achieve spontaneous aneurysm occlusion. The goal of this study was to assess the utility of computational fluid dynamics (CFD) in preoperative GIA treatment planning. We hypothesise that CFD simulations will predict treatment results. A fluid-structure interaction (FSI) CFD investigation was performed for the entire arterial brain circulation. The analyses were performed in three patient-specific CT angiogram models. The first served as the reference geometry with a C6 internal carotid artery (ICA) GIA, the second a proximal parent artery occlusion (PAO) and virtual bypass to the frontal M2 branch of the middle cerebral artery (MCA), and the third a proximal PAO in combination with a temporal M2 branch bypass. The volume of "old blood", flow residence time (FRT), dynamic viscosity and haemodynamic changes were also analysed. The "old blood" within the aneurysm in the bypass models reached 41 % after 20 cardiac cycles while in the reference model it was fully washed out. In Bypass 2 "old blood" was also observed in the main trunk of the MCA after 20 cardiac cycles. Extrapolation of the results yielded a duration of 4 years required to replace the "old blood" inside the aneurysm after bypass revascularization. In both bypass models a 7-fold increase in mean blood viscosity in the aneurysm region was noted. Bypass revascularization combined with proximal PAO favours thrombosis. Areas prone to thrombus formation, and subsequently the treatment outcomes, were accurately identified in the preoperative model. Virtual surgical operations can give a remarkable insight into haemodynamics that could support operative decision-making.

Fluid-Structure Interaction Simulations of the Initiation Process of Cerebral Aneurysms

Background: Hemodynamics during the growth process of cerebral aneurysms are incompletely understood. We developed a novel fluid-structure interaction analysis method for the identification of relevant scenarios of aneurysm onset. Method: This method integrates both fluid dynamics and structural mechanics, as well as their mutual interaction, for a comprehensive analysis. Patients with a single unruptured cerebral aneurysm were included. Results: Overall, three scenarios were identified. In scenario A, wall shear stress (WSS) was low, and the oscillatory shear index (OSI) was high in large areas within the region of aneurysm onset (RAO). In scenario B, the quantities indicated a reversed behavior, where WSS was high and OSI was low. In the last scenario C, a behavior in-between was found, with scenarios A and B coexisting simultaneously in the RAO. Structural mechanics demonstrated a similar but independent trend. Further, we analyzed the change in hemodynamics between the onset and a fully developed aneurysm. While scenarios A and C remained unchanged during aneurysm growth, 47% of aneurysms in scenario B changed into scenario A and 20% into scenario C. Conclusions: In conclusion, these findings suggest that WSS and the OSI are reciprocally regulated, and both low and high WSS/OSI conditions can lead to aneurysm onset.

Risk factors and predictive indicators of rupture in cerebral aneurysms

Cerebral aneurysms are abnormal dilations of blood vessels in the brain that have the potential to rupture, leading to subarachnoid hemorrhage and other serious complications. Early detection and prediction of aneurysm rupture are crucial for effective management and prevention of rupture-related morbidities and mortalities. This review aims to summarize the current knowledge on risk factors and predictive indicators of rupture in cerebral aneurysms. Morphological characteristics such as aneurysm size, shape, and location, as well as hemodynamic factors including blood flow patterns and wall shear stress, have been identified as important factors influencing aneurysm stability and rupture risk. In addition to these traditional factors, emerging evidence suggests that biological and genetic factors, such as inflammation, extracellular matrix remodeling, and genetic polymorphisms, may also play significant roles in aneurysm rupture. Furthermore, advancements in computational fluid dynamics and machine learning algorithms have enabled the development of novel predictive models for rupture risk assessment. However, challenges remain in accurately predicting aneurysm rupture, and further research is needed to validate these predictors and integrate them into clinical practice. By elucidating and identifying the various risk factors and predictive indicators associated with aneurysm rupture, we can enhance personalized risk assessment and optimize treatment strategies for patients with cerebral aneurysms.

Investigation of Relationship between Hemodynamic and Morphometric Characteristics of Aortas in Pediatric Patients

Background: The utilization of hemodynamic parameters, whose estimation is often cumbersome, can fasten diagnostics and decision-making related to congenital heart diseases. The main goal of this study is to investigate the relationship between hemodynamic and morphometric features of the thoracic aorta and to construct corresponding predictive models. Methods: Multi-slice spiral computed tomography images of the aortas of patients with coarctation diagnoses and patients without cardiac or vascular diseases were evaluated to obtain numerical models of the aorta and branches of the aortic arch. Hemodynamic characteristics were estimated in key subdomains of the aorta and three branches using computational fluid dynamics methods. The key morphometric features (diameters) were calculated at locations in proximity to the domains, where hemodynamic characteristics are evaluated. Results: The functional dependencies for velocities and pressure on the corresponding diameters have been fitted, and a metamodel has been constructed employing the predicted values from these models. Conclusions: The metamodel demonstrated high accuracy in classifying aortas into their respective types, thereby confirming the adequacy of the predicted hemodynamic characteristics by morphometric characteristics. The proposed methodology is applicable to other heart diseases without fundamental changes.

Description of the local hemodynamic environment in intracranial aneurysm wall subdivisions

Intracranial aneurysms (IAs) pose severe health risks influenced by hemodynamics. This study focuses on the intricate characterization of hemodynamic conditions within the IA walls and their influence on bleb development, aiming to enhance understanding of aneurysm stability and the risk of rupture. The methods emphasized utilizing a comprehensive dataset of 359 IAs and 213 IA blebs from 268 patients to reconstruct patient-specific vascular models, analyzing blood flow using finite element methods to solve the unsteady Navier-Stokes equations, the segmentation of aneurysm wall subregions and the hemodynamic metrics wall shear stress (WSS), its metrics, and the critical points in WSS fields were computed and analyzed across different aneurysm subregions defined by saccular, streamwise, and topographical divisions. The results revealed significant variations in these metrics, correlating distinct hemodynamic environments with wall features on the aneurysm walls, such as bleb formation. Critical findings indicated that regions with low WSS and high OSI, particularly in the body and central regions of aneurysms, are prone to conditions that promote bleb formation. Conversely, areas exposed to high WSS and positive divergence, like the aneurysm neck, inflow, and outflow regions, exhibited a different but substantial risk profile for bleb development, influenced by flow impingements and convergences. These insights highlight the complexity of aneurysm behavior, suggesting that both high and low-shear environments can contribute to aneurysm pathology through distinct mechanisms.

Measurement of Wall Shear Rate Across the Entire Vascular Wall Using Ultrasound Speckle Decorrelation

OBJECTIVE

The accurate measurement of the wall shear rate (WSR) plays a crucial role in the early diagnosis of cardiovascular disease progression and acute events such as aneurysms and atherosclerotic plaque ruptures. To address this need, the speckle decorrelation (SDC) technique has been used to measure WSR based on the 2-D out-of-plane blood flow speed. This technique is particularly advantageous because it enables the use of a 1-D array transducer to measure WSR over the entire luminal area. This study aims to develop a region-based singular value decomposition (SVD) filtering technique that selectively suppresses clutter noise in the vascular region to measure WSR using SDC.

METHOD

Ultrasound simulations, in-vitro flow experiments, and an in-vivo human study were conducted to evaluate the feasibility of this method's clinical application.

RESULTS

The results demonstrated that WSR can be effectively measured across entire vascular walls using a conventional 1-D array transducer along with the proposed methodology.

CONCLUSION

This study successfully demonstrates a noninvasive and accurate SDC-based method for measuring vital vascular WSR. This approach holds significant promise for assessing vascular WSR in both healthy individuals and high-risk cardiovascular disease patients.

A model with multiple intracranial aneurysms: possible hemodynamic mechanisms of aneurysmal initiation, rupture and recurrence

BACKGROUND

Hemodynamic factors play an important role in aneurysm initiation, growth, rupture, and recurrence, while the mechanism of the hemodynamic characteristics is still controversial. A unique model of multiple aneurysms (initiation, growth, rupture, and recurrence) is helpful to avoids the confounders and further explore the possible hemodynamic mechanisms of aneurysm in different states.

METHODS

We present a model with multiple aneurysms, and including the states of initiation, growth, rupture, and recurrence, discuss the proposed mechanisms, and describe computational fluid dynamic model that was used to evaluate the likely hemodynamic effect of different states of the aneurysms.

RESULTS

The hemodynamic analysis suggests that high flow impingement and high WSS distribution at normal parent artery was found before aneurysmal initiation. The WSS distribution and flow velocity were decreased in the new sac after aneurysmal growth. Low WSS was the risk hemodynamic factor for aneurysmal rupture. High flow concentration region on the neck plane after coil embolization still marked in recanalized aneurysm.

CONCLUSIONS

Associations have been identified between high flow impingement and aneurysm recanalization, while low WSS is linked to the rupture of aneurysms. High flow concentration and high WSS distribution at normal artery associated with aneurysm initiation and growth, while after growth, the high-risk hemodynamics of aneurysm rupture was occurred, which is low WSS at aneurysm dome.

Hemodynamic factors of spontaneous vertebral artery dissecting aneurysms assessed with numerical and deep learning algorithms: Role of blood pressure and asymmetry

BACKGROUND AND OBJECTIVES

The pathophysiology of spontaneous vertebral artery dissecting aneurysms (SVADA) is poorly understood. Our goal is to investigate the hemodynamic factors contributing to their formation using computational fluid dynamics (CFD) and deep learning algorithms.

METHODS

We have developed software that can use patient imagery as input to recreate the vertebrobasilar arterial system, both with and without SVADA, which we used in a series of three patients. To obtain the kinematic blood flow data before and after the aneurysm forms, we utilized numerical methods to solve the complex Navier-Stokes partial differential equations. This was accomplished through the application of a finite volume solver (OpenFoam/Helyx OS). Additionally, we trained a neural ordinary differential equation (NODE) to learn and replicate the dynamical streamlines obtained from the computational fluid dynamics (CFD) simulations.

RESULTS

In all three cases, we observed that the equilibrium of blood pressure distributions across the VAs, at a specific vertical level, accurately predicted the future SVADA location. In the two cases where there was a dominant VA, the dissection occurred on the dominant artery where blood pressure was lower compared to the contralateral side. The SVADA sac was characterized by reduced wall shear stress (WSS) and decreased velocity magnitude related to increased turbulence. The presence of a high WSS gradient at the boundary of the SVADA may explain its extension. Streamlines generated by CFD were learned with a neural ordinary differential equation (NODE) capable of capturing the system's dynamics to output meaningful predictions of the flow vector field upon aneurysm formation.

CONCLUSION

In our series, asymmetry in the vertebrobasilar blood pressure distributions at and proximal to the site of the future SVADA accurately predicted its location in all patients. Deep learning algorithms can be trained to model blood flow patterns within biological systems, offering an alternative to the computationally intensive CFD. This technology has the potential to find practical applications in clinical settings.

CFD analysis of the hyper-viscous effects on blood flow across abdominal aortic aneurysm in COVID patients: multiphysics approach

Recent research has shown that individuals suffering from COVID-19 are accommodating an elevated level of blood viscosity due to the morphological changes in blood cells. As viscosity is a major flow parameter influencing the flow across a stenosis or an aneurysm, the examination of the significance of hyperviscosity in COVID patients is imperative in arterial pathologies. In this research, we have considered a patient-specific case in which the aneurysm is located along the abdominal aortal walls. Recent research on the side effects of COVID-19 voiced out the various effects on the circulatory system of humans. Also, as abdominal aneurysms exist very often among individuals, causing the death of 150-200 million every year, the hyper-viscous effects of blood on the flow across the diseased aorta are explored by considering the elevated viscosity levels. In vitro explorations contribute considerably to the clinical methods and treatments to be regarded. The objective of the present inquest is to research the flow field in aneurysmatic-COVID-affected patients considering the elastic nature of vessel walls, using the arbitrary Lagrangian-Eulerian approach. The study supports the various clinical findings that voiced the detrimental effects associated with blood hyperviscosity. The simulation results obtained, by solving the fluid mechanics' equations coupled with the solid mechanics' equations, employing a FEM solver suggest that the elevated stress imparted by the hyper-viscous flows on the walls of the aneurysmal aorta can trigger the fastening of the aneurysmal sac enlargement or rupture.

Evaluation of predictive models of aneurysm focal growth and bleb development using machine learning techniques

BACKGROUND

The presence of blebs increases the rupture risk of intracranial aneurysms (IAs).

OBJECTIVE

To evaluate whether cross-sectional bleb formation models can identify aneurysms with focalized enlargement in longitudinal series.

METHODS

Hemodynamic, geometric, and anatomical variables derived from computational fluid dynamics models of 2265 IAs from a cross-sectional dataset were used to train machine learning (ML) models for bleb development. ML algorithms, including logistic regression, random forest, bagging method, support vector machine, and K-nearest neighbors, were validated using an independent cross-sectional dataset of 266 IAs. The models' ability to identify aneurysms with focalized enlargement was evaluated using a separate longitudinal dataset of 174 IAs. Model performance was quantified by the area under the receiving operating characteristic curve (AUC), the sensitivity and specificity, positive predictive value, negative predictive value, F1 score, balanced accuracy, and misclassification error.

RESULTS

The final model, with three hemodynamic and four geometrical variables, along with aneurysm location and morphology, identified strong inflow jets, non-uniform wall shear stress with high peaks, larger sizes, and elongated shapes as indicators of a higher risk of focal growth over time. The logistic regression model demonstrated the best performance on the longitudinal series, achieving an AUC of 0.9, sensitivity of 85%, specificity of 75%, balanced accuracy of 80%, and a misclassification error of 21%.

CONCLUSIONS

Models trained with cross-sectional data can identify aneurysms prone to future focalized growth with good accuracy. These models could potentially be used as early indicators of future risk in clinical practice.

Longitudinal wall shear stress evaluation using centerline projection approach in the numerical simulations of the patient-based carotid artery

In this numerical study, areas of the carotid bifurcation and of a distal stenosis in the internal carotid artery are closely observed to evaluate the patient's current risks of ischemic stroke. An indicator for the vessel wall defects is the stress exerted by blood on the vessel tissue, typically expressed by the amplitude of the wall shear stress vector (WSS) and its oscillatory shear index. To detect negative shear stresses corresponding with reversal flow, we perform orientation-based shear evaluation. We investigate the longitudinal component of the wall shear vector, where tangential vectors aligned longitudinally with the vessel are necessary. However, resulting from imaging segmentation resolution of patients' computed tomography angiography scans and stenotic regions, the geometry model's mesh is non-smooth on its surface areas and the automatically generated tangential vector field is discontinuous and multi-directional, making an interpretation of our orientation-based risk indicators unreliable. We improve the evaluation of longitudinal shear stress by applying the projection of the vessel's centerline to the surface to construct smooth tangential field aligned longitudinally with the vessel. We validate our approach for the longitudinal WSS component and the corresponding oscillatory index by comparing them to results obtained using automatically generated tangents in both rigid and elastic vessel modeling and to amplitude-based indicators. We present the major benefit of our longitudinal WSS evaluation based on its directionality for the cardiovascular risk assessment, which is the detection of negative WSS indicating persistent reversal or transverse flow. This is impossible in the case of the amplitude-based WSS.

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.

Computational Fluid Dynamic Simulations of Cerebral Aneurysms

Computational fluid dynamics (CFD) simulations have been introduced to enable individualized risk prognosis for patients with unruptured cerebral aneurysms. The present contribution provides an overview of the biomechanical and physiological principles of aneurysm formation and rupture. It describes the computational steps of the CFD and the evaluated parameters. The clinical value of CFD is then discussed based on a recent literature review. Finally, we discuss current methodological limitations and possible future developments to overcome the actual drawbacks of CFD.

Assessing the risk of intracranial aneurysm rupture using computational fluid dynamics: a pilot study

Objective

This study compared 2 representative cases with ruptured aneurysms to explore the role of hemodynamic and morphological parameters in evaluating the rupture risk of intracranial aneurysms (IAs).

Methods

CTA and 3-dimensional rotational angiography (3DRA) of 3 IAs in 2 patients were retrospectively analyzed in this study. Hemodynamics and morphological parameters were compared between a ruptured IA and an unruptured IA in case1, and between before and after aneurysm rupture in case 2.

Results

In case 1, the ruptured aneurysm had larger morphological parameters including size ratio (SR), aspect ratio (AR), aneurysm vessel angle (θF), Aneurysm inclination angle (θA), Undulation index (UI), Ellipticity index (EI), and Non-sphericity Index (NSI) than the unruptured aneurysm. And oscillatory shear index (OSI) is also larger. Higher rupture resemblance score (RRS) was shown in the ruptured aneurysm. In case 2, the aneurysm had one daughter sac after 2 years. Partial morphological and hemodynamic parameters including SR, AR, θF, θA, UI, EI, NSI, OSI, and relative residence time (RRT) increased, and normalized wall shear stress (NWSS) was significantly reduced. RRS increased during this period.

Conclusion

SR and OSI may have predictive values for the risk of intracranial aneurysm rupture. It is possible that WSS Changes before and after IA rupture, yet the influence of high or low WSS on growth and rupture of IA remains unclear. RRS is promising to be used in the clinical assessment of the rupture risk of IAs and to guide the formulation of treatment plans.

A matched-pair case control study identifying hemodynamic predictors of cerebral aneurysm growth using computational fluid dynamics

Introduction: Initiation and progression of cerebral aneurysms is known to be driven by complex interactions between biological and hemodynamic factors, but the hemodynamic mechanism which drives aneurysm growth is unclear. We employed robust modeling and computational methods, including temporal and spatial convergence studies, to study hemodynamic characteristics of cerebral aneurysms and identify differences in these characteristics between growing and stable aneurysms. Methods: Eleven pairs of growing and non-growing cerebral aneurysms, matched in both size and location, were modeled from MRA and CTA images, then simulated using computational fluid dynamics (CFD). Key hemodynamic characteristics, including wall shear stress (WSS), oscillatory shear index (OSI), and portion of the aneurysm under low shear, were evaluated. Statistical analysis was then performed using paired Wilcoxon rank sum tests. Results: The portion of the aneurysm dome under 70% of the parent artery mean wall shear stress was higher in growing aneurysms than in stable aneurysms and had the highest significance among the tested metrics (p = 0.08). Other metrics of area under low shear had similar levels of significance. Discussion: These results align with previously observed hemodynamic trends in cerebral aneurysms, indicating a promising direction for future study of low shear area and aneurysm growth. We also found that mesh resolution significantly affected simulated WSS in cerebral aneurysms. This establishes that robust computational modeling methods are necessary for high fidelity results. Together, this work demonstrates that complex hemodynamics are at play within cerebral aneurysms, and robust modeling and simulation methods are needed to further study this topic.

Reinforcement learning for patient-specific optimal stenting of intracranial aneurysms

Developing new capabilities to predict the risk of intracranial aneurysm rupture and to improve treatment outcomes in the follow-up of endovascular repair is of tremendous medical and societal interest, both to support decision-making and assessment of treatment options by medical doctors, and to improve the life quality and expectancy of patients. This study aims at identifying and characterizing novel flow-deviator stent devices through a high-fidelity computational framework that combines state-of-the-art numerical methods to accurately describe the mechanical exchanges between the blood flow, the aneurysm, and the flow-deviator and deep reinforcement learning algorithms to identify a new stent concepts enabling patient-specific treatment via accurate adjustment of the functional parameters in the implanted state.

Near-wall hemodynamic parameters quantification in in vitro intracranial aneurysms with 7 T PC-MRI

OBJECTIVE

Wall shear stress (WSS) and its derived spatiotemporal parameters have proven to play a major role on intracranial aneurysms (IAs) growth and rupture. This study aims to demonstrate how ultra-high field (UHF) 7 T phase contrast magnetic resonance imaging (PC-MRI) coupled with advanced image acceleration techniques allows a highly resolved visualization of near-wall hemodynamic parameters patterns in in vitro IAs, paving the way for more robust risk assessment of their growth and rupture.

MATERIALS AND METHODS

We performed pulsatile flow measurements inside three in vitro models of patient-specific IAs using 7 T PC-MRI. To this end, we built an MRI-compatible test bench, which faithfully reproduced a typical physiological intracranial flow rate in the models.

RESULTS

The ultra-high field 7 T images revealed WSS patterns with high spatiotemporal resolution. Interestingly, the high oscillatory shear index values were found in the core of low WSS vortical structures and in flow stream intersecting regions. In contrast, maxima of WSS occurred around the impinging jet sites.

CONCLUSIONS

We showed that the elevated signal-to-noise ratio arising from 7 T PC-MRI enabled to resolve high and low WSS patterns with a high degree of detail.

Aneurysm wall enhancement, hemodynamics, and morphology of intracranial fusiform aneurysms

Background and objective

Intracranial fusiform aneurysms (IFAs) are considered to have a complex pathophysiology process and poor natural history. The purpose of this study was to investigate the pathophysiological mechanisms of IFAs based on the characteristics of aneurysm wall enhancement (AWE), hemodynamics, and morphology.

Methods

A total of 21 patients with 21 IFAs (seven fusiform types, seven dolichoectatic types, and seven transitional types) were included in this study. Morphological parameters of IFAs were measured from the vascular model, including the maximum diameter (Dmax), maximum length (Lmax), and centerline curvature and torsion of fusiform aneurysms. The three-dimensional (3D) distribution of AWE in IFAs was obtained based on high-resolution magnetic resonance imaging (HR-MRI). Hemodynamic parameters including time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), gradient oscillatory number (GON), and relative residence time (RRT) were extracted by computational fluid dynamics (CFD) analysis of the vascular model, and the relationship between these parameters and AWE was investigated.

Results

The results showed that Dmax (p = 0.007), Lmax (p = 0.022), enhancement area (p = 0.002), and proportion of enhancement area (p = 0.006) were significantly different among three IFA types, and the transitional type had the largest Dmax, Lmax, and enhancement area. Compared with the non-enhanced regions of IFAs, the enhanced regions had lower TAWSS but higher OSI, GON, and RRT (p &lt; 0.001). Furthermore, Spearman’s correlation analysis showed that AWE was negatively correlated with TAWSS, but positively correlated with OSI, GON, and RRT.

Conclusion

There were significant differences in AWE distributions and morphological features among the three IFA types. Additionally, AWE was positively associated with the aneurysm size, OSI, GON, and RRT, while negatively correlated with TAWSS. However, the underlying pathological mechanism of the three fusiform aneurysm types needs to be further studied.

Aortic wall shear stress in bicuspid aortic valve disease-10-year follow-up

Background

Bicuspid aortic valve (BAV) disease leads to deviant helical flow patterns especially in the mid-ascending aorta (AAo), potentially causing wall alterations such as aortic dilation and dissection. Among others, wall shear stress (WSS) could contribute to the prediction of long-term outcome of patients with BAV. 4D flow in cardiovascular magnetic resonance (CMR) has been established as a valid method for flow visualization and WSS estimation. The aim of this study is to reevaluate flow patterns and WSS in patients with BAV 10 years after the initial evaluation.

Methods

Fifteen patients (median age 34.0 years) with BAV were re-evaluated 10 years after the initial study from 2008/2009 using 4D flow by CMR. Our particular patient cohort met the same inclusion criteria as in 2008/2009, all without enlargement of the aorta or valvular impairment at that time. Flow patterns, aortic diameters, WSS and distensibility were calculated in different aortic regions of interest (ROI) with dedicated software tools.

Results

Indexed aortic diameters in the descending aorta (DAo), but especially in the AAo did not change in the 10-year period. Median difference 0.05 cm/m² (95% CI: 0.01 to 0.22; P=0.06) for AAo and median difference -0.08 cm/m² (95% CI: -0.12 to 0.01; P=0.07) for DAo. WSS values were lower in 2018/2019 at all measured levels. Aortic distensibility decreased by median 25.6% in the AAo, while stiffness increased concordantly (median +23.6%).

Conclusions

After a ten years' follow-up of patients with isolated BAV disease, indexed aortic diameters did not change in this patient cohort. WSS was lower compared to values generated 10 years earlier. Possibly a drop of WSS in BAV could serve as a marker for a benign long-term course and implementation of more conservative treatment strategies.

Computational Fluid Dynamics Analysis of Blood Flow Changes during the Growth of Saccular Abdominal Aortic Aneurysm

Computational fluid dynamics analysis of the growth process of saccular abdominal aortic aneurysm was performed. A 3D model of aortic aneurysm was created based on CT images. Properties in terms of wall shear stress, mean flow velocity, mean pressure, energy loss, and pressure loss coefficient were calculated using thermal fluid analysis software "ANSYS CFX." As the aneurysm expanded, the mean flow velocity decreased and the wall shear stress, mean pressure, energy loss, and pressure loss coefficient increased. Wall shear stress increased when the aneurysm was small, suggesting that is related to the development and growth of the aneurysm. (This is secondary publication from J Jpn Coll Angiol 2021; 61: 3-10.).

Wall Shear Stress Estimation for 4D Flow MRI Using Navier-Stokes Equation Correction

This study introduces a novel wall shear stress (WSS) estimation method for 4D flow MRI. The method improves the WSS accuracy by using the reconstructed pressure gradient and the flow-physics constraints to correct velocity gradient estimation. The method was tested on synthetic 4D flow data of analytical Womersley flow and flow in cerebral aneurysms and applied to in vivo 4D flow data acquired in cerebral aneurysms and aortas. The proposed method's performance was compared to the state-of-the-art method based on smooth-spline fitting of velocity profile and the WSS calculated from uncorrected velocity gradient. The proposed method improved the WSS accuracy by as much as 100% for the Womersley flow and reduced the underestimation of mean WSS by 39 to 50% for the synthetic aneurysmal flow. The predicted mean WSS from the in vivo aneurysmal data using the proposed method was 31 to 50% higher than the other methods. The predicted aortic WSS using the proposed method was 3 to 6 times higher than the other methods and was consistent with previous CFD studies and the results from recently developed methods that take into account the limited spatial resolution of 4D flow MRI. The proposed method improves the accuracy of WSS estimation from 4D flow MRI, which can help predict blood vessel remodeling and progression of cardiovascular diseases.

Subarachnoid Hemorrhage from Ruptured Internal Carotid Artery Aneurysm: Association with Arterial Tortuosity

OBJECTIVE

Many researchers have found a correlation between tortuous arteries and development of aneurysms in cerebral arteries. We decided to determine whether tortuosity of the internal carotid artery can be related to its aneurysm rupture.

METHODS

We retrospectively analyzed the internal carotid artery anatomy of 149 patients with internal carotid artery aneurysms. For each patient, we calculated relative length (RL), sum of angle metrics (SOAM), triangular index (TI), product of angle distance (PAD), and inflection count metrics (ICM).

RESULTS

A total of 33 patients (22.15%) had subarachnoid hemorrhage. These patients had significantly lower SOAM (0.31 ± 0.17 vs. 0.42 ± 0.21; P < 0.01), TI (0.27 ± 0.09 vs. 0.31 ± 0.11; P = 0.03) and ICM (0.25 ± 0.11 vs. 0.31 ± 0.17; P = 0.04). In multivariate logistic regression analysis, higher SOAM (odds ratio, 0.780; 95% confidence interval, 0.619-0.961; P = 0.025) remained independently associated with lower risk of internal carotid artery aneurysm rupture. In addition, we found significant positive correlation of aneurysm dome size with SOAM (R = 0.224; P = 0.013) and PAD (0.269; P < 0.01). Our study also showed that age (R = 0.252; P = 0.036), Glasgow Coma Scale score (R = -0.706; P < 0.01), and TI (R = -0.249; P = 0.042) were independently correlated with modified Rankin Scale score on discharge.

CONCLUSIONS

Lower tortuosity might be a protective factor against internal carotid artery aneurysm rupture and poor outcome after subarachnoid hemorrhage. Higher tortuosity is correlated with internal carotid artery aneurysm growth.

Challenges in Modeling Hemodynamics in Cerebral Aneurysms Related to Arteriovenous Malformations

PURPOSE

The significantly higher incidence of aneurysms in patients with arteriovenous malformations (AVMs) suggests a strong hemodynamic relationship between these lesions. The presence of an AVM alters hemodynamics in proximal vessels by drastically changing the distal resistance, thus affecting intra-aneurysmal flow. This study discusses the challenges associated with patient-specific modeling of aneurysms in the presence of AVMs.

METHODS

We explore how the presence of a generic distal AVM affects upstream aneurysms by examining the relationship between distal resistance and aneurysmal wall shear stress using physiologically realistic estimates for the influence of the AVM on hemodynamics. Using image-based computational models of aneurysms and surrounding vasculature, aneurysmal wall-shear stress is calculated for a range of distal resistances corresponding to the presence of AVMs of various sizes and compared with a control case representing the absence of an AVM.

RESULTS

In the patient cases considered, the alteration in aneurysmal wall shear stress due to the presence of an AVM is considerable, as much as 19 times the base case wall shear stress. Furthermore, the relationship between aneurysmal wall shear stress and distal resistance is shown to be highly geometry-dependent and nonlinear. In most cases, the range of physiologically realistic possibilities for AVM-related distal resistance are so large that patient-specific flow measurements are necessary for meaningful predictions of wall shear stress.

CONCLUSIONS

The presented work offers insight on the impact of distal AVMs on aneurysmal wall shear stress using physiologically realistic computational models. Patient-specific modeling of hemodynamics in aneurysms and associated AVMs has great potential for understanding lesion pathogenesis, surgical planning, and assessing the effect of treatment of one lesion relative to another. However, we show that modeling approaches cannot usually meaningfully quantify the impact of AVMs if based solely on imaging data from CT and X-ray angiography, currently used in clinical practice. Based on recent studies, it appears that 4D flow MRI is one promising approach to obtaining meaningful patient-specific flow boundary conditions that improve modeling fidelity.

4D Flow MRI in the portal venous system: imaging and analysis methods, and clinical applications

Thus far, ultrasound, CT, and 2D cine phase-contrast MRI has been adopted to evaluate blood flow and vascular morphology in the portal venous system; however, all these techniques have some shortcomings, such as limited field of view and difficulty in accurately evaluating blood flow. A new imaging technique, namely 3D cine phase-contrast (4D Flow) MRI, can acquire blood flow data of the entire abdomen at once and in a time-resolved manner, allowing visual, quantitative, and comprehensive assessment of blood flow in the portal venous system. In addition, a retrospective blood flow analysis, i.e., "retrospective flowmetry," is possible. Although the development of 4D Flow MRI for the portal system has been delayed compared to that for the arterial system owing to the lower flow velocity of the portal venous system and the presence of respiratory artifacts, several useful reports have recently been published as the technology has advanced. In the first part of this narrative review article, technical considerations of image acquisition and analysis methods of 4D Flow MRI for the portal venous system and the validations of their results are described. In the second part, the current clinical application of 4D Flow MRI for the portal venous system is reviewed.

Minimum wall shear stress points and their underlying intra-aneurysmal flow structures of unruptured cerebral aneurysms on 4D flow MRI

BACKGROUND AND PURPOSE

Minimum wall shear stress (Min-WSS) points may be associated with wall instability of unruptured cerebral aneurysms. We aimed to investigate the relationship between the locations of Min-WSS points and their underlying intra-aneurysmal flow structure patterns in unruptured cerebral aneurysms using four-dimensional (4D) flow magnetic resonance imaging (MRI).

MATERIALS AND METHODS

Min-WSS points and the intra-aneurysmal flow structure patterns were identified in 50 unruptured aneurysms by 4D flow MRI.

RESULTS

The Min-WSS points were located around a vortex core tip in 31 (62.0%) aneurysms and on an intra-bleb vortex center in 7 (14.0%). Sixteen (32.0%) aneurysms had the Min-WSS points on the aneurysmal apex, and in 24 (48.0%) were on the neck. The Min-WSS values of aneurysms with the Min-WSS points on an intra-bleb flow were significantly lower than those of the other groups (P = 0.030). Aneurysms with the Min-WSS points on the neck had significantly higher Min-WSS values than the other aneurysms (P = 0.008).

CONCLUSIONS

The location of the Min-WSS point was corresponding to the vortex core or center in 76% of all aneurysms. The underlying intra-aneurysmal flow structure and location of the Min-WSS point affect the Min-WSS value. Further studies are needed to characterize Min-WSS points to identify aneurysms with a higher risk of wall instability.

CFD to Quantify Idealized Intra-Aneurysmal Blood Flow in Response to Regular and Flow Diverter Stent Treatment

Cerebral aneurysms are pathological dilatations of the vessels supplying the brain. They carry a certain risk of rupture, which in turn, results in a high risk of mortality and morbidity. Flow diverters (FDs) are high-density meshed stents which are implanted in the vessel segment harboring an intracranial aneurysm to cover the entrance of the aneurysm, thus reducing the blood flow into the aneurysm, promoting thrombosis formation and stable occlusion, which prevents rupture or growth of the aneurysm. In the present study, the blood flow in an idealized aneurysm, treated with an FD stent and a regular stent (RS), were modeled and analyzed considering their design, surface area porosity, and flow reduction to investigate the quantitative and qualitative effect of the stent on intra-aneurysmal hemodynamics. CFD simulations were conducted before and after treatment. Significant reductions were observed for most hemodynamic variables with the use of stents, during both the peak systolic and late diastolic cardiac cycles. FD reduces the intra-aneurysmal wall shear stress (WSS), inflow, and aneurysmal flow velocity, and increases the turnover time when compared to the RS; therefore, the possibility of aneurysm thrombotic occlusion is likely to increase, reducing the risk of rupture in cerebral aneurysms.

The Cerebral Arterial Wall in the Development and Growth of Intracranial Aneurysms

A considerable number of people harbor intracranial aneurysms (IA), which is a focal or segmental disease of the arterial wall. The pathophysiologic mechanisms of IAs formation, growth, and rupture are complex. The mechanism also differs with respect to the type of aneurysm. In broad aspects, aneurysms may be considered a disease of the vessel wall. In addition to the classic risk factors and the genetic/environmental conditions, altered structural and pathologic events along with the interaction of the surrounding environment and luminal flow dynamics contribute to the aneurysm’s development and growth. In this review, we have tried to simplify the complex interaction of a multitude of events in relation to vessel wall in the formation and growth of IAs.

Research on the effect of visceral artery Aneurysm's cardiac morphological variation on hemodynamic situation based on time-resolved CT-scan and computational fluid dynamics

BACKGROUND AND OBJECTIVE

Muscular arteries and related aneurysms keep deforming during the cardiac cycle. However, current patient-specific computational fluid dynamics (CFD) analyses of aneurysms are usually based on individual cardiac phase images. The cardiac deformation and displacement characteristics of muscle arteries and aneurysms, as well as their impact on CFD results, have not been adequately explored. The present study tried to illustrate the cardiac morphological variation of visceral muscular arteries (VMAs) & aneurysms (VAAs) and evaluate its influence on the hemodynamic situation at lesion locations.

METHODS

Four-dimensional computed tomography angiogram (4D-CTA) images of six patients with VAAs were acquired. Medical image registration is used to capture cardiac variations of VMAs. The steady-state CFD simulation is performed on twelve different time-phase geometries. Deformation, displacement, wall shear stress (WSS), velocity, and pressure values at pathological locations are compared to illustrate the deforming characteristics of VAAs and their influence on CFD simulation results.

RESULTS

The deformation and displacement characteristics of lesion locations for six specific patients show a pulsatile pattern. Maximum displacements are always less than 4 mm. The ratio fluctuations of endovascular cavity volume and vascular inner wall surface area, which were employed to depict cardiac deformation, are always less than 20%. According to CFD simulations based on deformed VMAs, WSS has a larger coefficient of variation (COV) than velocity and pressure. Except for one patient's WSS, the COVs of different hemodynamic parameters obtained from simulation results are always less than 10%.

CONCLUSIONS

Based on 4D-CTA images, we confirmed that cardiovascular circulation has a periodic impact on the morphologic characteristics of VMAs. A wave that has extended throughout the studied region is observed. It has a dominant influence on the displacement of VMAs. According to CFD results, the influence of the VMAs' deformation and displacement on different hemodynamic parameters is distinct. The variance in WSS is more prominent compared to pressure and velocity. On most occasions, the influence of the VMAs' periodic deformation and displacement on simulation results is insignificant. However, the variant simulation results induced by deforming VMAs cannot be simply ignored.

Hemodynamic Characteristic Analysis of Aneurysm Wall Enhancement in Unruptured Middle Cerebral Artery Aneurysm

Background and Purpose

Aneurysm wall enhancement (AWE) on vessel wall magnetic resonance imaging has been suggested as a marker of the unstable status of intracranial aneurysm (IA) and may predict IA rupture risk. However, the role of abnormal hemodynamics in unruptured IAs with AWE remains poorly understood. This study aimed to determine the association between abnormal hemodynamics and AWE in unruptured middle cerebral artery (MCA) aneurysms.

Methods

A total of 28 patients with 32 bifurcation aneurysms of the middle cerebral artery>3mm in size were retrospectively selected for this study. Vessel wall magnetic resonance images were reviewed, and the AWE pattern of each aneurysm was classified as no AWE, partial AWE, and circumferential AWE. Computational fluid dynamics were used to calculate the hemodynamic variables of each aneurysm. Univariate and multivariate analyses investigated the association between AWE and hemodynamic variables.

Results

AWE was present in 13 aneurysms (40.6%), with 7 (21.9%) showing partial AWE and 6 (18.7%) showing circumferential AWE. Kruskal-Wallis H analysis revealed that hemodynamic variables including wall shear stress (WSS), oscillatory shear index, aneurysm pressure (AP), relative residence time, and low shear area (LSA) were significantly associated with AWE (p < 0.05). Further ordinal logistic regression analysis found that WSS was the only factor with a significant association with AWE (p = 0.048); similar trends were identified for LSA (p = 0.055) and AP (p = 0.058). Spearman's correlation analysis showed that AWE was negatively correlated with WSS (rs = -0.622, p < 0.001) and AP (rs = -0.535, p = 0.002) but positively correlated with LSA (rs = 0.774, p < 0.001).

Conclusion

Low wall shear stress, low aneurysm pressure, and increased low shear area were associated with aneurysm wall enhancement on vessel wall magnetic resonance imaging in unruptured cerebral aneurysms. These abnormal hemodynamic parameters may induce inflammation and cause aneurysm wall enhancement. However, the association between these parameters and their underlying pathological mechanisms requires further investigation.

Aneurysm-on-a-Chip: Setting Flow Parameters for Microfluidic Endothelial Cultures Based on Computational Fluid Dynamics Modeling of Intracranial Aneurysms

Intracranial aneurysms are pouch-like extrusions from the vessels at the base of the brain which can rupture and cause a subarachnoid hemorrhage. The pathophysiological mechanism of aneurysm formation is thought to be a consequence of blood flow (hemodynamic) induced changes on the endothelium. In this study, the results of a personalized aneurysm-on-a-chip model using patient-specific flow parameters and patient-specific cells are presented. CT imaging was used to calculate CFD parameters using an immersed boundary method. A microfluidic device either cultured with human umbilical vein endothelial cells (HUVECs) or human induced pluripotent stem cell-derived endothelial cells (hiPSC-EC) was used. Both types of endothelial cells were exposed for 24 h to either 0.03 Pa or 1.5 Pa shear stress, corresponding to regions of low shear and high shear in the computational aneurysm model, respectively. As a control, both cell types were also cultured under static conditions for 24 h as a control. Both HUVEC and hiPSC-EC cultures presented as confluent monolayers with no particular cell alignment in static or low shear conditions. Under high shear conditions HUVEC elongated and aligned in the direction of the flow. HiPSC-EC exhibited reduced cell numbers, monolayer gap formation and cells with aberrant, spread-out morphology. Future research should focus on hiPSC-EC stabilization to allow personalized intracranial aneurysm models.

Numerical Assessment of the Risk of Abnormal Endothelialization for Diverter Devices: Clinical Data Driven Numerical Study

Numerical modeling is an effective tool for preoperative planning. The present work is devoted to a retrospective analysis of neurosurgical treatments for the occlusion of cerebral aneurysms using flow-diverters and hemodynamic factors affecting stent endothelization. Several different geometric approaches have been considered for virtual flow-diverters deployment. A comparative analysis of hemodynamic parameters as a result of computational modeling has been carried out basing on the four clinical cases: one successful treatment, one with no occlusion and two with in stent stenosis. For the first time, a quantitative assessment of both: the limiting magnitude of shear stresses that are necessary for the occurrence of in stent stenosis (MaxWSS > 1.23) and for conditions in which endothelialization is insufficiently active and occlusion of the cervical part of the aneurysm does not occur (MaxWSS < 1.68)—has been statistacally proven (p < 0.01).

EFFECTS OF THE ANEURYSM TREATMENT METHOD ON COIL PACKING DENSITY AND ITS RELATIONSHIP WITH THE OSTIUM AREA

PURPOSE

To test the following hypotheses: (1) balloon or stent assistance increases coil packing density (CPD) in the endovascular treatment of intracranial aneurysms and (2) CPD is related to the ostium area and volume of the aneurysm.

MATERIALS AND METHODS

This retrospective study included 60 ruptured and unruptured (54 and 6, respectively) aneurysms treated between August 2017 and December 2019 in the authors' institution with simple coiling (n=18), balloon-assisted coiling (n=7), or stent-assisted coiling (n=35). Aneurysm volume and ostium area measurements were obtained using commercial software from three-dimensional digital subtraction angiography images. Coil sizes were retrieved from patient files, and coil volume measurements were obtained using https://www.angiocalc.com/. Analysis of covariance, multivariable covariance analysis, and Pearson's correlation analyses were performed.

RESULTS

The median values for the aneurysm volume, coil volume, CPD, and ostium area were 63.4 (range 5.5-1771.4) mm³, 23.13 (2.03-296.95) mm³, 33.29% (13.41%-81.02%), and 10.7 (2.7-49.9) mm², respectively. Multivariate analysis showed that the CPD values were not significantly different among the treatment groups, although the ostium area differed significantly between the simple and stent-assisted coiling groups (P<.05). Pearson's correlations showed that similar to aneurysm volume, the ostium area was negatively correlated with CPD (r=-0.321, P<.05).

CONCLUSION

In cerebral aneurysms treated with balloon- or stent-assisted coiling, the CPD value does not differ from that in aneurysms treated with simple coiling.

Artificial Intelligence Applications in Intracranial Aneurysm: Achievements, Challenges and Opportunities

Intracranial aneurysms present in about 3% of the general population and the number of detected aneurysms is continuously rising with the advances in imaging techniques. Intracranial aneurysm rupture carries a high risk of death or permanent disabilities; therefore assessment of the intracranial aneurysm along the entire course is of great clinical importance. Given the outstanding performance of artificial intelligence (AI) in image-based tasks, many AI-based applications have emerged in recent years for the assessment of intracranial aneurysms. In this review we will summarize the state-of-the-art of AI applications in intracranial aneurysms, emphasizing the achievements, and exploring the challenges. We will also discuss the future prospects and potential opportunities. This article provides an updated view of the AI applications in intracranial aneurysms and may act as a basis for guiding the related future works.

4D flow MRI hemodynamic biomarkers for cerebrovascular diseases

Alterations in cerebral blood flow are common in several neurological diseases among the elderly including stroke, cerebral small vessel disease, vascular dementia, and Alzheimer's disease. 4D flow magnetic resonance imaging (MRI) is a relatively new technique to investigate cerebrovascular disease, and makes it possible to obtain time-resolved blood flow measurements of the entire cerebral arterial venous vasculature and can be used to derive a repertoire of hemodynamic biomarkers indicative of cerebrovascular health. The information that can be obtained from one single 4D flow MRI scan allows both the investigation of aberrant flow patterns at a focal location in the vasculature as well as estimations of brain-wide disturbances in blood flow. Such focal and global hemodynamic biomarkers show the potential of being sensitive to impending cerebrovascular disease and disease progression and can also become useful during planning and follow-up of interventions aiming to restore a normal cerebral circulation. Here, we describe 4D flow MRI approaches for analyzing the cerebral vasculature. We then survey key hemodynamic biomarkers that can be reliably assessed using the technique. Finally, we highlight cerebrovascular diseases where one or multiple hemodynamic biomarkers are of central interest.