Large eddy simulation as a fast and accurate engineering approach for the simulation of rotary blood pumps

Turbulent flow field in maglev centrifugal blood pumps of CH-VAD and HeartMate III: secondary flow and its effects on pump performance

Secondary flow path is one of the crucial aspects during the design of centrifugal blood pumps. Small clearance size increases stress level and blood damage, while large clearance size can improve blood washout and reduce stress level. Nonetheless, large clearance also leads to strong secondary flows, causing further blood damage. Maglev blood pumps rely on magnetic force to achieve rotor suspension and allow more design freedom of clearance size. This study aims to characterize turbulent flow field and secondary flow as well as its effects on the primary flow and pump performance, in two representative commercial maglev blood pumps of CH-VAD and HeartMate III, which feature distinct designs of secondary flow path. The narrow and long secondary flow path of CH-VAD resulted in low secondary flow rates and low disturbance to the primary flow. The flow loss and blood damage potential of the CH-VAD mainly occurred at the secondary flow path, as well as the blade clearances. By contrast, the wide clearances in HeartMate III induced significant disturbance to the primary flow, resulting in large incidence angle, strong secondary flows and high flow loss. At higher flow rates, the incidence angle was even larger, causing larger separation, leading to a significant decrease of efficiency and steeper performance curve compared with CH-VAD. This study shows that maglev bearings do not guarantee good blood compatibility, and more attention should be paid to the influence of secondary flows on pump performance when designing centrifugal blood pumps.

Effect of impeller rotational phase on the FDA blood pump velocity fields

BACKGROUND

The Food and Drug Administration (FDA) blood pump is an open-source benchmark cardiovascular device introduced for validating computational and experimental performance analysis tools. The time-resolved velocity field for the whole impeller has not been established, as is undertaken in this particle image velocimetry (PIV) study. The level of instantaneous velocity fluctuations is important, to assess the flow-induced rotor vibrations which may contribute to the total blood damage.

METHODS

To document these factors, time-resolved two-dimensional PIV experiments were performed that were precisely phase-locked with the impeller rotation angle. The velocity fields in the impeller and in the volute conformed with the previous single blade passage experiments of literature.

RESULTS

Depending on the impeller orientation, present experiments showed that volute outlet nozzle flow can fluctuate up to 34% during impeller rotation, with a maximum standard experimental uncertainty of 2.2%. Likewise, the flow fields in each impeller passage also altered in average 33.5%. Considerably different vortex patterns were observed for different blade passages, with the largest vortical structures reaching an average core radii of 7 mm. The constant volute area employed in the FDA pump design contributes to the observed velocity imbalance, as illustrated in our velocity measurements.

CONCLUSIONS

By introducing the impeller orientation parameter for the nozzle flow, this study considers the possible uncertainties influencing pump flow. Expanding the available literature data, analysis of inter-blade relative velocity fields is provided here for the first-time to the best of our knowledge. Consequently, our research fills a critical knowledge gap in the understanding of the flow dynamics of an important benchmark cardiovascular device. This study prompts the need for improved hydrodynamic designs and optimized devices to be used as benchmark test devices, to build more confidence and safety in future ventricular assist device performance assessment studies.

Resistance valves in circulatory loops have a significant impact on in vitro evaluation of blood damage caused by blood pumps: a computational study

Background: Hemolysis and its complications are major concerns during the clinical application of blood pumps. In-vitro circulatory testing loops have been employed as the key procedure to evaluate the hemolytic and thrombogenic performance of blood pumps during the development phase and before preclinical in-vivo animal studies. Except for the blood damage induced by the pump under test, blood damage induced by loop components such as the resistance valve may affect the accuracy, reproducibility, and intercomparability of test results. Methods: This study quantitatively investigated the impact of the resistance valve on in vitro evaluation of blood damage caused by blood pumps under different operating points. A series of idealized tubing models under the resistance valve with different openings were created. Three pumps - the FDA benchmark pump, the HeartMate 3 LVAD, and the CH-VAD - were involved in hypothetical tests. Eight operating points were chosen to cover a relatively wide spectrum of testing scenarios. Computational fluid dynamics (CFD) simulations of the tubing and pump models were conducted at the same operating points. Results and Conclusion: Overall, hemolysis and platelet activation induced by a typical resistance valve are equivalent to 17%-45% and 14%-60%, respectively, of those induced by the pump itself. Both ratios varied greatly with flow rate, valve opening and pump models. Differences in blood damage levels between different blood pumps or working conditions can be attenuated by up to 45%. Thus, hemolysis and platelet activation induced by the resistance valve significantly affect the accuracy of in-vitro hemocompatibility evaluations of blood pumps. A more accurate and credible method for hemocompatibility evaluations of blood pumps will benefit from these findings.

