Computational modelling of clot development in patient-specific cerebral aneurysm cases

Influence of vortical structures on fibrin clot formation in cerebral aneurysms: A two-dimensional computational study

Thrombosis is an important contributor to cerebral aneurysm growth and progression. A number of sophisticated multiscale and multiphase in silico models have been developed with a view towards interventional planning. Many of these models are able to account for clotting outcomes, but do not provide detailed insight into the role of flow during clot development. In this study, we present idealised, two-dimensional in silico cerebral fibrin clot model based on computational fluid dynamics (CFD), biochemical modelling and variable porosity, permeability, and diffusivity. The model captures fibrin clot growth in cerebral aneurysms over a period at least 1000 s in five different geometries. The fibrin clot growth results were compared to an experiment presented in literature. The biochemistry was found to be more sensitive to mesh size compared to the haemodynamics, while larger timesteps overpredicted clot size in pulsatile flow. When variable diffusivity was used, the predicted clot size was 25.4% lesser than that with constant diffusivity. The predicted clot size in pulsatile flow was 14.6% greater than in plug flow. Different vortex modes were observed in plug and pulsatile flow; the latter presented smaller intermediate modes where the main vortex was smaller and less likely to disrupt the growing fibrin clot. Furthermore, smaller vortex modes were seen to support fibrin clot propagation across geometries. The model clearly demonstrates how the growing fibrin clot alters vortical structures within the aneurysm sac and how this changing flow, in turn, shapes the growing fibrin clot.

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

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

Model reduction of coagulation cascade based on genetic algorithm

Fibrin is an important product of the coagulation cascade, and plays an eminent role in platelet stabilization. Since coagulation cascade models typically involve the reaction kinetics of dozens of proteins, which will incur burdensome computational costs when coupled to blood flow in complex geometries, researchers often ignore this process when constructing thrombosis models. However, previous studies have shown that fundamental aspects of coagulation can be reproduced with simpler models, which motivated us to obtain a reduced-order model of fibrin generation through a systematic approach. Therefore, we introduced a semi-automatic framework to perform model-reduction of cascade reactions in this study, which consisted of two processes. Specifically, the retained protein species and cascade reactions were determined based on published studies and simulation results from the full cascade model, while the optimal reaction rates for the new cascade network were determined using a genetic algorithm. The framework has been applied to a 19-species coagulation model that triggers fibrin generation in internal fields via reactive boundaries, and a 10-species reduced-order model was obtained to reproduce the kinetics of fibrinogenesis in the full cascade model at different boundary tissue factor concentrations. This reduced-order model of fibrinogenesis would be valuable for thrombosis modeling that considers both the coagulation cascade and platelet activity. Furthermore, the framework proposed herein can also be applied to the reductions of other cascade reaction models.

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.

Effect of Pulsatility on the Transport of Thrombin in an Idealized Cerebral Aneurysm Geometry

Computational models of cerebral aneurysm thrombosis are designed for use in research and clinical applications. A steady flow assumption is applied in many of these models. To explore the accuracy of this assumption a pulsatile-flow thrombin-transport computational fluid dynamics (CFD) model, which uses a symmetrical idealized aneurysm geometry, was developed. First, a steady-flow computational model was developed and validated using data from an in vitro experiment, based on particle image velocimetry (PIV). The experimental data revealed an asymmetric flow pattern in the aneurysm. The validated computational model was subsequently altered to incorporate pulsatility, by applying a data-derived flow function at the inlet boundary. For both the steady and pulsatile computational models, a scalar function simulating thrombin generation was applied at the aneurysm wall. To determine the influence of pulsatility on thrombin transport, the outputs of the steady model were compared to the outputs of the pulsatile model. The comparison revealed that in the pulsatile case, an average of 10.2% less thrombin accumulates within the aneurysm than the steady case for any given time, due to periodic losses of a significant amount of thrombin-concentrated blood from the aneurysm into the parent vessel’s bloodstream. These findings demonstrate that pulsatility may change clotting outcomes in cerebral aneurysms.

An integrated fluid-chemical model toward modeling the thrombus formation in an idealized model of aortic dissection

