Determination of a shear rate threshold for thrombus formation in intracranial aneurysms

Spontaneous disappearance of a small unruptured cerebral aneurysm in the clinoid segment of the internal carotid artery: A case report and literature review

Background

Various degrees of thrombosis have been reported in patients with giant aneurysms. However, small, unruptured aneurysms rarely resolve spontaneously. Herein, we report a case of a small unruptured aneurysm in the clinoid segment (C3) of the left internal carotid artery (ICA) that showed almost complete occlusion at the 1-year follow-up.

Case Description

A 66-year-old woman developed a subarachnoid hemorrhage on the left side of the perimesencephalic cistern. Cerebral angiography performed on admission revealed no evidence of hemorrhage. Subsequent cerebral angiography on day 12 revealed a dissecting aneurysm on a branch of the superior cerebellar artery (SCA), and the patient underwent parental artery occlusion with 25% n-butyl-2-cyanoacrylate. The postoperative course was uneventful, and the patient was discharged on day 22 with a modified Rankin Scale score of 1. The 1 year follow-up cerebral angiogram demonstrated that the dissecting aneurysm in the SCA branch remained occluded. Notably, a small 2-mm unruptured aneurysm in the clinoid segment (C3) of the left ICA, which was present at the onset of subarachnoid hemorrhage, was almost completely occluded without intervention. Magnetic resonance angiography 1 year after spontaneous resolution of the aneurysm showed no apparent recurrence.

Conclusion

This case highlights that even small, unruptured aneurysms can develop spontaneous occlusions.

Factors associated with false lumen changes in patients with superior mesenteric artery dissection

BACKGROUND

False lumen changes (FLCs) are the main reference for the prognosis judgment and treatment plan selection for type IIa superior mesenteric artery dissection (SMAD).

METHODS

For this retrospective study, 55 patients with symptomatic type IIa SMAD were included. Computational fluid dynamics (CFD) analysis was used to explore the hemodynamic basis of FLCs. Correlation and multiple linear regression analyses were performed to identify clinical, morphological and hemodynamic factors associated with FLCs.

RESULTS

The FLCs of patients with successful conservative treatment (n = 29) are significantly higher than those with failed conservative treatment (n = 26) (58.5 ± 21.1% vs 10.9 ± 17.4%, p < 0.0001). Positive correlations were seen between FLCs and the morphological parameters false lumen length (FLL)/dissection entrance length (DEL) and FLL. In terms of hemodynamic parameters, negative correlations were seen between FLCs and time-averaged wall shear stress (TAWSS), vorticity, and high areas of TAWSS and vorticity, whereas positive correlations were seen between FLCs and oscillatory shear index (OSI), relative residence time (RRT), and high areas of OSI and RRT. Multiple linear regression analysis identified symptom duration (odds ratio [OR], 0.93; 95% CI, 0.91-0.96; p < 0.0001), FLL/DEL (OR, 1.30; 95% CI, 1.01-1.67; p = 0.044), and high RRT area (OR, 2.03; 95% CI, 1.48-2.78; p < 0.0001) as predictors of FLCs.

CONCLUSION

The clinical predictor symptom duration, morphological factor FLL/DEL, and the hemodynamic factor high RRT area can serve as predictors of FLCs in patients with symptomatic type IIa SMAD.

Hemodynamics of thrombus formation in intracranial aneurysms: An in silico observational study

How prevalent is spontaneous thrombosis in a population containing all sizes of intracranial aneurysms? How can we calibrate computational models of thrombosis based on published data? How does spontaneous thrombosis differ in normo- and hypertensive subjects? We address the first question through a thorough analysis of published datasets that provide spontaneous thrombosis rates across different aneurysm characteristics. This analysis provides data for a subgroup of the general population of aneurysms, namely, those of large and giant size (>10 mm). Based on these observed spontaneous thrombosis rates, our computational modeling platform enables the first in silico observational study of spontaneous thrombosis prevalence across a broader set of aneurysm phenotypes. We generate 109 virtual patients and use a novel approach to calibrate two trigger thresholds: residence time and shear rate, thus addressing the second question. We then address the third question by utilizing this calibrated model to provide new insight into the effects of hypertension on spontaneous thrombosis. We demonstrate how a mechanistic thrombosis model calibrated on an intracranial aneurysm cohort can help estimate spontaneous thrombosis prevalence in a broader aneurysm population. This study is enabled through a fully automatic multi-scale modeling pipeline. We use the clinical spontaneous thrombosis data as an indirect population-level validation of a complex computational modeling framework. Furthermore, our framework allows exploration of the influence of hypertension in spontaneous thrombosis. This lays the foundation for in silico clinical trials of cerebrovascular devices in high-risk populations, e.g., assessing the performance of flow diverters in aneurysms for hypertensive patients.

