Numerical investigation of the effect of cannula placement on thrombosis

Smoothed particle hydrodynamics based FSI simulation of the native and mechanical heart valves in a patient-specific aortic model

The failure of the aortic heart valve is common, resulting in deterioration of the pumping function of the heart. For the end stage valve failure, bi-leaflet mechanical valve (most popular artificial valve) is implanted. However, due to its non-physiological behaviour, a significant alteration is observed in the normal haemodynamics of the aorta. While in-vivo experimentation of a human heart valve (native and artificial) is a formidable task, in-silico study using computational fluid dynamics (CFD) with fluid structure interaction (FSI) is an effective and economic tool for investigating the haemodynamics of natural and artificial heart valves. In the present work, a haemodynamic model of a natural and mechanical heart valve has been developed using meshless particle-based smoothed particle hydrodynamics (SPH). In order to further enhance its clinical relevance, this study employs a patient-specific vascular geometry and presents a successful validation against traditional finite volume method and 4D magnetic resonance imaging (MRI) data. The results have demonstrated that SPH is ideally suited to simulate the heart valve function due to its Lagrangian description of motion, which is a favourable feature for FSI. In addition, a novel methodology for the estimation of the wall shear stress (WSS) and other related haemodynamic parameters have been proposed from the SPH perspective. Finally, a detailed comparison of the haemodynamic parameters has been carried out for both native and mechanical aortic valve, with a particular emphasis on the clinical risks associated with the mechanical valve.

Left ventricular assist device and pump thrombosis: the importance of the inflow cannula position

Pump thrombosis is a devastating complication after left ventricular assist device implantation. This study aims to elucidate the relation between left ventricular assist device implantation angle and risk of pump thrombosis. Between November 2010 and March 2020, 53 left ventricular assist device-patients underwent a computed tomography scan. Using a 3-dimensional multiplanar reformation the left ventricular axis was reconstructed to measure the implantation angle of the inflow cannula. All patients were retrospectively analyzed for the occurrence of pump thrombosis. In 10 (91%) patients with a pump thrombosis, the implantation angle was towards the lateral wall of the left ventricle. In only 20 patients (49%) of the patients without a pump thrombosis the inflow cannula pointed towards the lateral wall of the left ventricle. The mean angle in patients with a pump thrombosis was 10.1 ± 11.9 degrees towards the lateral wall of the left ventricle compared to 4.1 ± 19.9 degrees towards the septum in non-pump thrombosis patients (P = 0.005). There was a trend towards a significant difference in time to first pump thrombosis between patients with a lateral or septal deviated left ventricular assist device (hazard ratio of 0.15, P = 0.07). This study demonstrates that left ventricular assist device implantation angle is associated with pump thrombosis. Almost all patients in whom a pump thrombosis occurred during follow-up had a left ventricular assist device implanted with the inflow-cannula pointing towards the lateral wall of the left ventricle.

Numerical study of the effect of LVAD inflow cannula positioning on thrombosis risk

Left ventricular assist devices (LVADs) have been increasingly used as a therapy for patients with end-stage heart failure. However, a growing number of clinical observations have shown that LVADs are associated with thromboembolic events, which are potentially related to the changes in intraventricular flow. Particularly, the flow fields around the inflow cannula (IC) of the LVAD. In this study, a fluid structure interaction (FSI) simulation was conducted to evaluate the hemodynamics of a patient specific left ventricle (LV) with varying LVAD IC orientations. The LV model was obtained from computed tomography scans and modeled to have contraction and relaxation during cardiac cycles following available experimental data of LV volume changes. The LV of the patient was assumed to have an end systolic volume of 223.7 mL and a stroke volume of 46.4 mL. Four different IC positions were considered: towards the (1) septum; (2) aortic valve (AV); (3) mitral valve (MV) and (4) inferior wall (IW). The potential thrombus growth around the IC was assumed to be caused by blood stagnation regions with low velocity (<5 mm/s) and low shear rate (<60/s) flow. Mean velocity magnitudes and low blood velocity regions around the IC were numerically obtained. To quantitatively compare the thrombosis risks of the four simulation cases, the time-averaged volumes of the low-velocity regions and the low shear rate regions were calculated. The intraventricular volumes of low velocity zones based on IC orientation are 1.42 mL toward the septum, 1.14 mL toward the AV, 0.93 mL toward the MV, and 1.24 mL toward the IW. The intraventricular volumes of low shear regions based on IC orientation are 11.54 mL toward the septum, 11.15 mL toward the AV, 9.24 mL toward the MV, and 10.7 mL toward the IW. IC orientation toward the MV results in lower volumetric regions of low flow and low shear within the ventricle, which consequently may lead to a reduced risk of thrombus formation.

