Studies of Turbulence Models in a Computational Fluid Dynamics Model of a Blood Pump

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.

Computational Fluid Dynamics Turbulence Model and Experimental Study for a Fontan Cavopulmonary Assist Device

Head-flow HQ curves for a Fontan cavopulmonary assist device (CPAD) were measured using a blood surrogate in a mock circulatory loop and simulated with various computational fluid dynamics (CFD) models. The tests benchmarked the CFD tools for further enhancement of the CPAD design. Recommended Reynolds-Averaged Navier-Stokes (RANS) CFD approaches for the development of conventional ventricular assist devices (VAD) were found to have shortcomings when applied to the Fontan CPAD, which is designed to neutralize off-condition obstruction risks that could contribute to a major adverse event. The no-obstruction condition is achieved with a von Karman pump, utilizing large clearances and small blade heights, which challenge conventional VAD RANS-based CFD hemodynamic simulations. High-fidelity large eddy simulation (LES) is always recommended; however, this may be cost-inhibitive for optimization studies in commercial settings, thus the reliance on RANS models. This study compares head and power predictions of various RANS turbulence models, employing experimental measurements and LES results as a basis for comparison. The models include standard k-ϵ, re-normalization group k-ϵ, realizable k-ϵ, shear stress transport (SST) k-ω, SST with transitional turbulence, and Generalized k-ω. For the pressure head predictions, it was observed that the standard k-ϵ model provided far better agreement with experiment. For the rotor torque, k-ϵ predictions were 30% lower than LES, while the SST and LES torque values were near identical. For the Fontan CPAD, the findings support using LES for the final design simulations, k-ϵ model for head and general flow simulation, and SST for power, shear stress, hemolysis, and thrombogenicity predictions.

CFD Study of the Effect of the Angle Pattern on Iliac Vein Compression Syndrome

Iliac vein compression syndrome (IVCS, or May-Thurner syndrome) occurs due to the compression of the left common iliac vein between the lumbar spine and right common iliac artery. Because most patients with compression are asymptomatic, the syndrome is difficult to diagnose based on the degree of anatomical compression. In this study, we investigated how the tilt angle of the left common iliac vein affects the flow patterns in the compressed blood vessel using three-dimensional computational fluid dynamic (CFD) simulations to determine the flow fields generated after compression sites. A patient-specific iliac venous CFD model was created to verify the boundary conditions and hemodynamic parameter set in this study. Thirty-one patient-specific CFD models with various iliac venous angles were developed using computed tomography (CT) angiograms. The angles between the right or left common iliac vein and inferior vena cava at the confluence level of the common iliac vein were defined as α1 and α2. Flow fields and vortex locations after compression were calculated and compared according to the tilt angle of the veins. Our results showed that α2 affected the incidence of flow field disturbance. At α2 angles greater than 60 degrees, the incidence rate of blood flow disturbance was 90%. In addition, when α2 and α1 + α2 angles were used as indicators, significant differences in tilt angle were found between veins with laminar, transitional, and turbulent flow (p < 0.05). Using this mathematical simulation, we concluded that the tilt angle of the left common iliac vein can be used as an auxiliary indicator to determine IVCS and its severity, and as a reference for clinical decision making.

Calibration of patient-specific boundary conditions for coupled CFD models of the aorta derived from 4D Flow-MRI

