Development of a computational model for macroscopic predictions of device-induced thrombosis

Exploring hemodynamic mechanisms and re-intervention strategies for partial false lumen thrombosis in Stanford type B aortic dissection after thoracic aortic endovascular repair

OBJECTIVES

False lumen (FL) thrombosis status for Stanford type B aortic dissection (TBAD) after thoracic endovascular aortic repair (TEVAR) is critical for evaluating aortic remodeling and long-term prognosis. This study aimed to monitor the morphology evolution of partial FL thrombosis (PFLT) and its hemodynamic conditions through an innovative approach, providing a re-intervention strategy from both morphologic and hemodynamic perspectives.

METHODS

Three-dimensional geometries are extracted from a five-year follow-up of CTA images for TBAD after TEVAR. The morphology and hemodynamics of PFLT are comprehensively analyzed based on patient-specific reconstructions and computational fluid dynamics (CFD). The impact of various strategies treating risk factors of PFLT, including proximal entry closure, left renal artery stenting, or accessory renal artery embolism on hemodynamics is assessed.

RESULTS

The introduced morphologic approaches appropriately reflected the evolution of PFLT. Gradual dilation of FL (surface area from 82.63cm2 to 98.84cm2, volume from 45.12 mL to 63.40 mL, increase in distal tear (from 3.72 cm to 4.32 cm), and fluctuation of thrombosis-blood lumen boundary are observed. For further surgical preparation in the absence of unanimously recognized re-intervention indicators, velocity and wall shear stress distributions reveal different simulated re-interventions have distinctly suppressive effects on hemodynamic conditions within the PFLT, providing valuable insights for further surgical preparation.

CONCLUSIONS

The present study demonstrates a re-intervention strategy for PFLT in TBAD patients after TEVAR utilizing morphologic and hemodynamic analyses. Acknowledging the deterioration of PFLT may result in adverse long-term outcomes, this strategy might offer an alternative approach for clinical monitoring and management of related patients.

Improving hydraulic performance of the left atrial assist device using computational fluid dynamics

BACKGROUND

The left atrial assist device (LAAD) is a novel continuous-flow pump designed to treat patients with heart failure with preserved ejection fraction, a growing type of heart failure, but with limited device-treatment options. The LAAD is implanted in the mitral plane and pumps blood from the left atrium into the left ventricle. The purpose of this study was to refine the initial design of the LAAD, using results from computational fluid dynamics (CFD) analyses to inform changes that could improve hydraulic performance and flow patterns within the LAAD.

METHODS

The initial design and three variations were simulated, exploring changes to the primary impeller blades, the housing shape, and the number, size, and curvature of the diffuser vanes. Several pump rotational speeds and flow rates spanning the intended range of use were modeled.

RESULTS

Guided by the insight gained from each design iteration, the final design incorporated impeller blades with improved alignment relative to the incoming flow and wider, more curved diffuser vanes that better aligned with the approaching flow from the volute. These design adjustments reduced flow separation within the impeller and diffuser regions. In vitro testing confirmed the CFD predicted improvement in the hydraulic performance of the revised LAAD flow path design.

CONCLUSIONS

The CFD results from this study provided insight into the key pump design-related parameters that can be adjusted to improve the LAAD's hydraulic performance and internal flow patterns. This work also provided a foundation for future studies assessing the LAAD's biocompatibility under clinical conditions.

A new method for scaling inlet flow waveform in hemodynamic analysis of aortic dissection

Computational fluid dynamics (CFD) simulations have shown great potentials in cardiovascular disease diagnosis and postoperative assessment. Patient-specific and well-tuned boundary conditions are key to obtaining accurate and reliable hemodynamic results. However, CFD simulations are usually performed under non-patient-specific flow conditions due to the absence of in vivo flow and pressure measurements. This study proposes a new method to overcome this challenge by tuning inlet boundary conditions using data extracted from electrocardiogram (ECG). Five patient-specific geometric models of type B aortic dissection were reconstructed from computed tomography (CT) images. Other available data included stoke volume (SV), ECG, and 4D-flow magnetic resonance imaging (MRI). ECG waveforms were processed to extract patient-specific systole to diastole ratio (SDR). Inlet boundary conditions were defined based on a generic aortic flow waveform tuned using (1) SV only, and (2) with ECG and SV (ECG + SV). 4D-flow MRI derived inlet boundary conditions were also used in patient-specific simulations to provide the gold standard for comparison and validation. Simulations using inlet flow waveform tuned with ECG + SV not only successfully reproduced flow distributions in the descending aorta but also provided accurate prediction of time-averaged wall shear stress (TAWSS) in the primary entry tear (PET) and abdominal regions, as well as maximum pressure difference, ∆Pmax, from the aortic root to the distal false lumen. Compared with simulations with inlet waveform tuned with SV alone, using ECG + SV in the tuning method significantly reduced the error in false lumen ejection fraction at the PET (from 149.1% to 6.2%), reduced errors in TAWSS at the PET (from 54.1% to 5.7%) and in the abdominal region (from 61.3% to 11.1%), and improved ∆Pmax prediction (from 283.1% to 18.8%) However, neither of these inlet waveforms could be used for accurate prediction of TAWSS in the ascending aorta. This study demonstrates the importance of SDR in tailoring inlet flow waveforms for patient-specific hemodynamic simulations. A well-tuned flow waveform is essential for ensuring that the simulation results are patient-specific, thereby enhancing the confidence and fidelity of computational tools in future clinical applications.

Simulation of the effect of hemolysis on thrombosis in blood-contacting medical devices

Heart failure, broadly characterized by the gradual decline of the ability of the heart to maintain adequate blood flow throughout the body's vascular network of veins and arteries, is one of the leading causes of death worldwide. Mechanical Circulatory Support is one of the few available alternative interventions for late-stage heart failure with reduced ejection fraction. A ventricular assist device is surgically implanted and connected to the left and or right heart ventricles to provide additional bloodflow, off-loading the work required by the heart to maintain circulation. Modern mechanical circulatory support devices generate non-physiological flow conditions that can lead to the damage and rupture of blood cells (hemolysis), and the formation of blood clots (thrombosis), which pose severe health risks to the patient. It is essential to improve prediction tools for blood damage to reduce the risk of hemolysis and thrombosis. A simulation-based approach examines the interaction between hemolysis and thrombosis. Incompressible finite-volume computational fluid dynamics simulations are executed on an open-hub axial flow ventricular assist device. A continuum model of thrombosis and the intrinsic coagulation process is extended to include the effect of hemolysis. The model accounts for the effect of activation of platelets by shear stress, paracrine signaling, adhesion, and hemoglobin and ADP released during hemolysis. The effect of hemolysis with thrombosis is modelled by accounting for the hyper-adhesivity of von-Willebrand Factor on extracellular hemoglobin, and the increased rate of platelet activation induced by ADP release. Thrombosis is assessed at varying inflow rates and rotor speeds, and cases are executed where thrombosis is affected by ADP release and Hb-induced hyper-adhesivity. It is found that there is a non-negligible effect from hemolysis on thrombosis across a range of rotor speeds, and that hyperadhesivity plays a dominant role in thrombus formation in the presence of hemolysis.

Characterizing thrombus adhesion strength on common cardiovascular device materials

Thrombus formation in blood-contacting medical devices is a major concern in the medical device industry, limiting the clinical efficacy of these devices. Further, a locally formed clot within the device has the potential to detach from the surface, posing a risk of embolization. Clot embolization from blood-contacting cardiovascular devices can result in serious complications like acute ischemic stroke and myocardial infarction. Therefore, clot embolization associated with device-induced thrombosis can be life-threatening and requires an enhanced fundamental understanding of embolization characteristics to come up with advanced intervention strategies. Therefore, this work aims to investigate the adhesive characteristics of blood clots on common biocompatible materials used in various cardiovascular devices. This study focuses on characterizing the adhesion strength of blood clots on materials such as polytetrafluoroethylene (PTFE), polyurethane (PU), polyether ether ketone (PEEK), nitinol, and titanium, frequently used in medical devices. In addition, the effect of incubation time on clot adhesion is explored. Results from this work demonstrated strongest clot adhesion to titanium with 3 h of incubation resulting in 1.06 ± 0.20 kPa detachment stresses. The clot adhesion strength on titanium was 51.5% higher than PEEK, 35.9% higher than PTFE, 63.1% higher than PU, and 35.4% higher than nitinol. Further, adhesion strength increases with incubation time for all materials. The percentage increase in detachment stress over incubation time (ranging from 30 min to 3 h) for polymers ranged from at least 108.75% (PEEK), 140.74% (PU), to 151.61% (PTFE). Whereas, for metallic surfaces, the percentage rise ranged from 70.21% (nitinol) to 89.28% (titanium). Confocal fluorescence imaging of clot remnants on the material surfaces revealed a well-bounded platelet-fibrin network at the residual region, representing a comparatively higher adhesive region than the non-residual zone of the surface.

Optimization design of semi-open impeller based on thrombogenicity in a blood pump

BACKGROUND

Blood clots are composed of aggregated fibrin and platelets, and thrombosis is the body's natural response to repairing injured blood vessels or stopping bleeding. However, when this process is activated abnormally, such as in a mechanical blood pump, it can lead to excessive thrombus formation. Therefore, how to avoid or reduce the probability of thrombus formation is an important indicator of the stable operation of a blood pump.

METHODS

In this paper, Lagrangian particle tracking trajectories are simulated to study platelet transport in a blood pump. The design of the thrombus blood pump was optimized using an orthogonal design method based on three factors: inlet angle, outlet angle, and blade number. The effect of blood pump pressure, rotational speed, impeller outlet angle, inlet angle, and number of blades on thrombus formation was analysed using Fluent software. The thrombogenic potential was derived by analyzing the trajectory and flow parameters of platelet particles in the blood pump, as well as the statistical parameters of residence time and stress accumulation thrombus in the platelet pump.

RESULTS

When the impeller inlet angle is 30°, the outlet angle is 20°, and the number of blades is 6, the probability of thrombus formation is minimized in the orthogonal design method, aligning with the requirements for blood pump performance.

CONCLUSIONS

These design parameters serve as a numerical guideline for optimizing the geometry of the semi-open impeller in blood pumps and provide a theoretical foundation for subsequent in vitro experiments.

Numerical Study on the Impact of Central Venous Catheter Placement on Blood Flow in the Cavo-Atrial Junction

An in silico study is performed to investigate fluid dynamic effects of central venous catheter (CVC) placement within patient-specific cavo-atrial junctions. Prior studies show the CVC infusing a liquid, but this study focuses on the placement without any liquid emerging from the CVC. A 7 or 15-French double-lumen CVC is placed virtually in two patient-specific models; the CVC tip location is altered to understand its effect on the venous flow field. Results show that the CVC impact is trivial on flow in the superior vena cava when the catheter-to-vein ratio ranges from 0.15 to 0.33. Results further demonstrate that when the CVC tip is directly in the right atrium, flow vortices in the right atrium result in elevated wall shear stress near the tip hole. A recirculation region characterizes a spatially variable flow field inside the CVC side hole. Furthermore, flow stagnation is present near the internal side hole corners but an elevated wall shear stress near the curvature of the side hole's exit. These results suggest that optimal CVC tip location is within the superior vena cava, so as to lower the potential for platelet activation due to elevated shear stresses and that CVC geometry and location depth in the central vein significantly influences the local CVC fluid dynamics. A thrombosis model also shows thrombus formation at the side hole and tip hole. After modifying the catheter design, the hemodynamics change, which alter thrombus formation. Future studies are warranted to study CVC design and placement location in an effort to minimize CVC-induced thrombosis incidence.

A multi-constituent model for assessing the effect of impeller shroud on the thrombosis potential of a centrifugal blood pump

Centrifugal blood pumps can be used for treating heart failure patients. However, pump thrombosis has remained one of the complications that trouble clinical treatment. This study analyzed the effect of impeller shroud on the thrombosis risk of the blood pump, and predicted areas prone to thrombosis. Multi-constituent transport equations were presented, considering mechanical activation and biochemical activation. It was found that activated platelets concentration can increase with shear stress and adenosine diphosphate(ADP) concentration increasing, and the highest risk of thrombosis inside the blood pump was under extracorporeal membrane oxygenation (ECMO) mode. Under the same condition, ADP concentration and thrombosis index of semi-shroud impeller can increase by 7.3% and 7.2% compared to the closed-shroud impeller. The main reason for the increase in thrombosis risk was owing to elevated scalar shear stress and more coagulation promoting factor-ADP released. The regions with higher thrombosis potential were in the center hole, top and bottom clearance. As a novelty, the findings revealed that impeller shroud can influence mechanical and biochemical activation factors. It is useful for identifying potential risk regions of thrombus formation based on relative comparisons.

Influence of Hematocrit Level and Integrin αIIbβIII Function on vWF-Mediated Platelet Adhesion and Shear-Induced Platelet Aggregation in a Sudden Expansion

Purpose

Shear-mediated thrombosis is a clinically relevant phenomenon that underlies excessive arterial thrombosis and device-induced thrombosis. Red blood cells are known to mechanically contribute to physiological hemostasis through margination of platelets and vWF, facilitating the unfurling of vWF multimers, and increasing the fraction of thrombus-contacting platelets. Shear also plays a role in this phenomenon, increasing both the degree of margination and the near-wall forces experienced by vWF and platelets leading to unfurling and activation. Despite this, the contribution of red blood cells in shear-induced platelet aggregation has not been fully investigated-specifically the effect of elevated hematocrit has not yet been demonstrated.

Methods

Here, a microfluidic model of a sudden expansion is presented as a platform for investigating platelet adhesion at hematocrits ranging from 0 to 60% and shear rates ranging from 1000 to 10,000 s-1. The sudden expansion geometry models nonphysiological flow separation characteristic to mechanical circulatory support devices, and the validatory framework of the FDA benchmark nozzle. PDMS microchannels were fabricated and coated with human collagen. Platelets were fluorescently tagged, and blood was reconstituted at variable hematocrit prior to perfusion experiments. Integrin function of selected blood samples was inhibited by a blocking antibody, and platelet adhesion and aggregation over the course of perfusion was monitored.

Results

Increasing shear rates at physiological and elevated hematocrit levels facilitate robust platelet adhesion and formation of large aggregates. Shear-induced platelet aggregation is demonstrated to be dependent on both αIIbβIII function and the presence of red blood cells. Inhibition of αIIbβIII results in an 86.4% reduction in overall platelet adhesion and an 85.7% reduction in thrombus size at 20-60% hematocrit. Hematocrit levels of 20% are inadequate for effective platelet margination and subsequent vWF tethering, resulting in notable decreases in platelet adhesion at 5000 and 10,000 s-1 compared to 40% and 60%. Inhibition of αIIbβIII triggered dramatic reductions in overall thrombus coverage and large aggregate formation. Stability of platelets tethered by vWF are demonstrated to be αIIbβIII-dependent, as adhesion of single platelets treated with A2A9, an anti-αIIbβIII blocking antibody, is transient and did not lead to sustained thrombus formation.

Conclusions

This study highlights driving factors in vWF-mediated platelet adhesion that are relevant to clinical suppression of shear-induced thrombosis and in vitro assays of platelet adhesion. Primarily, increasing hematocrit promotes platelet margination, permitting shear-induced platelet aggregation through αIIbβIII-mediated adhesion at supraphysiological shear rates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12195-024-00796-0.

A numerical study of the effect of thrombus breakdown on predicted thrombus formation and growth

Thrombosis is a complex biological process which involves many biochemical reactions and is influenced by blood flow. Various computational models have been developed to simulate natural thrombosis in diseases such as aortic dissection (AD), and device-induced thrombosis in blood-contacting biomedical devices. While most hemodynamics-based models consider the role of low shear stress in the initiation and growth of thrombus, they often ignore the effect of thrombus breakdown induced by elevated shear stress. In this study, a new shear stress-induced thrombus breakdown function is proposed and implemented in our previously published thrombosis model. The performance of the refined model is assessed by quantitative comparison with experimental data on thrombus formation in a backward-facing step geometry, and qualitative comparison with in vivo data obtained from an AD patient. Our results show that incorporating thrombus breakdown improves accuracy in predicted thrombus volume and captures the same pattern of thrombus evolution as measured experimentally and in vivo. In the backward-facing step geometry, thrombus breakdown impedes growth over the step and downstream, allowing a stable thrombus to be reached more quickly. Moreover, the predicted thrombus volume, height and length are in better agreement with the experimental measurements compared to the original model which does not consider thrombus breakdown. In the patient-specific AD, the refined model outperforms the original model in predicting the extent and location of thrombosis. In conclusion, the effect of thrombus breakdown is not negligible and should be included in computational models of thrombosis.

Investigating clot-flow interactions by integrating intravital imaging with in silico modeling for analysis of flow, transport, and hemodynamic forces

As a blood clot forms, grows, deforms, and embolizes following a vascular injury, local clot-flow interactions lead to a highly dynamic flow environment. The local flow influences transport of biochemical species relevant for clotting, and determines the forces on the clot that in turn lead to clot deformation and embolization. Despite this central role, quantitative characterization of this dynamic clot-flow interaction and flow environment in the clot neighborhood remains a major challenge. Here, we propose an approach that integrates dynamic intravital imaging with computer geometric modeling and computational flow and transport modeling to develop a unified in silico framework to quantify the dynamic clot-flow interactions. We outline the development of the methodology referred to as Intravital Integrated In Silico Modeling or IVISim, and then demonstrate the method on a sample set of simulations comprising clot formation following laser injury in two mouse cremaster arteriole injury model data: one wild-type mouse case, and one diYF knockout mouse case. Simulation predictions are verified against experimental observations of transport of caged fluorescent Albumin (cAlb) in both models. Through these simulations, we illustrate how the IVISim methodology can provide insights into hemostatic processes, the role of flow and clot-flow interactions, and enable further investigations comparing and contrasting different biological model scenarios and parameter variations.

Clot embolization studies and computational framework for embolization in a canonical tube model

Despite recent advances in the development of computational methods of modeling thrombosis, relatively little effort has been made in developing methods of modeling blood clot embolization. Such a model would provide substantially greater understanding of the mechanics of embolization, as in-vitro and in-vivo characterization of embolization is difficult. Here, a method of computationally simulating embolization is developed. Experiments are performed of blood clots formed in a polycarbonate tube, where phosphate-buffered saline is run through the tube at increasing flow rates until the clot embolizes. The experiments revealed embolization can be initiated by leading edge and trailing edge detachment or by non-uniform detachment. Stress-relaxation experiments are also performed to establish values of constitutive parameters for subsequent simulations. The embolization in the tube is reproduced in silico using a multiphase volume-of-fluid approach, where the clot is modeled as viscoelastic. By varying the constitutive parameters at the wall, embolization can be reproduced in-silico at varying flow rates, and a range of constitutive parameters fitting the experiments is reported. Here, the leading edge embolization is simulated at flow rates consistent with the experiments demonstrating excellent agreement in this specific behavior.

Finite Element Implementation of Computational Fluid Dynamics With Reactive Neutral and Charged Solute Transport in FEBio

The objective of this study was to implement a novel fluid-solutes solver into the open-source finite element software FEBio, that extended available modeling capabilities for biological fluids and fluid-solute mixtures. Using a reactive mixture framework, this solver accommodates diffusion, convection, chemical reactions, electrical charge effects, and external body forces, without requiring stabilization methods that were deemed necessary in previous computational implementations of the convection-diffusion-reaction equation at high Peclet numbers. Verification and validation problems demonstrated the ability of this solver to produce solutions for Peclet numbers as high as 1011, spanning the range of physiological conditions for convection-dominated solute transport. This outcome was facilitated by the use of a formulation that accommodates realistic values for solvent compressibility, and by expressing the solute mass balance such that it properly captured convective transport by the solvent and produced a natural boundary condition of zero diffusive solute flux at outflow boundaries. Since this numerical scheme was not necessarily foolproof, guidelines were included to achieve better outcomes that minimize or eliminate the potential occurrence of numerical artifacts. The fluid-solutes solver presented in this study represents an important and novel advancement in the modeling capabilities for biomechanics and biophysics as it allows modeling of mechanobiological processes via the incorporation of chemical reactions involving neutral or charged solutes within dynamic fluid flow. The incorporation of charged solutes in a reactive framework represents a significant novelty of this solver. This framework also applies to a broader range of nonbiological applications.

clotFoam: An open-source framework to simulate blood clot formation under arterial flow

Blood clotting involves the coupled processes of platelet aggregation and coagulation. Simulating clotting under flow in complex geometries is challenging due to multiple temporal and spatial scales and high computational cost. clotFoam is an open-source software developed in OpenFOAM that employs a continuum model of platelet advection, diffusion, and aggregation in a dynamic fluid environment and a simplified coagulation model with proteins that advect, diffuse, and react within the fluid and with wall-bound species through reactive boundary conditions. Our framework provides the foundation on which one can build more complex models and perform reliable simulations in almost any computational domain.

Performance of a Hydrogel Coated Nitinol with Oligonucleotide-Modified Nanoparticles Within Turbulent Conditions of Blood-Contacting Devices

INTRODUCTION

Hydrogels offer a wide range of applications in the antithrombotic modification of biomedical devices. The functionalization of these hydrogels with potentially drug-laden nanoparticles in the context of deviceassociated turbulence is critically under-studied. Thus, the purpose of this study was to use a hydrogel-coating nitinol surface as a model to understand the functions of hydrogels and the capture of nanoparticles under clinically relevant flow conditions.

METHODS

Nitinol was coated by an oligonucleotide (ON) functionalized hydrogel. Nanoparticles were functionalized with complementary oligonucleotides (CONs). The capture of CONfunctionalized nanoparticles by the ON-functionalized hydrogel surfaces was studied under both static and dynamic attachment conditions. Fluorescent-labelling of nanoparticles was utilized to assess capture efficacy and resistance to removal by device-relevant flow conditions.

RESULTS

The specificity of the ON-CON bond was verified, exhibiting a dose-dependent attachment response. The hydrogel coating was resistant to stripping by flow, retaining >95% after exposure to one hour of turbulent flow. Attachment of nanoparticles to the hydrogel was higher in the static condition than under laminar flow (p < 0.01), but comparable to that of attachment under turbulent flow. Modified nitinol samples underwent one hour of flow treatment under both laminar and turbulent regimes and demonstrated decreased nanoparticle loss following static conjugation rather than turbulent conjugation (36.1% vs 53.8%, p < 0.05). There was no significant difference in nanoparticle functionalization by upstream injection between laminar and turbulent flow.

CONCLUSION

The results demonstrate promising potential of hydrogelfunctionalized nitinol for capturing nanoparticles using nucleic acid hybridization. The hydrogel structure and ONCON bond integrity both demonstrated a resistance to mechanical damage and loss of biomolecular functionalization by exposure to turbulence. Further investigation is warranted to highlight drug delivery and antithrombogenic modification applications of nanoparticle-functionalized hydrogels.

Shear-driven modelling of thrombus formation in type B aortic dissection

Background: Type B aortic dissection (TBAD) is a dangerous pathological condition with a high mortality rate. TBAD is initiated by an intimal tear that allows blood to flow between the aortic wall layers, causing them to separate. As a result, alongside the original aorta (true lumen), a false lumen (FL) develops. TBAD compromises the whole cardiovascular system, in the worst case resulting in complete aortic rupture. Clinical studies have shown that dilation and rupture of the FL are related to the failure of the FL to thrombose. Complete FL thrombosis has been found to improve the clinical outcomes of patients with chronic TBAD and is the desired outcome of any treatment. Partial FL thrombosis has been associated with late dissection-related deaths and the requirement for re-intervention, thus the level of FL thrombosis is dominant in classifying the risk of TBAD patients. Therefore, it is important to investigate and understand under which conditions complete thrombosis of the FL occurs.

Method: Local FL hemodynamics play an essential role in thrombus formation and growth. In this study, we developed a simplified phenomenological model to predict FL thrombosis in TBAD under physiological flow conditions. Based on an existing shear-driven thrombosis model, a comprehensive model reduction study was performed to improve computational efficiency. The reduced model has been implemented in Ansys CFX and applied to a TBAD case following thoracic endovascular aortic repair (TEVAR) to test the model. Predicted thrombus formation based on post-TEVAR geometry at 1-month was compared to actual thrombus formation observed on a 3-year follow-up CT scan.

Results: The predicted FL status is in excellent agreement with the 3-year follow-up scan, both in terms of thrombus location and total volume, thus validating the new model. The computational cost of the new model is significantly lower than the previous thrombus model, with an approximate 65% reduction in computational time. Such improvement means the new model is a significant step towards clinical applicability.

Conclusion: The thrombosis model developed in this study is accurate and efficient at predicting FL thrombosis based on patient-specific data, and may assist clinicians in choosing individualized treatments in the future.

A fibrin enhanced thrombosis model for medical devices operating at low shear regimes or large surface areas

Over the past decade, much of the development of computational models of device-related thrombosis has focused on platelet activity. While those models have been successful in predicting thrombus formation in medical devices operating at high shear rates (> 5000 s⁻¹), they cannot be directly applied to low-shear devices, such as blood oxygenators and catheters, where emerging information suggest that fibrin formation is the predominant mechanism of clotting and platelet activity plays a secondary role. In the current work, we augment an existing platelet-based model of thrombosis with a partial model of the coagulation cascade that includes contact activation of factor XII and fibrin production. To calibrate the model, we simulate a backward-facing-step flow channel that has been extensively characterized in-vitro. Next, we perform blood perfusion experiments through a microfluidic chamber mimicking a hollow fiber membrane oxygenator and validate the model against these observations. The simulation results closely match the time evolution of the thrombus height and length in the backward-facing-step experiment. Application of the model to the microfluidic hollow fiber bundle chamber capture both gross features such as the increasing clotting trend towards the outlet of the chamber, as well as finer local features such as the structure of fibrin around individual hollow fibers. Our results are in line with recent findings that suggest fibrin production, through contact activation of factor XII, drives the thrombus formation in medical devices operating at low shear rates with large surface area to volume ratios.

Patients treated with blood-contacting medical devices suffer from clotting complications. Over the past decades, a great effort has been made to develop computational tools to predict and prevent clot formation in these devices. However, most models have focused on platelet activity and neglected other important parts of the problem such as the coagulation cascade reactions that lead to fibrin formation. In the current work, we incorporate this missing element into a well-established and validated model for platelet activity. We then use this novel approach to predict thrombus formation in two experimental configurations. Our results confirm that to accurately predict the clotting process in devices where surface area to volume ratios are large and flow velocity and shear stresses remain low, coagulation reactions and subsequent fibrin formation must be considered. This new model could have great implications for the design and optimization of medical devices such as blood oxygenators. In the long term, the model could evolve into a functional tool to inform anticoagulation therapies for these patients.

Toward modeling thrombosis and thromboembolism in laminar and turbulent flow regimes

Thrombosis and thromboembolism are deadly risk factors in blood-contacting biomedical devices, and in-silico models of thrombosis are attractive tools to understand the mechanics of these processes, though the simulation of thromboembolism remains underdeveloped. The purpose of this study is to modify an existing computational thrombosis model to allow for thromboembolism and to investigate the behavior of the modified model at a range of flow rates. The new and existing models are observed to lead to similar predictions of thrombosis in a canonical backward-facing step geometry across flow rates, and neither model predicts thrombosis in a turbulent flow. Simulations are performed by increasing flow rates in the case of a clot formed at lower flow to induce embolization. While embolization is observed, most of the clot breakdown is by shear rather than by breakup and subsequent transport of clotted material, and further work is required in the formulation and validation of embolization. This model provides a framework to further investigate thromboembolization.

A New Mathematical Numerical Model to Evaluate the Risk of Thrombosis in Three Clinical Ventricular Assist Devices

(1) Background: Thrombosis is the main complication in patients supported with ventricular assist devices (VAD). Models that accurately predict the risk of thrombus formation in VADs are still lacking. When VADs are clinically assisted, their complex geometric configuration and high rotating speed inevitably generate complex flow fields and high shear stress. These non-physiological factors can damage blood cells and proteins, release coagulant factors and trigger thrombosis. In this study, a more accurate model for thrombus assessment was constructed by integrating parameters such as shear stress, residence time and coagulant factors, so as to accurately assess the probability of thrombosis in three clinical VADs. (2) Methods: A mathematical model was constructed to assess platelet activation and thrombosis within VADs. By solving the transport equation, the influence of various factors such as shear stress, residence time and coagulation factors on platelet activation was considered. The diffusion equation was applied to determine the role of activated platelets and substance deposition on thrombus formation. The momentum equation was introduced to describe the obstruction to blood flow when thrombus is formed, and finally a more comprehensive and accurate model for thrombus assessment in patients with VAD was obtained. Numerical simulations of three clinically VADs (CH-VAD, HVAD and HMII) were performed using this model. The simulation results were compared with experimental data on platelet activation caused by the three VADs. The simulated thrombogenic potential in different regions of MHII was compared with the frequency of thrombosis occurring in the regions in clinic. The regions of high thrombotic risk for HVAD and HMII observed in experiments were compared with the regions predicted by simulation. (3) Results: It was found that the percentage of activated platelets within the VAD obtained by solving the thrombosis model developed in this study was in high agreement with the experimental data (r² = 0.984), the likelihood of thrombosis in the regions of the simulation showed excellent correlation with the clinical statistics (r² = 0.994), and the regions of high thrombotic risk predicted by the simulation were consistent with the experimental results. Further study revealed that the three clinical VADs (CH-VAD, HVAD and HMII) were prone to thrombus formation in the inner side of the secondary flow passage, the clearance between cone and impeller, and the corner region of the inlet pipe, respectively. The risk of platelet activation and thrombus formation for the three VADs was low to high for CH-VAD, HVAD, and HM II, respectively. (4) Conclusions: In this study, a more comprehensive and accurate thrombosis model was constructed by combining parameters such as shear stress, residence time, and coagulation factors. Simulation results of thrombotic risk received with this model showed excellent correlation with experimental and clinical data. It is important for determining the degree of platelet activation in VAD and identifying regions prone to thrombus formation, as well as guiding the optimal design of VAD and clinical treatment.

Computational Prediction of Thrombosis in Food and Drug Administration's Benchmark Nozzle

Thrombosis seriously threatens human cardiovascular health and the safe operation of medical devices. The Food and Drug Administration’s (FDA) benchmark nozzle model was designed to include the typical structure of medical devices. However, the thrombosis in the FDA nozzle has yet not been investigated. The objective of this study is to predict the thrombus formation process in the idealized medical device by coupling computational fluid dynamics and a macroscopic hemodynamic-based thrombus model. We developed the hemodynamic-based thrombus model by considering the effect of platelet consumption. The thrombus model was quantitatively validated by referring to the latest thrombosis experiment, which was performed in a backward-facing step with human blood flow. The same setup was applied in the FDA nozzle to simulate the thrombus formation process. The thrombus shaped like a ring was firstly observed in the FDA benchmark nozzle. Subsequently, the accuracy of the shear-stress transport turbulence model was confirmed in different turbulent flow conditions. Five scenarios with different Reynolds numbers were carried out. We found that turbulence could change the shape of centrosymmetric thrombus to axisymmetric and high Reynolds number blood flow would delay or even prevent thrombosis. Overall, the present study reports the thrombosis process in the FDA benchmark nozzle using the numerical simulation method, and the primary findings may shed light on the effect of turbulence on thrombosis.

Development of computational models for microtesla-level magnetic brain scanning: a novel avenue for device development

BACKGROUND

Detection of locally increased blood concentration and perfusion is critical for assessment of functional cortical activity as well as diagnosis of conditions such as intracerebral hemorrhage (ICH). Current paradigms for assessment of regional blood concentration in the brain rely on computed tomography (CT), magnetic resonance imaging (MRI), and perfusion blood oxygen level dependent functional magnetic resonance imaging (BOLD-fMRI).

RESULTS

In this study, we developed computational models to test the feasibility of novel magnetic sensors capable of detecting hemodynamic changes within the brain on a microtesla-level. We show that low-field magnetic sensors can accurately detect changes in magnetic flux density and eddy current damping signals resulting from increases in local blood concentration. These models predicted that blood volume changes as small as 1.26 mL may be resolved by the sensors, implying potential use for diagnosis of ICH and assessment of regional blood flow as a proxy for cerebral metabolism and neuronal activity. We then translated findings from our computational model to demonstrate feasibility of accurate detection of modeled ICH in a simulated human cadaver setting.

CONCLUSIONS

Overall, microtesla-level magnetic scanning is feasible, safe, and has distinct advantages compared to current standards of care. Computational modeling may facilitate rapid prototype development and testing of novel medical devices with minimal risk to human participants prior to device construction and clinical trials.

An Accelerated Thrombosis Model for Computational Fluid Dynamics Simulations in Rotary Blood Pumps

Purpose

Thrombosis ranks among the major complications in blood-carrying medical devices and a better understanding to influence the design related contribution to thrombosis is desirable. Over the past years many computational models of thrombosis have been developed. However, numerically cheap models able to predict localized thrombus risk in complex geometries are still lacking. The aim of the study was to develop and test a computationally efficient model for thrombus risk prediction in rotary blood pumps.

Methods

We used a two-stage approach to calculate thrombus risk. The first stage involves the computation of velocity and pressure fields by computational fluid dynamic simulations. At the second stage, platelet activation by mechanical and chemical stimuli was determined through species transport with an Eulerian approach. The model was compared with existing clinical data on thrombus deposition within the HeartMate II. Furthermore, an operating point and model parameter sensitivity analysis was performed.

Results

Our model shows good correlation (R² > 0.93) with clinical data and identifies the bearing and outlet stator region of the HeartMate II as the location most prone to thrombus formation. The calculation of thrombus risk requires an additional 10–20 core hours of computation time.

Conclusion

The concentration of activated platelets can be used as a surrogate and computationally low-cost marker to determine potential risk regions of thrombus deposition in a blood pump. Relative comparisons of thrombus risk are possible even considering the intrinsic uncertainty in model parameters and operating conditions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13239-021-00606-y.

Mathematical and Computational Modeling of Device-Induced Thrombosis

Given the extensive and routine use of cardiovascular devices, a major limiting factor to their success is the thrombotic rate that occurs. This both poses direct risk to the patient and requires counterbalancing with anticoagulation and other treatment strategies, contributing additional risks. Developing a better understanding of the mechanisms of device-induced thrombosis to aid in device design and medical management of patients is critical to advance the ubiquitous use and durability. Thus, mathematical and computational modelling of device-induced thrombosis has received significant attention recently, but challenges remain. Additional areas that need to be explored include microscopic/macroscopic approaches, reconciling physical and numerical timescales, immune/inflammatory responses, experimental validation, and incorporating pathologies and blood conditions. Addressing these areas will provide engineers and clinicians the tools to provide safe and effective cardiovascular devices.

Computational modeling of blood component transport related to coronary artery thrombosis in Kawasaki disease

Coronary artery thrombosis is the major risk associated with Kawasaki disease (KD). Long-term management of KD patients with persistent aneurysms requires a thrombotic risk assessment and clinical decisions regarding the administration of anticoagulation therapy. Computational fluid dynamics has demonstrated that abnormal KD coronary artery hemodynamics can be associated with thrombosis. However, the underlying mechanisms of clot formation are not yet fully understood. Here we present a new model incorporating data from patient-specific simulated velocity fields to track platelet activation and accumulation. We use a system of Reaction-Advection-Diffusion equations solved with a stabilized finite element method to describe the evolution of non-activated platelets and activated platelet concentrations [AP], local concentrations of adenosine diphosphate (ADP) and poly-phosphate (PolyP). The activation of platelets is modeled as a function of shear-rate exposure and local concentration of agonists. We compared the distribution of activated platelets in a healthy coronary case and six cases with coronary artery aneurysms caused by KD, including three with confirmed thrombosis. Results show spatial correlation between regions of higher concentration of activated platelets and the reported location of the clot, suggesting predictive capabilities of this model towards identifying regions at high risk for thrombosis. Also, the concentration levels of ADP and PolyP in cases with confirmed thrombosis are higher than the reported critical values associated with platelet aggregation (ADP) and activation of the intrinsic coagulation pathway (PolyP). These findings suggest the potential initiation of a coagulation pathway even in the absence of an extrinsic factor. Finally, computational simulations show that in regions of flow stagnation, biochemical activation, as a result of local agonist concentration, is dominant. Identifying the leading factors to a pro-coagulant environment in each case—mechanical or biochemical—could help define improved strategies for thrombosis prevention tailored for each patient.

Computational studies aiming to model thrombosis often rely on an arterial wall injury. Collagen and other extracellular matrix components are exposed to the bloodstream, which facilitates platelet adhesion to the wall and subsequent clot formation. However, these models are not adequate to explain thrombosis in other settings where even in the absence of a focal lesion, clots may still form under certain flow conditions. Coronary artery aneurysm thrombosis following KD is an example of the need to understand the mechanisms of thrombus initiation in the absence of an extrinsic factor. This study provides a new framework to investigate thrombus initiation in KD from a patient-specific perspective, which integrates fluid mechanics and biochemistry and which could help quantify the pro-coagulant environment induced by the aneurysm and become a predictive tool. The work presented here has broad relevance to other clinical situations where flow stagnation and transport are driving factors in thrombus formation.

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.

Hybrid modelling for stroke care: Review and suggestions of new approaches for risk assessment and simulation of scenarios

Stroke is an example of a complex and multi-factorial disease involving multiple organs, timescales, and disease mechanisms. To deal with this complexity, and to realize Precision Medicine of stroke, mathematical models are needed. Such approaches include: 1) machine learning, 2) bioinformatic network models, and 3) mechanistic models. Since these three approaches have complementary strengths and weaknesses, a hybrid modelling approach combining them would be the most beneficial. However, no concrete approach ready to be implemented for a specific disease has been presented to date. In this paper, we both review the strengths and weaknesses of the three approaches, and propose a roadmap for hybrid modelling in the case of stroke care. We focus on two main tasks needed for the clinical setting: a) For stroke risk calculation, we propose a new two-step approach, where non-linear mixed effects models and bioinformatic network models yield biomarkers which are used as input to a machine learning model and b) For simulation of care scenarios, we propose a new four-step approach, which revolves around iterations between simulations of the mechanistic models and imputations of non-modelled or non-measured variables. We illustrate and discuss the different approaches in the context of Precision Medicine for stroke.

