The design of blood-handling medical devices like implantable ventricular assist devices requires analysis of their hydraulic and hematologic properties. For this purpose, Computational Fluid Dynamics methods are employed to predict blood flow in prototypes. In order to assess biocompatibility, a model is applied to the flow field to estimate red blood cell damage (hemolysis) induced by fluid stress. We propose a new blood damage model formulation. In contrast to related models, it can readily be transformed to an Eulerian frame. This allows for a more direct analysis of critical flow regions like boundary layers and recirculation zones. The model is validated in select benchmark flows. Its discretization is discussed in a stabilized space-time finite element context. The appropriate technique for representing blood cells’ three-dimensional orientation within the flow is investigated with respect to numerical stability and computational efficiency.