In fact, if $f(x) = \frac 1 2\|x\|_p^2$ with $p\in (1,2]$, then we always have that $f$ is $(p-1)$-strongly convex wrt to $\|\cdot\|_p$, that is
$$D_f(x,y) \geq (p-1)\|x-y\|_p^2.$$
Proof.
The standard fact from duality says that $\mu$-strong convexity of $f$ wrt the primal norm $\|\cdot\|$ is equivalent to $\frac{1}{\mu}$-smoothness of the conjugate function $f^*$ wrt the dual norm $\|\cdot\|_*$. In other words, the following inequality must hold
\begin{equation}
f^*(u)\leq f^*(v) + \langle \nabla f^*(v),u-v\rangle + \frac{1}{2\mu}\|u-v\|_*^2.
\end{equation}
For twice differentiable $f^*$ it is equivalent to
\begin{equation}\label{eq:2}
\langle \nabla^2f^*(u)(v),v\rangle \leq \frac{1}{\mu} \|v\|_*^2
\end{equation}
This equivalence follows from the Taylor formula.
In our case, $f(x) = \frac 12 \|x\|_p^2$ with $p\in (1,2]$. Then $f^*(u) = \frac 1 2\|u\|_q^2$ with $q=\frac{p}{p-1}\geq 2$, and hence $f^*$ is twice differentiable. Thus, in order to establish the desired, we must prove
$$\langle \nabla^2f^*(u)v,v\rangle \leq \frac{1}{p-1} \|v\|_q^2.$$
This proof is a bit tedious, but overall it is just an application of the Hölder inequality. I refer to the paper A. Ben-Tal et al. "The ordered subsets mirror descent optimization method with applications to tomography" (see Appendix there).