Hemodynamic Parameters Predict In-stent Thrombosis After Multibranched Endovascular Repair of Complex Abdominal Aortic Aneurysms: A Retrospective Study of Branched Stent-Graft Thrombosis

Background: Branch vessel occlusion is reported in endovascular repair of aortic pathology. This study aimed to evaluate the hemodynamic indicators associated with in-stent thrombosis (IST) of a branched stent-graft (BSG) after endovascular aortic repair (EVAR) of a complex abdominal aortic aneurysm.

Methods: A retrospective evaluation was performed based on the computed tomography (CT) scans and clinical data of three patients who underwent multi-branched endovascular repair. Patient-specific 3-dimensional models were reconstructed, and hemodynamic analysis was performed for IST. Hemodynamics-related parameters including time-averaged wall shear stress (TAWSS), oscillatory shear stress index (OSI), and relative residence time (RRT) were compared among the individual patients.

Results: The flow velocity, TAWSS, OSI, and RRT were radically changed in the area of the IST. In BSGs, IST tended to occur in the regions of hemodynamic alteration near the bends in the device, where a decreased flow velocity (<0.6 m/s) and TAWSS (<0.8 Pa) and an elevated OSI (>0.2) and RRT (>5 s) were consistently observed.

Conclusions: Hemodynamic perturbations in BSGs cause a predisposition to IST, which can be predicted by a series of changes in the flow parameters. Early hemodynamic analysis might be useful for identifying and remediating IST after multibranched endovascular repair.