Numerical investigation of the effects of cannula geometry on hydraulic blood flow to prevent the risk of thrombosis

Despite significant advances in left ventricular assist devices and the cannula, unfavorable events leading to the death of patients, including bleeding, infection, neurological disorders, hemolysis, and thrombosis, are still being reported. Local parameters of blood flow, including static flow, vorticity and critical values of shear stress on the wall of ventricle and cannula, increase the risk of thrombosis. Therefore, the analysis of blood flow domains inside the ventricle and cannula is necessary to investigate the probability of forming thrombosis in the cannula of left ventricular assist devices. In this study, blood flow is investigated in a Medtronic DLP 16F clinical cannula by using computational fluid dynamics through three-dimensional modeling of the left ventricle and cannula based on real geometry. Apart from the fact that blood is considered non-Newtonian fluid, the effect of heart movement in the left ventricle is also applied. In this research, blood flow in the cannula has been examined and some problems resulting from the use of the cannulas have been investigated. The results indicate that changing the geometry of input holes, such as their number and size, on the tip of the cannula, alter the probability of forming thrombosis and the standard mode shows a better performance.

Small Left Ventricular Size Is an Independent Risk Factor for Ventricular Assist Device Thrombosis

The prevalence of ventricular assist device (VAD) therapy has continued to increase due to a stagnant donor supply and growing advanced heart failure (HF) population. We hypothesize that left ventricular (LV) size strongly influences biocompatibility and risk of thrombosis. Unsteady computational fluid dynamics (CFD) was used in conjunction with patient-derived computational modeling and virtual surgery with a standard, apically implanted inflow cannula. A dual-focus approach of evaluating thrombogenicity was employed: platelet-based metrics to characterize the platelet environment and flow-based metrics to investigate hemodynamics. Left ventricular end-diastolic dimensions (LVEDds) ranging from 4.5 to 6.5 cm were studied and ranked according to relative thrombogenic potential. Over 150,000 platelets were individually tracked in each LV model over 15 cardiac cycles. As LV size decreased, platelets experienced markedly increased shear stress histories (SHs), whereas platelet residence time (RT) in the LV increased with size. The complex interplay between increased SH and longer RT has profound implications on thrombogenicity, with a significantly higher proportion of platelets in small LVs having long RT times and being subjected to high SH, contributing to thrombus formation. Our data suggest that small LV size, rather than decreased VAD speed, is the primary pathologic mechanism responsible for the increased incidence of thrombosis observed in VAD patients with small LVs.

Magnetic Resonance Imaging Based Flow Field and Lagrangian Particle Tracking From a Left Ventricular Assist Device

This study explores the optimal left ventricular assist device (LVAD) cannula outflow configuration in a patient-specific replica of the aorta. The volumetric velocity field is measured using phase-contrast magnetic resonance imaging (PC-MRI) under a physiologically relevant steady flow. The effect of the LVAD outflow graft insertion site and anastomosis angle on the transport of embolic particles to cranial vessels is studied by solving the particle equation of motion for spheres in the range of 0.1-1.0 mm using the measured three-dimensional (3D) velocity field. Results show that for a given aorta anatomy, it is possible to design the cannula graft location and terminal curvature so that the probability of embolic transport to the cranial vessels is significantly minimized. This is particularly important since the complex flow pattern in each cannula case affects the embolic trajectories differently, and hence the common assumption that particles distribute by the volumetric flow division does not hold.

Left Ventricular Assist Device Inflow Cannula Angle and Thrombosis Risk

BACKGROUND

