Patient-specific arterial wall generation for intracranial aneurysms with a variable and a near realistic vessel wall thickness for FSI studies

BACKGROUND AND OBJECTIVE

Imaging methodologies such as, computed tomography (CT) aid in three-dimensional (3D) reconstruction of patient-specific aneurysms. The radiological data is useful in understanding their location, shape, size, and disease progression. However, there are serious impediments in discerning the blood vessel wall thickness due to limitations in the current imaging modalities. This further restricts the ability to perform high-fidelity fluid structure interaction (FSI) studies for an accurate assessment of rupture risk. FSI studies would require the arterial wall mesh to be generated to determine realistic maximum allowable wall stresses by performing coupled calculations for the hemodynamic forces with the arterial walls.

METHODS

In the present study, a novel methodology is developed to geometrically model variable vessel wall thickness for the lumen isosurface extracted from CT scan slices of patient-specific aneurysms based on clinical and histopathological inputs. FSI simulations are carried out with the reconstructed models to assess the importance of near realistic wall thickness model on rupture risk predictions.

RESULTS

During surgery, clinicians often observe translucent vessel walls, indicating the presence of thin regions. The need to generate variable vessel wall thickness model, that embodies the wall thickness gradation, is closer to such clinical observations. Hence, corresponding FSI simulations performed can improve clinical outcomes. Considerable differences in the magnitude of instantaneous wall shear stresses and von Mises stresses in the walls of the aneurysm was observed between a uniform wall thickness and a variable wall thickness model.

CONCLUSION

In the present study, a variable vessel wall thickness generation algorithm is implemented. It was shown that, a realistic wall thickness modeling is necessary for an accurate prediction of the shear stresses on the wall as well as von Mises stresses in the wall. FSI simulations are performed to demonstrate the utility of variable wall thickness modeling.

Morphometric and Hemodynamic Analysis of the Compressed Iliac Vein

PURPOSE

To investigate the relationship between the morphological structure and hemodynamic properties of the compressed iliac vein and explore the reason for the formation of thrombosis in the compressed iliac vein.

MATERIALS AND METHODS

A total of 11 patients with iliac vein compression syndrome (IVCS) were included in this study, and their iliac veins were reconstructed in 3 dimensions (3D). The morphological structures of the iliac veins (confluence angle, degree of stenosis) were analyzed based on the 3D model. Variations in the hemodynamic properties of the iliac vein were investigated at 4 typical moments in one cardiac cycle, and the relationship between the different morphological configurations and the pressure difference was investigated.

RESULTS

In the region of the compressed iliac vein, the blood flow velocity is accelerated and the pressure changes abruptly accompanied by the increase in pressure difference. Higher time averaged wall shear stress (TAWSS) and lower relative residence time (RRT) appeared in stenosis regions of compressed iliac vein, while TAWSS was low and RRT was large near the stenosis position. There was a strong positive correlation between the degree of stenosis and the pressure difference (r=0.894), and a positive correlation between the confluence angle of the iliac vein and the pressure difference (r=0.638).

CONCLUSION

The morphological structure of the compressed iliac vein has an obvious influence on the hemodynamic surroundings; the pressure difference becomes larger when the degree of stenosis and the confluence angle increase. The iliac vein luminal areas with low TAWSS and high RRT near the compressed location can impede blood flow and lead to accumulation of blood components, which may increase the risk of thrombosis formation and should be fully considered in the treatment of IVCS.

Impact of left atrial wall motion assumptions in fluid simulations on proposed predictors of thrombus formation

Atrial fibrillation (AF) poses a significant risk of stroke due to thrombus formation, which primarily occurs in the left atrial appendage (LAA). Medical image-based computational fluid dynamics (CFD) simulations can provide valuable insight into patient-specific hemodynamics and could potentially enhance personalized assessment of thrombus risk. However, the importance of accurately representing the left atrial (LA) wall dynamics has not been fully resolved. In this study, we compared four modeling scenarios; rigid walls, a generic wall motion based on a reference motion, a semi-generic wall motion based on patient-specific motion, and patient-specific wall motion based on medical images. We considered a LA geometry acquired from 4D computed tomography during AF, systematically performed convergence tests to assess the numerical accuracy of our solution strategy, and quantified the differences between the four approaches. The results revealed that wall motion had no discernible impact on LA cavity hemodynamics, nor on the markers that indicate thrombus formation. However, the flow patterns within the LAA deviated significantly in the rigid model, indicating that the assumption of rigid walls may lead to errors in the estimated risk factors. In contrast, the generic, semi-generic, and patient-specific cases were qualitatively similar. The results highlight the crucial role of wall motion on hemodynamics and predictors of thrombus formation, and also demonstrate the potential of using a generic motion model as a surrogate for the more complex patient-specific motion. While the present study considered a single case, the employed CFD framework is entirely open-source and designed for adaptability, allowing for integration of additional models and generic motions.

Endovascular Treatment of Intracranial Aneurysm: The Importance of the Rheological Model in Blood Flow Simulations

Hemodynamics in intracranial aneurysm strongly depends on the non-Newtonian blood behavior due to the large number of suspended cells and the ability of red blood cells to deform and aggregate. However, most numerical investigations on intracranial hemodynamics adopt the Newtonian hypothesis to model blood flow and predict aneurysm occlusion. The aim of this study was to analyze the effect of the blood rheological model on the hemodynamics of intracranial aneurysms in the presence or absence of endovascular treatment. A numerical investigation was performed under pulsatile flow conditions in a patient-specific aneurysm with and without the insertion of an appropriately reconstructed flow diverter stent (FDS). The numerical simulations were performed using Newtonian and non-Newtonian assumptions for blood rheology. In all cases, FDS placement reduced the intra-aneurysmal velocity and increased the relative residence time (RRT) on the aneurysmal wall, indicating progressive thrombus formation and aneurysm occlusion. However, the Newtonian model largely overestimated RRT values and consequent aneurysm healing with respect to the non-Newtonian models. Due to the non-Newtonian blood properties and the large discrepancy between Newtonian and non-Newtonian simulations, the Newtonian hypothesis should not be used in the study of the hemodynamics of intracranial aneurysm, especially in the presence of endovascular treatment.

Differences and Correlations of Morphological and Hemodynamic Parameters between Anterior Circulation Bifurcation and Side-wall Aneurysms

OBJECTIVE

The objective of this research was to explore the difference and correlation of the morphological and hemodynamic features between sidewall and bifurcation aneurysms in anterior circulation arteries, utilizing computational fluid dynamics as a tool for analysis.

METHODS

In line with the designated inclusion criteria, this study covered 160 aneurysms identified in 131 patients who received treatment at Union Hospital of Tongji Medical College, Huazhong University of Science and Technology, China, from January 2021 to September 2022. Utilizing follow-up digital subtraction angiography (DSA) data, these cases were classified into two distinct groups: the sidewall aneurysm group and the bifurcation aneurysm group. Morphological and hemodynamic parameters in the immediate preoperative period were meticulously calculated and examined in both groups using a three-dimensional DSA reconstruction model.

RESULTS

No significant differences were found in the morphological or hemodynamic parameters of bifurcation aneurysms at varied locations within the anterior circulation. However, pronounced differences were identified between sidewall and bifurcation aneurysms in terms of morphological parameters such as the diameter of the parent vessel (Dvessel), inflow angle (θF), and size ratio (SR), as well as the hemodynamic parameter of inflow concentration index (ICI) (P<0.001). Notably, only the SR exhibited a significant correlation with multiple hemodynamic parameters (P<0.001), while the ICI was closely related to several morphological parameters (R>0.5, P<0.001).

CONCLUSIONS

The significant differences in certain morphological and hemodynamic parameters between sidewall and bifurcation aneurysms emphasize the importance to contemplate variances in threshold values for these parameters when evaluating the risk of rupture in anterior circulation aneurysms. Whether it is a bifurcation or sidewall aneurysm, these disparities should be considered. The morphological parameter SR has the potential to be a valuable clinical tool for promptly distinguishing the distinct rupture risks associated with sidewall and bifurcation aneurysms.

Hemodynamic study on the therapeutic effects of varying diameter embolic coils in the treatment of intracranial aneurysms