Refining a numerical model for device-induced thrombosis and investigating the effects of non-Newtonian blood models

Thrombosis is one of the main causes of failure in device implantation. Computational thrombosis simulation is a convenient approach to evaluate the risk of thrombosis for a device. However, thrombosis is a complicated process involving multiple species and reactions. Application of a macroscopic, single-scale computational model for device-induced thrombosis is a cost-effective approach. The current study has refined an existing thrombosis model, which simulated thrombosis by tracing four species in blood: non-activated platelets, activated platelets, surface adherent platelets, and ADP. Platelets are activated mechanically by shear stress, and chemically by ADP. Platelet adhesion occurs on surfaces with low wall shear stress with platelet aggregation inhibited in regions of high shear stress. The study improves the existing thrombosis model by: 1) Modifying the chemical platelet activation function so that ADP activates platelets; 2) Modifying the function describing thrombus deposition and growth to distinguish between thrombus deposition on wall surfaces and thrombus growth on existing thrombus surfaces; 3) Modifying the thrombus breakdown function to allow for thrombus breakdown by shear stress; 4) Modeling blood flow as non-Newtonian. The results show that the inclusion of ADP and the use of a non-Newtonian model improve agreement with experimental data.

A two-phase flow approach for modeling blood stasis and estimating the thrombosis potential of a ventricular assist device

