# Numerical Simulation of LVAD Inflow Cannulas with Different Tip

Authors: Guang-Mao Liu; Hai-Bo Chen; Fu-Liang Luo; Yan Zhang; Han-Song Sun; Jian-Ye Zhou; Sheng-Shou Hu
Journal: International Journal of Chemical Engineering
 (2012)

---

## Abstract

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.


---

## Body

## 1. Introduction

LVADs have been widely applied to assist the circulation of heart failure (HF) patients. Some LVADs can prolong HF patients’ lives for more than two years [1–5]. Now, most LVADs are equipped with complete accessories such as designed inflow and outflow cannulas to ensure that the LVADs can work safely and longer. The inflow cannula is the connection between left ventricle and LVAD. The inflow cannulas are mainly made of good blood compatibility materials such as titanium alloy, silicon rubber, or polythene. The shapes of inflow cannulas include blunt, beveled, caged, trumpet, and beak [6]. For example, the inflow cannula of Berlin Heart INCOR is made of silicon rubber with beak-shaped tip [7, 8]. And the inflow cannula of Debakey LVAD is made of titanium alloy with trumpet shaped tip [9, 10]. The tip shape is a critical factor to fluid dynamic performance. The flow field and stress distribution in left ventricle and cannulas depend on the inflow cannulas tip structure at a certain flow rate and operative procedures such as insertion depth and insertion angle. Different tip structures of inflow cannulas can cause various degree of thrombosis and suction leading to obstruction that may further influence the LVAD’s service life and the survival of the patients [11, 12]. Thrombosis is mainly caused by static and/or vortical flow in the inflow cannulas.Mechanical hemolysis is another important aspect of the inflow cannulas. The hemolysis in the inflow cannulas is mainly caused by the high shear stress of flow in the inflow cannulas. The effective numerical prediction of mechanical hemolysis has been built [13–15]. In vitro and in vivo mechanical hemolysis tests have been studied previously [16–18].To obtain a good performance in blood compatibility, the tip shape is the crucial factor to reduce static and vortical flow area in an inflow cannula design. The objective of this paper is to find which kind of tip shape is better for LVAD inflow cannula. The fluid dynamics of the inflow cannula are studied with numerical simulation or test in a mock loop model. Numerical simulation can simulate and analyze the flow velocity, stress distribution, and wall shear stress of cannula insertion [19, 20]. In this paper, we built the numerical model of blunt, beveled, caged, and trumpet-tipped cannulas inserted into left ventricle, respectively. The velocity distribution and stress distribution around the cannula tips and in cannulas were calculated. The blood compatibility and risk of suction leading to obstruction [21] of the four insertions were compared. The analysis results can help us select tip shape and design an inflow cannula which is suitable for long-term LVAD.


## 2. Materials and Methods

### 2.1. Geometry Model

Figure1 illustrated the 3D model of four kinds of inflow cannulas with different tip shapes built by Pro-E software. The cannulas were made of titanium alloy. The four cannulas were assumed to be in the same size where the inner diameter was d=12 mm, and the wall thickness was s=1 mm. And they were considered as straight pipes to eliminate effects of curved portion. The cannulas’ length out of ventricle is l=46 mm.3D model of inflow cannula tip: (a) blunt, (b) beveled, (c) trumpet, and (d) caged.
(a)(b)(c)(d)Figure2 showed the geometry model of LV for calculation. According to Fraser et al. and Tsukiya et al. [19, 22], the geometry of the left ventricle was simplified without regarding the heart real geometry, insertion depth, and insertion angle. The insertion point was assumed the apex of left ventricle. The inlet surface diameter of LV was set to 80 mm. The LV wall thickness was 10 mm and uniform. The depth from the inlet surface to the apex was 70 mm.Figure 2
LV geometry model used to compare the inflow cannula tip.


### 2.2. Computational Model

