Computational fluid dynamics analysis of thrombosis potential in left ventricular assist device drainage cannulae

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.

Anti-Clogging Effect of Continuous Hemodiafiltration With Blood Recirculation

Pediatric dialysis requires low flow from the body, but greater flow is needed to prevent clogging. As a solution, we developed a new continuous hemodiafiltration system with blood recirculation (CHDF-R), which enables separate settings for blood flow from the body and to the hemofilter. We compared CHDF-R with conventional continuous hemodiafiltration (CHDF) of bovine plasma and blood by monitoring the transmembrane pressure (TMP) and observing the hemofilter membrane surface. When using bovine plasma, the postdialysis TMP with CHDF-R was significantly lower than with CHDF (median CHDF, 23.7; median CHDF-R, 18.1; p = 0.029). Likewise, when using bovine blood, the postdialysis TMP was also significantly lower with CHDF-R than with CHDF (median CHDF, 150; median CHDF-R, 100; p = 0.029). Moreover, the area of clogged membrane was significantly smaller with CHDF-R than with CHDF, and the inner membrane surface showed less material deposition with CHDF-R than CHDF. CHDF-R thus appears to suppress accumulation of clogging substances by producing higher shear stress within hollow fiber membranes.

Risk factors for cannula-associated arterial thrombosis following extracorporeal membrane oxygenation support: a retrospective study

BACKGROUND

Hemostatic dysfunction during extracorporeal membrane oxygenation (ECMO) due to blood-circuit interaction and the consequences of shear stress imposed by flow rates lead to rapid coagulation cascade and thrombus formation in the ECMO system and blood vessels. We aimed to identify the incidence and risk factors for cannula-associated arterial thrombosis (CaAT) post-decannulation.

METHODS

A retrospective study of patients undergoing arterial cannula removal following ECMO was performed. We evaluated the incidence of CaAT and compared the characteristics, ECMO machine parameters, cannula sizes, number of blood products transfused during ECMO, and daily hemostasis parameters in patients with and without CaAT. Multivariate analysis identified the risk factors for CaAT.

RESULTS

Forty-seven patients requiring venoarterial ECMO (VA-ECMO) or hybrid methods were recruited for thrombosis screening. The median Sequential Organ Failure Assessment score was 11 (interquartile range, 8-13). CaAT occurred in 29 patients (61.7%), with thrombosis in the superficial femoral artery accounting for 51.7% of cases. The rate of limb ischemia complications in the CaAT group was 17.2%. Multivariate analysis determined that the ECMO flow rate-body surface area (BSA) ratio (100 ml/min/m2) was an independent factor for CaAT, with an odds ratio of 0.79 (95% confidence interval, 0.66-0.95; P=0.014).

CONCLUSIONS

We found that the incidence of CaAT was 61.7% following successful decannulation from VA-ECMO or hybrid modes, and the ECMO flow rate-BSA ratio was an independent risk factor for CaAT. We suggest screening for arterial thrombosis following VA-ECMO, and further research is needed to determine the risks and benefits of such screening.

Numerical and experimental investigation of a lighthouse tip drainage cannula used in extracorporeal membrane oxygenation

BACKGROUND

Extracorporeal membrane oxygenation is a life-saving therapy used in case of acute respiratory/circulatory failure. Exposure of blood to non-physiological surfaces and high shear stresses is related to hemolytic damage and platelet activation. A detailed knowledge of the fluid dynamics of the components under different scenarios is thus paramount to assess the thrombogenicity of the circuit.

METHODS

An investigation of the flow structures developing in a conventional lighthouse tip (single-staged) drainage cannula was performed with cross-validated computational fluid dynamics and particle image velocimetry. The aim was to quantify the variation in drainage performance and stress levels induced by different fluid models, hematocrit and vessel-to-cannula flow rate ratios.

RESULTS

The results showed that the 90° bends of the flow through the side holes created a recirculation zone inside the cannula which increased residence time. Flow structures resembling a jet in a crossflow were also observed. The use of different hematocrits did not significantly affect drainage performances. The most proximal set of holes drained the largest fraction of fluid. However, different flow rate ratios altered the flow rate drained through the tip. The use of 2D data led to a 50% underestimation of shear rate levels. In the drainage zone the non-Newtonian behavior of blood was less relevant.

CONCLUSIONS

The most proximal holes drained the largest amount of fluid. The flow features and distribution of flow rates among the holes showed little dependence on the hematocrit. The non-Newtonian behavior of blood had a small influence on the dynamics of the flow.

Improved Flow Dynamics of Extracorporeal Membrane Oxygenation via Design Modification of Dual-Lumen Cannulas

