Write the functional derivative of the Dirichlet energy I have the following functional:
$J(u) = \int_{\Omega} \frac{|\nabla u|^p}{p}\,d\Omega$
and I want to compute its functional derivative along the direction of an arbitrary test function $v\in H_0^1$.
I tried applying the limit definition 
$J'(u)v = \lim_{\varepsilon \to 0} \frac{J(u+\varepsilon v)-J(u)}{\varepsilon}$
but I have some problems with the non linear part of the integral.
I know that the expected result should be
$J'(u)v=\int_{\Omega}|\nabla u|^{p-2}\nabla u \cdot \nabla v \,d\Omega$ 
since imposing $J'(u)v=0$ is equivalent to solve weakly the homogeneus p-Laplace equation, but I don't know how to proceed computing this derivative.
 A: In general when $f: \Bbb{R}^n \to \Bbb{R}$ is the euclidean norm, defined by
\begin{equation}
f(x) = |x| = \sqrt{\sum_{i=1}^n x_i^2}
\end{equation}
then away from the origin, $f$ is $C^{\infty}$ and its derivative is given by
\begin{equation}
f'(x)(\eta) = \left\langle \dfrac{x}{|x|}, \eta \right\rangle
\end{equation}
(for all $x \neq 0$, and for all $\eta \in \Bbb{R}^n$). 
Now, for convenince, define $g: \Bbb{R} \times \Bbb{R}^n \to \Bbb{R}^n$,
\begin{equation}
g(\varepsilon,x) = \nabla u(x) + \varepsilon \nabla v(x)
\end{equation}
The energy functional's derivative you wish to calculate is given by
\begin{align}
J'(u)(v) &= \dfrac{d}{d \varepsilon} \bigg|_{\varepsilon = 0} J(u+ \varepsilon v) \\
&= \dfrac{d}{d \varepsilon} \bigg|_{\varepsilon = 0} \int_{ x \in\Omega} \dfrac{\left[ (f \circ g)(\varepsilon, x) \right]^p}{p}  \\
&= \int_{ x \in\Omega}  \dfrac{\partial}{\partial \varepsilon} \bigg|_{\varepsilon = 0} \dfrac{\left[ (f \circ g)(\varepsilon, x) \right]^p}{p} ,
\end{align}
This last equality is by Leibniz's integral rule for differentiation under the integral sign. Now, inside, we just have to make use of the chain rule inside. Doing so yields
\begin{align}
J'(u)(v) &= \int_{ x \in\Omega}  \dfrac{p \cdot \left[ f(g(0,x)) \right]^{p-1}}{p} \cdot f'[g(0,x)] \left( \dfrac{\partial g}{\partial \varepsilon}(0,x) \right) \\
&= \int_{ x \in\Omega} |\nabla u(x)|^{p-1} \cdot \left \langle \dfrac{\nabla u(x)}{|\nabla u(x)|} , \nabla v(x) \right \rangle \\
&= \int_{x \in \Omega} |\nabla u(x)|^{p-2} \cdot \left \langle \nabla u(x), \nabla v(x) \right \rangle
\end{align}
Modulo some notation, this is precisely what you wanted to prove.

By the way, my entire answer assumes that all the functions are nice and differentiable, and in particular that if $ p< 2$, then $\nabla u(x)$ doesn't vanish anywhere, so that the norm is differentiable. If this isn't the case, then the argument will necessarily have to be more involved... I'm not too sure how one might proceed in this case though.
If $p\geq 2$, then we can write $|\nabla u|^p = \langle \nabla u, \nabla u\rangle^{p/2}$. In this case, the functions will be differentiable even regardless of whether or not $\nabla u$ vanishes; this is because the function $t \mapsto t^{p/2}$ is everywhere differentiable, and the inner product $\langle \cdot , \cdot \rangle$ is also everywhere differentiable, hence their composition will be as well. Of course, in this case, we will have to slightly modify the above argument, but the final answer is the same
