This paper numerically evaluates the accuracy and performance of a stabilized finite element Reduced Order Modelling (ROM) approach that is designed to simulate pulsatile blood flows. The method is able to estimate fluid flow parametric solutions of interest in hemodynamics by considering the hematocrit percentage (Hct%) and the heart rate (f(c)) as parameters. The method estimates off-trained parametric scenarios which have not been included in the training data set composing the ROM basis and that can adopt arbitrary values from specific patient conditions. The hematocrit percentages modifies the viscous properties of blood, which are incorporated into the problem physics through a power-law model. The hematocrit percentages are 5 = Hct% = 50, ranging from severe anemia to physiological values. The pulsatile flow condition is modeled via a Womersley function, with f(c)bpm (beats per minute) ranging between 60bpm = f(c) = 120bpm. These frequencies comprise physiological conditions of patients at rest and tachycardia or moderate exercise. Two-and-three-dimensional flows are simulated inside a representative geometry of a carotid artery, where the accuracy is tested by studying the effect of the number of components of the ROM approximating basis and the number of sampling points composing the training data set. The parametric calculation verifies that the proposed approach constitutes a valuable computational tool for simulating complex fluid flow hemodynamics and can be applied to clinical decision-making.