Veno-venous extracorporeal membrane oxygenation (VV-ECMO) supports patients with severe respiratory failure not responding to conventional treatments. Single-site jugular venous cannulation with dual-lumen cannulas (DLC) have several advantages over traditional single-lumen cannulas, however, bleeding and thrombosis are common, limiting their clinical utility. This study numerically investigated the effects of DLC side holes on blood flow dynamics since the maximum wall shear stress (WSS) occurs around the side holes. A DLC based on the Avalon Elite 27Fr model was implanted into an idealized 3D model of the vena cava and right atrium (RA). Eight DLCs were developed by changing the number, diameter, and spacing of side holes through an iterative design process. Physiologic flow at the inferior vena cava (IVC) and superior vena cava (SVC) were applied along with a partial ECMO support of 2 L/min. The SST k-ω turbulent model was solved for 6.4 seconds. WSS, washout, stagnation volume, and recirculation were compared. For all DLCs, no stasis region lasted more than one cardiac cycle and a complete washout was obtained in less than 4 seconds. Due to the IVC and SVC backflows, maximum WSS occurred around the DLC side holes at late systole and late diastole. A DLC with 16 and three side holes within the IVC and SVC, respectively, reduced the maximum WSS by up to 67% over the Avalon Elite 27Fr. Improved DLCs provided a more uniform WSS distribution with lower WSS around the side holes, potentially reducing the chance of thrombosis and bleeding.

Ventricular Flow Dynamics with an Intra-Ventricular Balloon Pump: An In Vitro Analysis

Due to the high treatment costs associated with durable ventricular assist devices, an intra-ventricular balloon pump (IVBP) was developed to provide low-cost, short-term support for patients suffering from severe heart failure. It is imperative that intraventricular flow dynamics are evaluated with an IVBP to ensure stagnation points, and potential regions for thrombus formation, are avoided. This study used particle image velocimetry to evaluate flow patterns within the left ventricle of a simulated severe heart failure patient with IVBP support to assess left ventricle pulsatility as an indicator of the likelihood of flow stasis. Two inflation timings were evaluated against the baseline severe heart failure condition: IVBP co-pulsation and IVBP counter-pulsation with respect to ventricular systole. IVBP co-pulsation was found to have a reduced velocity range compared to the severe heart failure condition (0.44 m/s compared to 0.54 m/s). IVBP co-pulsation demonstrated an increase in peak velocities (0.25 m/s directed toward the aortic valve during systole, as opposed to 0.2 m/s in severe heart failure), indicating constructive energy in systole and cardiac output (1.7 L/min increase with respect to severe heart failure baseline - 3.5 L/min) throughout the cardiac cycle. IVBP counter-pulsation, while exhibiting the greatest peak systolic velocity directed to the aortic valve (0.4 m/s) was found to counterasct the natural vortex flow pattern during ventricular filling, as well as inducing a secondary ventricular pulse during diastole and a 23% increase in left ventricle end-diastolic volume (indicative of dilation). Ideal IVBP actuation timing did not result in reduced intraventricular pulsatility, indicating promising blood washout.

Improved Drainage Cannula Design to Reduce Thrombosis in Veno-Arterial Extracorporeal Membrane Oxygenation

Thrombosis is a potentially life-threatening complication in veno-arterial extracorporeal membrane oxygenation (ECMO) circuits, which may originate from the drainage cannula due to unfavorable blood flow dynamics. This study aims to numerically investigate the effect of cannula design parameters on local fluid dynamics, and thus thrombosis potential, within ECMO drainage cannulas. A control cannula based on the geometry of a 17 Fr Medtronic drainage cannula concentrically placed in an idealized, rigid-walled geometry of the right atrium and superior and inferior vena cava was numerically modeled. Simulated flow dynamics in the control cannula were systematically compared with 10 unique cannula designs which incorporated changes to side hole diameter, the spacing between side holes, and side hole angles. Local blood velocities, maximum wall shear stress (WSS), and blood residence time were used to predict the risk of thrombosis. Numerical results were experimentally validated using particle image velocimetry. The control cannula exhibited low blood velocities (59 mm/s) at the cannula tip, which may promote thrombosis. Through a reduction in the side hole diameter (2 mm), the spacing between the side holes (3 mm) and alteration in the side hole angle (30° relative to the flow direction), WSS was reduced by 52%, and cannula tip blood velocity was increased by 560% compared to the control cannula. This study suggests that simple geometrical changes can significantly alter the risk of thrombosis in ECMO drainage cannulas.

Extended finite element method for fluid-structure interaction in wave membrane blood pump

