The pairwise velocity probability density function in models with local primordial non-Gaussianity