Endovascular coiling is the predominant method for treating cerebral aneurysms. Extensive reports on selecting coil length, hardness, and material are available. However, the impact of coil diameter on postoperative outcomes remains unclear. This study enrolled six personalized geometric models of intracranial aneurysms: three bifurcation aneurysms and three sidewall aneurysms. Four coil models were constructed by changing the coil diameter. Coil embolization was simulated using the finite element method. Computational fluid dynamics was used to characterize hemodynamics in the aneurysms after embolization. Evaluation parameters included velocity reduction, wall shear stress (WSS), low WSS (LWSS), oscillatory shear index (OSI), relative residence time (RRT), and residual flow volume in the aneurysms. At the peak time (t = 0.17 s), the proportion of LWSS area in bifurcation aneurysms increase with the rise in coil diameter: 0.8D, 71.28 ± 12.62% versus 1D, 74.97 ± 19.17% versus 1.2D, 78.88 ± 18.56% versus 1.4D, 84.00 ± 11.53% (mean ± SD). The proportion of high OSI area decreases as the coil diameter increases: 0.8D, 4.41% ± 2.82% versus 1.0D, 3.78 ± 3.33% versus 1.2D, 2.28% ± 1.77% versus 1.4D, 1.58% ± 1.11% (mean ± SD). The proportion of high RRT area increases as the coil diameter rises: 0.8D, 3.40% ± 1.68% versus 1.0D, 7.67 ± 4.12% versus 1.2D, 9.84% ± 9.50% versus 1.4D, 22.29% ± 14.28% (mean ± SD). Side wall aneurysms do not exhibit the aforementioned trend. Bifurcation aneurysms plugged with a coil of 1.4 times the diameter have the largest RFVs (<10 mm/s) within the group. Aforementioned patterns are not found in sidewall aneurysms. In the treatment of aneurysms with coiling, varying coil diameters can result in different hemodynamic environments within the aneurysm. Larger coil diameters have improved hemodynamic performance for bifurcation aneurysms. However, coil diameter and embolization effectiveness have no significant relationship for sidewall aneurysms.

Plasma levels of CD36 and glutathione as biomarkers for ruptured intracranial aneurysm

Evidence has proved that intracranial aneurysm (IA) formation and rupture might be closely related to inflammatory response and oxidative stress. Our objective was to evaluate the potential of CD36 and glutathione (GSH) as biomarkers for IA. In this study, the enzyme-linked immunosorbent assay was used to measure the plasma levels of CD36 and GSH in 30 IA patients and 30 healthy controls. Then, correlation analysis, receiver operating characteristic (ROC) curve, and logistic regression analysis were performed. The results showed that the plasma level of CD36 in IA patients was significantly higher than that in the control group (P < 0.0001), and plasma GSH was significantly lower compared with that in the control group (P < 0.0001). ROC analysis showed that CD36 and GSH had high sensitivity (90.0 and 96.6%) and specificity (96.6 and 86.6%) for IA diagnosis. The combined sensitivity and specificity achieved were 100 and 100%, respectively. The plasma levels of CD36 and GSH did not show a significant correlation with age, the Glasgow Coma Scale, Hunter-Hess score, aneurysm size, aneurysm height, aneurysm neck, and aspect ratio. The AUC of the logistic regression model based on CD36 and GSH was 0.505. Our results suggested that the combination of plasma CD36 and GSH could serve as potential biomarkers for IA rupture.

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.

Identification of miRNAs Involved in Intracranial Aneurysm Rupture in Cigarette-Smoking Patients

INTRODUCTION

Smoking is an independent risk factor for the formation and rupture of intracranial aneurysms (IA). However, the underlying mechanism remains unclear.

METHODS

In this study, we performed miRNA sequencing on plasma from 10 smoking patients with IA, 10 non-smoking patients with IA, and 10 healthy controls. The differentially expressed miRNAs (DE miRNAs) between smoking and non-smoking patients with IA were identified. Functional and pathway enrichment analysis is employed to investigate the potential functions of those DE miRNA target genes. The correlations with the clinical parameters were assessed using receiver operating characteristic curve (ROC) analysis.

RESULTS

In total, we identified 428 DE miRNAs. Functional enrichment analysis showed the target genes were significantly enriched in biological aspects related to cell characteristics, such as cell cycle, cell differentiation, and cell migration. Pathway analysis showed DE miRNAs mainly enriched in the PI3K-Akt signaling pathway, Focal adhesion, and JAK-STAT signaling pathway. The expressions of miR-574-5p, miR-151a-3p, and miR-652-3p correlated well with aneurysm parameters. The AUC of miR-574-5p, miR-151a-3p, and miR-652-3p were 97%, 92%, and 99%, respectively.

CONCLUSION

Our study indicated that smoking significantly altered the plasma miRNA profile in patients with IA. The expression of miR-574-5p, miR-151a-3p, and miR-652-3p correlated with aneurysm parameters, which may play a significant role in the formation and rupture of IA.

Predicting Morphological Changes to Vessel Walls Adjacent to Unruptured Cerebral Aneurysms Using Computational Fluid Dynamics

Objective  This study compared intraoperative findings with preoperative computed tomography angiography (CTA) and computational fluid dynamics (CFD) analysis of perianeurysmal findings for the indication of possible vessel wall thinning. Materials and Methods  Participants comprised 38 patients with unruptured middle cerebral artery aneurysms treated by surgical clipping at our hospital between May 2020 and April 2021. We defined parent artery radiation sign (PARS) as the presence of each of the following three findings in CFD analysis based on preoperative CTA: (1) impingement of the stream line on the outer parent vessel wall of the aneurysm; (2) radiation of wall shear stress vectors outwards from the same site; and (3) increased wall pressure compared with the surrounding area. CFD analysis showing PARS was compared with intraoperative findings. Results  In all nine cases with PARS, no morphological abnormalities were found in the same area on CTA. However, intraoperative findings showed thinning of the parent artery wall in one of the nine cases and formation of a very small mass in three cases, differing from CTA findings. All nine patients underwent additional clipping and/or wrapping and coating at the site of PARS. Conclusion  Detecting thinning of the vessel wall or the presence of a microaneurysm may be difficult in endovascular therapy, which is based on the visualization of the vessel lumen. CFD analysis suggests the necessity of confirming findings for the vessel wall around an aneurysm by direct manipulation, as the presence of PARS may indicate partial thinning of the vessel wall or formation of a microaneurysm.

Haemodynamic Effects on the Development and Stability of Atherosclerotic Plaques in Arterial Blood Vessel

Studying the formation and stability of atherosclerotic plaques in the hemodynamic field is essential for understanding the growth mechanism and preventive treatment of atherosclerotic plaques. In this paper, based on a multiplayer porous wall model, we established a two-way fluid-solid interaction with time-varying inlet flow. The lipid-rich necrotic core (LRNC) and stress in atherosclerotic plaque were described for analyzing the stability of atherosclerotic plaques during the plaque growth by solving advection-diffusion-reaction equations with finite-element method. It was found that LRNC appeared when the lipid levels of apoptotic materials (such as macrophages, foam cells) in the plaque reached a specified lower concentration, and increased with the plaque growth. LRNC was positively correlated with the blood pressure and was negatively correlated with the blood flow velocity. The maximum stress was mainly located at the necrotic core and gradually moved toward the left shoulder of the plaque with the plaque growth, which increases the plaque instability and the risk of the plaque shedding. The computational model may contribute to understanding the mechanisms of early atherosclerotic plaque growth and the risk of instability in the plaque growth.

Numerical investigation of unruptured middle cerebral artery bifurcation aneurysms: influence of aspect ratio

An aneurysm is a disease condition, which is due to the pathological weakening of an arterial wall. These aneurysms are often found in various branch points and bifurcations of an artery in the cerebral circulation. Most aneurysms come to medical attention, either due to brain hemorrhages caused by rupture or found unruptured. To consider surgically invasive treatment modalities, clinicians need scientific methods such as, hemodynamic analysis to assess rupture risk. The arterial wall loses its structural integrity when wall shear stress (WSS) and other hemodynamic parameters exceed a certain threshold. In the present study, numerical simulations are carried out for unruptured middle cerebral artery (MCA) aneurysms. Three distinct representative sizes are chosen from a larger patient pool of 26 MCA aneurysms. Logically, these aneurysms represent three growth stages of any patient with similar anatomical structure. Simulations are performed to compare the three growth phases (with different aspect ratios) of an aneurysm and correlate their hemodynamic parameters. Simulations with patient specific boundary conditions reveal that, aneurysms with a higher aspect ratio (AR) correspond to an attendant decrease in both time-averaged wall shear stress (TAWSS) and spatial wall shear stress gradients (WSSG). Smaller MCAs were observed to have higher positive wall shear stress divergence (WSSD), exemplifying the tensile nature of arterial wall stretching. Present study identifies positive wall shear stress divergence (PWSSD) to be a potential biomarker for evaluating the growth of an aneurysm.

Simulating Subject-Specific Aortic Hemodynamic Effects of Valvular Lesions in Rheumatic Heart Disease