Numerical simulations of cardiac blood pump systems are integral to the optimization of device design, hydraulic performance and hemocompatibility. In wave membrane blood pumps, blood propulsion arises from the wave propagation along an oscillating immersed membrane, which generates small pockets of fluid that are pushed towards the outlet against an adverse pressure gradient. We studied the Fluid-Structure Interaction between the oscillating membrane and the blood flow via three-dimensional simulations using the Extended Finite Element Method (XFEM), an unfitted numerical technique that avoids remeshing by using a fluid fixed mesh. Our three-dimensional numerical simulations in a realistic pump geometry highlighted, for the first time in this field of application, that XFEM is a reliable strategy to handle complex industrial problems. Moreover, they showed the role of the membrane deformation in promoting a blood flow towards the outlet despite an adverse pressure gradient. We also simulated the pump system at different pressure conditions and we validated the numerical results against in-vitro experimental data.

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.

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.

Tipless transseptal cannula concept combines improved hemodynamic properties and risk-reduced placement: An in silico proof-of-concept

As a leading cause of death worldwide, heart failure is a serious medical condition in which many critically ill patients require temporary mechanical circulatory support (MCS) as a bridge-to-recovery or bridge-to-decision. In many cases, the TandemHeart system is used to unload the left heart by draining blood from the left atrium (LA) to the femoral artery via a transseptal multistage cannula. However, even though the correct positioning of the cannula is crucial for a safe treatment, the long cannula tip currently used in transseptal cannulas complicates positioning, making the cannula vulnerable to displacement during MCS. To overcome these limitations, we propose the development of a new tipless transseptal cannula with improved hemodynamic properties. We discuss the tipless cannula concept by comparing it to the common multistage cannula concept using computational fluid dynamics simulations and assess the flow field in the LA, the wall shear stresses (WSS), and the pressure loss. Across the two distinct time points of end-systole and end-diastole and two drainage flow rates of 3.5 and 5.0 L/min, we find a more homogeneous inlet flow pattern for the tipless cannula concept, accompanied by a remarkably reduced area of platelet-activating WSS (up to 10-times smaller area compared to the multistage cannula). Moreover, pressure loss is up to 14.5% lower in the tipless cannula concept, confirming overall improved hemodynamic properties of the tipless cannula concept. Finally, a diameter-dependent study reveals that lower WSS and pressure losses can be further reduced by large-lumen designs for any simulation setting. Overall, our results suggest that a tipless cannula concept remedies the crucial disadvantages of a long-tip multistage cannula by reducing the risk of misplacement, and it furthermore promotes optimized hemodynamics. With this successful proof-of-concept, we underscore the potential for and encourage the realization of further experimental investigations regarding the development of a tipless transseptal cannula for MCS.

Berlin Heart EXCOR Paediatric Ventricular Assist Device: Does Weight Matter?

BACKGROUND

Berlin Heart EXCOR (BH) ventricular assist devices provide mechanical long-term circulatory support in children with end-stage heart failure, as a bridge to transplantation or to recovery. Most studies are from large-volume paediatric cardiac centres.

AIM

The aim of this study was to analyse the experiences of three French centres and to compare these with available published data.

METHOD

We performed a retrospective observational study of three paediatric cardiac intensive care units. All children supported with BH devices were included. Morbidity and mortality data were collected and risk factors analysed.

RESULTS

Fifty-four (54) patients (54% male) were included. Survival rate was 73% while on a BH device. Median age at BH device implantation was 17 months (range 2-180 months). The predominant indication was dilated cardiomyopathy (61%). Bi-ventricular assist device was used in 25 (46%) cases. The total length of long-term circulatory support was 3,373 days, with a mean length per patient of 62.5 days (range 5-267 days). Thirty-two (32) patients were transplanted (59%) and seven (13%) were successfully weaned. Type and length of support did not influence morbidity. Main complications were renal dysfunction (57%), bleeding (41%), and infection (39%). In multivariate analysis, a weight <5 kg was significantly associated with higher mortality.

CONCLUSIONS

The weight seems to be the most important risk factor of mortality in this precarious condition.

Influence of Hemodialysis Catheter Insertion on Hemodynamics in the Central Veins

