A new compressible hyperelastic model for the multi-axial deformation of blood clot occlusions in vessels

Biomechanical profiling of in vitro blood clots: sensitivity to sex, age, and blood composition in a healthy adult population

Blood clots' mechanical properties are important in both their physiological role and in the initiation and progression of thromboembolic diseases. Because studying blood clot properties in vivo is difficult, many prior studies have investigated the properties of in vitro clots instead. However, much remains to be understood about in vitro clots, especially those derived from human blood. For example, the association between subject-specific factors and clot mechanical properties is currently unknown. Our objective is to fill this knowledge gap and study the sensitivity of in vitro blood clots to subject-specific factors, including sex, age, and blood composition. We drew blood from healthy adults aged 19-46, coagulated clots into mechanical test specimens, and characterized their properties. Specifically, we quantified clot stiffness, fracture toughness, contractility, and hysteresis. We then quantified the relative dependence of those properties on subject-specific factors, including sex, age, and blood composition. We found that there is significant variation in clot properties within healthy subjects. Clots from female subjects' blood are stiffer, more resistant to fracture, and show more hysteresis compared to clots from male subjects. However, we found no association between clot properties and age and only a weak association with clot composition, e.g., hematocrit. Finally, even together, sex, age, and blood composition only moderately explain the observed variability in clot mechanical properties. Our work therefore suggests that in vitro clots may capture relevant information not reflected in standard clinical data. Future studies should investigate in vitro clots' potential as biomarkers for thrombotic risk and treatment response.

Rupture mechanics of blood clot fibrin fibers: A coarse-grained model study

Thrombosis, when occurring undesirably, disrupts normal blood flow and poses significant medical challenges. As the skeleton of blood clots, fibrin fibers play a vital role in the formation and fragmentation of blood clots. Thus, studying the deformation and fracture characteristics of fibrin fiber networks is the key factor to solve a series of health problems caused by thrombosis. This study employs a coarse-grained model of fibrin fibers to investigate the rupture dynamics of fibrin fiber networks. We propose a new method for generating biomimetic fibrin fiber networks to simulate their spatial geometry in blood clots. We examine the mechanical characteristics and rupture behaviors of fibrin fiber networks under various conditions, including fiber junction density, fiber tortuosity, fiber strength, and the strain limit of single fiber rupture in both tension and simple shear cases. Our findings indicate that the stress-strain relationship of the fibrin fiber network follows a similar pattern to that of individual fibers, characterized by a shortened entropy stretching phase and an extended transition phase. Fiber junction density, fiber strength, and single fiber rupture limit predominantly influence the stress of the network, while fiber tortuosity governs the strain behavior. The availability of more fibers in shear cases to bear the load results in delayed rupture compared to tension cases. With consideration of different factors of fibrin fibers in networks, this work provides a more realistic description of the mechanical deformation process in fibrin fiber networks, offering new insights into their rupture and failure mechanisms. These findings could inspire novel approaches and methodologies for understanding the fracture of fibrin networks during a surgical thrombectomy.

Impact of thrombus composition on virtual thrombectomy procedures using human clot analogues mechanical data

Endovascular thrombectomy (EVT) aims at restoring blood flow in case of acute ischemic stroke by removing the thrombus occluding a large cerebral artery. During the procedure with stent-retriever, the thrombus is captured within the device, which is then retrieved, subjecting the thrombus to several forces, potentially leading to its fragmentation. In silico studies, along with mechanical characterisation of thrombi, can enhance our understanding of the EVT, helping the development of new devices and interventional strategies. Our group previously validated a numerical approach to study EVT able to account for thrombus fragmentation. In this study, the same methodology was employed to explore the applicability of the chosen failure criterion to EVT simulations and the impact of thrombus composition on the outcome of the in silico procedure. For the first time, human clot analogues experimental data were applied to this methodology. Clot analogues of three different compositions were tested, and a material model incorporating failure was calibrated, followed by a verification analysis. Finally, the calibrated material model was used to perform EVT simulations, combining the three tested thrombus compositions with three different stent retriever models. The experimental tests confirmed a compression-tension asymmetry in the stress-strain curves, showing decreasing stiffness with increasing the red blood cell (RBC) content. Applying the resulting material models to EVT simulations demonstrated: (i) the dependency of the failure criterion on the thrombus mesh size, (ii) a greater tendency for RBC-rich thrombi to fragment, and (iii) increased difficulty in retrieving RBC-poor thrombi compared to RBC-rich thrombi.