Thrombosis and its related events have become a major concern during the development and optimization of ventricular assist devices (VADs, also called blood pumps), and limit their clinical use and economic benefits. Attempts have been made to model the thrombosis formation, considering hemodynamic and biochemical processes. However, the complexities and computational expenses are prohibitive. Blood stasis is one of the key factors which may lead to the formation of thrombosis and excessive thromboembolic risks for patients. This study proposed a novel approach for modeling blood stasis, based on a two-phase flow principle. The locations of blood residual can be tracked over time, so that regions of blood stasis can be identified. The blood stasis in an axial blood pump is simulated under various working conditions, the results agree well with the experimental results. In contrast, conventional hemodynamic metrics such as velocity, time-averaged wall shear stress (TAWSS), and relative residence time (RRT), were contradictory in judging risk of blood stasis and thrombosis, and inconsistent with experimental results. We also found that the pump operating at the designed rotational speed is less prone to blood stasis. The model provides an efficient and fast alternative for evaluating blood stasis and thrombosis potential in blood pumps, and will be a valuable addition to the tools to support the design and improvement of VADs.

In vitro real-time magnetic resonance imaging for quantification of thrombosis

OBJECTIVES

Thrombosis is a leading cause of failure for cardiovascular devices. While computational simulations are a powerful tool to predict thrombosis and evaluate the risk for medical devices, limited experimental data are available to validate the simulations. The aim of the current study is to provide experimental data of a growing thrombus for device-induced thrombosis.

