Numerical investigation of volute tongue design on hemodynamic characteristics and hemolysis of the centrifugal blood pump

In the design of rotary blood pumps, the optimization of design parameters plays an essential role in enhancing the hydrodynamic performance and hemocompatibility. This study investigates the influence of the volute tongue angle as a volute geometric parameter on the hemodynamic characteristics of a blood pump. A numerical investigation on five different versions of volute designs is carried out by utilizing a computational fluid dynamics (CFD) software ANSYS-FLUENT. The effect of volute tongue angle is evaluated regarding the hydrodynamic performance, circumferential pressure distribution, the radial force, and the blood damage potential. A series of volute configurations are constructed with a fixed radial gap (5%), but varying tongue angles ranging from 10 to 50°. The relative hemolysis is assessed with the Eulerian based empirical power-law blood damage model. The pressure-flow rate characteristics of the volute designs at a range of rotational speeds are obtained from the experimental measurements by using the blood analog fluid. The results indicate an inverse relationship between hydraulic performance and the tongue angle; at higher tongue angles, a decrease in performance was observed. However, a higher tongue angle improves the net radial force acting on the impeller. The pump achieves the optimized performance at 20° of the tongue angle with the relatively high hydrodynamic performance and minor blood damage risk.


Introduction
Continuous-flow cardiac assist pumps have been successfully used to treat end-stage heart failure patients as destination therapy, bridge-to-transplantation or bridgeto-recovery [1][2][3]. However, hemolysis reduction is still a significant design challenge for blood pumps [2,3]. Even though the impeller geometry plays an essential role in the overall performance of a centrifugal blood pump [4], the volute structure is equally critical for assuring steady flow lines and good kinetic energy to potential energy conversion for achieving greater hydraulic performance and a lower hemolysis risk [5,6]. Moreover, the uniformity of the pressure distribution around the impeller is a crucial parameter for designing non-contact bearings for centrifugal pumps. Non-axisymmetric nature of the volute geometry generates a non-uniform pressure distribution around the impeller, thus inducing a radial force critical for the magnetic suspension of pumps. Several researchers have studied the effect of volute geometry on the performance of centrifugal pumps at larger scales to address this problem [7,8]. The volute tongue region has been particularly studied in various engineering fields for improving the pump endurance, bearing life, noise reduction, magnetic levitation and the flow separation [7]. However, recent studies on rotary blood pumps have mostly focused on the impeller geometry and blade parameters to optimize the pump performance and minimize the blood damage risk [4,9]. To date, a few researchers have numerically investigated the hydraulic effect of volute parameters such as the volute tongue in rotary blood pumps [5], yet in these studies the effect of volute tongue position on the hemolytic characteristics is missing.
Computational fluid dynamics (CFD) has been broadly used for the design and optimization of rotary blood pumps [9][10][11]. This method is a reliable technique to investigate the overall flow characteristics, and to simplify the development period of the device, reducing the time and cost of experimentation [2,3]. As the blood damage potential has recently been incorporated into these CFD analyses, various numerical models have been adopted by several research groups to estimate the hemolysis and optimize the design of rotary blood pumps [10,11]. Although the researchers are still questioning the accuracy of the computational blood damage models for predicting the absolute level of hemolysis, it has become an accepted method for calculating the relative hemolysis level for the design optimization [10,12]. The ISO 14,708-5 standard also suggests using the CFD simulations for evaluating the relative design changes rather than determining the absolute quantities [13].
In this study, five volute configurations were constructed in a virtual environment with the tongue angles (TA) ranging from 10° to 50° in 10° increments and numerical simulations were conducted to evaluate the effect of tongue position. Firstly, the hydrodynamic performance and the relative hemolysis characteristic were numerically investigated by employing the Eulerian blood damage estimation approach in a CFD model. In the next step, the presented models were rapid-prototyped and a hydraulic experiment was conducted with a blood analog fluid to analyze their pressure-volumetric flow rate characteristics. Finally, the hemolytic and hydraulic performance of volute geometries were compared based on the in-silico findings and in-vitro measurements. The objective of this study is to investigate the influence of the tongue angle on the pump hemodynamics and optimizing the tongue design for a better matching volute and impeller region.