Collateral flow and pulsatility during large vessel occlusions: insights from a quantitative in vitro study

Acute ischemic stroke caused by large vessel occlusions is being increasingly treated with neurovascular interventions. The hemodynamics within the collateral system of the circle of Willis (CoW) hemodynamics play a fundamental role in therapy success. However, transient in vivo data on pathological collateral flow during large vessel occlusions are not available. Moreover, there are no flow models that accurately simulate the hemodynamic conditions in the CoW during large vessel occlusions. We used a circulatory loop to generate highly reproducible cerebrovascular-like flows and pressures and used non-invasive flow visualization and high-resolution flow and pressure measurements to acquire detailed, time-dependent hemodynamics inside an anatomical phantom of the CoW. After calibrating a physiological reference case, we induced occlusions in the 1. middle cerebral artery, 2. terminal carotid artery, and 3. basilar artery; and measured the left posterior communicating artery flow. Mean arterial pressure and pulse pressure remained unchanged in the different occlusion cases compared to the physiological reference case, while total cerebral flow decreased by up to 19%. In all three occlusion cases, reversed flow was found in the left posterior communicating artery compared to the reference case with different flow magnitudes and pulsatility index values. The experimental results were compared with clinical findings, demonstrating the capability of this realistic cerebrovascular flow setup. This novel cerebrovascular flow setup opens the possibility for investigating different topics of neurovascular interventions under various clinical conditions in controlled preclinical laboratory studies.

Mechanical properties of clot made from human and bovine whole blood differ significantly

Thromboembolism - that is, clot formation and the subsequent fragmentation of clot - is a leading cause of death worldwide. Clots' mechanical properties are critical determinants of both the embolization process and the pathophysiological consequences thereof. Thus, understanding and quantifying the mechanical properties of clots is important to our ability to treat and prevent thromboembolic disease. However, assessing these properties from in vivo clots is experimentally challenging. Therefore, we and others have turned to studying in vitro clot mimics instead. Unfortunately, there are significant discrepancies in the reported properties of these clot mimics, which have been hypothesized to arise from differences in experimental techniques and blood sources. The goal of our current work is therefore to compare the mechanical behavior of clots made from the two most common sources, human and bovine blood, using the same experimental techniques. To this end, we tested clots under pure shear with and without initial cracks, under cyclic loading, and under stress relaxation. Based on these data, we computed and compared stiffness, strength, work-to-rupture, fracture toughness, relaxation time constants, and prestrain. While clots from both sources behaved qualitatively similarly, they differed quantitatively in almost every metric. We also correlated each mechanical metric to measures of blood composition. Thereby, we traced this inter-species variability in clot mechanics back to significant differences in hematocrit, but not platelet count. Thus, our work suggests that the results of past studies that have used bovine blood to make in vitro mimics - without adjusting blood composition - should be interpreted carefully. Future studies about the mechanical properties of blood clots should focus on human blood alone.

In silico analysis of embolism in cerebral arteries using fluid-structure interaction method

