Optimization of a centrifugal blood pump designed using an industrial method through experimental and numerical study

With improved treatment of coronary artery disease, more patients are surviving until heart failure occurs. This leads to an increase in patients needing devices for struggling with heart failure. Ventricular assist devices are known as the mainstay of these devices. This study aimed to design a centrifugal pump as a ventricular assist device. In order to design the pump, firstly, the geometrical parameters of the pump, including the gap distance, blade height, and position of the outlet relative to the blade, were investigated. Finally, the selected configuration, which had all the appropriate characteristics, both hydraulically and physiologically, was used for the rest of the study. The study of the blade, as the main component in energy transfer to the blood, in a centrifugal pump, has been considered in the present study. In this regard, the point-to-point design method, which is used in industrial applications, was implemented. The designer chooses the relationship between the blade angles at each radius in the point-to-point method. The present study selected logarithmic and second-order relations for designing the blade's profile. In total, 58 blades were examined in this study, which differed regarding blade inlet and outlet angles and the relationship between angle and radial position. ANSYS CFX 17.0 software was utilized to simulate blades' performances, and a benchmark pump provided by the US Food and Drug Administration (FDA) was used to validate the numerical simulations. Then, the selected impeller from the numerical investigation was manufactured, and its performance was compared experimentally with the FDA benchmark pump. A hydraulic test rig was also developed for experimental studies. The results showed that among the blades designed in this study, the blade with an input angle of 45° and an output angle of 55°, which is designed to implement a logarithmic relationship, has the best performance. The selected impeller configuration can increase the total head (at least by 20%) at different flow rates compared to the FDA pump.

Design and optimisation of an Intra-Aortic Shrouded rotor axial pump

Undesirable side effects in patients with a LVAD (Left Ventricular Assist Device) pump fitted include blood damage, thrombosis, blood traumatisation, and End-Organ Disfunctions. These side effects have generally been attributed to the high wall shear stresses and the induced turbulent flow. In this study, we introduce a novel design to address these effects by lowering the rotational speed and providing an optimum flow path design to minimise blood damage. We present an initial scheme for a new Intra-Aortic Shrouded Rotary Axial Pump and develop a sequence of pump geometries, for which the Taguchi Design Optimisation Method has been applied. We apply CFD tools to simulate the pressure rise, pump performance, hydraulic efficiency, wall shear stress, exposure time and mass flow rate. A prototype pump has been tested in a mock cardiovascular circuit using a water-glycerol solution. The optimum design delivered the desired pressure/mass flow rate characteristics at a significantly low rpm (2900 rpm). As a result, the estimated blood damage index is low, matching the design requirements. The theoretical performance was matched by experimental results.

The influence of non-conformal grid interfaces on the results of large eddy simulation of centrifugal blood pumps

BACKGROUND

Computational fluid dynamics has been widely used to assist the design and evaluation of blood pumps. Discretization errors associated with computational grid may influence the credibility of numerical simulations. Non-conformal grid interfaces commonly exist in rotary machines, including rotary blood pumps. Should grid size across the interface differ greatly, large errors may occur.

METHODS

This study explored the effects of non-conformal grid interface on the prediction of the flow field and hemolysis in blood pumps using large eddy simulation (LES). Two benchmarks, a nozzle model and a centrifugal blood pump were chosen as test cases.

RESULTS

This study found that non-conformal grid interfaces with considerable change of grid sizes led to discontinuities of flow variables and brought errors to metrics such as pressure head (7%) and hemolysis (up to 14%).

CONCLUSIONS

The results on the full unstructured grid are more accurate with negligible changes of flow variables across the non-conformal grid interface. A full unstructured grid should be employed for centrifugal blood pumps to minimize the influence of non-conformal grid interfaces for LES simulations.

Distribution and history of extensional stresses on vWF surrogate molecules in turbulent flow

The configuration of proteins is critical for their biochemical behavior. Mechanical stresses that act on them can affect their behavior leading to the development of decease. The von Willebrand factor (vWF) protein circulating with the blood loses its efficacy when it undergoes non-physiological hemodynamic stresses. While often overlooked, extensional stresses can affect the structure of vWF at much lower stress levels than shear stresses. The statistical distribution of extensional stress as it applies on models of the vWF molecule within turbulent flow was examined here. The stress on the molecules of the protein was calculated with computations that utilized a Lagrangian approach for the determination of the molecule trajectories in the flow filed. The history of the stresses on the proteins was also calculated. Two different flow fields were considered as models of typical flows in cardiovascular mechanical devises, one was a Poiseuille flow and the other was a Poiseuille–Couette flow field. The data showed that the distribution of stresses is important for the design of blood flow devices because the average stress can be below the critical value for protein damage, but tails of the distribution can be outside the critical stress regime.