Rheumatic heart disease (RHD) is a neglected tropical disease despite the substantial global health burden. In this study, we aimed to develop a lower cost method of modeling aortic blood flow using subject-specific velocity profiles, aiding our understanding of RHD's consequences on the structure and function of the ascending aorta. Echocardiography and cardiovascular magnetic resonance (CMR) are often used for diagnosis, including valve dysfunction assessments. However, there is a need to further characterize aortic valve lesions to improve treatment options and timing for patients, while using accessible and affordable imaging strategies. Here, we simulated effects of RHD aortic valve lesions on the aorta using computational fluid dynamics (CFD). We hypothesized that inlet velocity distribution and wall shear stress (WSS) will differ between RHD and non-RHD individuals, as well as between subject-specific and standard Womersley velocity profiles. Phase-contrast CMR data from South Africa of six RHD subjects with aortic stenosis and/or regurgitation and six matched controls were used to estimate subject-specific velocity inlet profiles and the mean velocity for Womersley profiles. Our findings were twofold. First, we found WSS in subject-specific RHD was significantly higher (p < 0.05) than control subject simulations, while Womersley simulation groups did not differ. Second, evaluating spatial velocity differences (ΔSV) between simulation types revealed that simulations of RHD had significantly higher ΔSV than non-RHD (p < 0.05), these results highlight the need for implementing subject-specific input into RHD CFD, which we demonstrate how to accomplish through accessible methods.

Trends in the Age of Patients Treated for Unruptured Intracranial Aneurysms from 1990 to 2020

BACKGROUND

The decision for treatment for unruptured intracranial aneurysms (UIAs) is often difficult. Innovation in endovascular devices have improved the benefit-to-risk profile especially for elderly patients; however, the treatment guidelines from the past decade often recommend conservative management. It is unknown how these changes have affected the overall age of the patients selected for treatment. Herein, we aimed to study potential changes in the average age of the patients that are being treated over time.

METHODS

A systematic search of the literature was performed to identify all studies describing the age of the UIAs that were treated by any modality. Scatter diagrams with trend lines were used to plot the age of the patients treated over time and assess the presence of a potential significant trend via statistical correlation tests.

RESULTS

A total of 280 studies including 83,437 UIAs treated between 1987 and 2021 met all eligibility criteria and were entered in the analysis. Mean age of the patients was 55.5 years, and 70.7% were female. There was a significant increasing trend in the age of the treated patients over time (Spearman r: 0.250; P < 0.001), with a 1-year increase in the average age of the treated patients every 5 years since 1987.

CONCLUSIONS

The present study indicates that based on the treated UIA patient data published in the literature, older UIAs are being treated over time. This trend is likely driven by safer treatments while suggesting that re-evaluation of certain UIA treatment decision scores may be of great interest.

Automatic 4D Flow MRI Segmentation Using the Standardized Difference of Means Velocity

We present a method to automatically segment 4D flow magnetic resonance imaging (MRI) by identifying net flow effects using the standardized difference of means (SDM) velocity. The SDM velocity quantifies the ratio between the net flow and observed flow pulsatility in each voxel. Vessel segmentation is performed using an F-test, identifying voxels with significantly higher SDM velocity values than background voxels. We compare the SDM segmentation algorithm against pseudo-complex difference (PCD) intensity segmentation of 4D flow measurements in in vitro cerebral aneurysm models and 10 in vitro Circle of Willis (CoW) datasets. We also compared the SDM algorithm to convolutional neural network (CNN) segmentation in 5 thoracic vasculature datasets. The in vitro flow phantom geometry is known, while the ground truth geometries for the CoW and thoracic aortas are derived from high-resolution time-of-flight (TOF) magnetic resonance angiography and manual segmentation, respectively. The SDM algorithm demonstrates greater robustness than PCD and CNN approaches and can be applied to 4D flow data from other vascular territories. The SDM to PCD comparison demonstrated an approximate 48% increase in sensitivity in vitro and 70% increase in the CoW, respectively; the SDM and CNN sensitivities were similar. The vessel surface derived from the SDM method was 46% closer to the in vitro surfaces and 72% closer to the in vitro TOF surfaces than the PCD approach. The SDM and CNN approaches both accurately identify vessel surfaces. The SDM algorithm is a repeatable segmentation method, enabling reliable computation of hemodynamic metrics associated with cardiovascular disease.

Abnormal flow pattern of low wall shear stress and high oscillatory shear index in spontaneous vertebral artery dissection with vertebral artery hypoplasia

Introduction

Spontaneous vertebral artery dissection (sVAD) might tend to develop in vertebral artery hypoplasia (VAH) with hemodynamic dysfunction and it is crucial to assess hemodynamics in sVAD with VAH to investigate this hypothesis. This retrospective study aimed to quantify hemodynamic parameters in patients with sVAD with VAH.

Methods

Patients who had suffered ischemic stroke due to an sVAD of VAH were enrolled in this retrospective study. The geometries of 14 patients (28 vessels) were reconstructed using Mimics and Geomagic Studio software from CT angiography (CTA). ANSYS ICEM and ANSYS FLUENT were utilized for mesh generation, set boundary conditions, solve governing equations, and perform numerical simulations. Slices were obtained at the upstream area, dissection or midstream area and downstream area of each VA. The blood flow patterns were visualized through instantaneous streamline and pressure at peak systole and late diastole. The hemodynamic parameters included pressure, velocity, time-averaged blood flow, time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), endothelial cell action potential (ECAP), relative residence time (RRT) and time-averaged nitric oxide production rate (TAR_NO).

Results

Significant focal increased velocity was present in the dissection area of steno-occlusive sVAD with VAH compared to other nondissected areas (0.910 m/s vs. 0.449 vs. 0.566, p < 0.001), while focal slow flow velocity was observed in the dissection area of aneurysmal dilatative sVAD with VAH according to velocity streamlines. Steno-occlusive sVAD with VAH arteries had a lower time-averaged blood flow (0.499 cm³/s vs. 2.268, p < 0.001), lower TAWSS (1.115 Pa vs. 2.437, p = 0.001), higher OSI (0.248 vs. 0.173, p = 0.006), higher ECAP (0.328 Pa^-1 vs. 0.094, p = 0.002), higher RRT (3.519 Pa^-1 vs. 1.044, p = 0.001) and deceased TAR_NO (104.014 nM/s vs. 158.195, p < 0.001) than the contralateral VAs.

Conclusion

Steno-occlusive sVAD with VAH patients had abnormal blood flow patterns of focal increased velocity, low time-averaged blood flow, low TAWSS, high OSI, high ECAP, high RRT and decreased TAR_NO. These results provide a good basis for further investigation of sVAD hemodynamics and support the applicability of the CFD method in testing the hemodynamic hypothesis of sVAD. More detailed hemodynamic conditions with different stages of sVAD are warranted in the future.

Hemodynamics of anterior circulation intracranial aneurysms with daughter blebs: investigating the multidirectionality of blood flow fields

Recent advances in diagnostic neuroradiological imaging, allowed the detection of unruptured intracranial aneurysms (IAs). The shape - irregular or multilobular - of the aneurysmal dome, is considered as a possible rupture risk factor, independently of the size, the location and patient medical background. Disturbed blood flow fields in particular is thought to play a key role in IAs progression. However, there is an absence of widely-used hemodynamic indices to quantify the extent of a multi-directional disturbed flow. We simulated blood flow in twelve patient-specific anterior circulation unruptured intracranial aneurysms with daughter blebs utilizing the spectral/hp element framework Nektar++. We simulated three cardiac cycles using a volumetric flow rate waveform while we considered blood as a Newtonian fluid. To investigate the multidirectionality of the blood flow fields, besides the time-averaged wall shear stress (TAWSS), we calculated the oscillatory shear index (OSI), the relative residence time (RRT) and the time-averaged cross flow index (TACFI). Our CFD simulations suggest that in the majority of our vascular models there is a formation of complex intrasaccular flow patterns, resulting to low and highly oscillating WSS, especially in the area of the daughter blebs. The existence of disturbed multi-directional blood flow fields is also evident by the distributions of the RRT and the TACFI. These findings further support the theory that IAs with daughter blebs are linked to a potentially increased rupture risk.

Hemodynamics and bio-mechanics of morphologically distinct saccular intracranial aneurysms at bifurcations: Idealised vs Patient-specific geometries