Study and modeling of the thrombosis of small cerebral aneurysms, with and without flow diverter, by the lattice Boltzmann method

BACKGROUND AND OBJECTIVE

Small cerebral aneurysms are currently commonly treated non-invasively by flow diverter device. These stents lead to thrombotic occlusion of the aneurysm soon after their placement. The purpose of this work is to model clotting into intracranial aneurysms with and without stents, using a non-Newtonian of blood behavior, and to investigate the importance of stent to generate desired thrombus in intracranial aneurysms.

METHOD

The description of blood flow is made by the Boltzmann lattice equations, while thrombosis is modeled by the "fluid age" model. The lattice Boltzmann method is a computational technique for simulating fluid dynamics. The method is based on a mesoscopic approach, where the fluid is represented by a set of particles that move and interact on a grid. The model for blood coagulation is described by lattice Boltzmann Method, and it doesn't take into account the complicated coagulation pathway, this main idea is developed using the model of residence time of blood: all fluid in the domain is assumed to be capable of clotting, given enough time. The fluid age is measured by a passive scalar using a transport equation, and the node coagulates if the fluid age increases enough. Three small aneurysms of different sizes and shapes with three stents of various porosities were used to test the ability of the model to predict thrombosis. The "occlusion rate" parameter is used to assess the effectiveness of the flow diverter device.

RESULTS

For the large aspect ratio factor, the occlusion is: 91% for flow diverter devise with seven struts. For medium aspect ratio, a rate of 80% is achieved. An occlusion rate of slightly more than 30% is obtained for very small aneurysms with low aspect ratio. The Newtonian model underestimates the volume of thrombosis generated. The difference in the prediction of the thrombosis volume between the Newtonian and no-Newtonian Carreau-Yasuda models is approximately 10%.

CONCLUSION

The occlusion rate is proportional to the aspect ratio form factor. For the large and medium aspect ratio factors, the occlusion is satisfactory. Concerning very small aneurysms with low aspect ratio, aneurysm occlusion is low. This rate can be improved to almost complete occlusion if the flow diverter device is doubled. The generality of the model suggests its extensibility toward any other type of thrombosis (stenosis, thrombosis in aortic aneurysms).

A comprehensive assessment of self-reported symptoms among patients harboring an unruptured intracranial aneurysm

Background

Approximately 3.2%-6% of the general population harbor an unruptured intracranial aneurysm (UIA). Ruptured aneurysms represent a significant healthcare burden, and preventing rupture relies on early detection and treatment. Most patients with UIAs are asymptomatic, and many of the symptoms associated with UIAs are nonspecific, which makes diagnosis challenging. This study explored symptoms associated with UIAs, the rate of resolution of such symptoms after microsurgical treatment, and the likely pathophysiology.

Methods

A retrospective review of patients with UIAs who underwent microsurgical treatment from January 1, 2014, to December 31, 2020, at a single quaternary center were identified. Analyses included the prevalence of nonspecific symptoms upon clinical presentation and postoperative follow-up; comparisons of symptomatology by aneurysmal location; and comparisons of patient demographics, aneurysmal characteristics, and poor neurologic outcome at postoperative follow-up stratified by symptomatic versus asymptomatic presentation.

Results

The analysis included 454 patients; 350 (77%) were symptomatic. The most common presenting symptom among all 454 patients was headache (n = 211 [46%]), followed by vertigo (n = 94 [21%]), cognitive disturbance (n = 68[15%]), and visual disturbance (n = 64 [14%]). Among 328 patients assessed for postoperative symptoms, 258 (79%) experienced symptom resolution or improvement.

Conclusion