Elongational Stresses and Cells

Fluid forces and their effects on cells have been researched for quite some time, especially in the realm of biology and medicine. Shear forces have been the primary emphasis, often attributed as being the main source of cell deformation/damage in devices like prosthetic heart valves and artificial organs. Less well understood and studied are extensional stresses which are often found in such devices, in bioreactors, and in normal blood circulation. Several microfluidic channels utilizing hyperbolic, abrupt, or tapered constrictions and cross-flow geometries, have been used to isolate the effects of extensional flow. Under such flow cell deformations, erythrocytes, leukocytes, and a variety of other cell types have been examined. Results suggest that extensional stresses cause larger deformation than shear stresses of the same magnitude. This has further implications in assessing cell injury from mechanical forces in artificial organs and bioreactors. The cells' greater sensitivity to extensional stress has found utility in mechanophenotyping devices, which have been successfully used to identify pathologies that affect cell deformability. Further application outside of biology includes disrupting cells for increased food product stability and harvesting macromolecules for biofuel. The effects of extensional stresses on cells remains an area meriting further study.

Hemolysis estimation in turbulent flow for the FDA critical path initiative centrifugal blood pump

Hemolysis in medical devices and implants has been a primary concern over the past fifty years. Turbulent flow in particular can cause cell trauma and hemolysis in such devices. In this work, the effects of turbulence on red blood cell (RBC) damage are examined by simulating the flow field through a centrifugal blood pump that has been identified as a case study through the critical path initiative of the US Food and Drug Administration (FDA). In this study, a new model was employed to predict hemolysis in the turbulent flow environment in the pump selected for the FDA critical path initiative. The operating conditions for a centrifugal blood pump were specified by the FDA for rotational speeds of 2500 and 3500 rpm. The model is based on the analysis of the smaller eddies within the turbulent flow field, since it is assumed that turbulent flow eddies with sizes comparable to the dimensions of the RBCs lead to cell trauma. The Kolmogorov length scale of the velocity field is used to identify such small eddies. Using model parameters obtained in prior work through comparisons to capillary and jet flow, it is found that hemolysis for the 2500-rpm pump was predicted well, while hemolysis for the 3500-rpm pump was overpredicted. Results indicate refinement of the model and empirical constants with better experimental data is needed.

Equivalent Scalar Stress Formulation Taking into Account Non-Resolved Turbulent Scales

PURPOSE

Cardiovascular engineering includes flows with fluid-dynamical stresses as a parameter of interest. Mechanical stresses are high-risk factors for blood damage and can be assessed by computational fluid dynamics. By now, it is not described how to calculate an adequate scalar stress out of turbulent flow regimes when the whole share of turbulence is not resolved by the simulation method and how this impacts the stress calculation.

METHODS

We conducted direct numerical simulations (DNS) of test cases (a turbulent channel flow and the FDA nozzle) in order to access all scales of flow movement. After validation of both DNS with literature und experimental data using magnetic resonance imaging, the mechanical stress is calculated as a baseline. Afterwards, same flows are calculated using state-of-the-art turbulence models. The stresses are computed for every result using our definition of an equivalent scalar stress, which includes the influence from respective turbulence model, by using the parameter dissipation. Afterwards, the results are compared with the baseline data.

RESULTS

The results show a good agreement regarding the computed stress. Even when no turbulence is resolved by the simulation method, the results agree well with DNS data. When the influence of non-resolved motion is neglected in the stress calculation, it is underpredicted in all cases.

CONCLUSION

With the used scalar stress formulation, it is possible to include information about the turbulence of the flow into the mechanical stress calculation even when the used simulation method does not resolve any turbulence.

COMPUTATIONAL STUDY ON THE PERFORMANCE OF A CENTRIFUGAL LVAD WITH THE IMPELLER DESIGNED BY INDUSTRIAL METHOD: PROPOSING SIMPLE-TO-MANUFACTURE LVAD’S IMPELLERS

