Hemodynamic analysis of outflow grafting positions of a ventricular assist device using closed-loop multiscale CFD simulations: Preliminary results

Computational investigation of outflow graft variation impact on hemocompatibility profile in LVADs

BACKGROUND

Hemocompatibility-related adverse events (HRAE) occur commonly in patients with left ventricular assist devices (LVADs) and add to morbidity and mortality. It is unclear whether the outflow graft orientation can impact flow conditions leading to HRAE. This study presents a simulation-based approach using exact patient anatomy from medical images to investigate the influence of outflow cannula orientation in modulating flow conditions leading to HRAEs.

METHODS

A 3D model of a proximal aorta and outflow graft was reconstructed from a computed tomography (CT) scan of an LVAD patient and virtually modified to model multiple cannula orientations (n = 10) by varying polar (cranio-caudal) (n = 5) and off-set (anterior-posterior) (n = 2) angles. Time-dependent computational flow simulations were then performed for each anatomical orientation. Qualitative and quantitative hemodynamics metrics of thrombogenicity including time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), endothelial cell platelet activation potential (ECAP), particle residence time (PRT), and platelet activation potential (PLAP) were analyzed.

RESULTS

Within the simulations performed, endothelial cell activation potential (ECAP) and particle residence time (PRT) were found to be lowest with a polar angle of 85°, regardless of offset angle. However, polar angles that produced parameters at levels least associated with thrombosis varied when the offset angle was changed from 0° to 12°. For offset angles of 0° and 12° respectively, flow shear was lowest at 65° and 75°, time averaged wall shear stress (TAWSS) was highest at 85° and 35°, and platelet activation potential (PLAP) was lowest at 65° and 45°.

CONCLUSION

This study suggests that computational fluid dynamic modeling based on patient-specific anatomy can be a powerful analytical tool when identifying optimal positioning of an LVAD. Contrary to previous work, our findings suggest that there may be an "ideal" outflow cannula for each individual patient based on a CFD-based hemocompatibility profile.

Dynamic characteristic modeling of left ventricular assist devices based on hysteresis effects

OBJECTIVE

The purpose of this study is to develop a new model for the dynamic characteristics of left ventricular assist devices (LVADs) interacting with the cardiovascular system under constant-speed modes.

METHODS

A new hysteresis model is established on the basis of the hysteresis effect and turbomachinery principles. The simulation results from the hysteresis model were compared with the inertia model. The in-vitro experiment results of a centrifugal pump (from literature) and the unsteady computational fluid dynamics (CFD) simulation results of an axial pump were used as the benchmarks.

RESULTS

Compared with the inertia model, at the partial support mode, the relative estimation error of the time to the maximum and minimum pump flow (Q) in the hysteresis model decreased at least 16.3% cardiac cycle (T_c) in the centrifugal pump and at least 1.9% T_c in the axial pump, indicating its ability to simulate more realistic Q fluctuations. Moreover, the hysteresis model could predict an accurate time distribution of different Q.

CONCLUSION

The hysteresis model provides a general calculation method for simulating the dynamic characteristics of constant-speed LVADs under interaction with the cardiovascular system. It is more accurate than the inertia model.

SIGNIFICANCE

The hysteresis model is helpful for the rapid estimation of unsteady dynamic characteristics in absence of a physical pump prototype at the preliminary design stage.

Asymmetric distribution of pathogenic low wall shear stress of the bilateral subclavian arteries: two case reports

