I need to calculate how the interior product $\iota_{\xi}$ for some vector field $\xi$ acts on the inner product of two forms. That is, if $\alpha$ and $\beta$ are both $p$-forms, then how do I calculate
$$\iota_{\xi}\left(\alpha\wedge\star\beta\right)?$$
In particular, I am mostly interested in the case where $\alpha=\beta$, namely, I want to calculate
$$\iota_{\xi}\left(\alpha\wedge\star\alpha\right),$$
which I have reason to believe is simply $2\iota_{\xi}\alpha\wedge\star\alpha$, but this could be wrong.
My attempt so far has been to use the fact that
$$\iota_{\xi}\left(\alpha\wedge\star\beta\right)=\iota_{\xi}\alpha\wedge\star\beta+(-1)^p\alpha\wedge\iota_{\xi}(\star\beta),$$
but then have no idea how to deal with $\iota_{\xi}(\star\beta)$.
Proof: Both sides are linear in $\xi,\omega$. So it suffices to verify it for $\xi=e_1$ and $\omega=e^{i_1\dots i_p}$, $1\leq i_1<i_2<\dots<i_p\leq n$.
You can now applying this to $i_\xi(\alpha\wedge\star\beta)$.
There is a special simplification available with $\alpha=\beta$: $\alpha\wedge\star\alpha=\lvert\alpha\rvert^2\,(\star1)$ so $$ i_\xi(\alpha\wedge\star\alpha)=\lvert\alpha\rvert^2\,i_\xi(\star1)=\lvert\alpha\rvert^2\,\star(\xi^\flat). $$