As heart failure prevalence continues to increase in the setting of a static donor supply, left ventricular assist device (LVAD) therapy for end-stage heart failure continues to grow. Anecdotal evidence suggests that malalignment of the LVAD inflow cannula may increase thrombosis risk, but this effect has not been explored mechanistically or quantified statistically. Our objective is to elucidate the impact of surgical angulation of the inflow cannula on thrombogenicity.

METHODS AND RESULTS

Unsteady computational fluid dynamics is used in conjunction with computational modeling and virtual surgery to model flow through the left ventricle for 5 different inflow cannula angulations. We use a holistic approach to evaluate thrombogenicity: platelet-based (Lagrangian) metrics to evaluate the platelet mechanical environment, combined with flow-based (Eulerian) metrics to investigate intraventricular hemodynamics. The thrombogenic potential of each LVAD inflow cannula angulation is quantitatively evaluated based on platelet shear stress history and residence time. Intraventricular hemodynamics are strongly influenced by LVAD inflow cannula angulation. Platelet behavior indicates elevated thrombogenic potential for certain inflow cannula angles, potentially leading to platelet activation. Our analysis demonstrates that the optimal range of inflow angulation is within 0±7° of the left ventricular apical axis.

CONCLUSIONS

Angulation of the inflow cannula >7° from the apical axis (axis connecting mitral valve and ventricular apex) leads to markedly unfavorable hemodynamics as determined by computational fluid dynamics. Computational hemodynamic simulations incorporating Lagrangian and Eulerian metrics are a powerful tool for studying optimization of LVAD implantation strategies, with the long-term potential of improving outcomes.

Blood damage in Left Ventricular Assist Devices: Pump thrombosis or system thrombosis?

Introduction:

Despite significant technical advancements in the design and manufacture of Left Ventricular Assist Devices, post-implant thrombotic and thromboembolic complications continue to affect long-term outcomes. Previous efforts, aimed at optimizing pump design as a means of reducing supraphysiologic shear stresses generated within the pump and associated prothrombotic shear-mediated platelet injury, have only partially altered the device hemocompatibility.

Methods:

We examined hemodynamic mechanisms that synergize with hypershear within the pump to contribute to the thrombogenic potential of the overall Left Ventricular Assist Device system.

Results:

Numerical simulations of blood flow in differing regions of the Left Ventricular Assist Device system, that is the diseased native left ventricle, the pump inflow cannula, the impeller, the outflow graft and the anastomosed downstream aorta, reveal that prothrombotic hemodynamic conditions might occur at these specific sites. Furthermore, we show that beyond hypershear, additional hemodynamic abnormalities exist within the pump, which may elicit platelet activation, such as recirculation zones and stagnant platelet trajectories. We also provide evidences that particular Left Ventricular Assist Device implantation configurations and specific post-implant patient management strategies, such as those allowing aortic valve opening, are more hemodynamically favorable and reduce the thrombotic risk.

Conclusion:

We extend the perspective of pump thrombosis secondary to the supraphysiologic shear stress environment of the pump to one of Left Ventricular Assist Device system thrombosis, raising the importance of comprehensive characterization of the different prothrombotic risk factors of the total system as the target to achieve enhanced hemocompatibility and improved clinical outcomes.

A Multiphysics Biventricular Cardiac Model: Simulations With a Left-Ventricular Assist Device

Computational models have become essential in predicting medical device efficacy prior to clinical studies. To investigate the performance of a left-ventricular assist device (LVAD), a fully-coupled cardiac fluid-electromechanics finite element model was developed, incorporating electrical activation, passive and active myocardial mechanics, as well as blood hemodynamics solved simultaneously in an idealized biventricular geometry. Electrical activation was initiated using a simplified Purkinje network with one-way coupling to the surrounding myocardium. Phenomenological action potential and excitation-contraction equations were adapted to trigger myocardial contraction. Action potential propagation was formulated within a material frame to emulate gap junction-controlled propagation, such that the activation sequence was independent of myocardial deformation. Passive cardiac mechanics were governed by a transverse isotropic hyperelastic constitutive formulation. Blood velocity and pressure were determined by the incompressible Navier-Stokes formulations with a closed-loop Windkessel circuit governing the circulatory load. To investigate heart-LVAD interaction, we reduced the left ventricular (LV) contraction stress to mimic a failing heart, and inserted a LVAD cannula at the LV apex with continuous flow governing the outflow rate. A proportional controller was implemented to determine the pump motor voltage whilst maintaining pump motor speed. Following LVAD insertion, the model revealed a change in the LV pressure-volume loop shape from rectangular to triangular. At higher pump speeds, aortic ejection ceased and the LV decompressed to smaller end diastolic volumes. After multiple cycles, the LV cavity gradually collapsed along with a drop in pump motor current. The model was therefore able to predict ventricular collapse, indicating its utility for future development of control algorithms and pre-clinical testing of LVADs to avoid LV collapse in recipients.

