Practical implications of the erroneous treatment of exposure time in the Eulerian hemolysis power law model.

Background: Eulerian and Lagrangian power-law formulations are both widely used for computational fluid dynamics (CFD) to predict flow-induced hemolysis in blood-contacting medical devices. Both are based on the same empirical power-law correlation between hemolysis and the shear stress and exposure time. In the Lagrangian approach, blood damage is predicted by tracking both the stress and exposure time along a finite number of pathlines in the domain. In the Eulerian approach, a scalar transport equation is solved for a time-linearized damage index within the entire domain. Previous analytical work has demonstrated that there is a fundamental problem with the treatment of exposure time in the Eulerian model formulation such that the only condition under which the model correctly represents the true exposure time is in a flow field with no streamwise velocity variation. However, the practical implications of this limitation have yet to be thoroughly investigated.

Methods: In this study, we demonstrate the inaccuracy of Eulerian hemolysis power-law model predictions due to the erroneous treatment of exposure time by systematically considering four benchmark test cases with increasing degrees of flow acceleration: Poiseuille flow through a straight tube, inclined Couette flow, and flow through a converging tube with two different convergence ratios.

Results: Compared with Lagrangian predictions, we show that, as flow acceleration becomes more pronounced, the resultant inaccuracy in the Eulerian predictions increases. For the inclined Couette flow case, there is a small degree of flow acceleration that yields a discrepancy in the range of 10% between Lagrangian and Eulerian predictions. For flows with a larger degree of acceleration, as occurs in the converging tube flow cases, the discrepancy is considerably larger (up to 257%).

Conclusion: The inaccuracy of hemolysis predictions due to the erroneous treatment of exposure time in the Eulerian power-law model can be significant when there is streamwise velocity variation in the flow field. These results may partially explain the extremely large variability in CFD hemolysis predictions reported in the literature between Lagrangian and Eulerian models.