Weak harmonic map Let $\Omega$ be the unit sphere in two Or three dimensions. Why is then $(\nabla u, \nabla w)=(|\nabla u|^2u,w)$ for all testfunctions w in $C_0^\infty(\Omega, R^n)$? How to compute it?
 A: I think you mean that $\Omega $ is the unit ball $\{x:|x|<1\}$ rather than the unit sphere $\{x:|x|=1\}$. Indeed, the Laplace equation for $u$ takes the form $\Delta u+|\nabla u|^2 u=0$ (with the Euclidean Laplacian $\Delta$) only when the domain of $u$ is flat and the target is a unit sphere $S^k$ of any dimension.  
We can calculate the derivative of $u$ using the product rule and the chain rule: 
$$
\nabla (|x|^{-1}x) = \nabla (|x|^{-1})\otimes x + |x|^{-1}\nabla x = 
(-1)|x|^{-3} x\otimes x + |x|^{-1} I
$$
where $I$ is the identity matrix. In a simpler form, 
$$
\nabla u = \frac{|x|^2I-x\otimes x}{|x|^3} \tag{1}
$$
With a little algebra we get
$$
\begin{align}
|\nabla u|^2 &= \frac{1}{|x|^6} \sum_{i,j=1}^n (|x|^2\delta_{ij}-x_ix_j)^2
\\&=\frac{1}{|x|^6} \sum_{i=1}^n \left\{ (|x|^2-x_i^2)^2 + x_i^2 \sum_{j\ne  i} x_j^2 \right\} \\&=\frac{1}{|x|^6}\sum_{i=1}^n \left\{ (|x|^2-x_i^2)|x|^2\right\} 
=\frac{n-1}{|x|^2}
\end{align}\tag{2}$$
Note that in two dimensions the expression in (2) is not integrable near the origin. Therefore, in this case $u$ does not belong to $H^1$ and is not weakly harmonic. (Aside: weakly harmonic maps in two dimensions are smooth, by a theorem of Hélein.)
In dimensions $n\ge 3$ we have $u\in W^{2,1}(\Omega)$ because the second-order derivatives of $u$ are homogeneous of degree $-2$. Thus, we can show that the identity holds in the strong sense:  $-\Delta u=|\nabla u|^2u$ in $L^1$.  Indeed, for any $i=1,\dots,n$ we have 
$$
\begin{align}
\Delta (|x|^{-1}x_i)&=\operatorname{div}\nabla (|x|^{-1}x_i)=\operatorname{div}\left( |x|^{-1}e_i-x_i|x|^{-3}x\right) \\ &= \frac{\partial}{\partial x_i}|x|^{-1} -\sum_{j=1}^n \frac{\partial}{\partial x_j} (x_i|x|^{-3}x_j) \\ 
&=-\frac{x_i}{|x|^3}-\sum_{j=1}^n x_i\left\{|x|^{-3}-3x_j^2|x|^{-5}\right\}-|x|^{-3}x_i
 \\&= -(n-1)\frac{x_i}{|x|^3}
\end{align}
$$
which matches $|\nabla u|^2u$ component by component.
