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.

Asymptotically consistent numerical approximation of hemolysis.

In a previous communication, we have proposed a numerical framework for the prediction of in vitro hemolysis indices in the preselection and optimization of medical devices. This numerical methodology is based on a novel interpretation of Giersiepen-Wurzinger blood damage correlation as a volume integration of a damage function over the computational domain. We now propose an improvement of this approach based on a hyperbolic equation of blood damage that is asymptotically consistent. Consequently, while the proposed correction has yet to be proven experimentally, it has the potential to numerically predict more realistic red blood cell destruction in the case of in vitro experiments. We also investigate the appropriate computation of the shear stress scalar of the damage fraction model. Finally, we assess the validity of this consistent approach with an analytical example and with some 3D examples.