Ischemic stroke, particularly embolic stroke, stands as a significant global contributor to mortality and long-term disabilities. This paper presents a comprehensive simulation of emboli motion through the middle cerebral artery (MCA), a prevalent site for embolic stroke. Our patient-specific computational model integrates major branches of the middle cerebral artery reconstructed from magnetic resonance angiography images, pulsatile flow dynamics, and emboli of varying geometries, sizes, and material properties. The fluid-structure interactions method is employed to simulate deformable emboli motion through the middle cerebral artery, allowing observation of hemodynamic changes in artery branches upon embolus entry. We investigated the impact of embolus presence on shear stress magnitude on artery walls, analyzed the effects of embolus material properties and geometries on embolus trajectory and motion dynamics within the middle cerebral artery. Additionally, we evaluated the non-Newtonian behavior of blood, comparing it with Newtonian blood behavior. Our findings highlight that embolus geometry significantly influences trajectory, motion patterns, and hemodynamics within middle cerebral artery branches. Emboli with visco-hyperelastic material properties experienced higher stresses upon collision with artery walls compared to those with hyperelastic properties. Furthermore, considering blood as a non-Newtonian fluid had notable effects on emboli stresses and trajectories within the artery, particularly during collisions. Notably, the maximum von Mises stress experienced in our study was 21.83 kPa, suggesting a very low probability of emboli breaking during movement, impact, and after coming to a stop. However, in certain situations, the magnitude of shear stress on them exceeded 1 kPa, increasing the likelihood of cracking and disintegration. These results serve as an initial step in anticipating critical clinical conditions arising from arterial embolism in the middle cerebral artery. They provide insights into the biomechanical parameters influencing embolism, contributing to improved clinical decision-making for stroke management.

Coronary hemodynamic simulation study

In this paper, a two-way fluid-structure coupling model is developed to simulate and analyze the hemodynamic process based on dynamic coronary angiography, and examine the influence of different hemodynamic parameters on coronary arteries in typical coronary stenosis lesions. Using the measured FFR pressure data of a patient, the pressure-time function curve is fitted to ensure the accuracy of the boundary conditions. The average error of the simulation pressure results compared to the test data is 6.74%. In addition, the results related to blood flow, pressure contour and wall shear stress contour in a typical cardiac cycle are obtained by simulation analysis. These results are found to be in good agreement with the laws of the real cardiac cycle, which verifies the rationality of the simulation. In conclusion, based on the modeling and hemodynamic simulation analysis process of dynamic coronary angiography, this paper proposes a method to assist the analysis and evaluation of coronary hemodynamic and functional parameters, which has certain practical significance.

Numerical Study of a Thrombus Migration Risk in Aneurysm After Coil Embolization in Patient Cases: FSI Modelling

PURPOSE

There are still many challenges for modelling a thrombus migration process in aneurysms. The main novelty of the present research lies in the modelling of aneurysm clot migration process in a realistic cerebral aneurysm, and the analysis of forces suffered by clots inside an aneurysm, through transient FSI simulations.

METHODS

The blood flow has been modelled using a Womersley velocity profile, and following the Carreau viscosity model. Hyperelastic Ogden model has been used for clot and isotropic linear elastic model for the artery walls. The FSI coupled model was implemented in ANSYS® software. The hemodynamic forces suffered by the clot have been quantified using eight different clot sizes and positions inside a real aneurysm.

RESULTS

The obtained results have shown that it is almost impossible for clots adjacent to aneurysm walls, to leave the aneurysm. Nevertheless, in clots positioned in the centre of the aneurysm, there is a real risk of clot migration. The risk of migration of a typical post-coiling intervention clot in an aneurysm, in contact with the wall and occupying a significant percentage of its volume is very low in the case studied, even in the presence of abnormally intense events, associated with sneezes or impacts.

CONCLUSIONS

The proposed methodology allows evaluating the clot migration risk, vital for evaluating the progress after endovascular interventions, it is a step forward in the personalized medicine, patient follow-up, and helping the medical team deciding the optimal treatment.

The influence of blood composition and loading frequency on the behavior of embolus analogs

