Influence of shear rate and surface chemistry on thrombus formation in micro-crevice

In Vitro and In Silico Characterization of the Aggregation of Thrombi on Textured Ventricular Cannula

The unacceptably high stroke rate associated with HeartMate 3 ventricular assist device (VAD) without signs of adherent pump thrombosis is hypothesized to be the result of the emboli produced by the inflow cannula, that are ingested and ejected from the pump. This in vitro and numerical study aimed to emulate the surface features and supraphysiological shear of a ventricular cannula to provide insight into their effect on thrombogenesis. Human whole blood was perfused at calibrated flow rates in a microfluidic channel to achieve shear rates 1000-7500 s-1, comparable to that experienced on the cannula. The channel contained periodic teeth representative of the rough sintered surface of the HeartMate 3 cannula. The deposition of fluorescently labeled platelets was visualized in real time and analyzed with a custom entity tracking algorithm. Numerical simulations of a multi-constituent thrombosis model were performed to simulate laminar blood flow in the channel. The sustained growth of adherent platelets was observed in all shear conditions ( p <  0.05). However, the greatest deposition was observed at the lower shear rates. The location of deposition with respect to the microfluidic teeth was also found to vary with shear rate. This was confirmed by CFD simulation. The entity tracking algorithm revealed the spatial variation of instances of embolic events. This result suggests that the sintered surface of the ventricular cannula may engender unstable thrombi with a greater likelihood of embolization at supraphysiological shear rates.

Hemodynamics in nutcracker syndrome: implications for diagnosis

BACKGROUND

Nutcracker syndrome is a disease characterized by complex symptoms, making its diagnosis challenging and often delayed, often resulting in a painful experience for the patients.

OBJECTIVE

This study aimed to investigate the pathogenesis of nutcracker syndrome through the perspective of hemodynamics by simulating blood flow with varying compression degrees of the left renal vein.

METHODS

3D patient-specific vascular models of the abdominal aorta, superior mesenteric artery and left renal vein were constructed based on CT images of patients suspected of having nutcracker syndrome. A hemodynamic simulation was then conducted using computational fluid dynamics to identify the correlation between alterations in hemodynamic parameters and varying degrees of compression.

RESULTS

The study indicated the presence of an evident gradient in velocity distribution over the left renal vein with relatively high degrees of stenosis (α ≤ 50°), with maximum velocity in the central region of the stenosis. Additionally, when the compression degree of the left renal vein increases, the pressure distribution of the left renal vein presents an increasing number of gradient layers. Furthermore, the wall shear stress shows a correlation with the variation of blood flow velocity, i.e., the increase of wall shear stress correlates with the acceleration of the blood flow velocity.

CONCLUSIONS

Using computational fluid dynamics as a non-invasive instrument to obtain the hemodynamic characteristics of nutcracker syndrome is feasible and could provide insights into the pathological mechanisms of the nutcracker syndrome supporting clinicians in diagnosis.

Multi-indicator analysis of mechanical blood damage with five clinical ventricular assist devices

PURPOSE

Device-induced blood damage contributes the hemolysis, thrombosis and bleeding complications in patients supported with ventricular assist device (VAD). This study aims to use a multi-indicator method to understand how devices causes blood damage and identify the "hot spots" of blood trauma within VADs.

METHODS

Computational fluid dynamics (CFD) methods were chosen to investigate the hemodynamic features of five clinical VADs (Impella 5.0, UltraMag, CHVAD, HVAD, and HeartMate II) under the same clinical support condition (flow rate of 4.5L/min, pressure head around 75 mmHg). A comprehensive multi-indicator evaluation method including hemodynamic parameters, hemolysis model, thrombotic potential model and bleeding probability model was used to analyze blood damage and assess the hemodynamic performance and hemocompatibility of these VADs.

RESULTS

Simulation results show that shear stress from 50 Pa to 100 Pa plays a major role in blood damage in Impella 5.0, UltraMag and CHVAD, while blood damage in HVAD and HeartMate II is mainly caused by shear stress greater than 100 Pa. Residence time was not the main factor for blood damage in Impella 5.0, and also makes a limited contribution to blood trauma in UltraMag and CHVAD, while it takes a critical role in elevating thrombotic potential in HVAD and HeartMate II. The distribution of regions of high hemolysis risk and high bleeding probability was similar for all these VADs and partially overlapped for high thrombotic potential regions. For Impella 5.0, regions with high hemolysis and bleeding risk were found mainly in the blade tip clearance and diffuser domains, high thrombotic potential regions were almost absent. For UltraMag, regions with high hemolysis, bleeding and thrombosis potential were found in two corners of the inlet pipe, the secondary flow passage, and the impeller eye. For CHVAD, the high-risk regions for hemolysis, bleeding and thrombosis are mainly in the inner side of the secondary flow passage and the middle region of the impeller passage. The narrow hydrodynamic clearance and impeller passage had a high risk of hemolysis and bleeding, and the clearance between the rotor and guide cone and the hydrodynamic clearance had high thrombotic potential. For HeartMate II, regions of high hemolysis risk and bleeding probability were found in the near-wall region of the straightener, the blade tip clearance and the diffuser domain. The corners of the inlet and outlet pipe and the straightener and diffuser regions had high thrombotic potential.

CONCLUSION

The risk of hemolysis, bleeding and thrombosis for these five VADs, in increasing order, was Impella 5.0, UltraMag, CHVAD, HVAD, and HeartMate II. Flow losses caused by the rotor mechanical movement, chaotic flow and narrow clearances increase the blood damage for all these VADs. The multi-indicator analysis can comprehensively evaluate the VAD performance with improved assessment accuracy of CFD.

A Continuum Model for the Unfolding of von Willebrand Factor

von Willebrand Factor is a mechano-sensitive protein circulating in blood that mediates platelet adhesion to subendothelial collagen and platelet aggregation at high shear rates. Its hemostatic function and thrombogenic effect, as well as susceptibility to enzymatic cleavage, are regulated by a conformational change from a collapsed globular state to a stretched state. Therefore, it is essential to account for the conformation of the vWF multimers when modeling vWF-mediated thrombosis or vWF degradation. We introduce a continuum model of vWF unfolding that is developed within the framework of our multi-constituent model of platelet-mediated thrombosis. The model considers two interconvertible vWF species corresponding to the collapsed and stretched conformational states. vWF unfolding takes place via two regimes: tumbling in simple shear and strong unfolding in flows with dominant extensional component. These two regimes were demonstrated in a Couette flow between parallel plates and an extensional flow in a cross-slot geometry. The vWF unfolding model was then verified in several microfluidic systems designed for inducing high-shear vWF-mediated thrombosis and screening for von Willebrand Disease. The model predicted high concentration of stretched vWF in key regions where occlusive thrombosis was observed experimentally. Strong unfolding caused by the extensional flow was limited to the center axis or middle plane of the channels, whereas vWF unfolding near the channel walls relied upon the shear tumbling mechanism. The continuum model of vWF unfolding presented in this work can be employed in numerical simulations of vWF-mediated thrombosis or vWF degradation in complex geometries. However, extending the model to 3-D arbitrary flows and turbulent flows will pose considerable challenges.