This cohort demonstrates that the clinical presentation of patients with UIAs can be associated with vague and nonspecific symptoms. Early detection is crucial to prevent aneurysmal subarachnoid hemorrhage. It is imperative that physicians not rule out aneurysms in the setting of nonspecific neurologic symptoms.

Evaluation of stent effect and thrombosis generation with different blood rheology on an intracranial aneurysm by the Lattice Boltzmann method

BACKGROUND AND OBJECTIVE

Treatment of intracranial aneurysms with flow-diverting stents prevents rupture by reducing blood flow and creating thrombosis within the aneurysm. This paper aims to assess the hemodynamic effect of placing stents with different struts (0, 3, 5, 7 struts) on intracranial aneurysms and to propose a simple prediction model of thrombosis zone without any further computational cost.

METHOD

Lattice Boltzmann method with different rheological models (Newtonian, Carreau-Yasuda, KL) of blood are used to study the hemodynamic effect of flow-diverting stents in the aneurysm. Pulsatile flow boundary conditions were applied in the inlet of the artery. The average Reynolds number was resulting Re = 111. The Lagrangian tracking of the particle was developed to assess the intra-aneurysmal blood stagnation. To predict the probable thrombose zone induced by flow-diverting stents, the shear rate threshold is utilized to determine the nodes of fluid to clot.

RESULTS

The results show that the flow patterns into the aneurysmal sac develop a vortex, decreasing after stent placement until disappearance for the stent with seven struts (porosity 71.4%). Velocity, shear rate, shear stress, trajectory, path length, and occlusion rate are compared before and after stent placement. These parameters decrease inversely with the porosity of the stent. The three models yield a closes result of the (velocity, shear rate, occlusion rate). Tracking the fluid-particle trajectory shows that the length of the particle paths decreases with the number of struts causing fluid to slow down and increase, consequently, the residence time into the sac.

CONCLUSION

The flow-diverting stents placement cause the reduction of dynamic flow within aneurysm. The reduction effect is almost the same below five struts (80% of porosity). The results show that, if our objective is restricted to estimating the hemodynamic effect, measured by (velocity, shear rate, occlusion rate), the differences between rheological behavior models are, practically, not significant, and the models can be used indifferently.

The effect of stiffened diabetic red blood cells on wall shear stress in a reconstructed 3D microaneurysm

Blood flow within the vasculature of the retina has been found to influence the progression of diabetic retinopathy. In this research cell resolved blood flow simulations are used to study the pulsatile flow of whole blood through a segmented retinal microaneurysm. Images were collected using adaptive optics optical coherence tomography of the retina of a patient with diabetic retinopathy, and a sidewall (sacciform) microaneurysm was segmented from the volumetric data. The original microaneurysm neck width was varied to produce two additional aneurysm geometries in order to probe the influence of neck width on the transport of red blood cells and platelets into the aneurysm. Red blood cell membrane stiffness was also increased to resolve the impact of rigid red blood cells, as a result of diabetes, in blood flow. Wall shear stress and wall shear stress gradients were calculated throughout the aneurysm domains, and the quantification of the influence of the red blood cells is presented. Average wall shear stress and wall shear stress gradients increased due to the increase of red blood cell membrane stiffness. Stiffened red blood cells were also found to induce higher local wall shear stress and wall shear stress gradients as they passed through the leading and draining parental vessels. Stiffened red blood cells were found to penetrate the aneurysm sac more than healthy red blood cells, as well as decreasing the margination of platelets to the vessel walls of the parental vessel, which caused a decrease in platelet penetration into the aneurysm sac.

The Role of Patient-Specific Morphological Features of the Left Atrial Appendage on the Thromboembolic Risk Under Atrial Fibrillation

Background

A large majority of thrombi causing ischemic complications under atrial fibrillation (AF) originate in the left atrial appendage (LAA), an anatomical structure departing from the left atrium, characterized by a large morphological variability between individuals. This work analyses the hemodynamics simulated for different patient-specific models of LAA by means of computational fluid-structure interaction studies, modeling the effect of the changes in contractility and shape resulting from AF.

Methods

Three operating conditions were analyzed: sinus rhythm, acute atrial fibrillation, and chronic atrial fibrillation. These were simulated on four patient-specific LAA morphologies, each associated with one of the main morphological variants identified from the common classification: chicken wing, cactus, windsock, and cauliflower. Active contractility of the wall muscle was calibrated on the basis of clinical evaluations of the filling and emptying volumes, and boundary conditions were imposed on the fluid to replicate physiological and pathological atrial pressures, typical of the various operating conditions.