Pump geometry
The physical design of the pump is comprised of a rotating impeller with five blade-splitters, an inlet guided by three straight blades, a volute converting the kinetic energy to the potential energy and two ceramic bearings for contact regions. The complete pump assembly is shown in Fig. 1. The pump employs the concept of magnetic coupling with the magnets embedded inside the impeller. A miniature DC-motor drives the impeller through magnetic coupling. The pump used in this study is 25 mm in height, 51 mm in length and houses a 45 mm diameter impeller inside. The major specifications of the pump are summarized in Table 1. The pump is designed to operate at a rotational speed of 2000 rpm, targeting a 5 L/min pump flow rate and 100 mmHg head pressure according to the clinical use at the hospitals [14,15]

Volute design
The flow pattern inside a centrifugal pump is quite intricate, and the pump characteristics are critically affected by the geometrical constraints such as impeller-volute interactions [7]. The pump geometry was designed using a CAD software (Siemens NX 10.0, Siemens AG, Munich, Germany) with a configuration to create an adjustable impeller-volute region. For the present study, the constant mean velocity method was used to design the volute geometry according to the Stepanoff design rule [16]. The pump specific speed (Ns) was set to 790 (m 3 /h, m) for all models based on Eq. (1). The tongue shape and the circumferential angular position (tongue angle, θ) were shown in Fig. 2.
Radial gap (δ) between the impeller periphery and volute tongue was identified as a percentage of the impeller outer diameter based on the guideline [16].
A total of five volute models, ranging from 10° to 50° tongue angles, were constructed for the specified pump geometry. The tip of the volute tongue was round, and the lowest volute tongue angle was defined at 10°. Angles less than 10° were disregarded to avoid an excessive sharp corner on the tongue tip. As the tongue angle ranged from 10° to 50°, the tongue thickness was varied based on the diameter of the arc. The radial gap (δ) between the proposed volute designs was kept constant at 5% of the impeller outer diameter by using Eq. (2). The volute outlet had a circular cross-section with a 9 mm end diameter. The tongue position was characterized by only changing the tongue angle, and the relative positions of other geometrical components such as volute throat and impeller-volute regions were fixed.

CFD method and boundary conditions
In this study, we employed a previously developed CFD Model and this model was validated against previously published experimental study of the FDA centrifugal pump [4]. A commercial software (ANSYS Fluent 18.0, ANSYS Inc., Pittsburgh, PA, USA) was used to model the flow pattern inside the blood pump. The sliding mesh method was applied between the rotating impeller and the stator casing to simulate the flow dynamics. The blood was modeled as an incompressible Newtonian fluid with 1050 kg/m3 density and 0.0035 Pas dynamic viscosity.
Although the blood is known as non-Newtonian fluid, it only exhibits non-Newtonian properties at low shear rates (< 100/s) [17]. Hence, the Newtonian fluid characteristic is adopted for the blood flow due to the high shear rates inside the rotary blood pumps [10]. During the grid generation process, orthogonal quality, aspect ratio, and the skewness were assured to be in the required range. Turbulence was simulated with the k-ω model by considering the mean γ + < 3 value with appropriately sized cells near domain walls. The final meshes contained approximately 10 million cells. The convergence precision was set to 10 -4 for the scaled residuals. The inlet and outlet ports were extended to a specific length to prevent boundary condition effects on the convergence. The impeller was set to complete one rotation in 100 steps indicating 3.6° rotation for each step, and the CFD time step was adjusted to 3.00 × 10 -4 s accordingly. The inlet boundary condition was set to flow rate of 5 L/min, and the outlet pressure was set to 0 Pa as nominal values at 2000 rpm. The absolute pressure difference was monitored between the inlet and outlet port. Hydraulic efficiency and the net radial force acting on the impeller were numerically calculated. The physical quantities such as shear stress, pressure and the torque values on the rotation axis were monitored to detect the flow conditions after 3  revolutions. User-Defined Functions (UDF) were compiled and used for the quantification of the blood damage.