MATERIALS AND METHODS

Thrombosis within a backward-facing step (BFS), or sudden expansion was investigated, using bovine and human blood circulated through the BFS model for 30 min, with a constant inflow rate of 0.76 L/min. Real-time three-dimensional flow-compensated magnetic resonance imaging (MRI), supported with Magnevist, a contrast agent improving thrombus delineation, was applied to quantify thrombus deposition and growth within the model.

RESULTS

The study showed that the BFS model induced a flow recirculation region, which facilitated thrombosis. By 30 min, in comparison to bovine blood, human blood resulted in smaller thrombus formation, in terms of the length (13.3 ± 0.6 vs. 18.1 ± 1.3 mm), height (2.3 ± 0.1 vs. 2.6 ± 0.04 mm), surface area exposed to blood (0.67 ± 0.03 vs 1.05 ± 0.08 cm²), and volume (0.069 ± 0.004 vs. 0.093 ± 0.007 cm³), with p < 0.01. Normalization of the thrombus measurements, which excluded the flow recirculation effects, suggested that the thrombus sizes increased during the first 15 min and stabilized after 20 min. Blood properties, including viscosity, hematocrit, and platelet count affected thrombosis.

CONCLUSION

For the first time, contrast agent-supported real-time MRI was performed to investigate thrombus deposition and growth within a sudden expansion. This study provides experimental data for device-induced thrombosis, which is valuable for validation of computational thrombosis simulations.