Ventricular flow dynamics with varying LVAD inflow cannula lengths: In-silico evaluation in a multiscale model

Left ventricular assist devices are associated with thromboembolic events, which are potentially caused by altered intraventricular flow. Due to patient variability, differences in apical wall thickness affects cannula insertion lengths, potentially promoting unfavourable intraventricular flow patterns which are thought to be correlated to the risk of thrombosis. This study aimed to present a 3D multiscale computational fluid dynamic model of the left ventricle (LV) developed using a commercial software, Ansys, and evaluate the risk of thrombosis with varying inflow cannula insertion lengths in a severely dilated LV. Based on a HeartWare HVAD inflow cannula, insertion lengths of 5, 19, 24 and 50 mm represented cases of apical hypertrophy, typical ranges of apical thicknesses and an experimental length, respectively. The risk of thrombosis was evaluated based on blood washout, residence time, instantaneous blood stagnation and a pulsatility index. By introducing fresh blood to displace pre-existing blood in the LV, after 5 cardiac cycles, 46.7%, 45.7%, 45.1% and 41.8% of pre-existing blood remained for insertion lengths of 5, 19, 24 and 50 mm, respectively. Compared to the 50 mm insertion, blood residence time was at least 9%, 7% and 6% higher with the 5, 19 and 24 mm insertion lengths, respectively. No instantaneous stagnation at the apex was observed directly after the E-wave. Pulsatility indices adjacent to the cannula increased with shorter insertion lengths. For the specific scenario studied, a longer insertion length, relative to LV size, may be advantageous to minimise thrombosis by increasing LV washout and reducing blood residence time.

Intermittent Aortic Valve Opening and Risk of Thrombosis in Ventricular Assist Device Patients

The current study evaluates quantitatively the impact that intermittent aortic valve (AV) opening has on the thrombogenicity in the aortic arch region for patients under left ventricular assist device (LVAD) therapy. The influence of flow through the AV, opening once every five cardiac cycles, on the flow patterns in the ascending aortic is measured in a patient-derived computed tomography image-based model, after LVAD implantation. The mechanical environment of flowing platelets is investigated, by statistical treatment of outliers in Lagrangian particle tracking, and thrombogenesis metrics (platelet residence times and activation state characterized by shear stress accumulation) are compared for the cases of closed AV versus intermittent AV opening. All hemodynamics metrics are improved by AV opening, even at a reduced frequency and flow rate. Residence times of platelets or microthrombi are reduced significantly by transvalvular flow, as are the shear stress history experienced and the shear stress magnitude and gradients on the aortic root endothelium. The findings of this device-neutral study support the multiple advantages of management that enables AV opening, providing a rationale for establishing this as a standard in long-term treatment and care for advanced heart failure patients.

LVAD Outflow Graft Angle and Thrombosis Risk

This study quantifies thrombogenic potential (TP) of a wide range of left ventricular assist device (LVAD) outflow graft anastomosis angles through state-of-the-art techniques: 3D imaged-based patient-specific models created via virtual surgery and unsteady computational fluid dynamics with Lagrangian particle tracking. This study aims at clarifying the influence of a single parameter (outflow graft angle) on the thrombogenesis associated with flow patterns in the aortic root after LVAD implantation. This is an important and poorly-understood aspect of LVAD therapy, because several studies have shown strong inter and intrapatient thrombogenic variability and current LVAD implantation strategies do not incorporate outflow graft angle optimization. Accurate platelet-level investigation, enabled by statistical treatment of outliers in Lagrangian particle tracking, demonstrates a strong influence of outflow graft anastomoses angle on thrombogenicity (platelet residence times and activation state characterized by shear stress accumulation) with significantly reduced TP for acutely-angled anastomosed outflow grafts. The methodology presented in this study provides a device-neutral platform for conducting comprehensive thrombogenicity evaluation of LVAD surgical configurations, empowering optimal patient-focused surgical strategies for long-term treatment and care for advanced heart failure patients.