Introduction: Patient-specific computational fluid dynamics (CFD) models permit analysis of complex intra-aortic hemodynamics in patients with aortic dissection (AD), where vessel morphology and disease severity are highly individualized. The simulated blood flow regime within these models is sensitive to the prescribed boundary conditions (BCs), so accurate BC selection is fundamental to achieve clinically relevant results. Methods: This study presents a novel reduced-order computational framework for the iterative flow-based calibration of 3-Element Windkessel Model (3EWM) parameters to generate patient-specific BCs. These parameters were calibrated using time-resolved flow information derived from retrospective four-dimensional flow magnetic resonance imaging (4D Flow-MRI). For a healthy and dissected case, blood flow was then investigated numerically in a fully coupled zero dimensional-three dimensional (0D-3D) numerical framework, where the vessel geometries were reconstructed from medical images. Calibration of the 3EWM parameters was automated and required ~3.5 min per branch. Results: With prescription of the calibrated BCs, the computed near-wall hemodynamics (time-averaged wall shear stress, oscillatory shear index) and perfusion distribution were consistent with clinical measurements and previous literature, yielding physiologically relevant results. BC calibration was particularly important in the AD case, where the complex flow regime was captured only after BC calibration. Discussion: This calibration methodology can therefore be applied in clinical cases where branch flow rates are known, for example, via 4D Flow-MRI or ultrasound, to generate patient-specific BCs for CFD models. It is then possible to elucidate, on a case-by-case basis, the highly individualized hemodynamics which occur due to geometric variations in aortic pathology high spatiotemporal resolution through CFD.

Computational Modeling of the Penn State Fontan Circulation Assist Device

To address the increasing number of failing Fontan patients, Penn State University and the Penn State Hershey Medical Center are developing a centrifugal blood pump for long-term mechanical support. Computational fluid dynamics (CFD) modeling of the Penn State Fontan Circulatory Assist Device (FCAD) was performed to understand hemodynamics within the pump and its potential for hemolysis and thrombosis. CFD velocity and pressure results were first validated against experimental data and found to be within the standard deviations of the velocities and within 5% of the pressures. Further simulations performed with a human blood model found that most of the fluid domain was subjected to low shear stress (<50 Pa), with areas of highest stress around the rotor blade tips that increased with pump flow rate and rotor speed (138-178 Pa). However, the stresses compared well to previous CFD studies of commercial blood pumps and remained mostly below common thresholds of hemolysis and platelet activation. Additionally, few regions of low shear rate were observed within the FCAD, signifying minimal potential for platelet adhesion. These results further emphasize the FCAD's potential that has been observed previously in experimental and animal studies.

Numerical assessment of hemodynamic perspectives of a left ventricular assist device and subsequent proposal for improvisation

Due to the unavailability of donors, the use of left ventricular assist devices has emerged to be a reliable line of alternative treatment for heart failure. However, ventricular assist devices (VAD) have been associated with several postoperative complications such as thrombosis, hemolysis, etc. Despite considerable improvements in technology, blood trauma due to high shear stress generation has been a major concern that is largely related to the geometrical feature of the VAD. This study aims to establish the design process of a centrifugal pump by considering several variations in the geometrical feature of a base design using the commercial solver ANSYS-CFX. To capture the uncertain behavior of blood as fluid, Newtonian, as well as non-Newtonian (Bird-Carreau model), models are used for flow field prediction. To assess the possibility of blood damage maximum wall shear stress and hemolysis index have been estimated for each operating point. The results of the simulations yield an optimized design of the pump based on parameters like pressure head generation, maximum shear stress, hydraulic efficiency, and hemolysis index. Further, the design methodology and the steps of development discussed in the paper can serve as a guideline for developing small centrifugal pumps handling blood.

Effect of Geometric Configuration of the Impeller on the Performance of Liquivac Pump: Single Phase Flow (Water)

Liquivac pumps, with their unique shaped twin start helical rotor, have found utility in various sectors but the major drawback limiting in their global exploitation is their low performance. This paper investigates the study of performance of the Liquivac pump produced by Tomlinson Hall Ltd. Experimental data was used to validate a numerical model developed in Ansys Fluent 20.2 for the Liquivac pump. Four different geometric models of the rotor were tested numerically to find the optimum design using blade number and pitch length as the criteria to achieve improved efficiency. The choice of turbulence model is an important factor in the most accurate prediction with computational fluid dynamics (CFD) simulation. Four different turbulence models were validated with experimental measurements. The realizable K-ε model gave the most accurate performance predictions with a relative deviation of 3.8%. So, the realizable K-ε model was employed for further parametric optimization of the rotor. The results indicate a reasonable improvement in the head and efficiency of the Liquivac pump with a new rotor geometry of four equidistant blades in the front, back and four flights with 30 mm pitch. This is attributed to the most favourable balance between the different losses and most guided and uniform flow inside the rotor channels.