Type B aortic dissection is a major aortic catastrophe that can be acutely complicated by rapid expansion, rupture, and malperfusion syndromes. The separation of the intima from aortic walls will form a second blood-filled lumen defined as "false lumen (FL)", where the thrombus is more likely to form due to the local stasis hemodynamic conditions. Complete thrombosis of FL is associated with a beneficial outcome while patency and partial thrombosis will lead to later complications. However, the thrombosis mechanism is still unclear and little is known about the impact of chemical species transported by blood flow on this process. The proteins involved in the coagulation cascade (CC) may play an important role in the process of thrombosis, especially in the activation and stabilization of platelets. Based on this hypothesis, a reduced-order fluid-chemical model was established to simulate CC in an aortic dissection phantom with two tears. A high level of fibrin is continuously observed at the top of the FL and some time-varying areas between two tears, indicating a high likelihood of thrombus formation there. This finding is consistent with the clinical observation. The time evolution of coagulation factors is greatly affected by local hemodynamics, especially in the high disturbance zone where the evolution has characteristics of periodic changes consistent with the flow field. The ability of the proposed model to reproduce the CC response provides a potential application to integrate with a model that can simulate platelet activities, forming a biochemical-based model which would help unveil the mechanisms of thrombosis in FL and the clinical decision of appropriate treatment.

Thrombin-Fibrinogen In Vitro Flow Model of Thrombus Growth in Cerebral Aneurysms

Cerebral aneurysms are balloon-like structures that develop on weakened areas of cerebral artery walls, with a significant risk of rupture. Thrombi formation is closely associated with cerebral aneurysms and has been observed both before and after intervention, leading to a wide variability of outcomes in patients with the condition. The attempt to manage the outcomes has led to the development of various computational models of cerebral aneurysm thrombosis. In the current study, we developed a simplified thrombin–fibrinogen flow system, based on commercially available purified human-derived plasma proteins, which enables thrombus growth and tracking in an idealized cerebral aneurysm geometry. A three-dimensional printed geometry of an idealized cerebral aneurysm and parent vessel configuration was developed. An unexpected outcome was that this phantom-based flow model allowed us to track clot growth over a period of time, by using optical imaging to record the progression of the growing clot into the flow field. Image processing techniques were subsequently used to extract important quantitative metrics from the imaging dataset, such as end point intracranial thrombus volume. The model clearly demonstrates that clot formation, in cerebral aneurysms, is a complex interplay between mechanics and biochemistry. This system is beneficial for verifying computational models of cerebral aneurysm thrombosis, particularly those focusing on initial angiographic occlusion outcomes, and will also assist manufacturers in optimizing interventional device designs.

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.

Microfluidic and computational study of structural properties and resistance to flow of blood clots under arterial shear

The ability of a blood clot to modulate blood flow is determined by the clot’s resistance, which depends on its structural features. For a flow with arterial shear, we investigated the characteristic patterns relating to clot shape, size, and composition on the one hand, and its viscous resistance, intraclot axial flow velocity, and shear distributions on the other. We used microfluidic technology to measure the kinetics of platelet, thrombin, and fibrin accumulation at a thrombogenic surface coated with collagen and tissue factor (TF), the key clot-formation trigger. We subsequently utilized the obtained data to perform additional calibration and validation of a detailed computational fluid dynamics model of spatial clot growth under flow. We then ran model simulations to gain insights into the resistance of clots formed under our experimental conditions. We found that increased thrombogenic surface length and TF surface density enhanced the bulk thrombin and fibrin generation in a nonadditive, synergistic way. The height of the platelet deposition domain—and, therefore, clot occlusivity—was rather robust to thrombogenic surface length and TF density variations, but consistently increased with time. Clot viscous resistance was non-uniform and tended to be higher in the fibrin-rich, inner “core” region of the clot. Interestingly, despite intraclot structure and viscous resistance variations, intraclot flow velocity variations were minor compared to the abrupt decrease in flow velocity around the platelet deposition region. Our results shed new light on the connection between the structure of clots under arterial shear and spatiotemporal variations in their resistance to flow.

A Virtual Comparison of the eCLIPs Device and Conventional Flow-Diverters as Treatment for Cerebral Bifurcation Aneurysms

PURPOSE

Effective, consistent, and complication-free treatment of cerebral bifurcation aneurysms remains elusive despite a pressing need, with the majority of lesions presenting in such locations. Current treatment options focus either on aneurysm coil retention, supported by a stent-like device positioned in the parent vessel lumen, or intrasaccular devices that disrupt flow within the aneurysm dome. A third alternative, i.e., the use of conventional (intraluminal) flow-diverters to treat such bifurcation aneurysms raises the problem that at least one daughter vessel needs to be jailed in such a deployment. The eCLIPs is a stent-like device that offers the possibility of flow-diversion at the aneurysm neck, without the drawbacks of daughter vessel occlusion or those of intrasaccular deployment.

METHODS

In this study the eCLIPs device was virtually deployed in five cerebral bifurcation aneurysms and compared with a conventional tubular flow-diverter device. Computational fluid dynamics (CFD) simulations of the aneurysm haemodynamic environment pre- and post-implantation were conducted, and focussed on metrics associated with successful aneurysm occlusion. Absolute and relative reductions in aneurysm inflow rate (Q) and time-averaged wall shear stress (TAWSS) were recorded.