Central venous catheter (CVC) related thrombosis is a major cause of CVC dysfunction in patients under hemodialysis. The aim of our study was to investigate the impact of CVC insertion on hemodynamics in the central veins and to examine the changes in hemodynamic environments that may be related to thrombus formation due to the implantation of CVC. Patient-specific models of the central veins with and without CVC were reconstructed based on computed tomography images. Flow patterns in the veins were numerically simulated to obtain hemodynamic parameters such as time-averaged wall shear stress (TAWSS), oscillating shear index (OSI), relative residence time (RRT), and normalized transverse wall shear stress (transWSS) under pulsatile flow. The non-Newtonian effects of blood flow were also analyzed using the Casson model. The insertion of CVC caused significant changes in the hemodynamic environment in the central veins. A greater disturbance and increase of velocity were observed in the central veins after the insertion of CVC. As a result, TAWSS and transWSS were markedly increased, but most parts of OSI and RRT decreased. Newtonian assumption of blood flow would overestimate the increase in TAWSS after CVC insertion. High wall shear stress (WSS) and flow disturbance, especially the multidirectionality of the flow, induced by the CVC may be a key factor in initiating thrombosis after CVC insertion. Accordingly, approaches to decrease the flow disturbance during CVC insertion may help restrain the occurrence of thrombosis. More case studies with pre-operative and postoperative modeling and clinical follow-up need to be performed to verify these findings. Non-Newtonian blood flow assumption is recommended in computational fluid dynamics (CFD) simulations of veins with CVCs.

Effect of inflow cannula side-hole number on drainage flow characteristics: flow dynamic analysis using numerical simulation

Background:

Venous drainage in cardiopulmonary bypass is a very important factor for safe cardiac surgery. However, the ideal shape of venous drainage cannula has not been determined. In the present study, we evaluated the effect of side-hole number under fixed total area and venous drainage flow to elucidate the effect of increasing the side-hole numbers.

Method:

Computed simulation of venous drainage was performed. Cannulas were divided into six models: an end-hole model (EH) and models containing four (4SH), six (6SH), eight (8SH), 10 (10SH) or 12 side-holes (12SH). Total orifice area of the side-holes was fixed to 120 mm2 on each side-hole cannula. The end-hole orifice area was 36.3 mm2. The total area of the side-holes was kept constant when the number of side-holes was increased.

Result:

The mean venous drainage flow rate of the EH, 4SH, 6SH, 8SH, 10SH and 12SH was 2.57, 2.52, 2.51, 2.50, 2.49, 2.41 L/min, respectively. The mean flow rate decreased in accordance with the increased number of side-holes.

Conclusion:

We speculate that flow separation at the most proximal site of the side-hole induces stagnation of flow and induces energy loss. This flow separation may hamper the main stream from the end-hole inlet, which is most effective with low shear stress. The EH cannula was associated with the best flow rate and flow profile. However, by increasing side-hole numbers, flow separation occurs on each side-hole, resulting in more energy loss than the EH cannula and flow rate reduction.

Experimental measurement and numerical modelling of dye washout for investigation of blood residence time in ventricular assist devices

Ventricular assist devices have become the standard therapy for end-stage heart failure. However, their use is still associated with severe adverse events related to the damage done to the blood by fluid dynamic stresses. This damage relates to both the stress magnitude and the length of time the blood is exposed to that stress. We created a dye washout technique which combines experimental and numerical approaches to measure the washout times of ventricular assist devices. The technique was used to investigate washout characteristics of three commercially available and clinically used ventricular assist devices: the CentriMag, HVAD and HeartMate II. The time taken to reach 5% dye concentration at the outlet (T05) was used as an indicator of the total residence time. At a typical level of cardiac support, 5 L/min and 100 mmHg, T05 was 0.93, 0.28 and 0.16 s for CentriMag, HVAD and HeartMate II, respectively, and increased to 5.06, 1.64 and 0.96 s for reduced cardiac support of 1 L/min. Regional variations in washout characteristics are described in this article. While the volume of the flow domain plays a large role in the differences in T05 between the ventricular assist devices, after standardising for ventricular assist device volume, the secondary flow path was found to increase T05 by 35%. The results explain quantitatively, for the first time, why the CentriMag, which exerts low shear stress magnitude, has still been found to cause acquired von Willebrand Syndrome in patients.

Effects of Cone-Shaped Bend Inlet Cannulas of an Axial Blood Pump on Thrombus Formation: An Experiment and Simulation Study

Background

Cannula shape and connection style influence the risk of thrombus formation in the blood pump by varying the blood flow characteristics inside the pump. Inlet cannulas should be designed based on the need for anatomical fit and reducing the risk of thrombus generation in the blood pump. The effects on thrombus formation of the cone-shaped bend inlet cannulas of axial blood pumps should be studied.

Material/Methods

The cannulas were designed as cone-shaped, with 1 bent section connecting 2 straight sections. Both the silicone tube and novel cone-shaped cannula were simulated for comparison. The flow fields of a blood pump with inlet cannula were simulated by computational fluid dynamics (CFD) at flows of 2.0, 2.5, and 3.0 liters per minute (lpm), with pump rotational speeds of 7500, 8000, and 8500 rpm, respectively. Then, 6 two-dimensional (2D) particle image velocimetry (PIV) tests were conducted and the velocity distributions were analyzed.

