Let $(M, g)$ be a Riemannian manifold and $X, Y$ vector fields on a compact subset $D$ of $M$ such that $Y$ is divergence free. I would like to show that $$\int_D g(\nabla_Y X, X) dV = 0 \tag{1}$$ where $\nabla$ is the Levi-Civita connection on $M$ and $dV$ the Riemannian volume form.
In order to show this I am thinking of somehow integrating by parts. I know on a Riemannian manifold there exists the formula $$\int_Mg(\text{grad}(f), X) dV = -\int_Mf \text{div}(X) dV + \int_{\partial M} fg(X,N) d\tilde{V}$$ for any $f \in C^\infty(M)$ and where $N$ is the outward-pointing unit normal vector field along $\partial M$ and $\tilde{g}$ is the induced Riemannian metric on $\partial M$. Is there a similar formula where instead of the gradient we have the covariant derivative as in (1)?
Alternatively, I know that for any $f \in C^\infty(D)$ we have $$\int_D L_Yf dV = 0$$ where $L_Y$ is the Lie derivative taken along $Y$. Is the Lie derivative here in anyway related to the integrand in (1)?