Results

The LAA volume and shear strain rates were analyzed over time and space for the different models. Globally, under AF conditions, all models were well aligned in terms of shear strain rate values and predicted levels of risk. Regions of low shear rate, typically associated with a higher risk of a clot, appeared to be promoted by sudden bends and focused at the trabecule and the lobes. These become substantially more pronounced and extended with AF, especially under acute conditions.

Conclusion

This work clarifies the role of active and passive contraction on the healthy hemodynamics in the LAA, analyzing the hemodynamic effect of AF that promotes clot formation. The study indicates that local LAA topological features are more directly associated with a thromboembolic risk than the global shape of the appendage, suggesting that more effective classification criteria should be identified.

Non-Newtonian Blood Modeling in Intracranial Aneurysm Hemodynamics: Impact on the Wall Shear Stress and Oscillatory Shear Index Metrics for Ruptured and Unruptured Cases

When simulating blood flow in intracranial aneurysms (IAs), the Newtonian model seems to be ubiquitous. However, analyzing the results from the few studies on this subject, the doubt remains on whether it is necessary to use non-Newtonian models in computational fluid dynamics (CFD) simulations of cerebral vascular flows. The objective of this study is to investigate whether different rheology models would influence the hemodynamic parameters related to the wall shear stress (WSS) for ruptured and unruptured IA cases, especially because ruptured aneurysms normally have morphological features, such as lobular regions and blebs, that could trigger non-Newtonian phenomena in the blood flow due to low shear rates. Using CFD in an open-source framework, we simulated four ruptured and four unruptured patient-specific aneurysms to assess the influence of the blood modeling on the main hemodynamic variables associated with aneurysm formation, growth, and rupture. Results for WSS and oscillatory shear index (OSI) and their metrics were obtained using Casson and Carreau-Yasuda non-Newtonian models and were compared with those obtained using the Newtonian model. We found that all differences between non-Newtonian and the Newtonian models were consistent among all cases irrespective of their rupture status. We further found that the WSS at peak systole is overestimated by more than 50% by using the non-Newtonian models, but its metrics based on time and surface averaged values are less affected-the maximum relative difference among the cases is 7% for the Casson model. On the other hand, the surface-averaged OSI is underestimated by more than 30% by the non-Newtonian models. These results suggest that it is recommended to investigate different blood rheology models in IAs simulations when specific parameters to characterize the flow are needed, such as peak-systole WSS and OSI.

In-silico trial of intracranial flow diverters replicates and expands insights from conventional clinical trials

The cost of clinical trials is ever-increasing. In-silico trials rely on virtual populations and interventions simulated using patient-specific models and may offer a solution to lower these costs. We present the flow diverter performance assessment (FD-PASS) in-silico trial, which models the treatment of intracranial aneurysms in 164 virtual patients with 82 distinct anatomies with a flow-diverting stent, using computational fluid dynamics to quantify post-treatment flow reduction. The predicted FD-PASS flow-diversion success rates replicate the values previously reported in three clinical trials. The in-silico approach allows broader investigation of factors associated with insufficient flow reduction than feasible in a conventional trial. Our findings demonstrate that in-silico trials of endovascular medical devices can: (i) replicate findings of conventional clinical trials, and (ii) perform virtual experiments and sub-group analyses that are difficult or impossible in conventional trials to discover new insights on treatment failure, e.g. in the presence of side-branches or hypertension.

In-silico trials rely on virtual populations and interventions simulated using patient-specific models and may offer a solution to lower costs. Here, the authors present the flow diverter performance assessment in-silico trial, which models the treatment of intracranial aneurysms with a flow-diverting stent.

Association Between Aneurysmal Haemodynamics and Device Microstructural Characteristics After Flow-Diversion Treatments With Dual Stents of Different Sizes: A Numerical Study

Objectives

Treating intracranial aneurysms with flow-diverting stents sometimes requires deployment of a second device. Herein we quantify the sizing effects of devices in dual-stent treatments upon the final stent microstructure and the post-treatment aneurysmal haemodynamics.

Methods

