Let $\Omega$ be a smooth bounded domain in Euclidean space. Let $u \in C_c^1(\Omega)$. The subscript indicates compact support. Let $1 < p < \infty$.
Can a $v \in C_c^1(\Omega)$ (or in its closure in the Sobolev space $W^{1,q}(\Omega)$, where $1/p + 1/q = 1$) be found such that $\nabla u \cdot \nabla v = |\nabla u|^p$ (a.e.) in $\Omega$?
Or merely $\int_\Omega \nabla u \cdot \nabla v = \int_\Omega |\nabla u|^p$?
Paraphrasing V. Barbu, Nonlinear differential equations of monotone types in Banach spaces, 2010, Springer, p.4:
Let $1 < p < \infty$ and $\Omega$ a bounded and open subset of $\mathbb{R}^n$. Write $X := W_0^{1,p}(\Omega)$. For $u \in X$ define $J(u) \in X'$ by $$J(u) = -\|u\|_X^{2-p} \sum_{i=1}^n \partial_i (|\partial_i u|^{p-2} \partial_i u).$$ Then $J : X \to X'$ is the duality mapping (it is indeed single-value by uniform convexity of $X$), for which $$\langle J(u), u \rangle = \|u\|_X^2 = \|J(u)\|_{X'}^2.$$ This (almost) answers the second part of the question.
As an aside, the duality mapping is, in general, a set-valued mapping $J : X \to 2^{X'}$ on a Banach space defined by $$J(x) = \{ f \in X' : f(x) = \|x\|_X^2 = \|f\|_{X'}^2 \},$$ and this is always non-empty by the Hahn-Banach theorem (analytic version).
End of paraphrase.
Let $u \in C_c^1(\Omega)$ be nonzero. Then $J(u)$ is a bounded linear functional on $H_0^1(\Omega)$. Let $v \in H_0^1(\Omega)$ solve $\int_\Omega \nabla v \cdot \nabla \varphi = \langle J(u), \varphi \rangle$ for all $\varphi \in C_c^1(\Omega)$. Then $\|u\|_X^2 = \langle J(u), u \rangle = \int_\Omega \nabla v \cdot \nabla u$.