Understanding the factors that influence the rupture of aneurysms is of primary concern to the clinicians, who are grappled with patient management. It is important to know how the relation between morphological features of the cerebral aneurysm, and the mechanical stresses on the containing arterial walls are influenced by the hemodynamic forces. Present study investigates three different shapes, which have been identified correspondingly in patient-specific scenarios as well. The primary objective is to categorize the bifurcation aneurysms into standard shapes such as, spherical, beehive and pear-shaped, based on patient-specific clinical studies and further compare and contrast the model aneurysms with the patient specific configurations, for their hemodynamic factors as well as the attendant stresses on the wall. MethodsComputational fluid dynamic simulations are performed accounting for the fluid-structure interaction (FSI) effects between the flowing fluid and the containing vessel wall. Blood is assumed to be Newtonian, while the arterial walls are assumed to be linearly elastic. A commercial solver is used for performing detailed calculations. Hemodynamic and bio-mechanical rupture predictions are carried out for the three different shapes. Observations derived from the idealised simulations are compared and contrasted against their patient-specific counterparts. ResultsFrom detailed numerical simulations, it was observed that pear-shaped aneurysms exhibit large re-circulation bubble and flow stagnation zone, with higher residence time for the particles, which may lead to atherosclerotic lesions. Beehive shape allows for maximum flow into the aneurysmal sac with concentrated jet impinging on the dome, leading to high values of maximum WSS (MWSS) resulting in great propensity to form a secondary bleb. However, flow field inside a spherical aneurysm is found to be stable with fewer vortices, and nearly uniform distribution of wall stresses are observed though-out the sac, which perhaps signifies hemodynamically and bio-mechanically stable condition. ConclusionCategorizing patient-specific intracranial aneurysms into standard shapes viz, spherical, beehive and pear could generalize the process of prediction of hemodynamic and bio-mechanical rupture indicators. Comparative assessment of the flow field and stresses reported from the simulations on idealised models, with corresponding patient-specific simulations reveal that, these studies could aid in understanding the generalised shape dependence of hemodynamic and bio-mechanical behaviour of aneurysms.

Hemodynamic Analysis Shows High Wall Shear Stress Is Associated with Intraoperatively Observed Thin Wall Regions of Intracranial Aneurysms

BACKGROUND

Studying the relationship between hemodynamics and local intracranial aneurysm (IA) pathobiology can help us understand the natural history of IA. We characterized the relationship between the IA wall appearance, using intraoperative imaging, and the hemodynamics from CFD simulations.

METHODS

Three-dimensional geometries of 15 IAs were constructed and used for CFD. Two-dimensional intraoperative images were subjected to wall classification using a machine learning approach, after which the wall type was mapped onto the 3D surface. IA wall regions included thick (white), normal (purple-crimson), and thin/translucent (red) regions. IA-wide and local statistical analyses were performed to assess the relationship between hemodynamics and wall type.

RESULTS

Thin regions of the IA sac had significantly higher WSS, Normalized WSS, WSS Divergence and Transverse WSS, compared to both normal and thick regions. Thicker regions tended to co-locate with significantly higher RRT than thin regions. These trends were observed on a local scale as well. Regression analysis showed a significant positive correlation between WSS and thin regions and a significant negative correlation between WSSD and thick regions.

CONCLUSION

Hemodynamic simulation results were associated with the intraoperatively observed IA wall type. We consistently found that elevated WSS and WSS_Norm were associated with thin regions of the IA wall rather than thick and normal regions.

Tortuosity of parent artery predicts in-stent stenosis after pipeline flow-diverter stenting for internal carotid artery aneurysms

Background and purpose

The relationship between the tortuosity of the parent artery and treatment outcomes is not well established. We investigate the association between parent artery tortuosity and flow diverter (FD) treatment outcomes in patients with internal carotid artery aneurysms in this study.

Methods

A retrospective review study was conducted to identify all patients with internal carotid artery aneurysms who were implanted with Pipeline embolization device (PED) between 2016 and 2020. The relationship between parent artery tortuosity and aneurysm complete occlusion (CO) and in-stent stenosis (ISS) was analyzed. The mathematical parameters “Curvature”, “torsion”, and “DM” extracted from the parent artery were utilized to quantify the parent artery tortuosity. A vascular narrowing of greater than 25% was categorized as ISS. Logistic regression analysis was used to identify significant independent predictors. Furthermore, we compared the performance of four machine learning algorithms and Logistic Regression model in predicting ISS.

Results

This research included 62 patients who with internal carotid artery aneurysms. In 49 (79%) cases, follow-up angiography (mean follow-up duration 11.7 ±7.3 months) revealed CO of the aneurysm. ISS was detected in 22 (35.5%) cases. According to univariate analysis, parent artery tortuosity and other variables were not associated with CO (p > 0.1). Maximum curvature (OR = 1.084; 95% CI, 1.008–1.165; p = 0.03) and DM (OR = 0.01; 95% CI, 0–0.488; p = 0.02) exhibited strong independent associations with ISS in multivariate analysis. The SVM model is superior to the conventional Logistic Regression model and the other models in predicting ISS.

Conclusions

The tortuosity of the parent artery may affect the treatment outcome of FD stenting. We found that parent artery tortuosity was associated with ISS, but not with aneurysm complete occlusion following PED stenting for internal carotid artery aneurysms in this study. Parent arteries with higher maximum curvature and lower DM were more likely to develop ISS.

Fluid-structure interaction (FSI) simulation for studying the impact of atherosclerosis on hemodynamics, arterial tissue remodeling, and initiation risk of intracranial aneurysms

The biomechanical and hemodynamic effects of atherosclerosis on the initiation of intracranial aneurysms (IA) are not yet clearly discovered. Also, studies for the observation of hemodynamic variation due to atherosclerotic stenosis and its impact on arterial remodeling and aneurysm genesis remain a controversial field of vascular engineering. The majority of studies performed are relevant to computational fluid dynamic (CFD) simulations. CFD studies are limited in consideration of blood and arterial tissue interactions. In this work, the interaction of the blood and vessel tissue because of atherosclerotic occlusions is studied by developing a fluid and structure interaction (FSI) analysis for the first time. The FSI presents a semi-realistic simulation environment to observe how the blood and vessels' structural interactions can increase the accuracy of the biomechanical study results. In the first step, many different intracranial vessels are modeled for an investigation of the biomechanical and hemodynamic effects of atherosclerosis in arterial tissue remodeling. Three physiological conditions of an intact artery, the artery with intracranial atherosclerosis (ICAS), and an atherosclerotic aneurysm (ACA) are employed in the models with required assumptions. Finally, the obtained outputs are studied with comparative and statistical analyses according to the intact model in a normal physiological condition. The results show that existing occlusions in the cross-sectional area of the arteries play a determinative role in changing the hemodynamic behavior of the arterial segments. The undesirable variations in blood velocity and pressure throughout the vessels increase the risk of arterial tissue remodeling and aneurysm formation.

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.

Increased Concentrations of Atherogenic Proteins in Aneurysm Sac Are Associated with Wall Enhancement of Unruptured Intracranial Aneurysm

Current MR-vessel wall imaging (VWI) of unruptured intracranial aneurysms (UIAs) permits the visualization of wall structures. Aneurysm wall enhancement (AWE) was associated with atherosclerotic remodeling of the aneurysm wall accompanied by infiltration of inflammatory cells, potentially contributing to rupture. This study sought to investigate whether the luminal concentrations of atherosclerotic proteins in the aneurysm sac were associated with increased wall enhancement of UIAs in VWI. Subjects undergoing endovascular treatment for UIAs were prospectively recruited. All subjects underwent evaluation using 3 T-MRI including pre/post contrast VWI of the UIAs. Blood samples were collected from the aneurysm sac and the parent artery during endovascular procedures. The presence of AWE was correlated with the delta difference in concentration between the aneurysm sac and the parent artery for each atherosclerotic protein. A total of consecutive 45 patients with 50 UIAs were enrolled. The delta differences of anti-oxidized low-density lipoprotein (LDL) antibody, small dense LDL, and lipoprotein(a) [Lp(a)] were significantly higher in UIAs with AWE compared with those without AWE (767.6 ± 1957.1 versus - 442.4 ± 1676.3 mIU/mL, p = 0.02, 114.8 ± 397.7 versus - 518.5 ± 1344.4 μg/mL, p = 0.04, and - 5.6 ± 11.3 versus - 28.7 ± 38.5 μg/mL, p = 0.01, respectively). In multivariate logistic regression analysis, the delta Lp(a) was significantly associated with AWE (p = 0.04). Increased concentrations of atherogenic proteins in the aneurysm sac were significantly associated with wall enhancement of UIAs. Future studies examining the effect of medications for atherosclerosis on the atherogenic proteins within the aneurysm sac and hence the wall enhancement are warranted.

Quantitative assessment of microstructural evolution of intracranial aneurysm wall by vessel wall imaging

PURPOSE

This study aimed to evaluate new quantitative parameters of aneurysm wall enhancement (AWE) on magnetic resonance vessel wall imaging (VWI) in differentiating between the stable and evolving unruptured intracranial aneurysms (UIAs).

METHODS

Thirty-eight consecutive patients with UIAs (27 stable and 11 evolving) underwent VWI with contrast-enhanced 3D T1 volume isotropic turbo spin echo acquisition. The voxel-based enhancement maps were created using pre- and post-contrast images. The aneurysmal lumen with signal suppression by black-blood method was segmented. Then, one voxel outer and inner layers of the lumen contour were automatically segmented. The shape features of the aneurysms and AWE of the two layers were compared between stable and evolving groups.