In cases of acute ischemic stroke (AIS), mechanical thrombectomy (MT) can be used to directly remove lodged thromboemboli. Despite improvements in patient outcomes, one of the key factors affecting MT success is the mechanical properties of the occlusive thrombus. Therefore, the goal of this study was to investigate the viscoelastic properties of embolus analogs (EAs) and determine the influence of EA hematocrit and loading frequency. Bovine blood EAs were created over a range of physiological hematocrits (0-60%) and cyclic uniaxial compression testing was performed at three loading frequencies to mimic in vivo loading conditions, followed by stress-relaxation testing. It was found that EAs exhibited behaviors typical of hyper-viscoelastic materials and that EA hematocrit played a large role in both EA stiffness and relaxation, with both parameters decreasing as hematocrit increased from 0 to 60%. The viscoelastic behavior of the EAs was also affected by the frequency at which they were loaded, with significant increases in peak stresses between the 0.5 and 2 Hz loaded EAs. Lower hematocrit EAs had very dense fibrin networks while the higher hematocrit EAs consisted of closely packed RBCs with little fibrin present. These results suggest that fibrin contributes to EA stiffness and relaxation behaviors while RBCs play a role in decreasing the overall viscous response and strain-rate dependency. An Ogden hyperelastic model was found to best reproduce the EA loading data while a 3-term Prony series was fit to the stress relaxation data. A hyper-viscoelastic modeling framework was then implemented combining the loading and stress-relaxation fits and the results could match the full cyclic loading data for EAs of varying hematocrit and loading frequency. The results of the experimental mechanical characterization and hyper-viscoelastic curve fitting can be incorporated in future modeling efforts to optimize mechanical thrombectomy for AIS patients.

An introduction to the Ogden model in biomechanics: benefits, implementation tools and limitations

Constitutive models are important to biomechanics for two key reasons. First, constitutive modelling is an essential component of characterizing tissues' mechanical properties for informing theoretical and computational models of biomechanical systems. Second, constitutive models can be used as a theoretical framework for extracting and comparing key quantities of interest from material characterization experiments. Over the past five decades, the Ogden model has emerged as a popular constitutive model in soft tissue biomechanics with relevance to both informing theoretical and computational models and to comparing material characterization experiments. The goal of this short review is threefold. First, we will discuss the broad relevance of the Ogden model to soft tissue biomechanics and the general characteristics of soft tissues that are suitable for approximating with the Ogden model. Second, we will highlight exemplary uses of the Ogden model in brain tissue, blood clot and other tissues. Finally, we offer a tutorial on fitting the one-term Ogden model to pure shear experimental data via both an analytical approximation of homogeneous deformation and a finite-element model of the tissue domain. Overall, we anticipate that this short review will serve as a practical introduction to the use of the Ogden model in biomechanics. This article is part of the theme issue 'The Ogden model of rubber mechanics: Fifty years of impact on nonlinear elasticity'.

An introduction to the Ogden model in biomechanics: benefits, implementation tools and limitations

Constitutive models are important to biomechanics for two key reasons. First, constitutive modelling is an essential component of characterizing tissues' mechanical properties for informing theoretical and computational models of biomechanical systems. Second, constitutive models can be used as a theoretical framework for extracting and comparing key quantities of interest from material characterization experiments. Over the past five decades, the Ogden model has emerged as a popular constitutive model in soft tissue biomechanics with relevance to both informing theoretical and computational models and to comparing material characterization experiments. The goal of this short review is threefold. First, we will discuss the broad relevance of the Ogden model to soft tissue biomechanics and the general characteristics of soft tissues that are suitable for approximating with the Ogden model. Second, we will highlight exemplary uses of the Ogden model in brain tissue, blood clot and other tissues. Finally, we offer a tutorial on fitting the one-term Ogden model to pure shear experimental data via both an analytical approximation of homogeneous deformation and a finite-element model of the tissue domain. Overall, we anticipate that this short review will serve as a practical introduction to the use of the Ogden model in biomechanics. This article is part of the theme issue 'The Ogden model of rubber mechanics: Fifty years of impact on nonlinear elasticity'.

