There is significant interest in the development and application of variable speed impeller-pump type ventricular assist devices designed to generate pulsatile blood flow. However, no study has so far been carried out to investigate the systemic cardiovascular response to various aspects of pump motion. In this article, a numerical model is constructed for the simulation of the cardiovascular response in the heart failure condition under representative cases of pulsatile impeller pump support. The native cardiovascular model is based on a previously validated model, and the impeller pump is modeled by directly fitting the pressure-flow curves that describe the pump characteristics. The model developed is applied to study circulatory dynamics under different degrees of phase shift and pulsation ratio in the pump motion profile. The characteristic variables are discussed as criteria for the evaluation of system response for comparison of the pulsatile flows. Simulation results show that a constant pump speed is the most efficient work mode for the rotary pump, and with the application of either a phase shift of 75% and a pulsation ratio of 0.5, or a phase shift of 42% and a pulsation ratio of 0.55, it is possible to generate arterial pulse pressure with the maximal magnitude of about 28 mmHg. However, this is achieved at the cost of reduced cardiac output and pump efficiency.