RESULTS

The shape features, including aneurysm volume, surface, and compacity were significantly different between the stable and evolving groups (P = 0.024, 0.028, and 0.033, respectively). Stable and evolving groups also differed significantly in the AWE at the union of outer and inner layers of the aneurysm wall (P = 0.0082) but not in that of the outer or inner layer alone. Multivariate logistic regression analysis revealed significant differences in aneurysm volume, surface, and AWE at the union of outer and inner layers between the two groups (P = 0.0029, 0.0092, and 0.0033, respectively). Receiver operating characteristics curve analysis revealed that the area under the curve of the logistic regression model was 0.89.

CONCLUSION

Quantitative combined analysis of aneurysm shape features and AWE of the union of outer and inner layers were effective for differentiating between stable and evolving UIAs.

Geometry and Flow Properties Affect the Phase Shift between Pressure and Shear Stress Waves in Blood Vessels

The phase shift between pressure and wall shear stress (WSS) has been associated with vascular diseases such as atherosclerosis and aneurysms. The present study aims to understand the effects of geometry and flow properties on the phase shift under the stiff wall assumption, using an immersed-boundary-lattice-Boltzmann method. For pulsatile flow in a straight pipe, the phase shift is known to increase with the Womersley number, but is independent of the flow speed (or the Reynolds number). For a complex geometry, such as a curved pipe, however, we find that the phase shift develops a strong dependence on the geometry and Reynolds number. We observed that the phase shift at the inner bend of the curved vessel and in the aneurysm dome is larger than that in a straight pipe. Moreover, the geometry affects the connection between the phase shift and other WSS-related metrics, such as time-averaged WSS (TAWSS). For straight and curved blood vessels, the phase shift behaves qualitatively similarly to and can thus be represented by the TAWSS, which is a widely used hemodynamic index. However, these observables significantly differ in other geometries, such as in aneurysms. In such cases, one needs to consider the phase shift as an independent quantity that may carry additional valuable information compared to well-established metrics.

Prediction of atherosclerotic changes in cavernous carotid aneurysms based on computational fluid dynamics analysis: a proof-of-concept study

PURPOSE

Recent computational fluid dynamics (CFD) studies have demonstrated the concurrence of atherosclerotic changes in regions exposed to prolonged blood residence. In this proof-of-concept study, we investigated a small but homogeneous cohort of large, cavernous carotid aneurysms (CCAs) to establish the clinical feasibility of CFD analysis in treatment planning, based on the association between pathophysiology and hemodynamics.

METHODS

This study included 15 patients with individual large CCAs. We identified calcifications, which indicated atherosclerotic changes, using the masking data of digital subtraction angiography. We conducted a CFD simulation under patient-specific inlet flow rates measured using magnetic resonance (MR) velocimetry. In the post-CFD analysis, we calculated the blood residence time ([Formula: see text]) and segmented the surface exposed to blood residence time over 1 s ([Formula: see text]). We measured the decrease in volume after flow diversion using the original time-of-flight MR angiography data.

RESULTS

Calcifications were observed in the region with [Formula: see text]. In addition, the ratio of [Formula: see text] to the surface of the aneurysmal domain exhibited a negative relationship with the rate of volume reduction at the 6- and 12-month follow-ups. Post-CFD visualization demonstrated that intra-aneurysmal swirling flow prolonged blood residence time under the condition of a small inlet flow rate, when compared to the aneurysmal volume.

CONCLUSION

The results of this study suggest the usefulness of CFD analysis for the diagnosis of atherosclerotic changes in large CCAs that may affect the therapeutic response after flow diversion.

Significance of Hemodynamics Biomarkers, Tissue Biomechanics and Numerical Simulations in the Pathogenesis of Ascending Thoracic Aortic Aneurysms

Guidelines for the treatment of aortic wall diseases are based on measurements of maximum aortic diameter. However, aortic rupture or dissections do occur for small aortic diameters. Growing scientific evidence underlines the importance of biomechanics and hemodynamics in aortic disease development and progression. Wall shear stress (WWS) is an important hemodynamics marker that depends on aortic wall morphology and on the aortic valve function. WSS could be helpful to interpret aortic wall remodeling and define personalized risk criteria. The complementarity of Computational Fluid Dynamics and 4D Magnetic Resonance Imaging as tools for WSS assessment is a promising reality. The potentiality of these innovative technologies will provide maps or atlases of hemodynamics biomarkers to predict aortic tissue dysfunction. Ongoing efforts should focus on the correlation between these non-invasive imaging biomarkers and clinico-pathologic situations for the implementation of personalized medicine in current clinical practice.

The association between hemodynamics and wall characteristics in human intracranial aneurysms: a review

Hemodynamics plays a key role in the natural history of intracranial aneurysms (IAs). However, studies exploring the association between aneurysmal hemodynamics and the biological and mechanical characteristics of the IA wall in humans are sparse. In this review, we survey the current body of literature, summarize the studies' methodologies and findings, and assess the degree of consensus among them. We used PubMed to perform a systematic review of studies that explored the association between hemodynamics and human IA wall features using different sources. We identified 28 publications characterizing aneurysmal flow and the IA wall: 4 using resected tissues, 17 using intraoperative images, and 7 using vessel wall magnetic resonance imaging (MRI). Based on correlation to IA tissue, higher flow conditions, such as high wall shear stress (WSS) with complex pattern and elevated pressure, were associated with degenerated walls and collagens with unphysiological orientation and faster synthesis. MRI studies strongly supported that low flow, characterized by low WSS and high blood residence time, was associated with thicker walls and post-contrast enhancement. While significant discrepancies were found among those utilized intraoperative images, they generally supported that thicker walls coexist at regions with prolonged residence time and that thinner regions are mainly exposed to higher pressure with complex WSS patterns. The current body of literature supports a theory of two general hemodynamic-biologic mechanisms for IA development. One, where low flow conditions are associated with thickening and atherosclerotic-like remodeling, and the other where high and impinging flow conditions are related to wall degeneration, thinning, and collagen remodeling.

A comprehensive comparison of various patient-specific CFD models of the left atrium for atrial fibrillation patients

BACKGROUND

Recently, advances in medical imaging, segmentation techniques, and high-performance computing have supported the use of patient-specific computational fluid dynamics (CFD) simulations. At present, CFD-compatible atrium geometries can be easily reconstructed from atrium images, providing important insight into the atrial fibrillation (AF) phenomenon, and assistance during therapy selection and surgical procedures. However, the hypothesis assumed for such CFD models should be adequately validated.

AIM

This work aims to perform an extensive study of the different hypotheses that are commonly assumed when performing atrial simulations for AF patients, as well as to evaluate and compare the range of indices that are usually applied to assess thrombus formation within the left atrium appendage (LAA).

METHODS

The atrial geometries of two AF patients have been segmented. The resulting geometries have been registered and interpolated to construct a dynamic mesh, which has been employed to compare the rigid and flexible models. Two families of hemodynamic indices have been calculated and compared: wall shear-based and blood age distribution-based.

RESULTS

The findings of this study illustrate the importance of validating the rigid atrium hypothesis when utilizing an AF CFD model. In particular, the absence of the A-wave contraction does not avoid a certain degree of passive atrial contraction, making the rigid model a poor approximation in some cases. Moreover, a new thrombosis predicting index has been proposed, i.e., M4, which has been shown to predict stasis more effectively than other indicators.

Correlation of Inflow Velocity Ratio Detected by Phase Contrast Magnetic Resonance Angiography with the Bleb Color of Unruptured Intracranial Aneurysms

BACKGROUND

Intraoperative rupture is the most fatal and catastrophic complication of surgery for unruptured intracranial aneurysms (UIAs); thus, it is extremely useful to predict reddish and thin-walled regions of the UIA before surgery. Although several studies have reported a relationship between the hemodynamic characteristics and intracranial aneurysm wall thickness, a consistent opinion is lacking. We aimed to investigate the relationship between objectively and quantitatively evaluated bleb wall color and hemodynamic characteristics using phase-contrast magnetic resonance angiography (PC-MRA).

METHODS

Ten patients diagnosed with UIA who underwent surgical clipping and preoperative magnetic resonance imaging along with PC-MRA were included in this study. Bleb wall color was evaluated from an intraoperative video. Based on the Red (R), Green, and Blue values, bleb wall redness (modified R value; mR) was calculated and compared with the hemodynamic characteristics obtained from PC-MRA.

RESULTS

The wall redness distribution of 18 blebs in 11 UIAs in 10 patients was analyzed. Bleb/neck inflow velocity ratio (Vb/Va: r = 0.66, P = 0.003) strongly correlated with mR, whereas bleb/neck inflow rate ratio (r = 0.58, P = 0.012) correlated moderately. Multivariate regression analysis revealed that only Vb/Va (P = 0.017) significantly correlated with mR. There was no correlation between wall shear stress and mR.