Four computational models of cannula inserted into LV were built by computational fluid dynamics software Workbench (ANSYS, Inc, Canonsburg, Pennsylvania, United States) according to the geometry models in Figures1 and 2. Figure 3 demonstrated the typical computational grid for the four inflow cannulas. We adopted unstructured numeric grid. Furthermore, blood was assumed to be Newtonian with density of ρ=1059 kg/m3 and viscosity of μ=0.0035 kg/m·s. The average flow rate was set to Q=5 L/min. Thus the inlet mass flow was 0.0088 kg/s. The static outlet pressure was set to 0 Mpa. The LV wall and cannula wall were considered to be stiffness boundary and immobile. The fluid-structure interaction simulation was complicated, so the interaction was not involved in this study [19, 20, 22].Numerical models of inflow cannulas insertion: (a) blunt, (b) beveled, (c) trumpet, and (d) caged.
(a)(b)(c)(d)The Reynolds number was used as the index to estimate the flow to be laminar or turbulent. When the index is greater than the critical Reynolds number, the flow is regarded as turbulent. The Reynolds number definition is(1)Re=VLν,
where V is the mean velocity of the section (m/s), L is the characteristic length (m), and ν is fluid dynamic viscosity (m2/s). For circular tube flow, the circular tube diameter d is the characteristic length L. The Reynolds number of inflow cannula blood flow is(2)Re=VLν=Q/(60*1000⋅π(d/2000)2)⋅d/1000μ/ρ=Qρ15πμd=2658.It is generally acknowledged that the critical Reynolds number of tube flow is 2320 [23]. So the flow in inflow cannula was turbulent flow. Therefore, K-ε standard turbulence model was employed in numerical computation. The computational grid was generated by Workbench Fluent software. The unstructured numeric grid was generated from cannula surface.


## 2.1. Geometry Model

Figure1 illustrated the 3D model of four kinds of inflow cannulas with different tip shapes built by Pro-E software. The cannulas were made of titanium alloy. The four cannulas were assumed to be in the same size where the inner diameter was d=12 mm, and the wall thickness was s=1 mm. And they were considered as straight pipes to eliminate effects of curved portion. The cannulas’ length out of ventricle is l=46 mm.3D model of inflow cannula tip: (a) blunt, (b) beveled, (c) trumpet, and (d) caged.
(a)(b)(c)(d)Figure2 showed the geometry model of LV for calculation. According to Fraser et al. and Tsukiya et al. [19, 22], the geometry of the left ventricle was simplified without regarding the heart real geometry, insertion depth, and insertion angle. The insertion point was assumed the apex of left ventricle. The inlet surface diameter of LV was set to 80 mm. The LV wall thickness was 10 mm and uniform. The depth from the inlet surface to the apex was 70 mm.Figure 2
LV geometry model used to compare the inflow cannula tip.


## 2.2. Computational Model

Four computational models of cannula inserted into LV were built by computational fluid dynamics software Workbench (ANSYS, Inc, Canonsburg, Pennsylvania, United States) according to the geometry models in Figures1 and 2. Figure 3 demonstrated the typical computational grid for the four inflow cannulas. We adopted unstructured numeric grid. Furthermore, blood was assumed to be Newtonian with density of ρ=1059 kg/m3 and viscosity of μ=0.0035 kg/m·s. The average flow rate was set to Q=5 L/min. Thus the inlet mass flow was 0.0088 kg/s. The static outlet pressure was set to 0 Mpa. The LV wall and cannula wall were considered to be stiffness boundary and immobile. The fluid-structure interaction simulation was complicated, so the interaction was not involved in this study [19, 20, 22].Numerical models of inflow cannulas insertion: (a) blunt, (b) beveled, (c) trumpet, and (d) caged.
(a)(b)(c)(d)The Reynolds number was used as the index to estimate the flow to be laminar or turbulent. When the index is greater than the critical Reynolds number, the flow is regarded as turbulent. The Reynolds number definition is(1)Re=VLν,
where V is the mean velocity of the section (m/s), L is the characteristic length (m), and ν is fluid dynamic viscosity (m2/s). For circular tube flow, the circular tube diameter d is the characteristic length L. The Reynolds number of inflow cannula blood flow is(2)Re=VLν=Q/(60*1000⋅π(d/2000)2)⋅d/1000μ/ρ=Qρ15πμd=2658.It is generally acknowledged that the critical Reynolds number of tube flow is 2320 [23]. So the flow in inflow cannula was turbulent flow. Therefore, K-ε standard turbulence model was employed in numerical computation. The computational grid was generated by Workbench Fluent software. The unstructured numeric grid was generated from cannula surface.