Proceedings From the 2019 Stanford Single Ventricle Scientific Summit: Advancing Science for Single Ventricle Patients: From Discovery to Clinical Applications

Because of remarkable advances in survival over the past 40 years, the worldwide population of individuals with single ventricle heart disease living with Fontan circulation has grown to ≈70 000, with nearly half aged >18 years. Survival to at least 30 years of age is now achievable for 75% of Fontan patients. On the other hand, single ventricle patients account for the largest group of the 6000 to 8000 children hospitalized with circulation failure, with or without heart failure annually in the United States, with the highest in‐hospital mortality. Because there is little understanding of the underlying mechanisms of heart failure, arrhythmias, pulmonary and lymphatic vascular abnormalities, and other morbidities, there are no specific treatments to maintain long‐term myocardial performance or to optimize overall patient outcomes.

Platelet deposition estimation: A novel method for emulating the pump thrombosis potential of blood pumps

Pump thrombosis potential exists in most blood pumps and limits their clinical use. To improve the pump thrombosis performance of blood pumps, a method for emulating the platelet deposition on the flow passage component surfaces inside blood pumps was presented and tested. The method emulates the blood platelet deposition, employing laser-induced fluorescence tracing technology. The blood pump was rotated in a mock circulation loop with deionized water filled with fluorescent particles. The component surfaces were then explored via laser. The fluorescent particles were induced by laser and imaged in a charge-coupled device (CCD) camera to show the distribution of fluorescent particles gathering on the blood pump component surfaces. The activated platelet deposition was emulated by fluorescent particle gathering. The experiment showed obvious particle gathering on the interface surfaces and cross-sectional surface (perpendicular to the flow). This platelet deposition estimation (PDE) method can be easily incorporated in the in vitro testing phase to analyze and decrease a pump's thrombosis potential before animal experimentation, thereby reducing the cost of blood pump development. This methodology of emulating blood platelet deposition indicates its potential for improving flow passage component structure and reducing device thrombosis of blood pumps.