Hemolysis estimation
An empirical power-law model has mostly been used to compute the hemolysis in blood pumps [10]. The hemolysis potential of the pump is quantified as a blood damage index (D) which indicates the increase in plasma-free hemoglobin. The power law model predicts the blood damage by correlating the plasma-free hemoglobin level with the shear stress and exposure time [18], where C, α and β are the constants obtained by postprocessing the experimental data, τ is the scalar stress, t is the exposure time of blood cell, dHb represents the damaged hemoglobin content, and Hb indicates the hemoglobin content. As these constants vary based on the conducted experiment, several researchers have proposed different values [10]. The constants Heuser et al. C = 1.8 × 10 −6 , α = 0.765, β = 1.991 [19] suggest (C = 1.8 × 10-6, α = 0.765, β = 1.991) were utilized in this numerical study due to the reported accuracy over experimental data [10].
Scalar shear stress (SSS) parameter can be found by the following expression [17]; The power-law equation was linearized to obtain the concentration of damaged blood [20,21]. This method directly computes the damage index from the velocity and the shear stress field by the volume integration over the entire flow domain. The time dependency of the power law was eliminated by the linearization [20,21]. The linearized blood damage index was then expressed as: A hyperbolic equation was also obtained by taking the time derivative of Eq. (6) and applying the Eulerian transform [20,21].

In vitro hydraulic test
A hydraulic test loop was constructed to validate and compare the hydraulic performance of the models. The loop contains a pump prototype, two piezoelectric pressure transducers, an ultrasonic flowmeter, a DC motor, a motor controller box, a DAQ system, a total of 2 m length flexible polyvinyl chloride tubing (3/8-inch ID) and a reservoir. The pump prototypes were manufactured from ABS plastic with a 3D printer (Ultra 3SP Envisiontec GMBH, Gladbeck, Germany). The brushless DC motor (Maxon Motor AG, Sachseln, Switzerland) was magnetically coupled to the impeller. A motor controller was implemented for controlling and reading the speed of the impeller. A water-glycerol mixture was prepared as a blood analog fluid (65% water, 35% glycerol) at the room temperature (25 °C) in accordance with the CFD simulations. As the fluid traveled through the loop, a valve was used for adjusting the resistance to the flow. The flow rate was evaluated via the ultrasonic flowmeter (FD-Q 10C, Keyence Corp., Osaka, Japan) which was placed on the outlet tubing. Two piezoelectric pressure transducers (RAG25R0, Kistler Holding AG, Winterthur, Switzerland) were placed through the inlet and the outlet tubing. A DAQ card (LabView National Instruments, Austin, TX, USA) and the interface were used to collect pressure readings with a 10 kHz sampling rate.

Results
The hydraulic and the hemolytic performances of five volute models were compared at the design point flow rate of 5 L/min. Table 2 illustrates the overall hydraulic performance of the models. The tongue angle was varied from 10° to 50° in 10° increments. The pressure head and the hydraulic efficiency were observed to decrease correspondingly as the tongue angle increase. The poorest hydraulic performance was realized in the 50° tongue angle with an efficiency of 44.8% and a pressure head of 93.2 mmHg, but the radial force was lowest at 0.2 N among all cases.
Pressure distributions of five different cases were compared in Fig. 3. The absolute pressure contours of the midvolute layer were plotted for all cases, respectively. A more uniform pressure field was observed around the impeller region at the 50° tongue angle. Moreover, an increasing tongue angle reduced the net radial force acting on the impeller. The velocity vector fields were displayed in Fig. 4. The relative velocity field presents the flow interaction between the outlet and the impeller outflow. At a low tongue angle, the flow was observed as well-guided along the volute region through the outlet. However, the flow was separated through to the suction side of the blades as the tongue angle increased. The findings of hydraulic experiments were also reported in Fig. 5. The pressure-flow performance was measured at 2000 rpm. The performance curves revealed that the configuration of 50° tongue angle had the poorest hydraulic output, and the hydraulic performance was observed to be negatively influenced by a larger tongue angle. The SSS was also computed based on Eq. (5). High shear stress zones were mostly concentrated around the trailing edges of the blades and the volute tongue region. The maximum and averaged SSS in the entire flow domain were shown in Fig. 6.
The most extensive shear stress was observed at the 50° of tongue angle, for which the maximum SSS was found to be 354.4 Pa, and the average SSS was calculated as 13.9 Pa. By comparing the relative blood damage of various tongue angles in Fig. 7, the minimum Relative Hemolysis Index (RHI), indicated to be optimal for the 20° tongue angle. The highest RHI was observed when the tongue angle was 50° which means a greater possibility of blood damage for larger tongue angles.