## 3. Results

### 3.1. Velocity Distribution

Figure4 showed velocity contour from inflow cannulas tip to outlet of four insertions. The flow velocity altered around the cannula tips and in cannulas. The flow velocity in ventricle was well distributed in the model of blunt and trumpet-tipped cannula. The maximal flow velocity was lower than 1 m/s. The velocity distribution was uniform and symmetrical along the z-axis after the blood flowed into the cannulas (Figures 4(a) and 4(c)). The uniform distribution velocity in cannulas could improve the flow field in LVADs because the blood flows into LVAD from inflow cannulas directly. That meant the blood compatibility of LVADs was advanced. The beveled and caged tips disturbed the flow in left ventricle. The maximal flow velocity reached to 1.2 m/s. The velocity distribution was not even or symmetrical along the z-axis in the cannulas flow (Figures 4(b) and 4(d)). The uneven velocity distribution in these two cannulas adversely effected flow in LVAD.Velocity contour from tip to outlet: (a) blunt, (b) beveled, (c) trumpet, and (d) caged.
(a)(b)(c)(d)Figure5 showed velocity vectors around inflow cannulas tips. There was obvious vortex region between the blunt-tipped cannula outer wall and ventricular inner wall near the insertion point (Figure 5 (a)). The flow velocity in the region was slower than 0.1 m/s. So this area could be potential of thrombosis. Nevertheless, the area in blunt-tipped cannula was not likely to form thrombus because the flow in it was smooth. There were six obvious vortex regions near the caged tip (Figure 5(b)). Four of the six vortex regions were in caged tipped cannula and the minimum velocity was lower than 0.07 m/s. Additionally, turbulent flow was obvious near the inlet holes areas. Therefore these regions were likely to form clot. Figure 5(c) showed the velocity vector of beveled tipped inflow cannula. There were two vortex regions near insertion. One region was in left ventricle, and the other was in beveled tipped cannula. Both of these two regions were likely to form clot because the velocity was lower than 0.08 m/s. There was not vortex or flow stasis field in the trumpet tipped cannula (Figure 5(d)). The even velocity was 0.7 m/s and well distributed. Thus trumpet-tipped cannula had good blood compatibility and was not likely to form thrombosis.Velocity vector near tip: (a) blunt, (b) caged, (c) beveled, and (d) trumpet.
(a)(b)(c)(d)


### 3.2. Static Pressure Distribution

Ventricular collapse and suction may occur when negative or low-pressure fields appeared. Furthermore, obstruction may happen if the tip could not resist the suction. Figure6 showed static pressure distribution of four inflow cannulas insertion. The blunt tipped cannula did not create low-pressure field near the ventricular inner wall (Figure 6(a)). The beveled tipped cannula created one low-pressure field near the inflow bevel which was likely to ventricular collapsing (Figure 6(b)). The trumpet-tipped cannula had uniform low-pressure distribution around inlet (Figure 6(c)) which may cause ventricular collapse. Nevertheless, obstruction may not happen because the trumpet structure can resist suction. The caged tipped cannula created serious low-pressure field near inlet holes (Figure 6(d)). Ventricular collapse and suction were likely to occur. But the inlet could not be obstructed completely because not all the nine inlet holes could be obstructed at the same time.Contour of static pressure: (a) blunt, (b) beveled, (c) trumpet, and (d) caged.
(a)(b)(c)(d)


