Prediction of mechanical hemolysis in medical devices via a Lagrangian strain-based multiscale model.

This work introduces a new Lagrangian strain-based model to predict the shear-induced hemolysis in biomedical devices. Current computational models for device-induced hemolysis usually utilize empirical fitting of the released free hemoglobin (Hb) in plasma from damaged red blood cells (RBCs). These empirical correlations contain parameters that depend on specific device and operating conditions, thus cannot be used to predict hemolysis in a general device. The proposed algorithm does not have any empirical parameters, thus can presumably be used for hemolysis prediction in various blood-wetting medical devices. In contrast to empirical correlations in which the Hb release is related to the shear stress and exposure time without considering the physical processes, the proposed model links flow-induced deformation of the RBC membrane to membrane permeabilization and Hb release. In this approach, once the steady-state numerical solution of blood flow in the device is obtained under a prescribed operating condition, sample path lines are traced from the inlet of the device to the outlet to calculate the history of the shear stress tensor. In solving the fluid flow, it is assumed that RBCs do not have any influence on the flow pattern. Along each path line, shear stress tensor will be input into a coarse-grained (CG) RBC model to calculate the RBC deformation. Then the correlations obtained from molecular dynamics (MD) simulations are applied to relate the local areal RBC deformation to the perforated area on the RBC membrane. Finally, Hb released out of transient pores is calculated over each path line via a diffusion equation considering the effects of the steric hindrance and increased hydrodynamic drag due to the size of the Hb molecule. The total index of hemolysis (IH) is calculated by integration of released Hb over all the path lines in the computational domain. Hemolysis generated in the Food and Drug Administration (FDA) nozzle and two blood pumps, that is, a CentriMag blood pump (a centrifugal pump) and HeartMate II (an axial pump), for different flow regimes including the laminar and turbulent flows are calculated via the proposed algorithm. In all the simulations, the numerical predicted IH is close to the range of experimental data. The results promisingly indicate that this multiscale approach can be used as a tool for predicting hemolysis and optimizing the hematologic design of other types of blood-wetting devices.