RESULTS

The eCLIPs device was found to perform in a similar qualitative fashion to tubular flow-diverters, with overall reduction of metrics being somewhat more modest however, when compared to such devices. Aneurysm inflow reduction and TAWSS reduction were typically 10-20% lower for the eCLIPs, when compared to a generic flow diverter (FD_BRAIDED) similar to devices currently in clinical use. The eCLIPs was less effective at diffusing inflow jets and at reducing the overall velocity of the flow, when compared to these devices. This result is likely due to the larger device pore size in the eCLIPs. Notably, it was found that the eCLIPs provided approximately equal resistance to flow entering and exiting the aneurysm, which was not true for the FD_BRAIDED device, where high-speed concentrations of outflow were seen at the aneurysm neck along with local TAWSS elevation. The clinical implications of such behaviour are not examined in detail here but could be significant.

CONCLUSIONS

Our findings indicate that the eCLIPs device acts as a flow-diverter for bifurcation aneurysms, with somewhat diminished occlusion properties comparing to tubular flow diverters but without the jailing and diminished flow evident in a daughter vessel associated with use of conventional devices.

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.

Introducing the pro-coagulant contact system in the numerical assessment of device-related thrombosis

Thrombosis is a major concern in blood-coated medical devices. Contact activation, which is the initial part of the coagulation cascade in device-related thrombosis, is not considered in current thrombus formation models. In the present study, pro-coagulant reactions including the contact activation system are coupled with a fluid solver in order to evaluate the potential of the contact system to initiate thrombin production. The biochemical/fluid model is applied to a backward-facing step configuration, a flow configuration that frequently appears in medical devices. In contrast to the in vivo thrombosis models in which a specific thrombotic zone (injury region) is set a priori by the user to initiate the coagulation reaction, a reactive surface boundary condition is applied to the whole device wall. Simulation results show large thrombin concentration in regions related to recirculation zones without the need of an a priori knowledge of the thrombus location. The numerical results align well with the regions prone to thrombosis observed in experimental results reported in the literature. This approach could complement thrombus formation models that take into account platelet activity and thrombus growth to optimize a wide range of medical devices.

Effect of different stent configurations using Lattice Boltzmann method and particles image velocimetry on artery bifurcation aneurysm problem

Proper design of stent for application at specific aneurysm effect arteries could help to reduce the issues with thrombosis and aneurysm. In this paper, four types of stent configuration namely half-Y (6 mm), half-Y (4 mm), cross-bar, and full-Y configuration will implanted on real 3D artery bifurcation aneurysm effected arteries. Comparisons were then conducted based on the flow patterns after stent placement using both LBM-based solver and PIV experimental findings. According to the data obtained from all 4 stent designs, the flow profiles and the computed velocity from both methods were in agreement with each other. Both methods found that half-Y (6 mm) stent configuration is by far the best configuration in reducing the blood velocity at the vicinity of the aneurysm sac. The analysis also show that the half-Y (6 mm) stent configuration recorded the highest percentage of velocity reduction and managed to substantially reduce the pressure at the bifurcation region. This high flow velocity reduction through the use of half-Y stent could consequently promote the formation of thrombus thereby reducing the risk of rupture in the aneurysm sac.

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.

Virtual endovascular treatment of intracranial aneurysms: models and uncertainty

Virtual endovascular treatment models (VETMs) have been developed with the view to aid interventional neuroradiologists and neurosurgeons to pre-operatively analyze the comparative efficacy and safety of endovascular treatments for intracranial aneurysms. Based on the current state of VETMs in aneurysm rupture risk stratification and in patient-specific prediction of treatment outcomes, we argue there is a need to go beyond personalized biomechanical flow modeling assuming deterministic parameters and error-free measurements. The mechanobiological effects associated with blood clot formation are important factors in therapeutic decision making and models of post-treatment intra-aneurysmal biology and biochemistry should be linked to the purely hemodynamic models to improve the predictive power of current VETMs. The influence of model and parameter uncertainties associated to each component of a VETM is, where feasible, quantified via a random-effects meta-analysis of the literature. This allows estimating the pooled effect size of these uncertainties on aneurysmal wall shear stress. From such meta-analyses, two main sources of uncertainty emerge where research efforts have so far been limited: (1) vascular wall distensibility, and (2) intra/intersubject systemic flow variations. In the future, we suggest that current deterministic computational simulations need to be extended with strategies for uncertainty mitigation, uncertainty exploration, and sensitivity reduction techniques. WIREs Syst Biol Med 2017, 9:e1385. doi: 10.1002/wsbm.1385 For further resources related to this article, please visit the WIREs website.

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.