Although the demand of donor hearts for patients with end-stage heart failure is growing, its supply has remained constant. Ventricular assist devices (VADs) provide a chance of finding donor heart by increasing waiting period. In this study, the main goal is to employ an industrial method (point-by-point method) for designing blades profile with a simplified geometry which can be produced by conventional manufacturing methods. In this study, a centrifugal continuous-flow rotary pump is designed and the effects of components’ different geometries on the left ventricular assist devices (LVADs) function are investigated. Moreover, both hydraulic performance and blood damages (hemolysis index (HI)) caused by the pump are considered as design criteria. ANSYS CFX 17 is used to analyze the performance of the designed LVAD. Additionally, the geometry of components are investigated based on fulfilling the required performance of the LVAD while reducing the blood damage level. Comparing the designed VAD with the commercial ones shows that the designed blade further improves the performance of the centrifugal LVAD. Therefore, designing the impeller’s blade profile with point-by-point method seems to be promising. Simplicity in manufacturing is considered to be a big advantage for a design which also leads to lower manufacturing costs. This study demonstrates how industrial design methods can be employed to design simple-to-manufacture impellers which are suitable for LVADs.

Enhanced discrete phase model for multiphase flow simulation of blood flow with high shear stress

Shear stress is often present in the blood flow within blood-contacting devices, which is the leading cause of hemolysis. However, the simulation method for blood flow with shear stress is still not perfect, especially the multiphase flow model and experimental verification. In this regard, this study proposes an enhanced discrete phase model for multiphase flow simulation of blood flow with shear stress. This simulation is based on the discrete phase model (DPM). According to the multiphase flow characteristics of blood, a virtual mass force model and a pressure gradient influence model are added to the calculation of cell particle motion. In the experimental verification, nozzle models were designed to simulate the flow with shear stress, varying the degree of shear stress through different nozzle sizes. The microscopic flow was measured by the Particle Image Velocimetry (PIV) experimental method. The comparison of the turbulence models and the verification of the simulation accuracy were carried out based on the experimental results. The result demonstrates that the simulation effect of the SST k-ω model is better than other standard turbulence models. Accuracy analysis proves that the simulation results are accurate and can capture the movement of cell-level particles in the flow with shear stress. The results of the research are conducive to obtaining accurate and comprehensive analysis results in the equipment development phase.

Computational modeling of the Food and Drug Administration's benchmark centrifugal blood pump

In order to simulate hemodynamics within centrifugal blood pumps and to predict pump hemolysis, CFD simulations must be thoroughly validated against experimental data. They must also account for and accurately model the specific working fluid in the pump, whether that is a blood-analog solution to match an experimental PIV study or animal blood in a hemolysis experiment. Therefore, the Food and Drug Administration (FDA) benchmark centrifugal blood pump and its database of experimental PIV and hemolysis data were used to thoroughly validate CFD simulations of the same blood pump. A Newtonian blood model was first used to compare to the PIV data with a blood analog fluid while hemolysis data were compared using a power-law hemolysis model fit to porcine blood data. A viscoelastic blood model was then incorporated into the CFD solver to investigate the importance of modeling blood's viscoelasticity in centrifugal pumps. The established computational framework, including a dynamic rotating mesh, animal blood-specific fluid properties and hemolysis modeling, and a k-ω SST turbulence model, was shown to more accurately predict pump pressure heads, velocity fields, and hemolysis compared to previously published CFD studies of the FDA centrifugal pump. The CFD simulations were able to match the FDA pressure and hemolysis data for multiple pump operating conditions, with the CFD results being within the standard deviations of the experimental results. While CFD radial velocity profiles between the impeller blades also compared well to the PIV velocity results, more work is still needed to address the large variability among both experimental and computational predictions of velocity in the diffuser outlet jet. Small differences were observed between the Newtonian and viscoelastic blood models in pressure head and hemolysis at the higher flow rate cases (FDA Conditions 4 and 5) but were more significant at lower flow rate and pump impeller speeds (FDA Condition 1). These results suggest that the importance of accounting for blood's viscoelasticity may be dependent on the specific blood pump operating conditions. This detailed computational framework with improved modeling techniques and an extensive validation procedure will be used in future CFD studies of centrifugal blood pumps to aid in device design and predictions of their biological responses.

Deformation of human red blood cells in extensional flow through a hyperbolic contraction

Flow-induced damage to red blood cells has been an issue of considerable importance since the introduction of the first cardiovascular devices. Early blood damage prediction models were based on measurements of damage by shear stress only. Subsequently, these models were extrapolated to include other components of the fluid stress tensor. However, the expanded models were not validated by measurements of damage in response to the added types of stress. Recent investigations have proposed that extensional stress might be more damaging to red cells than shear stress. In this study, experiments were conducted to compare human red cell deformation under laminar extensional stress versus laminar shear stress. It was found that the deformation caused by shear stress is matched by that produced by an extensional stress that is approximately 34 times smaller. Assuming that blood damage scales directly with cell deformation, this result indicates that mechanistic blood damage prediction models should weigh extensional stress more than shear stress.