Platelet adhesion emulation: A novel method for estimating the device thrombosis potential of a ventricular assist device

Device thrombosis inside ventricular assist devices remains a limitation to their long-term clinical use. Thrombosis potential exists in almost all ventricular assist devices because the device-induced high shear stress and vortices can activate platelets, which then aggregate and adhere to the surfaces inside the ventricular assist device. To decrease the device thrombosis potential of long-term use of ventricular assist devices, a methodology entitled platelet adhesion emulation for predicting the thrombosis potential and thrombosis position inside the ventricular assist devices is developed. The platelet adhesion emulation methodology combines numerical simulations with in vitro experiments by correlating the structure of the flow passage components within the ventricular assist device with the platelet adhesion to estimate the thrombosis potential and location, with the goal of developing ventricular assist devices with optimized antithrombotic performance. Platelet adhesion emulation is aimed at decreasing the device thrombus potential of ventricular assist devices. The platelet adhesion emulation effectiveness is validated by simulating and testing an axial left ventricular assist device. The blood velocity relative to the surfaces of the flow passage components is calculated to estimate the platelet adhesion potential, indicating the probability of thrombus formation on the surfaces. Platelet adhesion emulation experiments conducted in a mock circulation loop with pump prototypes show the distribution of platelet adhesion on the surfaces. This methodology of emulating the device thrombosis distribution indicates the potential for improving the component structure and reducing the device thrombosis of ventricular assist devices.

Automated reduction of blood coagulation models

Mathematical modeling of thrombosis typically involves modeling the coagulation cascade. Models of coagulation generally involve the reaction kinetics for dozens of proteins. The resulting system of equations is difficult to parameterize, and its numerical solution is challenging when coupled to blood flow or other physics important to clotting. Prior research suggests that essential aspects of coagulation may be reproduced by simpler models. This evidence motivates a systematic approach to model reduction. We herein introduce an automated framework to generate reduced-order models of blood coagulation. The framework consists of nested optimizations, where an outer optimization selects the optimal species for the reduced-order model and an inner optimization selects the optimal reaction rates for the new coagulation network. The framework was tested on an established 34-species coagulation model to rigorously consider what level of model fidelity is necessary to capture essential coagulation dynamics. The results indicate that a nine-species reduced-order model is sufficient to reproduce the thrombin dynamics of the benchmark 34-species model for a range of tissue factor concentrations, including those not included in the optimization process. Further model reduction begins to compromise the ability to capture the thrombin generation process. The framework proposed herein enables automated development of reduced-order models of coagulation that maintain essential dynamics used to model thrombosis.

In vivo measurement of blood clot mechanics from computational fluid dynamics based on intravital microscopy images

Ischemia which leads to heart attacks and strokes is one of the major causes of death in the world. Whether an occlusion occurs or not depends on the ability of a growing thrombus to resist flow forces exerted on its structure. This manuscript provides the first known in vivo measurement of how much stress a clot can withstand, before yielding to the surrounding blood flow. Namely, Lattice-Boltzmann Method flow simulations are performed based on 3D clot geometries, which are estimated from intravital microscopy images of laser-induced injuries in cremaster microvasculature of live mice. In addition to reporting the blood clot yield stresses, we also show that the thrombus "core" does not experience significant deformation, while its "shell" does. This indicates that the shell is more prone to embolization. Therefore, drugs should be designed to target the shell selectively, while leaving the core intact to minimize excessive bleeding. Finally, we laid down a foundation for a nondimensionalization procedure which unraveled a relationship between clot mechanics and biology. Hence, the proposed framework could ultimately lead to a unified theory of thrombogenesis, capable of explaining all clotting events. Thus, the findings presented herein will be beneficial to the understanding and treatment of heart attacks, strokes and hemophilia.

A CFD-based Kriging surrogate modeling approach for predicting device-specific hemolysis power law coefficients in blood-contacting medical devices

Most stress-based hemolysis models used in computational fluid dynamics (CFD) are based on an empirical power law correlation between hemolysis generation and the flow-induced stress and exposure time. Empirical model coefficients are typically determined by fitting global hemolysis measurements in simplified blood shearing devices under uniform shear conditions and with well-defined exposure times. CFD simulations using these idealized global empirical coefficients are then performed to predict hemolysis in a medical device with complex hemodynamics. The applicability, however, of this traditional approach of using idealized coefficients for a real device with varying exposure times and non-uniform shear is currently unknown. In this study, we propose a new approach for determining device- and species-specific hemolysis power law coefficients (C, a, and b). The approach consists of calculating multiple hemolysis solutions using different sets of coefficients to map the hemolysis response field in three-dimensional (C, a, b) parameter space. The resultant response field is then compared with experimental data in the same device to determine the coefficients that when incorporated into the locally defined power law model yield correct global hemolysis predictions. We first develop the generalized approach by deriving analytical solutions for simple uniform and non-uniform shear flows (planar Couette flow and circular Poiseuille flow, respectively) that allow us to continuously map the hemolysis solution in (C, a, b) parameter space. We then extend our approach to more practical cases relevant to blood-contacting medical devices by replacing the requirement for an analytical solution in our generalized approach with CFD and Kriging surrogate modeling. Finally, we apply our verified CFD-based Kriging surrogate modeling approach to predict the device- and species-specific power law coefficients for developing laminar flow in a small capillary tube. We show that the resultant coefficients are much different than traditional idealized coefficients obtained from simplified uniform shear experiments and that using such idealized coefficients yields a highly inaccurate prediction of hemolysis that is in error by more than 2000% compared to experiments. Our approach and surrogate modeling framework may be applied to more complex medical devices and readily extended to determine empirical coefficients for other continuum-based models of hemolysis and other forms of flow-induced blood damage (e.g., platelet activation and thrombosis).

Hemodynamic variables in aneurysms are associated with thrombotic risk in children with Kawasaki disease

BACKGROUND

Thrombosis is a major adverse outcome associated with coronary artery aneurysms (CAAs) resulting from Kawasaki disease (KD). Clinical guidelines recommend initiation of anticoagulation therapy with maximum CAA diameter (D_max) ≥8 mm or Z-score ≥ 10. Here, we investigate the role of aneurysm hemodynamics as a superior method for thrombotic risk stratification in KD patients.

METHODS AND RESULTS

