Mechanical hemolysis in blood flow: user-independent predictions with the solution of a partial differential equation

This paper presents for the first time numerical predictions of mechanical blood hemolysis obtained by solving a hyperbolic partial differential equation (PDE) modelling the hemolysis in a Eulerian frame of reference. This provides hemolysis predictions over the entire computational domain as an alternative to the Lagrangian approach consisting in evaluating cell hemolysis along their trajectories. The solution of a PDE over a computational domain, such as in the approach presented herein, yields a unique solution. This is a clear advantage over the Lagrangian approach, which requires the human-made choice of a limited number of trajectories for integration and inevitably results in the incomplete coverage of the computational domain. The hyperbolic hemolysis model is solved with a Discontinuous Galerkin finite element method. The solution algorithm also includes adaptive remeshing to provide high accuracy simulations. Predictions of the modified index of hemolysis (MIH) are presented for flows in dialysis cannulae and sudden contractions. MIH predictions for cannulae differ significantly from those obtained by other authors using the Lagrangian approach. The predictions for flows in sudden contractions are used, along with our own experimental measurements, to assess the value of the threshold shear stress required for hemolysis that is included in the hemolysis model.