CONCLUSIONS

The bleb redness of UIAs and Vb/Va, calculated using PC-MRA, showed a significantly greater correlation. Thus, it is possible to predict bleb thickness noninvasively before surgery. This will facilitate more detailed pre- and intraoperative strategies for clipping and coiling for safe surgery.

Contribution of atherosclerotic plaque location and severity to the near-wall hemodynamics of the carotid bifurcation: an experimental study and FSI modeling

Atherosclerosis is initiated by endothelial injury that is related to abnormal values of hemodynamic parameters such as wall shear stress (WSS), oscillatory shear index (OSI) and stress phase angle (SPA), which are more common in arterial bifurcations due to the complex structure. An experimental model of human carotid bifurcation with accurate geometrical and mechanical features was set up, and using realistic pulsatile flow rates, the inlet and outlet pressure pulses were measured for normal and stenosed models with 40% and 80% severities at common carotid (CCA), internal carotid (ICA) and external carotid (ECA) arteries. Based on the obtained experimental data, fluid-structure models were developed to obtain WSS, OSI, and SPA and evaluate pathological consequences at different locations. Mild severity had minor impact, however, inducing severe 80% stenosis in each branch led to considerable localized changes of hemodynamic parameters both in the stenosis site and other locations. This included sharp increases in WSS values accompanied by very low values close to zero before and after the peaks. Severe stenosis not only caused significant changes in the local artery, but also in other branches. OSI and SPA were less sensitive to stenosis, although high peaks were observed on bifurcation site for the stenosis at ECA. The interconnection of arteries at carotid bifurcation results in altered pressure/flow patterns in all branches when a stenosis is applied in any site. Such effect confirms pathological findings that atherosclerotic plaques are observed simultaneously in different carotid branches, although with different degrees of plaque growth and severity.

Study on the association of wall shear stress and vessel structural stress with atherosclerosis: An experimental animal study

BACKGROUND AND AIMS

Artery is subject to wall shear stress (WSS) and vessel structural stress (VSS) simultaneously. This study is designed to explore the role of VSS in development of atherosclerosis.

METHODS

Silastic collars were deployed on the carotid to create two constrictions on 13 rabbits for a distinct mechanical environment at the constriction. MRI was performed to visualize arteries' configuration. Animals with high fat (n = 9; Model-group) and normal diet (n = 4; Control-group) were sacrificed after 16 weeks. 3D fluid-structure interaction analysis was performed to quantify WSS and VSS simultaneously.

RESULTS

Twenty plaques were found in Model-group and 3 in Control-group. In Model-group, 8 plaques located proximally to the first constriction (Region-1, close to the heart) and 7 distally to the second (Region-2, close to the head) and 5 plaques were found on the contralateral side of 3 rabbits. Plaques at Region-1 tended to be bigger than those at Region-2 and the macrophage density at these locations was comparable. Minimum time-averaged WSS (TAWSS) in Region-1 was significantly higher than that in Region-2, and both maximum oscillatory shear index (OSI) and particle relative residence time (RRT) were significantly lower. Peak and mean VSS in Region-1 were significantly higher than those in Region-2. Correlation analyses indicated that low TAWSS, high OSI and RRT were only associated with plaque in Region-2, while lesions in Region-1 were only associated with high VSS. Moreover, only VSS was associated with wall thickness of plaque-free regions in both regions.

CONCLUSIONS

VSS might contribute to the initialization and development of atherosclerosis solely or in combination with WSS.

Computational Fluid Dynamics for Cerebral Aneurysms in Clinical Settings

Hemodynamics is thought to play an important role in the pathogenesis of cerebral aneurysms and recent development of computer technology makes it possible to simulate blood flow using high-resolution 3D images within several hours. A lot of studies of computational fluid dynamics (CFD) for cerebral aneurysms were reported; therefore, application of CFD for cerebral aneurysms in clinical settings is reviewed in this article.CFD for cerebral aneurysms using a patient-specific geometry model was first reported in 2003 and it has been revealing that hemodynamics brings a certain contribution to understanding aneurysm pathology, including initiation, growth and rupture. Based on the knowledge of the state-of-the-art techniques, this review treats the decision-making process for using CFD in several clinical settings. We introduce our CFD procedure using digital imaging and communication in medicine (DICOM) datasets of 3D CT angiography or 3D rotational angiography. In addition, we review rupture status, hyperplastic remodeling of aneurysm wall, and recurrence of coiled aneurysms using the hemodynamic parameters such as wall shear stress (WSS), oscillatory shear index (OSI), aneurysmal inflow rate coefficient (AIRC), and residual flow volume (RFV).

Increased tortuosity of basilar artery might be associated with higher risk of aneurysm development

Objectives

We analysed tortuosity of basilar artery (BA) to determine its relationship with the presence of aneurysm.

Methods

We retrospectively analysed 71 patients with BA aneurysms along with 71 age- and risk factors-matched control patients without BA aneurysm. From patients’ medical records, we obtained their history including previous and current diseases and medications. For each patient, we calculated relative length (RL), sum of angle metrics (SOAM), triangular index (TI), product of angle distance (PAD) and inflexion count metrics (ICM). We used t-test and Mann-Whitney U test for continuous variables and χ² test for dichotomised variables. To find independent predictors of BA aneurysm, we employed logistic regression analysis.

Results

We found significant positive correlation between age and SOAM (R = 0.195, p = 0.02) and PAD (R = 0.199, p = 0.018). Our study also showed that patients with BA aneurysm had significantly higher SOAM (0.21 ± 0.16 vs. 0.11 ± 0.08; p < 0.01), PAD (0.30 ± 0.19 vs. 0.18 ± 0.11; p < 0.01), TI (0.23 ± 0.23 vs. 0.10 ± 0.16; p < 0.01) and ICM (0.20 ± 0.16 vs. 0.15 ± 0.11; p = 0.045). In multivariate logistic regression analysis, after adjustment for all possible confounders, SOAM (OR = 1.086; 95% CI 1.046–1.136; p < 0.01) and TI (OR = 1.004; 95%C: 1.002–1.006; p < 0.01) remained independently associated with higher risk of BA aneurysm.

Conclusions

Increased tortuosity of BA is associated with higher risk of its aneurysm development.

Key Points

• Basilar artery sum of angle metrics and product of angle distance are correlated with age.

• Basilar artery tortuosity is independently associated with higher risk of its aneurysm development.

• Basilar artery tortuosity is positively correlated with its diameter and bifurcation angle.

Computational Hemodynamic Modeling of Arterial Aneurysms: A Mini-Review

Arterial aneurysms are pathological dilations of blood vessels, which can be of clinical concern due to thrombosis, dissection, or rupture. Aneurysms can form throughout the arterial system, including intracranial, thoracic, abdominal, visceral, peripheral, or coronary arteries. Currently, aneurysm diameter and expansion rates are the most commonly used metrics to assess rupture risk. Surgical or endovascular interventions are clinical treatment options, but are invasive and associated with risk for the patient. For aneurysms in locations where thrombosis is the primary concern, diameter is also used to determine the level of therapeutic anticoagulation, a treatment that increases the possibility of internal bleeding. Since simple diameter is often insufficient to reliably determine rupture and thrombosis risk, computational hemodynamic simulations are being developed to help assess when an intervention is warranted. Created from subject-specific data, computational models have the potential to be used to predict growth, dissection, rupture, and thrombus-formation risk based on hemodynamic parameters, including wall shear stress, oscillatory shear index, residence time, and anomalous blood flow patterns. Generally, endothelial damage and flow stagnation within aneurysms can lead to coagulation, inflammation, and the release of proteases, which alter extracellular matrix composition, increasing risk of rupture. In this review, we highlight recent work that investigates aneurysm geometry, model parameter assumptions, and other specific considerations that influence computational aneurysm simulations. By highlighting modeling validation and verification approaches, we hope to inspire future computational efforts aimed at improving our understanding of aneurysm pathology and treatment risk stratification.

What does computational fluid dynamics tell us about intracranial aneurysms? A meta-analysis and critical review

Despite the plethora of published studies on intracranial aneurysms (IAs) hemodynamic using computational fluid dynamics (CFD), limited progress has been made towards understanding the complex physics and biology underlying IA pathophysiology. Guided by 1733 published papers, we review and discuss the contemporary IA hemodynamics paradigm established through two decades of IA CFD simulations. We have traced the historical origins of simplified CFD models which impede the progress of comprehending IA pathology. We also delve into the debate concerning the Newtonian fluid assumption used to represent blood flow computationally. We evidently demonstrate that the Newtonian assumption, used in almost 90% of studies, might be insufficient to describe IA hemodynamics. In addition, some fundamental properties of the Navier-Stokes equation are revisited in supplementary material to highlight some widely spread misconceptions regarding wall shear stress (WSS) and its derivatives. Conclusively, our study draws a roadmap for next-generation IA CFD models to help researchers investigate the pathophysiology of IAs.