Discussion
In general, the comparison of five models revealed that the major differences might be observed in the pump performance depending on the volute profile. Pressure contours indicated that increasing the tongue angle promotes the decrease in the corresponding net radial force on the impeller. However, the hydraulic efficiency and the head pressure decreased at higher tongue angles. Thus, the hydraulic experiments confirmed the in-silico characteristics of the pumps.
It was observed that the SSS distributions were negatively impacted by the larger tongue angle. Nevertheless, the lowest tongue angle (10°) generated a relatively high SSS value of 347.8 Pa. This observation might be explained by the fact that too small of a tongue angle leads to an unusually sharp corner on the volute tongue region where the fluid is subjected to intense shear. Consequently, the intense shearing motion results in a high shear rate and excessive acceleration in a small area. The tongue angles of 10° and 20° were observed to be superior among the models regarding the hydraulic output and the blood damage.
Although the hydraulic output was slightly better for the 10° tongue angle, 20° tongue angle represents 13% lower hemolysis than that of 10°. However, the higher net radial force acting on the impeller was observed in the lower tongue angles. Near the volute tongue region and the trailing edge of blades, the blood damage risk is usually higher due to shear stress. Velocity vector fields showed that a large tongue angle resulted in flow separation. A well-matching volute tongue and impeller outflow should minimize the flow separation and the net radial force acting on the impeller. Eventually, a lower blood damage can be achieved. Henceforth, factors such as the hydraulic output, the net radial force and the amount of blood damage are critical indicators for pump designers.
The presented study reveals encouraging outputs, but there are also some limitations in the experimental study. the tolerances of pump geometry were around 100μ and removal of the support material was achieved by hand polishing. Maximum measurement errors of flow meter and pressure sensors were ± 0.15 L/min and ± 1.875 mmHg, respectively, corresponds to the 0.5% of the full scale. Nevertheless, the hydraulic performance of the physical prototypes at the nominal condition was in good agreement with the in-silico analysis. Overall, the findings of this study do not necessarily put forth a definitive tongue angle in the pump design, but the reported effects of the tongue angle can serve as a guideline for the pump optimization. The tongue angle should be adjusted based on the design requirements such as hemolysis, hydraulic efficiency, magnetic suspension force, etc. Moreover, the reported radial force data is limited to highlighting the impact of volute tongue on the non-uniform pressure distribution around the impeller region. Further unsteady analyses at different flow rates might be needed to assess the magnetic suspension requirements for active levitation control. It should also be noted that the evaluation of the accuracy and credibility of the presented hemolysis estimation approach is beyond the scope of this study. This estimation method was adopted from the existing studies and implemented as a CFD tool for ranking the volute geometries.

Conclusion
In summary, this study explores impact of the volute tongue position on the hydraulic and hemolytic performance of a centrifugal blood pump. The tongue angle of 20° demonstrated the minimum blood damage risk among five models (10° to 50°) with an adequate hydraulic outputs. The blood damage potential and radial imbalance increases gradually as the tongue angle grows larger than 20°. Although the 10° tongue angle represents the best hydraulic output and efficiency, yet it indicates slightly higher hemolysis index compared to 20° tongue angle. Our findings showed a maximum increase in the hydraulic efficiency and pressure head up to %9.1 and %8.8 respectively when the volute tongue was modified. However, the same modifications resulted a substantial improvement of in-silico hemolysis estimation up to 77.8% at the design point.
This observation proves valuable for further research into this subject by repeating this study by considering the limitations. In our research we observed this by working on blood analogue fluid is used for the physical tests and the hydraulic efficiency was estimated from the in-silico findings. Also, depending on the impeller-volute interaction type and dimensional tolerance, the optimal volute tongue angle could vary significantly between different blood pump designs and further research is required considering this variation. Acknowledging these limitations, this study does not claim to establish a gold standard method for designing a volute geometry, rather it reveals the effect of the tongue position on a fixed impeller-varying volute tongue couple. Future studies with in-vitro blood tests and different impeller-volute examples may reveal further insights about the volute tongue position in LVAD design. Finally, this study provides an original insight to blood pump designers for finding the best volute tongue position to reduce the blood damage risk, increase the efficiency, minimize the power needs of a magnetic bearing or decrease the possible wear of the mechanical bearings.
Data availability and material All authors confirm that all data and materials as well as software application or custom code support their published claims and comply with field standards

Compliance with ethical standards
Conflicts of interest The authors declare that they have no conflict of interest. Consent for publication All authors agreed with the content and that all gave explicit consent to submit and that they obtained consent from the responsible authorities at the institute/organization where the work has been carried out.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.