Results

A low-velocity region was located inside the pump entrance when a soft silicone tube was used. At 8500 rpm and 3.0 lpm working condition, the minimum velocity inside the pump with cone-shaped cannulas was 2.5×10⁻¹ m/s. The cone-shaped cannulas eliminated the low-velocity region inside the pump. Both CFD and PIV results showed that the low-velocity region did not spread to the entrance of the blood pump within the flow range from 2.0 lpm to 7.0 lpm.

Conclusions

The designed cone-shaped bent cannulas can eliminate the low-velocity region inside the blood pump and reduce the risk of thrombus formation in the blood pump.

Device thrombosis and pre-clinical blood flow models for assessing antithrombogenic efficacy of drug-device combinations

Thrombosis associated with blood-contacting devices is a complex process involving several component interactions that have eluded precise definition. Extensive investigations of individual biological modules such as protein adsorption, coagulation cascade activation and platelet activation/adhesion/aggregation have provided an initial foundation for developing biomaterials for blood-contacting devices, but a material that is intrinsically non-thrombogenic is yet to be developed. The well-recognized association between fluid dynamics parameters such as shear stress, vortices, stagnation and thrombotic processes such as platelet aggregation and coagulation aggravate thrombosis on most device geometries that elicit these flow disturbances. Thus, antithrombotic drugs that were developed to treat thrombosis associated with vascular diseases such as atherosclerosis have also been adapted to mitigate the risk of device thrombosis. However, balancing the risk of bleeding with the antithrombotic efficacy of these drugs continues to be a challenge, and surface modification of devices with these drug molecules to mitigate device thrombosis locally has been explored. Pre-clinical blood flow models to test the effectiveness of these drug-device combinations have also evolved and several in-vitro, ex-vivo, and in-vivo test configurations are available with their attendant merits and limitations. Despite considerable efforts toward iterative design and testing of blood contacting devices and antithrombogenic surface modifications, device thrombosis remains an unsolved problem.

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.

Percutaneous Double Lumen Cannula for Right Ventricle Assist Device System: A Computational Fluid Dynamics Study

OBJECTIVES

Our goal is to develop a double lumen cannula (DLC) for a percutaneous right ventricular assist device (pRVAD) in order to eliminate two open chest surgeries for RVAD installation and removal. The objective of this study was to evaluate the performance, flow pattern, blood hemolysis, and thrombosis potential of the pRVAD DLC.

METHODS

Computational fluid dynamics (CFD), using the finite volume method, was performed on the pRVAD DLC. For Reynolds numbers <4000, the laminar model was used to describe the blood flow behavior, while shear-stress transport k-ω model was used for Reynolds numbers >4000. Bench testing with a 27 Fr prototype was performed to validate the CFD calculations.

RESULTS

There was <1.3% difference between the CFD and experimental pressure drop results. The Lagrangian approach revealed a low index of hemolysis (0.012% in drainage lumen and 0.0073% in infusion lumen) at 5 l/min flow rate. Blood stagnancy and recirculation regions were found in the CFD analysis, indicating a potential risk for thrombosis.

CONCLUSIONS

The pRVAD DLC can handle up to 5 l/min flow with limited potential hemolysis. Further modification of the pRVAD DLC is needed to address blood stagnancy and recirculation.

In Vivo Evaluation of Active and Passive Physiological Control Systems for Rotary Left and Right Ventricular Assist Devices

Preventing ventricular suction and venous congestion through balancing flow rates and circulatory volumes with dual rotary ventricular assist devices (VADs) configured for biventricular support is clinically challenging due to their low preload and high afterload sensitivities relative to the natural heart. This study presents the in vivo evaluation of several physiological control systems, which aim to prevent ventricular suction and venous congestion. The control systems included a sensor-based, master/slave (MS) controller that altered left and right VAD speed based on pressure and flow; a sensor-less compliant inflow cannula (IC), which altered inlet resistance and, therefore, pump flow based on preload; a sensor-less compliant outflow cannula (OC) on the right VAD, which altered outlet resistance and thus pump flow based on afterload; and a combined controller, which incorporated the MS controller, compliant IC, and compliant OC. Each control system was evaluated in vivo under step increases in systemic (SVR ∼1400-2400 dyne/s/cm(5) ) and pulmonary (PVR ∼200-1000 dyne/s/cm(5) ) vascular resistances in four sheep supported by dual rotary VADs in a biventricular assist configuration. Constant speed support was also evaluated for comparison and resulted in suction events during all resistance increases and pulmonary congestion during SVR increases. The MS controller reduced suction events and prevented congestion through an initial sharp reduction in pump flow followed by a gradual return to baseline (5.0 L/min). The compliant IC prevented suction events; however, reduced pump flows and pulmonary congestion were noted during the SVR increase. The compliant OC maintained pump flow close to baseline (5.0 L/min) and prevented suction and congestion during PVR increases. The combined controller responded similarly to the MS controller to prevent suction and congestion events in all cases while providing a backup system in the event of single controller failure.

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.

