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.
BellumkondaLBondeP.Ventricular assist device therapy for heart failure-past, present, and future. Int Anesthesiol Clin2012; 50: 123–145.
2.
KawaboriMKuriharaCSugiuraT, et al. Continuousflow left ventricular assist device implantation in patients with a small left ventricle. Ann Thorac Surg2018; 105(3): 799–806.
3.
MagruderJTGrimmJCCrawfordTC, et al. Survival after orthotopic heart transplantation in patients undergoing bridge to transplantation with the HeartWare HVAD versus the Heartmate II. Ann Thorac Surg2016; 103: 1505–1511.
4.
CuiGAkersWJScottMJ, et al. Diagnosis of LVAD thrombus using a high-avidity fibrin-specific 99m Tc probe. Theranostics2018; 8(4): 1168–1179.
5.
GaffeyACChenCWChungJJ, et al. Is there a difference in bleeding after left ventricular assist device implant: centrifugal versus axial?J Cardiothorac Surg2018; 13(1): 22.
6.
ShahreyarMBob-ManuelTKhouzamRN, et al. Trends, predictors and outcomes of ischemic stroke and intracranial hemorrhage in patients with a left ventricular assist device. Ann Transl Med2018; 6(1): 5.
7.
MolinaEJBoyceSW.Current status of left ventricular assist device technology. Semin Thorac Cardiovasc Surg2013; 25(1): 56–63.
8.
MaltaisSKilicANathanS, et al. PREVENtion of HeartMate II pump thrombosis through clinical management: the PREVENT multi-center study. J Heart Lung Transplant2016; 36(1): 1–12.
9.
ShahPTantryUSBlidenKP, et al. Bleeding and thrombosis associated with ventricular assist device therapy. J Heart Lung Transplant2017; 36(11): 1164–1173.
10.
BirschmannIDittrichMEllerT, et al. Ambient hemolysis and activation of coagulation is different between HeartMate II and HeartWare left ventricular assist devices. J Heart Lung Transplant2014; 33(1): 80–87.
11.
NetukaIKvasnickaTKvasnickaJ, et al. Evaluation of von Willebrand factor with a fully magnetically levitated centrifugal continuous-flow left ventricular assist device in advanced heart failure. J Heart Lung Transplant2016; 35(7): 860–867.
12.
MehraMRNakaYUrielN, et al. A fully magnetically levitated circulatory pump for advanced heart failure. N Engl J Med2017; 376(5): 440–450.
13.
KormosRLCowgerJPaganiFD, et al. The Society of Thoracic Surgeons Intermacs database annual report: evolving indications, outcomes, and scientific partnerships. J Heart Lung Transplant2019; 38(2): 114–126.
14.
SonntagS.Numerical washout study of a pulsatile total. Int J Artif Organs2014; 37(3): 241–252.
15.
GoubergritsL.Numerical modeling of blood damage current status, challenges and prospects. J Med Devices2006; 3(5): 527–531.
16.
TaylorJOMeyerRSDeutschS, et al. Development of a computational model for macroscopic predictions of device-induced thrombosis. Biomech Model Mechanobiol2016; 15(6): 1713–1731.
17.
GoodmanPDBarlowETCrapoPM, et al. Computational model of device-induced thrombosis and thromboembolism. Ann Biomed Eng2005; 33(6): 780–797.
18.
WuWTJamiolkowskiMAWagnerWR, et al. Multiconstituent simulation of thrombus deposition. Sci Rep2017; 7: 42720.
19.
WuWTYangFWuJ, et al. High fidelity computational simulation of thrombus formation in Thoratec HeartMate II continuous flow ventricular assist device. Sci Rep2016; 6: 38025.
20.
TseKMChiuPLeeHP, et al. Investigation of hemodynamics in the development of dissecting aneurysm within patient-specific dissecting aneurismal aortas using computational fluid dynamics (CFD) simulations. J Biomech2011; 44(5): 827–836.
21.
Dillon-MurphyDNooraniANordslettenD, et al. Multi-modality image-based computational analysis of haemodynamics in aortic dissection. Biomech Model Mechanobiol2016: 15(4): 857–876.
22.
NautaFJLauKDArthursCJ, et al. Computational fluid dynamics and aortic thrombus formation following thoracic endovascular aortic repair. Ann Thorac Surg2017; 103(6): 1914–1921.
23.
MenichiniCChengZGibbsRG, et al. A computational model for false lumen thrombosis in type B aortic dissection following thoracic endovascular repair. J Biomech2018; 66: 36–43.
24.
GalloDSteinmanDAMorbiducciU.Insights into the co-localization of magnitude-based versus direction-based indicators of disturbed shear at the Carotid Bifurcation. J Biomech2016; 49(12): 2413–2419.
25.
BarkDLParaANKuDN.Correlation of thrombosis growth rate to pathological wall shear rate during platelet accumulation. Biotechnol Bioeng2012; 109(10): 2642–2650.
26.
MalaspinasOTurjmanARibeirodeSousaD, et al. A spatio-temporal model for spontaneous thrombus formation in cerebral aneurysms. J Theor Biol2016; 394: 68–76.
27.
XuHPiccinelliMLeshnowerBG, et al. Coupled morphological-hemodynamic computational analysis of type B aortic dissection: a longitudinal study. Ann Biomed Eng2018; 46(7): 927–939.
28.
GoubergritsLTimmelTAffeldK, et al. Characterization of an artificial valve flow using the numerical dye washout visualization technique: application to the monoleaflet valve with purged flow. Artif Organs2006; 30(8): 642–650.
29.
WuPBoehningFGross-HardtS, et al. On the accuracy of hemolysis models in Couette-type blood shearing devices. Artif Organs2018; 42: E290–E303.
30.
WiegmannLThamsenBZélicourtD, et al. Fluid dynamics in the HeartMate 3: influence of the artificial pulse feature and residual cardiac pulsation. Artif Organs2018; 43(4): 363–376.
31.
HenryFSShortlandAPIudicelloF, et al. Flow in a simple model skeletal muscle ventricle: comparison between numerical and physical simulations. J Biomech Eng1997; 19: 13–9.
32.
GoubergritsLKertzscherUAffeldK, Numerical dye washout method as a tool for characterizing the heart valve flow: study of three standard mechanical heart valves. ASAIO J2008; 54(1): 50–57.
33.
MolteniAMasriZLowK, et al. Experimental measurement and numerical modelling of dye washout for investigation of blood residence time in ventricular assist devices. Int J Artif Organs2018; 41(4): 201–212.
34.
MartorellJSantomaPKolandaiveluK, et al. Extent of flow recirculation governs expression of atherosclerotic and thrombotic biomarkers in arterial bifurcations. Cardiovasc Res2014; 103: 37–46.
35.
LiuGMChenHBHouJF, et al. Platelet adhesion emulation: a novel method for estimating the device thrombosis potential of a ventricular assist device. Int J Artif Organs2019; 43(4): 252–257.
36.
FerzigerJHPericM.Computational methods for fluid dynamics. Berlin: Springer, 2002.
37.
MerrillEW.Rheology of blood. Physiol Rev1969; 49: 836–886.