Tetralogy of Fallot Surgical Repair: Shunt Configurations, Ductus Arteriosus and the Circle of Willis

In this study, hemodynamic performance of three novel shunt configurations that are considered for the surgical repair of tetralogy of Fallot (TOF) disease are investigated in detail. Clinical experience suggests that the shunt location, connecting angle, and its diameter can influence the post-operative physiology and the neurodevelopment of the neonatal patient. An experimentally validated second order computational fluid dynamics (CFD) solver and a parametric neonatal diseased great artery model that incorporates the ductus arteriosus (DA) and the full patient-specific circle of Willis (CoW) are employed. Standard truncated resistance CFD boundary conditions are compared with the full cerebral arterial system, which resulted 21, -13, and 37% difference in flow rate at the brachiocephalic, left carotid, and subclavian arteries, respectively. Flow splits at the aortic arch and cerebral arteries are calculated and found to change with shunt configuration significantly for TOF disease. The central direct shunt (direct shunt) has pulmonary flow 5% higher than central oblique shunt (oblique shunt) and 23% higher than modified Blalock Taussig shunt (RPA shunt) while the DA is closed. Maximum wall shear stress (WSS) in the direct shunt configuration is 9 and 60% higher than that of the oblique and RPA shunts, respectively. Patent DA, significantly eliminated the pulmonary flow control function of the shunt repair. These results suggests that, due to the higher flow rates at the pulmonary arteries, the direct shunt, rather than the central oblique, or right pulmonary artery shunts could be preferred by the surgeon. This extended model introduced new hemodynamic performance indices for the cerebral circulation that can correlate with the post-operative neurodevelopment quality of the patient.

Numerical prediction of thrombus risk in an anatomically dilated left ventricle: the effect of inflow cannula designs

BACKGROUND

Implantation of a rotary blood pump (RBP) can cause non-physiological flow fields in the left ventricle (LV) which may trigger thrombosis. Different inflow cannula geometry can affect LV flow fields. The aim of this study was to determine the effect of inflow cannula geometry on intraventricular flow under full LV support in a patient specific model.

METHODS

Computed tomography angiography imaging of the LV was performed on a RBP candidate to develop a patient-specific model. Five inflow cannulae were evaluated, which were modelled on those used clinically or under development. The inflow cannulae are described as a crown like tip, thin walled tubular tip, large filleted tip, trumpet like tip and an inferiorly flared cannula. Placement of the inflow cannula was at the LV apex with the central axis intersecting the centre of the mitral valve. Full support was simulated by prescribing 5 l/min across the mitral valve. Thrombus risk was evaluated by identifying regions of stagnation. Rate of LV washout was assessed using a volume of fluid model. Relative haemolysis index and blood residence time was calculated using an Eulerian approach.

RESULTS

The inferiorly flared inflow cannula had the lowest thrombus risk due to low stagnation volumes. All cannulae had similar rates of LV washout and blood residence time. The crown like tip and thin walled tubular tip resulted in relatively higher blood damage indices within the LV.

CONCLUSION

Changes in intraventricular flow due to variances in cannula geometry resulted in different stagnation volumes. Cannula geometry does not appreciably affect LV washout rates and blood residence time. The patient specific, full support computational fluid dynamic model provided a repeatable platform to investigate the effects of inflow cannula geometry on intraventricular flow.

Virtual implantation of a novel LVAD: toward computer-assisted surgery for heart failure

BACKGROUND

Mechanical and hemodynamic factors are among the determinants of patient-device interaction and early-term and long-term outcomes in left ventricular assist device (LVAD) recipients.

MATERIAL AND METHODS