### 3.3. Wall Shear Stress Distribution

Shear stress distribution was analyzed to evaluate hemolytic potential. Hemolysis is related to wall shear stress (WSS) value and blood exposure time [13–15, 24–26]. The computational fluid dynamics model with Eulerian formulation to estimate hemolysis in prosthesis device had been studied by Lacasse et al. [13], Song et al. [14], and Arvand et al. [15].We estimated the damage to red cells with (3) [14]:(3)D=∫outletinlet1.8×10-6⋅τ1.991⋅dT0.765=∑outletinlet1.8×10-6⋅τ1.991⋅ΔT0.765,
where D symbolizes the blood damage index, τ is shear stress (Pa), and T is blood exposure time (s). The particle displacement is calculated using forward Euler integration of the particle velocity over time step (δt) as (4):(4)xin=xio+dxiodtδt,
where x denotes the displacement and the superscripts o and n refer to old and new values, respectively.The blood damage index of the four inflow cannulas was shown in Figure7. The estimation value of blood damage was between 5 × 10−5% to 20 × 10−5%. The acceptable estimation value of blood damage for LVAD was 10−3% [27, 28]. So the blood damage in inflow cannula was negligible. The shear stress in inflow cannula was not able to cause hemolysis.Figure 7
Blood damage index of inflow cannula.


### 3.4. Limitations

This study eliminated the factors of real ventricle geometry, insertion depth, and insertion angle. But actually, these factors may change the flow distribution in LV and inflow cannula and eventually cause the suction, and obstruction. So the effect of insertion depth and angle on thrombosis potential, suction and obstruction should be further studied in the future.The other limitation in this study is the simplified boundary conditions with rigid cavity and uniform flow inlet assumption. In reality, the flow in LV is much more complicated. The flow and pressure in ventricle change all the time with the beating. As the purpose of this study is to compare the tip shapes of the inflow cannulas, we think that the simplified boundary conditions are acceptable. Of course, the final outcome should be verified in animal test. And based on the results, profound study with more exact conditions is need in the future.


## 3.1. Velocity Distribution

Figure4 showed velocity contour from inflow cannulas tip to outlet of four insertions. The flow velocity altered around the cannula tips and in cannulas. The flow velocity in ventricle was well distributed in the model of blunt and trumpet-tipped cannula. The maximal flow velocity was lower than 1 m/s. The velocity distribution was uniform and symmetrical along the z-axis after the blood flowed into the cannulas (Figures 4(a) and 4(c)). The uniform distribution velocity in cannulas could improve the flow field in LVADs because the blood flows into LVAD from inflow cannulas directly. That meant the blood compatibility of LVADs was advanced. The beveled and caged tips disturbed the flow in left ventricle. The maximal flow velocity reached to 1.2 m/s. The velocity distribution was not even or symmetrical along the z-axis in the cannulas flow (Figures 4(b) and 4(d)). The uneven velocity distribution in these two cannulas adversely effected flow in LVAD.Velocity contour from tip to outlet: (a) blunt, (b) beveled, (c) trumpet, and (d) caged.
(a)(b)(c)(d)Figure5 showed velocity vectors around inflow cannulas tips. There was obvious vortex region between the blunt-tipped cannula outer wall and ventricular inner wall near the insertion point (Figure 5 (a)). The flow velocity in the region was slower than 0.1 m/s. So this area could be potential of thrombosis. Nevertheless, the area in blunt-tipped cannula was not likely to form thrombus because the flow in it was smooth. There were six obvious vortex regions near the caged tip (Figure 5(b)). Four of the six vortex regions were in caged tipped cannula and the minimum velocity was lower than 0.07 m/s. Additionally, turbulent flow was obvious near the inlet holes areas. Therefore these regions were likely to form clot. Figure 5(c) showed the velocity vector of beveled tipped inflow cannula. There were two vortex regions near insertion. One region was in left ventricle, and the other was in beveled tipped cannula. Both of these two regions were likely to form clot because the velocity was lower than 0.08 m/s. There was not vortex or flow stasis field in the trumpet tipped cannula (Figure 5(d)). The even velocity was 0.7 m/s and well distributed. Thus trumpet-tipped cannula had good blood compatibility and was not likely to form thrombosis.Velocity vector near tip: (a) blunt, (b) caged, (c) beveled, and (d) trumpet.
(a)(b)(c)(d)


