We know the following:
$$\boldsymbol\nabla = \hat r \frac{\partial }{\partial r} + \hat \theta \frac1r \frac{\partial }{\partial \theta} + \hat\phi\frac{1}{r\sin\theta} \frac{\partial }{\partial \phi}$$
It's clear that: $ \hat{r}·\vec{p} \; \phi =-i\hbar \frac{\partial \phi}{\partial r}$
But:
$$\begin{align}\vec p\cdot \hat r \Psi & = -~i\hbar \boldsymbol \nabla\cdot (\hat r \Psi)\\&= -~i\hbar ((\boldsymbol \nabla \Psi)\cdot \hat r +(\boldsymbol\nabla \cdot \hat{r})\Psi)\\ &=-~i\hbar \left(\frac{\partial}{\partial r}\Psi +\frac2r \Psi\right)\\ &\ne \hat r\cdot \vec p \cdot \Psi \end{align}$$
My question is, why $\nabla · \hat{r} = \frac{2}{r}$ ?