Outside-In Hemofiltration for Prolonged Operation without Clogging

Hemofiltration (HF) is used extensively for continuous renal replacement therapy, but long-term treatment is limited by thrombosis leading to fiber clogging. Maximum filter life is typically less than 20 hours. We have achieved for the first time continuous and consistent hemofiltration for more than 100 hours using outside-in hemofiltration with the blood flow into the inter-fiber space (IFS). Although thrombi do deposit in the IFS, they have minimal affect on the blood flow and filtrate flux due to the three-dimensional system of interconnected hydrodynamic flow channels in the IFS. Microscopic examination of sections of the fiber bundle showed that deposited thrombi have dimensions about the size of the gaps between the hollow fibers and remain isolated from each other. A simple mathematical model is developed to describe the effect of thrombus deposition on the fluid flow that accounts for the enhanced performance arising from the interconnected flow. The hydrodynamic advantage of outside-in HF decreases at low anticoagulant concentration due to the instability in the blood and the very high volume fraction of thrombi that deposit in the entrance zone of the filter. These results clearly demonstrate the significant potential advantages of using outside-in hemofiltration for long-term renal replacement therapy.

A compliant, banded outflow cannula for decreased afterload sensitivity of rotary right ventricular assist devices

Biventricular support with dual rotary ventricular assist devices (VADs) has been implemented clinically with restriction of the right VAD (RVAD) outflow cannula to artificially increase afterload and, therefore, operate within recommended design speed ranges. However, the low preload and high afterload sensitivity of these devices increase the susceptibility of suction events. Active control systems are prone to sensor drift or inaccurate inferred (sensor-less) data, therefore an alternative solution may be of benefit. This study presents the in vitro evaluation of a compliant outflow cannula designed to passively decrease the afterload sensitivity of rotary RVADs and minimize left-sided suction events. A one-way fluid-structure interaction model was initially used to produce a design with suitable flow dynamics and radial deformation. The resultant geometry was cast with different initial cross-sectional restrictions and concentrations of a softening diluent before evaluation in a mock circulation loop. Pulmonary vascular resistance (PVR) was increased from 50 dyne s/cm(5) until left-sided suction events occurred with each compliant cannula and a rigid, 4.5 mm diameter outflow cannula for comparison. Early suction events (PVR ∼ 300 dyne s/cm(5) ) were observed with the rigid outflow cannula. Addition of the compliant section with an initial 3 mm diameter restriction and 10% diluent expanded the outflow restriction as PVR increased, thus increasing RVAD flow rate and preventing left-sided suction events at PVR levels beyond 1000 dyne s/cm(5) . Therefore, the compliant, restricted outflow cannula provided a passive control system to assist in the prevention of suction events with rotary biventricular support while maintaining pump speeds within normal ranges of operation.

Numerical investigation of three patterns of motion in an electromagnetic pulsatile VAD

