The following identity is used in the treatment of the multivariate wave equation (Evans on p. 74)
Lemma 2 (Some useful identities). Let $\phi:\Bbb R\rightarrow \Bbb R$ be $C^{k+1}$. Then for $k=1,2,\dots$
$$\left(\frac{d^2}{dr^2}\right)\left(\frac 1r\frac{d}{dr}\right)^{k-1}\left(r^{2k-1}\phi(r)\right)=\left(\frac 1r\frac{d}{dr}\right)^{k}\left(r^{2k}\phi'(r)\right)$$
For $k=1$, both sides equal $2\phi'+r\phi''$; for $k=2$, they equal $8\phi'+7r\phi''+r^2\phi'''$. Abbreviating $\frac{d}{dr}$ with $d$, I think that the following identity will be useful in a proof by induction
\begin{aligned} d \frac 1r d &= \frac 1r d^{-1}d^2\\ &=\left(\frac 1r d + \frac 1{r^2}\right)d^{-1}d^2 \\ &=\frac 1r d^2 + \frac 1{r^2}d \end{aligned}
From this, it also follows that $d^2=rd\frac 1r d-\frac 1r d$. However, I tried to insert these identities in multiple places without success.
I suppose $r\ne 0$ here. Fix $r_0\ne 0$. Let $P$ be the $(k+1)$th Taylor polynomial of $\phi$ at $r_0$, i.e., $P^{(j)}(r_0) = \phi^{(j)}(r_0)$ for $j=0, \dots, k+1$. Note that if we expand both sides of the claimed identity, they will have derivatives of $\phi$ only up to order $k+1$. Therefore, both sides, evaluated at $r_0$, remain the same if we replace $\phi$ by $P$.
So if suffices to prove the identity for polynomials. For that, it suffices to prove it for monomials (by linearity). For a monomial $r^n$, the statement becomes $$\left(\frac{d^2}{dr^2}\right)\left(\frac 1r\frac{d}{dr}\right)^{k-1}\left(r^{2k+n-1}\right) = n\left(\frac 1r\frac{d}{dr}\right)^{k}\left(r^{2k+n-1}\right)\tag1$$ Let $$Q(r) = \left(\frac{d}{dr}\right)\left(\frac 1r\frac{d}{dr}\right)^{k-1}\left(r^{2k+n-1}\right)$$ This is a monomial of degree $2k+n-1 - 2(k-1) -1 =n$. That is, $Q$ is $cr^n$ for some constant $c$. In terms of $Q$, the identity (1) becomes $$ \frac{d}{dr}Q = n\frac{1}{r}Q $$ which is obvious: both sides are $ncr^{n-1}$.