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.