Influence of Inlet Boundary Conditions on the Prediction of Flow Field and Hemolysis in Blood Pumps Using Large-Eddy Simulation

Inlet boundary conditions (BC) are one of the uncertainties which may influence the prediction of flow field and hemolysis in blood pumps. This study investigated the influence of inlet BC, including the length of inlet pipe, type of inlet BC (mass flow rate or experimental velocity profile) and turbulent intensity (no perturbation, 5%, 10%, 20%) on the prediction of flow field and hemolysis of a benchmark centrifugal blood pump (the FDA blood pump) and a commercial axial blood pump (Heartmate II), using large-eddy simulation. The results show that the influence of boundary conditions on integral pump performance metrics, including pressure head and hemolysis, is negligible. The influence on local flow structures, such as velocity distributions, mainly existed in the inlet. For the centrifugal FDA blood pump, the influence of type of inlet BC and inlet position on velocity distributions can also be observed at the diffuser. Overall, the effects of position of inlet and type of inlet BC need to be considered if local flow structures are the focus, while the influence of turbulent intensity is negligible and need not be accounted for during numerical simulations of blood pumps.

A Hemodynamic Analysis of the Thrombosis Within Occluded Coronary Arterial Fistulas With Terminal Aneurysms Using a Blood Stasis Model

Objective: The aim of this study is to numerically evaluate thrombosis risk within occluded coronary arterial fistulas (CAF) with terminal aneurysms, and provide guidance in choosing occlusion positions, with clinical observations as reference. Method: Four patients with CAF were studied, with different occlusion positions in actual treatments. Hemodynamics simulations were conducted, with blood residue predicted using the blood stasis model. Three types of models (untreated model, aneurysm-reserved model and aneurysm-removed model) were studeid for each patient. Four metrics, i.e., proportion of high oscillatory shear index (OSI), area of high OSI, old blood volume fraction (OBVF)) and old blood volume (OBV) was obtained to distinguish the thrombosis risk of different treatments (proximal or distal occlusion), comparing with the follow-up CTA. Results: For all the postopertive models, the high OBVF, high OSI(>0.3) and low time-averaged wall shear stress (TAWSS) regions were mainly at the distal fistula, indicating these regions were prone to thrombosis. The regions where blood residue remains are roughly regions of high OSI, corresponding well with clinical observations. In contrast, TAWSS failed to distinguish the difference in thrombosis risk. Absolute values (area of high OSI, OBV) can better reflect the degree of thrombosis risk between treatment types compared with percentage values (proportion of high OSI, OBVF). By comparing with the actual clinical treatments and observations, the OBV is superior to the area of high OSI in determining treatment type. Conclusion: The OBV, a volumetric parameter for blood stasis, can better account for the CAF thrombosis and reflect the degree of blood stasis compared with OSI or TAWSS, is a more appropriate metric for thrombosis in the fistula. Together with morphological parameters, the OBV could guide clinicians to formulate more appropriate surgical plans, which is of great significance for the preoperative evaluation and treatment prognosis of CAF patients.

Effect of Conical Spiral Flow Channel and Impeller Parameters on Flow Field and Hemolysis Performance of an Axial Magnetic Blood Pump

For a blood pump, the blood flow channel and impeller parameters directly affect the performance of the pump and the resulting blood circulation. The flow channel in particular has a great impact on the hydraulic performance of the pump (e.g., flow and pressure), which directly determines the overall performance of the blood pump. Traditional bearing-supported blood pumps can cause mechanical damage to blood cells, leading to hemolysis and thrombosis. In this study, therefore, we designed a conical spiral axial blood pump with magnetic levitation. The blood pump was supported by electrodynamic bearings in the radial direction and electromagnetic bearings in the axial direction. The impeller and the front and rear hubs were integrated to minimize blood stagnation and reduce the formation of thrombosis. The hub had a conical spiral flow channel design, which not only reduced the size of the impeller but also increased blood flow and pressure while meeting the design requirements. Computational fluid dynamics (CFD) analysis was used to analyze the flow field of the axial blood pump, a power function model was used to establish a hemolysis prediction model, and the particle tracking method was used to obtain the flow trajectories of individual blood cells, thereby predicting hemolysis-related performance of the blood pump. The simulation results showed that the main high shear stress area in the blood pump was located in the impeller inlet and the clearance between the top of the impeller and the inner chamber of the blood pump. When the hub taper angle of the blood pump was 0.72° and the clearance was 0.3 mm, the average hemolysis prediction value was 0.00216. This prediction value was smaller than that of traditional axial blood pumps. These findings can provide an important reference for the structural design of axial blood pumps and for reducing the hemolysis prediction value.

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.