Elevated Lipid Infiltration Is Associated With Cerebral Aneurysm Rupture

Background: Intracranial aneurysm wall degradation can be associated with lipid infiltration. However, the relationship between lipid infiltration and aneurysm rupture has not been explored quantitatively. To investigate the correlation between lipid infiltration and aneurysm rupture, we utilized patient-specific simulation of low-density lipoprotein (LDL) transport to analyze lipid infiltration in the cerebral aneurysm wall.

Methods: Sixty-two aneurysms were analyzed. Patient blood pressure, plasma LDL concentration, and three-dimensional angiographic images were obtained to simulate LDL transport in aneurysms. Morphological, hemodynamic, and lipid accumulation parameters were compared between ruptures and unruptured groups. Multivariate logistic regression was also performed to determine parameters that are independently associated with rupture.

Results: Size ratio, wall shear stress, low shear area, relative residence time, area-averaged LDL infiltration rate, and maximum LDL infiltration rate were significant parameters in univariate analysis (P < 0.05). Multivariate analysis revealed that only average LDL infiltration remained as a significant variable (P < 0.05). The prediction model derived showed good performance for rupture prediction (AUC, 0.885; 95% CI, 0.794–0.976).

Conclusions: Ruptured aneurysms showed significantly higher LDL infiltration compared to unruptured ones. Our results suggested that lipid infiltration may promote aneurysm rupture. Lipid infiltration characteristics should be considered when assessing aneurysm rupture risk.

Circumferential wall enhancement in evolving intracranial aneurysms on magnetic resonance vessel wall imaging

OBJECTIVE

Recent MR vessel wall imaging studies have indicated intracranial aneurysms in the active state could show circumferential enhancement along the aneurysm wall (CEAW). While ruptured aneurysms frequently show CEAW, CEAW in unruptured aneurysms at the evolving state (i.e., growing or symptomatic) has not been studied in detail. The authors quantitatively assessed the degree of CEAW in evolving unruptured aneurysms by comparing it separately to that in stable unruptured and ruptured aneurysms.

METHODS

A quantitative analysis of CEAW was performed in 26 consecutive evolving aneurysms using MR vessel wall imaging. Three-dimensional T1-weighted fast spin echo sequences were obtained before and after contrast media injection, and the contrast ratio of the aneurysm wall against the pituitary stalk (CRstalk) was calculated as the indicator of CEAW. Aneurysm characteristics of evolving aneurysms were compared with those of 69 stable unruptured and 67 ruptured aneurysms.

RESULTS

The CRstalk values in evolving aneurysms were significantly higher than those in stable aneurysms (0.54 vs 0.34, p < 0.0001), and lower than those in ruptured aneurysms (0.54 vs 0.83, p < 0.0002). In multivariable analysis, CRstalk remained significant when comparing evolving with stable aneurysms (odds ratio [OR] 12.23, 95% confidence interval [CI] 3.53-42.41), and with ruptured aneurysms (OR 0.083, 95% CI 0.022-0.310).

CONCLUSIONS

The CEAW in evolving aneurysms was higher than those in stable aneurysms, and lower than those in ruptured aneurysms. The degree of CEAW may indicate the process leading to rupture of intracranial aneurysms, which can be useful additional information to determine an indication for surgical treatment of unruptured aneurysms.

Influence of Relative Residence Time on Side-Wall Aneurysm Inception

BACKGROUND

Relative residence time (RRT) is a marker of disturbed blood flow, marked by low magnitude and high oscillatory wall shear stress (WSS). The relation between solute residence time in proximity to the vascular endothelium and the atherosclerotic process is well appreciated in the literature.

OBJECTIVE

To assess the influence of RRT on side-wall aneurysm inception to better understand the role of atherosclerosis in aneurysm formation.

METHODS

Fourteen side-wall internal carotid artery aneurysms from the Aneurisk repository which met criteria for parent vessel reconstruction were reconstructed with Vascular Modeling Toolkit. Computational fluid dynamics analysis was carried out in Fluent. RRT was calculated in MATLAB (The MathWorks Inc, Natick, Massachusetts). We analyzed the results for correlations, defined as presence or absence of local elevations in RRT in specific regions of vasculature.

RESULTS

RRT was concluded to be negatively correlated with aneurysm inception in this study of side-wall internal carotid artery aneurysms, with 12/14 cases yielding the absence of local RRT elevations within or in close proximity of the removed ostium. Subsequent analysis of WSS showed that 11 of 14 aneurysms were formed in an atheroprotective environment, with only 1 of 14 formed in an atherogenic environment. Two models were found to be of indeterminate environment.

CONCLUSION

Atherogenesis and atherosclerosis have long been thought to be a major inciting factor responsible for the formation of aneurysms in the cerebral vasculature. We propose that inception of side-wall aneurysms occurs in hemodynamic environments that promote an atheroprotective endothelial phenotype and that the atheroprotective phenotype is therefore aneurysmogenic.

Calcification in Human Intracranial Aneurysms Is Highly Prevalent and Displays Both Atherosclerotic and Nonatherosclerotic Types

OBJECTIVE

Although the clinical and biological importance of calcification is well recognized for the extracerebral vasculature, its role in cerebral vascular disease, particularly, intracranial aneurysms (IAs), remains poorly understood. Extracerebrally, 2 distinct mechanisms drive calcification, a nonatherosclerotic, rapid mineralization in the media and a slower, inflammation driven, atherosclerotic mechanism in the intima. This study aims to determine the prevalence, distribution, and type (atherosclerotic, nonatherosclerotic) of calcification in IAs and assess differences in occurrence between ruptured and unruptured IAs. Approach and Results: Sixty-five 65 IA specimens (48 unruptured, 17 ruptured) were resected perioperatively. Calcification and lipid pools were analyzed nondestructively in intact samples using high resolution (0.35 μm) microcomputed tomography. Calcification is highly prevalent (78%) appearing as micro (<500 µm), meso (500 µm-1 mm), and macro (>1 mm) calcifications. Calcification manifests in IAs as both nonatherosclerotic (calcification distinct from lipid pools) and atherosclerotic (calcification in the presence of lipid pools) with 3 wall types: Type I-only calcification, no lipid pools (20/51, 39%), Type II-calcification and lipid pools, not colocalized (19/51, 37%), Type III-calcification colocalized with lipid pools (12/51, 24%). Ruptured IAs either had no calcifications or had nonatherosclerotic micro- or meso-calcifications (Type I or II), without macro-calcifications.

CONCLUSIONS

Calcification in IAs is substantially more prevalent than previously reported and presents as both nonatherosclerotic and atherosclerotic types. Notably, ruptured aneurysms had only nonatherosclerotic calcification, had significantly lower calcification fraction, and did not contain macrocalcifications. Improved understanding of the role of calcification in IA pathology should lead to new therapeutic targets.

Vessel Wall Imaging Predicts the Presence of Atherosclerotic Lesions in Unruptured Intracranial Aneurysms

BACKGROUND

Recent studies have suggested that magnetic resonance vessel wall imaging (VWI) can visualize thickened intracranial aneurysm wall. We aimed to investigate correlations between VWI findings and intraoperative aneurysm wall features based on the hypothesis that VWI can visualize atherosclerotic changes in unruptured intracranial aneurysm (UIA) walls.

METHODS

A total of 36 microsurgically treated UIAs were retrospectively reviewed. All aneurysms underwent VWI before microsurgical clipping, and fusion images with time-of-flight magnetic resonance angiography were created to localize aneurysm wall enhancement (AWE) lesions. Intraoperatively, 2 neurosurgeons who were blinded to the VWI findings evaluated the aneurysm wall features, giving each aneurysm an atherosclerosis score on a 5-point scale (5: yellowish, 4: whitish, 3: normal, 2: slightly reddish, 1: reddish). We defined atherosclerotic lesions as those having average scores ≥4. We evaluated the rate of correspondence between AWE lesions and atherosclerotic lesions, and the factors associated with AWE.

RESULTS

Sixteen of the 36 UIAs (44%) were identified as AWE. The sensitivity, specificity, positive predictive value, and negative predictive value of correspondence between AWE lesions and atherosclerotic lesions were 79%, 94%, 94%, and 80%, respectively. The average atherosclerosis scores (4.2 ± 0.5 vs. 2.7 ± 0.9; P < 0.001) were significantly higher in aneurysms with AWE. Twelve of 16 UIAs with wall enhancement had wall thinning adjacent to the part with AWE.

CONCLUSIONS