Hemolysis and thrombus formation which are critical concerns in designing a long-term implantable ventricular assist device (VAD) have impeded the widespread use of VADs. In this study, thus, the three-dimensional fluid domain of blood flow in a small bichamber positive displacement VAD (25 ml) with a magnetically levitated moving pusher plate was simulated by the means of a finite element package called ADINA. To optimize the function of the pump for minimizing shear stress induced blood damage, three different driver patterns (linear, sinusoidal, and Guyton's pulse) were investigated. The first pattern produced a constant flow, whereas the two others created pulsatile flows. The flow pattern and the distribution of shear stress of each pattern were observed for comparison. It was revealed that the three types of motions may induce less than 0.06% red blood cell damage. Moreover, in comparison to the other patterns not only did the sinusoidal motion of the pusher plate cause less risk of hemolysis, but in comparison to the linear pattern, it produced a pulsatile flow which reduced the stagnation areas in chambers, lowering the probability of thrombosis. In addition, this motion eliminates the probability of cavitations as compared with the Guyton's pulse pattern.

Local hemodynamic disturbance accelerates early thrombosis of small-caliber expanded polytetrafluoroethylene grafts

OBJECTIVE

Small-caliber expanded polytetrafluoroethylene (ePTFE) grafts cannot be used widely in the clinical situation because of early thrombosis and occlusion. This unsolved and under-recognized problem warrants further investigation.

METHODS

Grafts of uncoated ePTFE (n = 6) and anti-CD34-coated ePTFE (n = 6) were implanted unilaterally into the carotid artery in 12 domestic pigs. Ultrasonography was used to test the proximal and distal anastomotic stoma morphology, diameters and blood velocities. A thrombosis instrument was used to examine the blood coagulation state. After seven days, the pigs were sacrificed and the implanted grafts were excised for general and histological analysis. Computational fluid dynamics was used to investigate the blood flow fields of the implanted grafts and to calculate parameters that might be indicative of thrombosis.

RESULTS

Thrombosis was detected in 10 of 12 (83.3%) implanted ePTFE grafts, 5 in uncoated grafts and 5 in anti-CD34-coated grafts. Endothelial cell coverage was observed in both uncoated and anti-CD34-coated grafts. No obvious abnormalities in anastomotic stoma or blood coagulation state were observed. Computer-based local hemodynamic simulation showed the low flexibility of synthetic ePTFE grafts caused obvious coarctation. Local wall pressure, velocity and wall shear stress were much higher than in the contralateral normal vessel.

CONCLUSIONS

The patency of small-caliber ePTFE grafts for clinical use is impaired by early thrombosis due to mixed causes. Local hemodynamic disturbance was the most powerful predictor of early thrombosis. Decreasing local hemodynamic disturbance, improving the quality of anastomotic stoma, selecting reasonable anticoagulation strategies and promoting rapid endothelialization may increase the long-term patency of small-caliber vascular grafts.

Numerical investigation of the effect of cannula placement on thrombosis

Despite the rapid advancement of left ventricular assist devices (LVADs), adverse events leading to deaths have been frequently reported in patients implanted with LVADs, including bleeding, infection, thromboembolism, neurological dysfunction and hemolysis. Cannulation forms an important component with regards to thrombus formation in assisted patients by varying the intraventricular flow distribution in the left ventricle (LV). To investigate the correlation between LVAD cannula placement and potential for thrombus formation, detailed analysis of the intraventricular flow field was carried out in the present study using a two way fluid structure interaction (FSI), axisymmetric model of a passive LV incorporating an inflow cannula. Three different cannula placements were simulated, with device insertion near the LV apex, penetrating one-fourth and mid-way into the LV long axis. The risk of thrombus formation is assessed by analyzing the intraventricular vorticity distribution and its associated vortex intensity, amount of stagnation flow in the ventricle as well as the level of wall shear stress. Our results show that the one-fourth placement of the cannula into the LV achieves the best performance in reducing the risk of thrombus formation. Compared to cannula placement near the apex, higher vortex intensity is achieved at the one-fourth placement, thus increasing wash out of platelets at the ventricular wall. One-fourth LV penetration produced negligible stagnation flow region near the apical wall region, helping to reduce platelet deposition on the surface of the cannula and the ventricular wall.

A quantitative comparison of mechanical blood damage parameters in rotary ventricular assist devices: shear stress, exposure time and hemolysis index

Ventricular assist devices (VADs) have already helped many patients with heart failure but have the potential to assist more patients if current problems with blood damage (hemolysis, platelet activation, thrombosis and emboli, and destruction of the von Willebrand factor (vWf)) can be eliminated. A step towards this goal is better understanding of the relationships between shear stress, exposure time, and blood damage and, from there, the development of numerical models for the different types of blood damage to enable the design of improved VADs. In this study, computational fluid dynamics (CFD) was used to calculate the hemodynamics in three clinical VADs and two investigational VADs and the shear stress, residence time, and hemolysis were investigated. A new scalar transport model for hemolysis was developed. The results were compared with in vitro measurements of the pressure head in each VAD and the hemolysis index in two VADs. A comparative analysis of the blood damage related fluid dynamic parameters and hemolysis index was performed among the VADs. Compared to the centrifugal VADs, the axial VADs had: higher mean scalar shear stress (sss); a wider range of sss, with larger maxima and larger percentage volumes at both low and high sss; and longer residence times at very high sss. The hemolysis predictions were in agreement with the experiments and showed that the axial VADs had a higher hemolysis index. The increased hemolysis in axial VADs compared to centrifugal VADs is a direct result of their higher shear stresses and longer residence times. Since platelet activation and destruction of the vWf also require high shear stresses, the flow conditions inside axial VADs are likely to result in more of these types of blood damage compared with centrifugal VADs.

Numerical Simulation of LVAD Inflow Cannulas with Different Tip

The tip structure of LVAD inflow cannula is one of major factors to lead adverse events such as thrombosis and suction leading to obstruction. In this research, four kinds of tips that had been used in inflow cannulas were selected and designed. The flow field of the four inflow cannulas inserted into the apex of left ventricle (LV) was numerically computed by computational fluid dynamics. The flow behavior was analyzed in order to compare the blood compatibility and suction in left ventricle and cannulas after the inflow cannulas with different tips were inserted to the apex of LV. The results showed that the cannula tip structure affected the LVAD performance. Among these four cannulas, the trumpet-tipped inflow cannula owned the best performance in smooth flow velocity distribution without backflow or low-velocity flow so that it was the best in blood compatibility. Nevertheless, the caged tipped cannula was the worst in blood compatibility. And the blunt-tipped and beveled tipped inflow cannulas may obstruct more easily than trumpet and caged tipped inflow cannulas because of their shape. The study indicated that the trumpet tip was the most preferable for the inflow cannula of long-term LVAD.

Fluid structure interaction simulation of left ventricular flow dynamics under left ventricular assist device support

For patient's receiving mechanical circulatory support, malfunction of the left ventricular assist device (LVADs) as well as mal-positioning of the cannula imposes serious threats to their life. It is therefore important to characterize the flow pattern and pressure distribution within the ventricle in the presence of an LVAD. In this paper, we present a 2D axisymmetric fluid structure interaction model of the passive left ventricle (LV) incorporating an LVAD cannula to simulate the effect of the LVAD cannula placement on the vortex dynamics. Results showed that larger recirculation area was formed at the cannula tip with increasing cannula insertion depth, and this is believed to reduce the risk of thrombus formation. Furthermore, we also simulated suction events (collapse of the LV) by closing the inlet. Vortex patterns were significantly altered under this condition, and the greatest LV wall displacement was observed at the part of the myocardium closest to the cannula tip.

Numerical characterization of hemodynamics conditions near aortic valve after implantation of Left Ventricular Assist Device

Left Ventricular Assist Devices (LVADs) are implantable mechanical pumps that temporarily aid the function of the left ventricle. The use of LVADs has been associated with thrombus formation next to the aortic valve and close to the anastomosis region, especially in patients in which the native cardiac function is negligible and the aortic valve remains closed. Stagnation points and recirculation zones have been implicated as the main fluid dynamics factors contributing to thrombus formation. The purpose of the present study was to develop and use computer simulations based on a fluid-structure interaction (FSI) solver to study flow conditions corresponding to different strategies in LVAD ascending aortic anastomosis providing a scenario with the lowest likelihood of thrombus formation. A novel FSI algorithm was developed to deal with the presence of multiple structures corresponding to different elastic properties of the native aorta and of the LVAD cannula. A sensitivity analysis of different variables was performed to assess their impact of flow conditions potentially leading to thrombus formation. It was found that the location of the anastomosis closest to the aortic valve (within 4 cm away from the valve) and at the angle of 30 minimizes the likelihood of thrombus formation. Furthermore, it was shown that the rigidity of the dacron anastomosis cannula plays almost no role in generating pathological conditions downstream from the anastomosis. Additionally, the flow analysis presented in this manuscript indicates that compliance of the cardiovascular tissue acts as a natural inhibitor of pathological flow conditions conducive to thrombus formation and should not be neglected in computer simulations.

The use of computational fluid dynamics in the development of ventricular assist devices

Progress in the field of prosthetic cardiovascular devices has significantly contributed to the rapid advancements in cardiac therapy during the last four decades. The concept of mechanical circulatory assistance was established with the first successful clinical use of heart-lung machines for cardiopulmonary bypass. Since then a variety of devices have been developed to replace or assist diseased components of the cardiovascular system. Ventricular assist devices (VADs) are basically mechanical pumps designed to augment or replace the function of one or more chambers of the failing heart. Computational Fluid Dynamics (CFD) is an attractive tool in the development process of VADs, allowing numerous different designs to be characterized for their functional performance virtually, for a wide range of operating conditions, without the physical device being fabricated. However, VADs operate in a flow regime which is traditionally difficult to simulate; the transitional region at the boundary of laminar and turbulent flow. Hence different methods have been used and the best approach is debatable. In addition to these fundamental fluid dynamic issues, blood consists of biological cells. Device-induced biological complications are a serious consequence of VAD use. The complications include blood damage (haemolysis, blood cell activation), thrombosis and emboli. Patients are required to take anticoagulation medication constantly which may cause bleeding. Despite many efforts blood damage models have still not been implemented satisfactorily into numerical analysis of VADs, which severely undermines the full potential of CFD. This paper reviews the current state of the art CFD for analysis of blood pumps, including a practical critical review of the studies to date, which should help device designers choose the most appropriate methods; a summary of blood damage models and the difficulties in implementing them into CFD; and current gaps in knowledge and areas for future work.