Combined stent-retriever and aspiration intra-arterial thrombectomy performance for fragmentable blood clots: A proof-of-concept computational study

Mechanical thrombectomy (MT) treatment of acute ischemic stroke (AIS) patients typically involves use of stent retrievers or aspiration catheters alone or in combination. For in silico trials of AIS patients, it is crucial to incorporate the possibility of thrombus fragmentation during the intervention. This study focuses on two aspects of the thrombectomy simulation: i) Thrombus fragmentation on the basis of a failure model calibrated with experimental tests on clot analogs; ii) the combined stent-retriever and aspiration catheter MT procedure is modeled by adding both the proximal balloon guide catheter and the distal access catheter. The adopted failure criterion is based on maximum principal stress threshold value. If elements of the thrombus exceed this criterion during the retrieval simulation, then they are deleted from the calculation. Comparison with in-vitro tests indicates that the simulation correctly reproduces the procedures predicting thrombus fragmentation in the case of red blood cells rich thrombi, whereas non-fragmentation is predicted for fibrin-rich thrombi. Modeling of balloon guide catheter prevents clot fragments' embolization to further distal territories during MT procedure.

Preclinical modeling of mechanical thrombectomy

Mechanical thrombectomy to treat large vessel occlusions (LVO) causing a stroke is one of the most effective treatments in medicine, with a number needed to treat to improve clinical outcomes as low as 2.6. As the name implies, it is a mechanical solution to a blocked artery and modeling these mechanics preclinically for device design, regulatory clearance and high-fidelity physician training made clinical applications possible. In vitro simulation of LVO is extensively used to characterize device performance in representative vascular anatomies with physiologically accurate hemodynamics. Embolus analogues, validated against clots extracted from patients, provide a realistic simulated use experience. In vitro experimentation produces quantitative results such as particle analysis of distal emboli generated during the procedure, as well as pressure and flow throughout the experiment. Animal modeling, used mostly for regulatory review, allows estimation of device safety. Other than one recent development, nearly all animal modeling does not incorporate the desired target organ, the brain, but rather is performed in the extracranial circulation. Computational modeling of the procedure remains at the earliest stages but represents an enormous opportunity to rapidly characterize and iterate new thrombectomy concepts as well as optimize procedure workflow. No preclinical model is a perfect surrogate; however, models available can answer important questions during device development and have to date been successful in delivering efficacious and safe devices producing excellent clinical outcomes. This review reflects on the developments of preclinical modeling of mechanical thrombectomy with particular focus on clinical translation, as well as articulate existing gaps requiring additional research.

Realistic computer modelling of stent retriever thrombectomy: a hybrid finite-element analysis-smoothed particle hydrodynamics model

Stent retriever thrombectomy is a pre-eminent treatment modality for large vessel ischaemic stroke. Simulation of thrombectomy could help understand stent and clot mechanics in failed cases and provide a digital testbed for the development of new, safer devices. Here, we present a novel, in silico thrombectomy method using a hybrid finite-element analysis (FEA) and smoothed particle hydrodynamics (SPH). Inspired by its biological structure and components, the blood clot was modelled with the hybrid FEA-SPH method. The Solitaire self-expanding stent was parametrically reconstructed from micro-CT imaging and was modelled as three-dimensional finite beam elements. Our simulation encompassed all steps of mechanical thrombectomy, including stent packaging, delivery and self-expansion into the clot, and clot extraction. To test the feasibility of our method, we simulated clot extraction in simple straight vessels. This was compared against in vitro thrombectomies using the same stent, vessel geometry, and clot size and composition. Comparisons with benchtop tests indicated that our model was able to accurately simulate clot deflection and penetration of stent wires into the clot, the relative movement of the clot and stent during extraction, and clot fragmentation/embolus formation. In this study, we demonstrated that coupling FEA and SPH techniques could realistically model stent retriever thrombectomy.