Fifteen sidewall ICA aneurysm geometries were included. Using a virtual stenting technique, we implanted either one or two stents for each aneurysm treatment considered, with each stent specified as one of two different sizes, yielding a total of two single-stent and fouir dual-stent treatment scenarios for each aneurysm. Three stent microstructural parameters and nine aneurysmal haemodynamic parameters were quantified and systematically compared across the 90 treatment scenarios.

Results

Deployment of a second stent further reduced the aneurysmal inflow rate (IR) and energy loss (EL) by, respectively, 14 ± 11% (p = 0.001) and 9 ± 12% (p = 0.056), relative to the untreated condition. Sizing effects of the earlier-deployed stent led to largest differences of 6.9% for the final IR reduction and 11.1% for the EL, whereas sizing effects from the later-deployed stent were minor (≤2.1%). The change in stent pore size was the only microstructural parameter demonstrating a strong correlation with the reduction in the post-treatment aneurysmal haemodynamics, in terms of the IR (r = 0.50, p < 0.001) and pressure drop (r = 0.63, p < 0.001).

Conclusion

Size of the earlier-deployed stent has substantial effects on the final haemodynamic outcomes after dual-stent treatment. The average pore size of stent wires at the aneurysm orifice shows promise as a potential index for predicting the efficacy of flow-diversion treatments.

Flow diverter modeled as heterogeneous and anisotropic porous medium: Simulation, experimental validation and case analysis

Simulation of flow diverter (FD) treated aneurysm can evaluate treatment efficacy and aid treatment planning. However, explicit modeling of thin wires of FD impose extremely high demand of computational resources and time, which limit its use in time-sensitive presurgical planning. One alternative approach is to model FD as homogenous porous medium, which saves time but with compromise in accuracy. We proposed a new method to model FD as heterogeneous and anisotropic porous medium whose properties were determined from local porosity. The new method was validated by comparing with PIV measurement from an in-vitro phantom. Simulation result was in good agreement with experimental measurement. Four patient cases were further analyzed to compare the new method with the homogenous porous media method. Results showed that in patient cases with curved artery, new method was preferred over the homogenous method, as the assumption of homogenous porosity led to overpredicted flow reduction effect by as much as 87.9%, which may lead to overoptimistic decision making and poor prognosis. Our new method can provide timely and accurate simulation to aid in the treatment planning of aneurysms.

Complete spontaneous thrombosis in unruptured non-giant intracranial aneurysms: A case report and systematic review

BACKGROUND AND AIM

Spontaneous partial or complete thrombosis of saccular unruptured intracranial aneurysm (UIAs) is a known occurrence in giant aneurysms. However, spontaneous complete thrombosis of non-giant aneurysms is a rare event in the natural history of UIAs. The aim of this paper is to report on the cases from literature of complete spontaneous thrombosis with a view to identify possible factors associated with this phenomenon.

MATERIAL AND METHODS

We performed a systematic review of the current literature on spontaneous complete thrombosis of saccular, non-giant, unruptured UIAs, including a case that we treated at our institution. We analysed the possible risk factors for thrombosis, association with ischemic events, rupture and recanalization. We reviewed the possible management's strategies for this group of patients described in literature to date.

RESULTS

We identified 26 patients for a total of 27 thrombosed aneurysms from the literature review (including our case). Thrombosis was prevalent in women, in the anterior circulation and in larger aneurysms. Endovascular events in the parent artery, either spontaneous or iatrogenic, were associated with spontaneous thrombosis in 4 cases. In 47 % of cases an antiplatelet treatment (AP) was introduced. Rupture and recanalization of the aneurysm were observed in 14 % and 33 % respectively. A larger size was the only factor statistically associated with rupture (P = 0041). AP was not statistically associated with recanalization or rupture of the aneurysm.

CONCLUSION

Complete spontaneous thrombosis is not a curative event. Its natural history is associated with recanalization, rupture and ischemic stroke. Conservative treatment with a clinical-radiological follow up and treatment with AP is a safe option for small aneurysms. Definitive aneurysmal exclusion should be considered in medium and large aneurysms due to the significant risks associated with untreated aneurysms.

Evaluation of a Desktop 3D Printed Rigid Refractive-Indexed-Matched Flow Phantom for PIV Measurements on Cerebral Aneurysms

PURPOSE