We are currently developing computer simulation tools aimed at (1) analyze the intrathoracic and intracavitary positioning of LVADs after implantation and establish correlation with clinical outcomes; (2) assist surgeons in the choice of device and of left ventricular coring site for optimized intrathoracic placement and function; and (3) facilitate the planning of less-invasive LVAD implantation. A virtual representation of LVAD (mesh device component) was created through cone-beam computed tomography and semiautomatic segmentation. A modular framework software (CamiTK, Grenoble, France) was used to create a three-dimensional representation of patients' computed tomography (CT) scan and incorporate the mesh device component for virtual implantation.

RESULTS

Device reconstruction was included into a dedicated software with the purposes of virtual implantation, based on the preoperative CT scan of surgical candidates.

CONCLUSIONS

We present herein the first digital reconstruction of the novel HeartMate 3 LVAD. Virtual implantation on the basis of preoperative CT scan is feasible within a user-friendly interactive software. Future studies will be focused on correlation with clinical variables.

Effect of Rotary Blood Pump Pulsatility on Potential Parameters of Blood Compatibility and Thrombosis in Inflow Cannula Tips

Purpose

Rotary Blood Pump (RBP) pulsatile strategies relative to the native cardiac cycle have been widely studied because of their benefits to hemodynamics. However, the effects that inducing pulses has on the blood compatibility of ventricular assist device (VAD) support have not yet been understood. Inflow cannulae have been found to be associated with thrombosis under conventional constant speed support of RBPs. To prevent further risks to blood compatibility, it is necessary to understand the relationship between cannula tip design and the induced pulsatility. The purpose of this study was to evaluate the flow field of 5 different tip geometries under RBP pulsatile support using stereo-particle image velocimetry (PIV).

Methods

Inflow cannulae with conventional tip geometries (blunt, blunt with 4 side ports, beveled with 3 side ports, and cage) and a custom designed crown tip were studied. All cannulae were interposed between a mixed-flow RBP and a silicone left ventricle. The contractile function and hemodynamics were reproduced in a mock circulation loop (MCL). The RBP was configured to induce synchronous and counter-synchronous pulses relative to cardiac cycles while supporting the failing ventricle.

Results

Between both pulsing strategies, low shear volume ([Formula: see text], potential parameter of thrombus formation) showed no significant difference. However, counter-synchronous pulsatile mode induced less increase of both high shear volume ([Formula: see text], potential parameter of platelet activation) and recirculation volume (Vz>0, potential parameter of thrombus formation).

Conclusions

Although the clinical relationship cannot be inferred from this measurement, when considering the inflow tips only, a necessary trade-off should be made between adverse effects on blood compatibility and benefits for hemodynamics during RBP pulsatile mode.

Tortuosity of coronary bifurcation as a potential local risk factor for atherosclerosis: CFD steady state study based on in vivo dynamic CT measurements

The purpose of the present study was to determine whether in vivo bifurcation geometric factors would permit prediction of the risk of atherosclerosis. It is worldwide accepted that low or oscillatory wall shear stress (WSS) is a robust hemodynamic factor in the development of atherosclerotic plaque and has a strong correlation with the local site of plaque deposition. However, it still remains unclear how coronary bifurcation geometries are correlated with such hemodynamic forces. Computational fluid dynamics simulations were performed on left main (LM) coronary bifurcation geometries derived from CT of eight patients without significant atherosclerosis. WSS amplitudes were accurately quantified at two high risk zones of atherosclerosis, namely at proximal left anterior descending artery (LAD) and at proximal left circumflex artery (LCx), and also at three high WSS concentration sites near the bifurcation. Statistical analysis was used to highlight relationships between WSS amplitudes calculated at these five zones of interest and various geometric factors. The tortuosity index of the LM-LAD segment appears to be an emergent geometric factor in determining the low WSS amplitude at proximal LAD. Strong correlations were found between the high WSS amplitudes calculated at the endothelial regions close to the flow divider. This study not only demonstrated that CT imaging studies of local risk factor for atherosclerosis could be clinically performed, but also showed that tortuosity of LM-LAD coronary branch could be used as a surrogate marker for the onset of atherosclerosis.