The effects of non-Newtonian blood modeling and pulsatility on hemodynamics in the food and drug administration's benchmark nozzle model

BACKGROUND

Computational fluid dynamics (CFD) is an important tool for predicting cardiovascular device performance. The FDA developed a benchmark nozzle model in which experimental and CFD data were compared, however, the studies were limited by steady flows and Newtonian models.

OBJECTIVE

Newtonian and non-Newtonian blood models will be compared under steady and pulsatile flows to evaluate their influence on hemodynamics in the FDA nozzle.

METHODS

CFD simulations were validated against the FDA data for steady flow with a Newtonian model. Further simulations were performed using Newtonian and non-Newtonian models under both steady and pulsatile flows.

RESULTS

CFD results were within the experimental standard deviations at nearly all locations and Reynolds numbers. The model differences were most evident at Re = 500, in the recirculation regions, and during diastole. The non-Newtonian model predicted blunter upstream velocity profiles, higher velocities in the throat, and differences in the recirculation flow patterns. The non-Newtonian model also predicted a greater pressure drop at Re = 500 with minimal differences observed at higher Reynolds numbers.

CONCLUSIONS

An improved modeling framework and validation procedure were used to further investigate hemodynamics in geometries relevant to cardiovascular devices and found that accounting for blood's non-Newtonian and pulsatile behavior can lead to large differences in predictions in hemodynamic parameters.

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.

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.

Mesh sensitivity analysis for quantitative shear stress assessment in blood pumps using computational fluid dynamics

The reduction of excessive, nonphysiologic shear stresses leading to blood trauma can be the key to overcome many of the associated complications in blood recirculating devices. In that regard, Computational Fluid Dynamics (CFD) are gaining in importance for the hydraulic and hemocompatibility assessment. Still, direct hemolysis assessments with CFD remain inaccurate and limited to qualitative comparisons rather than quantitative predictions. An underestimated quantity for improved blood damage prediction accuracy is the influence of near-wall mesh resolution on shear stress quantification in regions of complex flows. This study investigated the necessary mesh refinement to quantify shear stress for two selected, meshing sensitive hotspots within a rotary centrifugal blood pump. The non-dimensional mesh characteristic number y+, which is known in the context of turbulence modelling, underestimated the maximum wall shear stress by 60% on average with the recommended value of 1, but was found to be exact below 0.1. To evaluate the meshing related error on the numerical hemolysis prediction, three-dimensional simulations of a generic centrifugal pump were performed with mesh sizes from 3 to 30 million elements. The respective hemolysis was calculated using an Eulerian scalar transport model. Mesh insensitivity was found below a maximum y+ of 0.2 necessitating 18 million mesh elements. A meshing related error of up to 25% was found for the coarser meshes. Further investigations need to address: 1) the transferability to other geometries and 2) potential adaptions on blood damage estimation models to allow better quantitative predictions.

HEMOLYSIS ANALYSIS OF AXIAL BLOOD PUMPS WITH VARIOUS STRUCTURE IMPELLERS

Low hemolysis is an important factor for axial blood pumps that has been used in patients with heart failure. The structure of impellers plays a key role in the hemolytic properties of axial blood pumps. Axial blood pumps with various structure impellers exhibit different hemolytic characteristic. In the present study, we aimed to investigate the type of impellers structures in axial blood pumps that contain the best low hemolytic properties. Also, it is expensive and time-consuming to validate the axial blood pump's hemolytic property by in vivo experiments. Therefore, in the present study, the numerical method was applied to analyze the hemolytic property in a blood pump. Specifically, the hemolysis of the pump was calculated by using a forward Euler approach based on the changes in shear stress and related exposure times along the particle trace lines. The different vane structures and rotational speed that affect hemolysis were analyzed and compared. The results showed that long–short alternant vanes exhibited the best hemolytic property which could be utilized in the optimization design of axial blood pumps.