Fabrication of a suitable flow model or phantom is critical to the study of biomedical fluid dynamics using optical flow visualization and measurement methods. The main difficulties arise from the optical properties of the model material, accuracy of the geometry and ease of fabrication.

METHODS

Conventionally an investment casting method has been used, but recently advancements in additive manufacturing techniques such as 3D printing have allowed the flow model to be printed directly with minimal post-processing steps. This study presents results of an investigation into the feasibility of fabrication of such models suitable for particle image velocimetry (PIV) using a common 3D printing Stereolithography process and photopolymer resin.

RESULTS

An idealised geometry of a cerebral aneurysm was printed to demonstrate its applicability for PIV experimentation. The material was shown to have a refractive index of 1.51, which can be refractive matched with a mixture of de-ionised water with ammonium thiocyanate (NH₄SCN). The images were of a quality that after applying common PIV pre-processing techniques and a PIV cross-correlation algorithm, the results produced were consistent within the aneurysm when compared to previous studies.

CONCLUSIONS

This study presents an alternative low-cost option for 3D printing of a flow phantom suitable for flow visualization simulations. The use of 3D printed flow phantoms reduces the complexity, time and effort required compared to conventional investment casting methods by removing the necessity of a multi-part process required with investment casting techniques.

A computational model for prediction of clot platelet content in flow-diverted intracranial aneurysms

Treatment of intracranial aneurysms with flow-diverting stents is a safe and minimally invasive technique. The goal is stable embolisation that facilitates stent endothelialisation, and elimination of the aneurysm. However, it is not fully understood why some aneurysms fail to develop a stable clot even with sufficient levels of flow reduction. Computational prediction of thrombus formation dynamics can help predict the post-operative response in such challenging cases. In this work, we propose a new model of thrombus formation and platelet dynamics inside intracranial aneurysms. Our novel contribution combines platelet activation and transport with fibrin generation, which is key to characterising stable and unstable thrombus. The model is based on two types of thrombus inside aneurysms: red thrombus (fibrin- and erythrocyte-rich) can be found in unstable clots, while white thrombus (fibrin- and platelet-rich) can be found in stable clots. The thrombus generation model is coupled to a CFD model and the flow-induced platelet index (FiPi) is defined as a quantitative measure of clot stability. Our model is validated against an in vitro phantom study of two flow-diverting stents with different sizing. We demonstrate that our model accurately predicts the lower thrombus stability in the oversized stent scenario. This opens possibilities for using computational simulations to improve endovascular treatment planning and reduce adverse events, such as delayed haemorrhage of flow-diverted aneurysms.

Cell-resolved blood flow simulations of saccular aneurysms: effects of pulsatility and aspect ratio

We study the effect of pulsatile flow on the transport of red blood cells (RBCs) and platelets into aneurysm geometries with varying dome-to-neck aspect ratios (AR). We use a validated two-dimensional lattice Boltzmann model for blood plasma with a discrete element method for both RBCs and platelets coupled by the immersed boundary method. Flow velocities and vessel diameters were matched with measurements of cerebral perforating arteries and flow was driven by a synthetic heartbeat curve typical for such vessel sizes. We observe a flow regime change as the aspect ratio increases from a momentum-driven regime in the small aspect ratio to a shear-driven regime in the larger aspect ratios. In the small aspect ratio case, we see the development of a re-circulation zone that exhibits a layering of high (greater than or equal to 7 s) and low (less than 7 s) residence cells. In the shear-driven regime, we see high and low residence cells well mixed, with an increasing population of cells that are trapped inside the aneurysm as the aspect ratio increases. In all cases, we observe aneurysms that are platelet-rich and red blood cell-poor when compared with their respective parental vessel populations. Pulsatility also plays a role in the small aspect ratio as we observe a smaller population of older trapped cells along the aneurysm wall in the pulsatile case when compared with a steady flow case. Pulsatility does not have a significant effect in shear-driven regime aspect ratios.

Estimation of Diabetic Retinal Microaneurysm Perfusion Parameters Based on Computational Fluid Dynamics Modeling of Adaptive Optics Scanning Laser Ophthalmoscopy