The effects of increasing blood flow on the pathogenic wall shear stress (pWSS) of subclavian arteries (SAs) are currently unclear. Patient-specific models of the SA were constructed based on computed tomographic images from two patients. Using the Ansys Fluent 19.0 transient laminar flow solver, the finite volume method was chosen to solve the Navier–Stokes equation governing fluid behavior. The time-averaged wall shear stress, ratio of risk area, cumulative ratio of risk area (P¯), ratio of risk time, and ratio contour of risk time were calculated to describe the temporal and spatial distributions of pWSS. Virtually all pWSS occurred during the diastolic phase. The P¯ was 2.3 and 1.29 times higher on the left than on the right in Patients 1 (P1) and 2 (P2), respectively. Increasing the blood flow volume of the left SA by 20%, 40%, and 60% led to a 9.27%, 15.10%, and 20.99% decrease in P¯ for P1 and a 5.74%, 11.55%, and 17.14% decrease in P¯ for P2, respectively, compared with baseline values. In conclusion, the left SA showed greater diastolic pWSS than the right SA, and increasing the blood flow volume reduced the pWSS in the left SA.

Understanding the influence of left ventricular assist device inflow cannula alignment and the risk of intraventricular thrombosis

Background

Adverse neurological events associated with left ventricular assist devices (LVADs) have been suspected to be related to thrombosis. This study aimed to understand the risks of thrombosis with variations in the implanted device orientation. A severely dilated pulsatile patient-specific left ventricle, modelled with computational fluid dynamics, was utilised to identify the risk of thrombosis for five cannulation angles. With respect to the inflow cannula axis directed towards the mitral valve, the other angles were 25° and 20° towards the septum and 20° and 30° towards the free wall.

Results

Inflow cannula angulation towards the free wall resulted in longer blood residence time within the ventricle, slower ventricular washout and reduced pulsatility indices along the septal wall. Based on the model, the ideal inflow cannula alignment to reduce the risk of thrombosis was angulation towards the mitral valve and up to parallel to the septum, avoiding the premature clearance of incoming blood.

Conclusions

This study indicates the potential effects of inflow cannulation angles and may guide optimised implantation configurations; however, the ideal approach will be influenced by other patient factors and is suspected to change over the course of support.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12938-021-00884-6.

Computational analyses of aortic blood flow under varying speed CF-LVAD support

Continuous Flow Left Ventricular Assist Devices (CF-LVADs) generally operate at a constant speed whilst supporting a failing heart. However, constant speed CF-LVAD support may cause complications and increase the morbidity rates in the patients. Therefore, different varying speed operating modes for CF-LVADs have been proposed to generate more physiological blood flow, which may reduce complication rates under constant speed CF-LVAD support. The proposed varying speed CF-LVAD algorithms simulate time-dependant dynamics and three dimensional blood flow patterns in aorta under varying speed CF-LVAD support remain unclear. The aim of this study is to evaluate three dimensional blood flow patterns in a patient-specific aorta model under co-pulsating and counter-pulsating CF-LVAD support modes driven by speed and flow rate control algorithms using numerical simulations. Aortic blood flow was evaluated for 10,000 rpm constant speed CF-LVAD support generating 4.71 L/min mean flow rate over a cardiac cycle. Co-pulsating and counter-pulsating CF-LVAD speed control operated the pump at the same average speed over a cardiac cycle and co-pulsating and counter-pulsating CF-LVAD flow rate control generated the same average flow rate over cardiac cycle as in the constant speed pump support. Simulation results show that the utilised counter-pulsating pump flow rate control may decrease the haemolysis to a third compared to the most commonly employed constant speed pump operating mode. Moreover, CF-LVAD support utilising counter-pulsating pump flow rate control generated the most favourable hemodynamic characteristics, i.e. low Dean number, least wall shear stress and least haemolysis values among the investigated cases.

How to evaluate the outflow tract of LVAD after minimally invasive implantation by 3D CT-scan

During a minimally invasive implantation technique, the outflow graft of left ventricular assist device (LVAD) is tunnelled blindly through the pericardium or left pleura, with an inability to assess for twisting or malposition. Three-dimensional computed tomography scan (CT-scan) has a role in qualitative evaluation of the different outflow tract configurations. The different surgical minimally invasive approaches include: (a) mini-sternotomy and left mini-thoracotomy, (b) right mini-thoracotomy and left mini-thoracotomy, (c) subclavian artery access and left mini-thoracotomy. The outflow graft could be anastomosed to the left axillary artery or the ascending aorta. CT-scan reconstruction using syngo InSpace4D (Siemens, Muenchen, Germany) was used to provide fast segmentation and high-resolution images. The 3D reconstructions permit an evaluation of different anastomosis configurations and to assess the route of outflow graft.

