Grid-Induced Numerical Errors for Shear Stresses and Essential Flow Variables in a Ventricular Assist Device: Crucial for Blood Damage Prediction?

Adverse events due to flow-induced blood damage remain a serious problem for blood pumps as cardiac support systems. The numerical prediction of blood damage via computational fluid dynamics (CFD) is a helpful tool for the design and optimization of reliable pumps. Blood damage prediction models primarily are based on the acting shear stresses, which are calculated by solving the Navier–Stokes equations on computational grids. The purpose of this paper is to analyze the influence of the spatial discretization and the associated discretization error on the shear stress calculation in a blood pump in comparison to other important flow quantities like the pressure head of the pump. Therefore, CFD analysis using seven unsteady Reynolds-averaged Navier–Stokes (URANS) simulations was performed. Two simple stress calculation indicators were applied to estimate the influence of the discretization on the results using an approach to calculate numerical uncertainties, which indicates discretization errors. For the finest grid with 19 × 106 elements, numerical uncertainties up to 20% for shear stresses were determined, while the pressure heads show smaller uncertainties with a maximum of 4.8%. No grid-independent solution for velocity gradient-dependent variables could be obtained on a grid size that is comparable to mesh sizes in state-of-the-art blood pump studies. It can be concluded that the grid size has a major influence on the shear stress calculation, and therefore, the potential blood damage prediction, and that the quantification of this error should always be taken into account.