Diabetic retinopathy (DR) is a leading cause of vision loss worldwide. Microaneurysms (MAs), which are abnormal outpouchings of the retinal vessels, are early and hallmark lesions of DR. The presence and severity of MAs are utilized to determine overall DR severity. In addition, MAs can directly contribute to retinal neural pathology by leaking fluid into the surrounding retina, causing abnormal central retinal thickening and thereby frequently leading to vision loss. Vascular perfusion parameters such as shear rate (SR) or wall shear stress (WSS) have been linked to blood clotting and endothelial cell dysfunction, respectively in non-retinal vasculature. However, despite the importance of MAs as a key aspect of diabetic retinal pathology, much remains unknown as to how structural characteristics of individual MAs are associated with these perfusion attributes. MA structural information obtained on high resolution adaptive optics scanning laser ophthalmoscopy (AOSLO) was utilized to estimate perfusion parameters through Computational Fluid Dynamics (CFD) analysis of the AOSLO images. The HemeLB flow solver was used to simulate steady-state and time-dependent fluid flow using both commodity hospital-based and high performance computing resources, depending on the degree of detail required in the simulations. Our results indicate that WSS is lowest in MA regions furthest away from the feeding vessels. Furthermore, areas of low SR are associated with clot location in saccular MAs. These findings suggest that morphology and CFD estimation of perfusion parameters may be useful tools for determining the likelihood of clot presence in individual diabetic MAs.

Thrombosis in Cerebral Aneurysms and the Computational Modeling Thereof: A Review

Thrombosis is a condition closely related to cerebral aneurysms and controlled thrombosis is the main purpose of endovascular embolization treatment. The mechanisms governing thrombus initiation and evolution in cerebral aneurysms have not been fully elucidated and this presents challenges for interventional planning. Significant effort has been directed towards developing computational methods aimed at streamlining the interventional planning process for unruptured cerebral aneurysm treatment. Included in these methods are computational models of thrombus development following endovascular device placement. The main challenge with developing computational models for thrombosis in disease cases is that there exists a wide body of literature that addresses various aspects of the clotting process, but it may not be obvious what information is of direct consequence for what modeling purpose (e.g., for understanding the effect of endovascular therapies). The aim of this review is to present the information so it will be of benefit to the community attempting to model cerebral aneurysm thrombosis for interventional planning purposes, in a simplified yet appropriate manner. The paper begins by explaining current understanding of physiological coagulation and highlights the documented distinctions between the physiological process and cerebral aneurysm thrombosis. Clinical observations of thrombosis following endovascular device placement are then presented. This is followed by a section detailing the demands placed on computational models developed for interventional planning. Finally, existing computational models of thrombosis are presented. This last section begins with description and discussion of physiological computational clotting models, as they are of immense value in understanding how to construct a general computational model of clotting. This is then followed by a review of computational models of clotting in cerebral aneurysms, specifically. Even though some progress has been made towards computational predictions of thrombosis following device placement in cerebral aneurysms, many gaps still remain. Answering the key questions will require the combined efforts of the clinical, experimental and computational communities.

Predicting false lumen thrombosis in patient-specific models of aortic dissection

Aortic dissection causes splitting of the aortic wall layers, allowing blood to enter a ‘false lumen’ (FL). For type B dissection, a significant predictor of patient outcomes is patency or thrombosis of the FL. Yet, no methods are currently available to assess the chances of FL thrombosis. In this study, we present a new computational model that is capable of predicting thrombus formation, growth and its effects on blood flow under physiological conditions. Predictions of thrombus formation and growth are based on fluid shear rate, residence time and platelet distribution, which are evaluated through convection–diffusion–reaction transport equations. The model is applied to a patient-specific type B dissection for which multiple follow-up scans are available. The predicted thrombus formation and growth patterns are in good qualitative agreement with clinical data, demonstrating the potential applicability of the model in predicting FL thrombosis for individual patients. Our results show that the extent and location of thrombosis are strongly influenced by aortic dissection geometry that may change over time. The high computational efficiency of our model makes it feasible for clinical applications. By predicting which aortic dissection patient is more likely to develop FL thrombosis, the model has great potential to be used as part of a clinical decision-making tool to assess the need for early endovascular intervention for individual dissection patients.

A physical description of the adhesion and aggregation of platelets