Innovative Modeling Techniques and 3D Printing in Patients with Left Ventricular Assist Devices: A Bridge from Bench to Clinical Practice

Left ventricular assist devices (LVAD) cause altered flow dynamics that may result in complications such as stroke, pump thrombosis, bleeding, or aortic regurgitation. Understanding altered flow dynamics is important in order to develop more efficient and durable pump configurations. In patients with LVAD, hemodynamic assessment is limited to imaging techniques such as echocardiography which precludes detailed assessment of fluid dynamics. In this review article, we present some innovative modeling techniques that are often used in device development or for research purposes, but have not been utilized clinically. Computational fluid dynamic (CFD) modeling is based on computer simulations and particle image velocimetry (PIV) employs ex vivo models that helps study fluid characteristics such as pressure, shear stress, and velocity. Both techniques may help elaborate our understanding of complications that occur with LVAD and could be potentially used in the future to troubleshoot LVAD-related alarms. These techniques coupled with 3D printing may also allow for patient-specific device implants, lowering the risk of complications increasing device durability.

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.

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.

Cannula and Pump Positions Are Associated With Left Ventricular Unloading and Clinical Outcome in Patients With HeartWare Left Ventricular Assist Device

BACKGROUND

Cannula and pump positions are associated with clinical outcomes such as device thrombosis in patients with HeartMate II; however, clinical implications of HVAD (HeartWare International, Framingham, Massachusetts) cannula position are unknown. This study aims to assess the relationship among cannula position, left ventricular (LV) unloading, and patient prognosis.

METHODS AND RESULTS

Twenty-seven HVAD patients (60.0 ± 12.6 years of age and 19 males [70%]) underwent ramp test. Device position was quantified from chest X-ray parameters obtained at the time of the hemodyamic ramp test: (1) cannula coronal angle, (2) pump depth, (3) cannula sagittal angle, and (4) pump area. Lower cannula coronal angle was associated with LV unloading (as measured by smaller LV diastolic dimension and lower pulmonary capillary wedge pressure). Smaller pump area was associated with LV dynamic unloading, as assessed by steeper negative slopes of LV diastolic dimension and pulmonary capillary wedge pressure during incremental rotational speed change. Cannula coronal angle ≤65° was associated with reduced heart failure readmission rate (hazard ratio, 10.33; P = .007 by log-rank test).

CONCLUSION

HVAD cannula and pump positions are associated with LV unloading and improved clinical outcomes. Prospective studies evaluating surgical techniques to ensure optimal device positioning and its effects on clinical outcomes are warranted.

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.

Influence of IABP-Induced Abdominal Occlusions on Aortic Hemodynamics: A Patient-Specific Computational Evaluation

Intraaortic balloon pump (IABP) is used as temporary mechanical assistance in case of cardiovascular diseases, even if different hemodynamic problems and, thus, clinical complications may happen, such as the decrease of visceral perfusion. A computational fluid dynamic (CFD) study was carried out to investigate the effects of different IABP-induced abdominal occlusions on patient-specific aortic flow. Two possible sizes (25 and 34 cm) and two locations (2 and 3 cm) of the balloon were compared, modeling four abdominal occlusions and numerically reproducing IAB inflation/deflation behavior. The results highlighted that the perfusion in renal, mesenteric, and iliac arteries decreases when the abdominal occlusion increases with balloon inflation. The study illustrates also how the balloon size affects the flow in aorta vessels in both locations, and that the positioning is of little relevance for the 34 cm balloon, whereas it influences the aortic flow very much in case of 25 cm IAB. This analysis demonstrates how the IAB-induced occlusion may vary the abdominal circulation; therefore, the correct size and positioning are emphasized for patient's outcome.