A new constitutive model for permanent deformation of blood clots with application to simulation of aspiration thrombectomy

As a first line option in the treatment of acute ischemic stroke (AIS), direct aspiration is a fast and effective technique with promising outcomes. In silico models are widely used for design and preclinical assessment of new developed devices and therapeutic methods. Accurate modelling of the mechanical behaviour of blood clot is a key factor in the design and simulation of aspiration devices. In this study we develop a new constitutive model which incorporates the unrecoverable plastic deformation of clots. The model is developed based on the deformation-induced microstructural changes in fibrin network, including the formation and dissociation of the cross-links between fibrin fibres. The model is calibrated using previously reported experimentally measured permanent clot deformation following uniaxial stretching. The calibrated plasticity model is then used to simulate aspiration thrombectomy. Results reveal that inclusion of permanent plastic deformation results in ∼ 15 % increase in clot aspiration length at an applied aspiration pressure of 100 mmHg. The constitutive law developed in this study provides a basis for improved design and evaluation of novel aspiration catheters leading to increased first-pass revascularization rate.

In vitro and in silico modeling of endovascular stroke treatments for acute ischemic stroke

Acute ischemic stroke occurs when a thrombus obstructs a cerebral artery, leading to sub-optimal blood perfusion to brain tissue. A recently developed, preventive treatment is the endovascular stroke treatment (EVT), which is a minimally invasive procedure, involving the use of stent-retrievers and/or aspiration catheters. Despite its increasing use, many critical factors of EVT are not well understood. In this respect, in vitro, and in silico studies have the great potential to help us deepen our understanding of the procedure, perform further device and procedural optimization, and help in clinical training. This review paper provides an overview of the previous in vitro and in silico evaluations of EVT treatments, with a special emphasis on the four main aspects of the adopted experimental and numerical set-ups: vessel, thrombus, device, and procedural settings.

Blood clot fracture properties are dependent on red blood cell and fibrin content

Thrombus fragmentation during endovascular stroke treatment, such as mechanical thrombectomy, leads to downstream emboli, resulting in poor clinical outcomes. Clinical studies suggest that fragmentation risk is dependent on clot composition. This current study presents the first experimental characterization of the composition-dependent fracture properties of blood clots, in addition to the development of a predictive model for blood clot fragmentation. A bespoke experimental test-rig and compact tension specimen fabrication has been developed to measure fracture toughness of thrombus material. Fracture tests are performed on three physiologically relevant clot compositions: a high-fibrin clot made from a 5% haematocrit (H) blood mixture, a medium-fibrin clot made form a 20% H blood mixture, a low-fibrin clot made from a 40% H blood mixture. Fracture toughness is observed to significantly increase with increasing fibrin content, i.e. red blood cell-rich clots are more prone to tear during loading compared to the fibrin-rich clots. Results also reveal that the mechanical behaviour of clot analogues is significantly different in compression and tension. Finite element cohesive zone modelling of clot fracture experiments show that fibrin fibres become highly aligned in the direction perpendicular to crack propagation, providing a significant toughening mechanism. The results presented in this study provide the first characterization of the composition-dependent fracture behaviour of blood clots and are of key importance for development of next-generation thrombectomy devices and clinical strategies. STATEMENT OF SIGNIFICANCE: This study provides a characterisation of the composition-dependent fracture toughness of blood clots. This entails the development of novel experimental techniques for fabrication and testing of blood clot compact tension fracture specimens. The study also develops cohesive zone models of fracture initiation and propagation in blood clots. Results reveal that the fracture resistance of fibrin-rich clots is significantly higher than red blood cell rich clots. Simulations also reveal that stretching and realignment of the fibrin network should be included in blood clot material models in order to accurately replicate compression-tension asymmetry and fibrin enhanced fracture toughness. The results of this study have potentially important clinical implications in terms of clot fracture risk and secondary embolization during mechanical thrombectomy procedures.