This work proposes a lumped parameter model (0D) for biological prosthetic heart valves which represents the flow regurgitation. This model is used to obtain a more precise model of a pulsatile pediatric ventricular assist device (pVAD). The 0D model consists of the cannula-valve inlet and cannula-valve outlet configurations, which represent the discrete transfer functions related to the conduction and non-conduction states of the input and output valves, respectively. The flow and pressure curves in the inlet and outlet cannulas of the heart valves are measured from a hydraulic simulator connected to a pVAD. To find this model, a sequence of steps is performed, which includes the data pre-processing, the structure model choice and the refinement of model parameters. From the simulation of the refined model, the regurgitation phenomenon was represented with fidelity of 90 % of the value measured in the nominal condition test. The great importance of pVAD modeling is the possibility of coupling to the 0D model of the human cardiovascular system so that its interactions can be studied and automatic control strategies developed.