NUMERICAL STUDIES OF AN AXIAL FLOW BLOOD PUMP WITH DIFFERENT DIFFUSER DESIGNS

Most axial flow blood pumps basically consist of a straightener, an impeller, and a diffuser. The diffuser plays a very important role in the performance of the pump to provide an adequate pressure head and to increase the hydraulic efficiency. During the development of an axial flow blood pump, irregular flow field near the diffuser hub is not desirable as it may induce thrombosis. In order to avoid this phenomenon, two approaches were adopted. In the first approach, the number of the diffuser blades was increased from three (B3, baseline model) to five (B5 model). It was observed that the flow field was improved, but the irregular flow patterns were not completely eliminated. In the second approach, we detached the blades from the diffuser hub (B3C2 model), which was integrated and rotated with the impeller hub. It was found that the rotary diffuser hub significantly improved the flow field, especially near the diffuser hub. Besides the detailed flow fields, the hydraulic and hematologic performances at various flow conditions were also estimated using computational fluid dynamics (CFD). Although each design has its own advantages and disadvantages, the B5 model was superior based on a comparative overview.

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.

Shape optimization of the diffuser blade of an axial blood pump by computational fluid dynamics

Computational fluid dynamics (CFD) has been a viable and effective way to predict hydraulic performance, flow field, and shear stress distribution within a blood pump. We developed an axial blood pump with CFD and carried out a CFD-based shape optimization of the diffuser blade to enhance pressure output and diminish backflow in the impeller-diffuser connecting region at a fixed design point. Our optimization combined a computer-aided design package, a mesh generator, and a CFD solver in an automation environment with process integration and optimization software. A genetic optimization algorithm was employed to find the pareto-optimal designs from which we could make trade-off decisions. Finally, a set of representative designs was analyzed and compared on the basis of the energy equation. The role of the inlet angle of the diffuser blade was analyzed, accompanied by its relationship with pressure output and backflow in the impeller-diffuser connecting region.

Computational fluid dynamics investigation of a centrifugal blood pump

In the development of a ventricular assist device, computational fluid dynamics (CFD) analysis is an efficient tool to obtain the best design before making the final prototype. In this study, different designs of a centrifugal blood pump were developed to investigate flow characteristics and performance. This study assumed the blood flow as being an incompressible homogeneous Newtonian fluid. A constant velocity was applied at the inlet; no slip boundary conditions were applied at device wall; and pressure boundary conditions were applied at the outlet. The CFD code used in this work was based on the finite volume method. In the future, the results of CFD analysis can be compared with flow visualization and hemolysis tests.

A Validated Computational Fluid Dynamics Model to Estimate Hemolysis in a Rotary Blood Pump

A major part of developing rotary blood pumps requires the optimization of hemolytic properties of the entire pump. Application of a suited computational fluid dynamics (CFD)-based hemolysis model allows approximation of blood damage in an early phase of the design process. Thus, a drastic reduction of time- and cost- intensive hemolysis experiments can be achieved. For the MicroDiagonal Pump (MDP), still under development at Helmholtz-Institute in Aachen, Germany, different pump configurations have been analyzed, both numerically and experimentally. The CFD model of the pump has been successfully validated based on the comparison of the pressure head curves (H-Q curves), as discussed in a prior publication. In the present study, the authors focus on the development of a semiempiric blood damage model using the CFD and in vitro hemolysis data. On the one hand, mean key characteristics (shear stress and exposure time) and other characteristics affecting blood damage have been calculated based on numerical data. On the other hand, in vitro hemolysis tests have been accomplished in order to determine the hemolytic curves of two different pump configurations (with the same impeller but different tip clearances). Finally, a new function based on a general power law has been defined by means of the mean key characteristics. The unknown constants of the function have been determined by multidimensional regression analysis using the hemolytic curves. For the final validation of this new blood damage model, the calculated and the in vitro obtained hemolysis indices at the specific VAD operating point have been compared for all pump configurations. The comparison showed an excellent agreement, both qualitatively and quantitatively.