## 3.2. Static Pressure Distribution

Ventricular collapse and suction may occur when negative or low-pressure fields appeared. Furthermore, obstruction may happen if the tip could not resist the suction. Figure6 showed static pressure distribution of four inflow cannulas insertion. The blunt tipped cannula did not create low-pressure field near the ventricular inner wall (Figure 6(a)). The beveled tipped cannula created one low-pressure field near the inflow bevel which was likely to ventricular collapsing (Figure 6(b)). The trumpet-tipped cannula had uniform low-pressure distribution around inlet (Figure 6(c)) which may cause ventricular collapse. Nevertheless, obstruction may not happen because the trumpet structure can resist suction. The caged tipped cannula created serious low-pressure field near inlet holes (Figure 6(d)). Ventricular collapse and suction were likely to occur. But the inlet could not be obstructed completely because not all the nine inlet holes could be obstructed at the same time.Contour of static pressure: (a) blunt, (b) beveled, (c) trumpet, and (d) caged.
(a)(b)(c)(d)


## 3.3. Wall Shear Stress Distribution

Shear stress distribution was analyzed to evaluate hemolytic potential. Hemolysis is related to wall shear stress (WSS) value and blood exposure time [13–15, 24–26]. The computational fluid dynamics model with Eulerian formulation to estimate hemolysis in prosthesis device had been studied by Lacasse et al. [13], Song et al. [14], and Arvand et al. [15].We estimated the damage to red cells with (3) [14]:(3)D=∫outletinlet1.8×10-6⋅τ1.991⋅dT0.765=∑outletinlet1.8×10-6⋅τ1.991⋅ΔT0.765,
where D symbolizes the blood damage index, τ is shear stress (Pa), and T is blood exposure time (s). The particle displacement is calculated using forward Euler integration of the particle velocity over time step (δt) as (4):(4)xin=xio+dxiodtδt,
where x denotes the displacement and the superscripts o and n refer to old and new values, respectively.The blood damage index of the four inflow cannulas was shown in Figure7. The estimation value of blood damage was between 5 × 10−5% to 20 × 10−5%. The acceptable estimation value of blood damage for LVAD was 10−3% [27, 28]. So the blood damage in inflow cannula was negligible. The shear stress in inflow cannula was not able to cause hemolysis.Figure 7
Blood damage index of inflow cannula.


## 3.4. Limitations

This study eliminated the factors of real ventricle geometry, insertion depth, and insertion angle. But actually, these factors may change the flow distribution in LV and inflow cannula and eventually cause the suction, and obstruction. So the effect of insertion depth and angle on thrombosis potential, suction and obstruction should be further studied in the future.The other limitation in this study is the simplified boundary conditions with rigid cavity and uniform flow inlet assumption. In reality, the flow in LV is much more complicated. The flow and pressure in ventricle change all the time with the beating. As the purpose of this study is to compare the tip shapes of the inflow cannulas, we think that the simplified boundary conditions are acceptable. Of course, the final outcome should be verified in animal test. And based on the results, profound study with more exact conditions is need in the future.


## 4. Conclusions

Despite the limitation, trumpet tip structure had the best effect on blood compatibility and was hard to generate suction. So trumpet-tipped inflow cannula was the best inflow cannula for LVAD. The caged tip structure had the worst impact on blood compatibility. The thrombosis potential of caged tipped cannula for long-term use was substantial.



---

*Source: 596960-2012-05-10.xml*