The early stages of clot formation in blood vessels involve platelet adhesion-aggregation. Although these mechanisms have been extensively studied, gaps in their understanding still persist. We have performed detailed in vitro experiments, using the well-known Impact-R device, and developed a numerical model to better describe and understand this phenomenon. Unlike previous studies, we took into account the differential role of pre-activated and non-activated platelets, as well as the three-dimensional nature of the aggregation process. Our investigation reveals that blood albumin is a major parameter limiting platelet aggregate formation in our experiment. Simulations are in very good agreement with observations and provide quantitative estimates of the adhesion and aggregation rates that are hard to measure experimentally. They also provide a value of the effective diffusion of platelets in blood subject to the shear rate produced by the Impact-R.

Virtual flow-diverter treatment planning: The effect of device placement on bifurcation aneurysm haemodynamics

Bifurcation aneurysms account for a large fraction of cerebral aneurysms and often present morphologies that render traditional endovascular treatments, such as coiling, challenging and problematic. Flow-diverter stents offer a potentially elegant treatment option for such aneurysms, but clinical use of these devices remains controversial. Specifically, the deployment of a flow-diverter device in a bifurcation entails jailing one or more potentially vital vessels with a low-porosity mesh designed to restrict the flow. When multiple device placement configurations exist, the most appropriate clinical decision becomes increasingly opaque. In this study, three bifurcation aneurysm geometries were virtually treated by flow-diverter device. Each aneurysm was selected to offer two possible device deployment positions. Flow-diverters similar to commercially available designs were deployed with a fast-deployment algorithm before transient and steady state computational fluid dynamics simulations were performed. Reductions in aneurysm inflow, mean wall shear stress and maximum wall shear stress, all factors often linked with aneurysm treatment outcome, were compared for different device configurations in each aneurysm. In each of the three aneurysms modelled, a particular preferential device placement was shown to offer superior performance with the greatest reduction in the flow metrics considered. In all the three aneurysm geometries, substantial variations in inflow reduction (up to 25.3%), mean wall shear stress reduction (up to 14.6%) and maximum wall shear stress reduction (up to 12.1%) were seen, which were all attributed to device placement alone. Optimal device placement was found to be non-trivial and highly aneurysm specific; in only one-third of the simulated geometries, the best overall performance was achieved by deploying a device in the daughter vessel with the highest flow rate. Good correspondence was seen between transient results and steady state computations that offered a significant reduction in simulation run time. If accurate steady state computations are combined with the fast-deployment algorithm used, the modest run time and corresponding hardware make a virtual treatment pipeline in the clinical setting a meaningful possibility.

Biology of Saccular Cerebral Aneurysms: A Review of Current Understanding and Future Directions

Understanding the biology of intracranial aneurysms is a clinical quandary. How these aneurysms form, progress, and rupture is poorly understood. Evidence indicates that well-established risk factors play a critical role, along with immunologic factors, in their development and clinical outcomes. Much of the expanding knowledge of the inception, progression, and rupture of intracranial aneurysms implicates inflammation as a critical mediator of aneurysm pathogenesis. Thus, therapeutic targets exploiting this arm of aneurysm pathogenesis have been implemented, often with promising outcomes.

A computational model based on fibrin accumulation for the prediction of stasis thrombosis following flow-diverting treatment in cerebral aneurysms

Flow diverters, the specially designed low porosity stents, have been used to redirect blood flow from entering aneurysm, which induces flow stasis in aneurysm and promote thrombosis for repairing aneurysm. However, it is not clear how thrombus develops following flow-diversion treatment. Our objective was to develop a computation model for the prediction of stasis-induced thrombosis following flow-diversion treatment in cerebral aneurysms. We proposed a hypothesis to initiate coagulation following flow-diversion treatment. An experimental model was used by ligating rat's right common carotid artery (RCCA) to create flow-stasis environment. Thrombus formed in RCCA as a result of flow stasis. The fibrin distributions in different sections along the axial length of RCCA were measured. The fibrin distribution predicted by our computational model displayed a trend of increase from the proximal neck to the distal tip, consistent with the experimental results on rats. The model was applied on a saccular aneurysm treated with flow diverter to investigate thrombus development following flow diversion. Thrombus was predicted to form inside the sac, and the aneurysm was occluded with only a small remnant neck remained. Our model can serve as a tool to evaluate flow-diversion treatment outcome and optimize the design of flow diverters.