We retrospectively studied ten KD patients with CAAs, including five patients who developed thrombosis. We constructed patient-specific anatomic models from cardiac magnetic resonance images and performed computational hemodynamic simulations using SimVascular. Our simulations incorporated pulsatile flow, deformable arterial walls and boundary conditions automatically tuned to match patient-specific arterial pressure and cardiac output. From simulation results, we derived local hemodynamic variables including time-averaged wall shear stress (TAWSS), low wall shear stress exposure, and oscillatory shear index (OSI). Local TAWSS was significantly lower in CAAs that developed thrombosis (1.2 ± 0.94 vs. 7.28 ± 9.77 dynes/cm², p = 0.006) and the fraction of CAA surface area exposed to low wall shear stress was larger (0.69 ± 0.17 vs. 0.25 ± 0.26%, p = 0.005). Similarly, longer residence times were obtained in branches where thrombosis was confirmed (9.07 ± 6.26 vs. 2.05 ± 2.91 cycles, p = 0.004). No significant differences were found for OSI or anatomical measurements such us D_max and Z-score. Assessment of thrombotic risk according to hemodynamic variables had higher sensitivity and specificity compared to standard clinical metrics (D_max, Z-score).

CONCLUSIONS

Hemodynamic variables can be obtained non-invasively via simulation and may provide improved thrombotic risk stratification compared to current diameter-based metrics, facilitating long-term clinical management of KD patients with persistent CAAs.

Numerical investigation of the influence of a bearing/shaft structure in an axial blood pump on the potential for device thrombosis

Adverse events caused by flow-induced thrombus formation around the bearing/shaft of an axial blood pump remain a serious problem for axial blood pumps. Moreover, excessive anticoagulation with thrombosis around the bearing potentially increases the risk of postoperative gastrointestinal bleeding. The purpose of this study is to analyze the influence of the bearing structure on the thrombosis potential of an axial blood pump. The bearing/shaft structure was embedded into an axial blood pump numerical model. The numerical simulation and analysis are focused on the low wall shear stresses, recirculation, and residence time close to the bearing region to evaluate the potential for thrombosis around the bearing. Then, the flow field near the blood pump bearing was tested via in vitro particle image velocimetry experiments to verify the numerical results. The simulation results showed that after embedding the bearing/shaft structure a recirculation zone appeared in the outlet guide vane bearing/shaft region, the residence time increased 11-fold in comparison to the pump without the bearing/shaft structure, the scalar shear stress in the shaft surface was less than 7.8 Pa, and the stress accumulation was less than 0.10 Pa s. The numerical results showed that platelets that flow through the bearing region are exposed to significantly lower wall shear stress and a longer residence time, leading to activated platelet adhesion. The reduced stress accumulation and increased time in the bearing region lead to increased platelet activation.

Abdominal aortic aneurysm endovascular repair: profiling post-implantation morphometry and hemodynamics with image-based computational fluid dynamics

Endovascular aneurysm repair (EVAR) has disseminated rapidly as an alternative to open surgical repair for the treatment of abdominal aortic aneurysms (AAAs), because of its reduced invasiveness, low mortality and morbidity rate. The effectiveness of the endovascular devices used in EVAR is always at question as postoperative adverse events can lead to re-intervention or to a possible fatal scenario for the circulatory system. Motivated by the assessment of the risks related to thrombus formation, here the impact of two different commercial endovascular grafts on local hemodynamics is explored through 20 image-based computational hemodynamic models of EVAR-treated patients (N=10 per each endograft model). Hemodynamic features, susceptible to promote thrombus formation, such as flow separation and recirculation, are quantitatively assessed and compared with the local hemodynamics established in image-based infrarenal abdominal aortic models of healthy subjects (N=10). The hemodynamic analysis is complemented by a geometrical characterization of the EVAR-induced reshaping of the infrarenal abdominal aortic vascular region. The findings of this study indicate that: (1) the clinically observed propensity to thrombus formation in devices used in EVAR strategies can be explained in terms of local hemodynamics by means of image-based computational hemodynamics approach; (2) reportedly pro-thrombotic hemodynamic structures are strongly correlated with the geometry of the aortoiliac tract postoperatively. In perspective, our study suggests that future clinical follow up studies could include a geometric analysis of the region of the implant, monitoring shape variations that can lead to hemodynamic disturbances of clinical significance.

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.

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.

High Frequency Components of Hemodynamic Shear Stress Profiles are a Major Determinant of Shear-Mediated Platelet Activation in Therapeutic Blood Recirculating Devices

We systematically analyzed the relative contributions of frequency component elements of hemodynamic shear stress waveforms encountered in cardiovascular blood recirculating devices as to overall platelet activation over time. We demonstrated that high frequency oscillations are the major determinants for priming, triggering and yielding activated “prothrombotic behavior” for stimulated platelets, even if the imparted shear stress has low magnitude and brief exposure time. Conversely, the low frequency components of the stress signal, with limited oscillations over time, did not induce significant activation, despite being of high magnitude and/or exposure time. In vitro data were compared with numerical predictions computed according to a recently proposed numerical model of shear-mediated platelet activation. The numerical model effectively resolved the correlation between platelet activation and the various frequency components examined. However, numerical predictions exhibited a different activation trend compared to experimental results for different time points of a stress activation sequence. With this study we provide a more fundamental understanding for the mechanobiological responsiveness of circulating platelets to the hemodynamic environment of cardiovascular devices, and the importance of these environments in mediating life-threatening thromboembolic complications associated with shear-mediated platelet activation. Experimental data will guide further optimization of the thromboresistance of cardiovascular implantable therapeutic devices.

Prediction of Thrombus Growth: Effect of Stenosis and Reynolds Number

Shear stresses play a major role in platelet-substrate interactions and thrombus formation and growth in blood flow, where under both pathological and physiological conditions platelet adhesion and accumulation occur. In this study, a shear-dependent continuum model for platelet activation, adhesion and aggregation is presented. The model was first verified under three different shear conditions and at two heparin levels. Three-dimensional simulations were then carried out to evaluate the performance of the model for severely damaged (stripped) aortas with mild and severe stenosis degrees in laminar flow regime. For these cases, linear shear-dependent functions were developed for platelet-surface and platelet-platelet adhesion rates. It was confirmed that the platelet adhesion rate is not only a function of Reynolds number (or wall shear rate) but also the stenosis severity of the vessel. General correlations for adhesion rates of platelets as functions of stenosis and Reynolds number were obtained based on these cases. Finally using the new platelet adhesion rates, the model was applied to different experimental systems and shown to agree well with measured platelet deposition.

Multi-Constituent Simulation of Thrombus Deposition

In this paper, we present a spatio-temporal mathematical model for simulating the formation and growth of a thrombus. Blood is treated as a multi-constituent mixture comprised of a linear fluid phase and a thrombus (solid) phase. The transport and reactions of 10 chemical and biological species are incorporated using a system of coupled convection-reaction-diffusion (CRD) equations to represent three processes in thrombus formation: initiation, propagation and stabilization. Computational fluid dynamic (CFD) simulations using the libraries of OpenFOAM were performed for two illustrative benchmark problems: in vivo thrombus growth in an injured blood vessel and in vitro thrombus deposition in micro-channels (1.5 mm × 1.6 mm × 0.1 mm) with small crevices (125 μm × 75 μm and 125 μm × 137 μm). For both problems, the simulated thrombus deposition agreed very well with experimental observations, both spatially and temporally. Based on the success with these two benchmark problems, which have very different flow conditions and biological environments, we believe that the current model will provide useful insight into the genesis of thrombosis in blood-wetted devices, and provide a tool for the design of less thrombogenic devices.