AWE lesions corresponded with intraoperatively confirmed atherosclerotic lesions of UIAs. Detecting these lesions would be valuable in exploring UIAs with wall degeneration.

Impact of Intracranial Aneurysm Morphology and Rupture Status on the Particle Residence Time

BACKGROUND AND PURPOSE

Aneurysm hemodynamics play an important role in aneurysm growth and subsequent rupture. Within the available hemodynamic characteristics, particle residence time (PRT) is relatively unexplored. However, some studies have shown that PRT is related to thrombus formation and inflammation. The goal of this study is to evaluate the association between PRT and aneurysm rupture and morphology.

METHODS

We determined the PRT for 113 aneurysms (61 unruptured, 53 ruptured) based on computational fluid dynamic models. Virtual particles were injected into the parent vessel and followed during multiple cardiac cycles. PRT was defined as the time needed for 99% of the particles that entered an aneurysm to leave the aneurysm. Subsequently, we evaluated the association between PRT, rupture, and morphology (aneurysm type, presence of blebs, or multiple lobulations).

RESULTS

PRT showed no significant difference between unruptured (1.1 seconds interquartile range [IQR .39-2.0 seconds]) and ruptured aneurysms (1.2 seconds [IQR .47-2.3 seconds]). PRT was influenced by aneurysm morphology. Longer PRTs were seen in bifurcation aneurysms (1.3 seconds [IQR .54-2.4 seconds], P = .01) and aneurysms with blebs or multiple lobulations (1.92 seconds [IQR .94-2.8 seconds], P < .001). Four of five partially thrombosed aneurysms had a long residence time (>1.9 seconds).

CONCLUSIONS

Our study shows an influence of aneurysm morphology on PRT. Nevertheless, it suggests that PRT cannot be used to differentiate unruptured and ruptured aneurysms.

Genetic correlates of wall shear stress in a patient-specific 3D-printed cerebral aneurysm model

OBJECTIVES

To study the correlation between wall shear stress and endothelial cell expression in a patient-specific, three-dimensional (3D)-printed model of a cerebral aneurysm.

MATERIALS AND METHODS

A 3D-printed model of a cerebral aneurysm was created from a patient's angiogram. After populating the model with human endothelial cells, it was exposed to media under flow for 24 hours. Endothelial cell morphology was characterized in five regions of the 3D-printed model using confocal microscopy. Endothelial cells were then harvested from distinct regions of the 3D-printed model for mRNA collection and gene analysis via quantitative polymerase chain reaction (qPCR.) Cell morphology and mRNA measurement were correlated with computational fluid dynamics simulations.

RESULTS

The model was successfully populated with endothelial cells, which survived under flow for 24 hours. Endothelial morphology showed alignment with flow in the proximal and distal parent vessel and aneurysm neck, but disorganization in the aneurysm dome. Genetic analysis of endothelial mRNA expression in the aneurysm dome and distal parent vessel was compared with the proximal parent vessels. ADAMTS-1 and NOS3 were downregulated in the aneurysm dome, while GJA4 was upregulated in the distal parent vessel. Disorganized morphology and decreased ADAMTS-1 and NOS3 expression correlated with areas of substantially lower wall shear stress and wall shear stress gradient in computational fluid dynamics simulations.

CONCLUSIONS

Creating 3D-printed models of patient-specific cerebral aneurysms populated with human endothelial cells is feasible. Analysis of these cells after exposure to flow demonstrates differences in both cell morphology and genetic expression, which correlate with areas of differential hemodynamic stress.

Increased tortuosity of ACA might be associated with increased risk of ACoA aneurysm development and less aneurysm dome size: a computer-aided analysis

Objectives

We decided to perform computer-aided analysis of the anterior cerebral artery (ACA) to check for a potential correlation with anterior communicating artery (ACoA) aneurysm presence and growth.

Methods

We retrospectively analyzed the ACA anatomy of 121 patients with ACoA aneurysms along with 121 age, risk factors, and vessel side-matched control patients without an ACoA aneurysm. We obtained their medical history and digital subtraction angiography (DSA) data from their medical records. For each patient’s DSA, we extracted curve representing the course of their ACA and calculated its relative length (RL), sum of angle metrics (SOAM), triangular index (TI), product of angle distance (PAD), and inflection count metrics (ICM).

Results

Patients with ACoA aneurysm had significantly higher RL (0.64 ± 0.23 vs. 0.56 ± 0.22; p < 0.01), SOAM (0.27 ± 0.19 vs. 0.18 ± 0.15; p < 0.01), PAD (0.12 ± 0.13 vs. 0.09 ± 0.11; p = 0.02), and TI (0.57 ± 0.14 vs. 0.44 ± 0.15; p < 0.01). In multivariate logistic regression analysis, after adjustment for possible confounders, SOAM (OR, 1.34; 95% CI, 1.12–1.63; p < 0.01) and TI (OR, 1.84; 95% CI, 1.47–2.35; p < 0.01) remained independently associated with higher risk of ACoA aneurysm. Additionally, we found significant negative correlations between TI and aneurysm dome size (R = − 0.194; p = 0.047).

Conclusions

Increased tortuosity of ACA might increase the risk of ACoA aneurysm development and decrease the risk of aneurysm growth.

Key Points

• Anterior cerebral artery’s sum of angle metrics is associated with hypertension as well as with history of ischemic stroke and myocardial infarction.

• Increased tortuosity of anterior cerebral artery might be associated with anterior communicating artery aneurysm development.

• Tortuosity of anterior cerebral artery is negatively correlated with anterior communicating artery aneurysm dome size.

Tortuosity of the Internal Carotid Artery and Its Clinical Significance in the Development of Aneurysms

Tortuosity of blood vessels is a common angiographic finding that may indicate systemic disease and can be correlated with vascular pathologies. In this work, we determined whether patients with and without internal carotid artery (ICA) aneurysm presented with differences in its tortuosity descriptors. We retrospectively analysed data of 298 patients hospitalized between January 2014 and June 2018. For each patient's imaging data, we extracted a curve representing the ICA course and measured its Relative Length (RL), Sum of Angle Metrics (SOAM), Product of Angle Distance (PAD), Triangular Index (TI), and Inflection Count Metrics (ICM). We found that patients with an ICA aneurysm had significantly lower RL (0.46 ± 0.19 vs. 0.51 ± 0.17; p = 0.023) and significantly higher SOAM (0.39 ± 0.21 vs. 0.32 ± 0.21 p = 0.003), PAD (0.38 ± 0.19 vs. 0.32 ± 0.21; p = 0.011), TI (0.30 ± 0.11 vs. 0.27 ± 0.14; p = 0.034), and ICM (0.30 ± 0.16 vs. 0.22 ± 0.12; p < 0.001). We found that that patients who presented with a subarachnoid hemorrhage had significantly higher PAD (0.46 ± 0.22 vs. 0.35 ± 0.20; p = 0.024). In conclusion, higher tortuosity of ICA is associated with ICA aneurysm presence.

Local Hemodynamic Conditions Associated with Focal Changes in the Intracranial Aneurysm Wall

BACKGROUND AND PURPOSE

Aneurysm hemodynamics has been associated with wall histology and inflammation. We investigated associations between local hemodynamics and focal wall changes visible intraoperatively.

MATERIALS AND METHODS

Computational fluid dynamics models were constructed from 3D images of 65 aneurysms treated surgically. Aneurysm regions with different visual appearances were identified in intraoperative videos: 1) "atherosclerotic" (yellow), 2) "hyperplastic" (white), 3) "thin" (red), 4) rupture site, and 5) "normal" (similar to parent artery), They were marked on 3D reconstructions. Regional hemodynamics was characterized by the following: wall shear stress, oscillatory shear index, relative residence time, wall shear stress gradient and divergence, gradient oscillatory number, and dynamic pressure; these were compared using the Mann-Whitney test.

RESULTS

Hyperplastic regions had lower average wall shear stress (P = .005) and pressure (P = .009) than normal regions. Flow conditions in atherosclerotic and hyperplastic regions were similar but had higher average relative residence time (P = .03) and oscillatory shear index (P = .04) than thin regions. Hyperplastic regions also had a higher average gradient oscillatory number (P = .002) than thin regions. Thin regions had lower average relative residence time (P < .001), oscillatory shear index (P = .006), and gradient oscillatory number (P < .001) than normal regions, and higher average wall shear stress (P = .006) and pressure (P = .009) than hyperplastic regions. Thin regions tended to be aligned with the flow stream, while atherosclerotic and hyperplastic regions tended to be aligned with recirculation zones.

CONCLUSIONS

Local hemodynamics is associated with visible focal wall changes. Slow swirling flow with low and oscillatory wall shear stress was associated with atherosclerotic and hyperplastic changes. High flow conditions prevalent in regions near the flow impingement site characterized by higher and less oscillatory wall shear stress were associated with local "thinning" of the wall.