Let $u:\mathbb{R}^3 \to \mathbb{R}$ a harmionic function, and let $S_r = \{x^2+y^2+z^2 = r^2 \}$
I want to show that if $0 < r_1 <r_2$, then:
$$ \frac{1}{\mathrm{area}(S_{r_1})} \int_{S_{r_1}} u^2 dS \le \frac{1}{\mathrm{area}(S_{r_2})} \int_{S_{r_2}} u^2 dS$$
I know that: $u^2(0) \le \frac{1}{\mathrm{area}(S_r)} \int_{S_r} u^2 dS$
So I can instead try to show $$\frac{1}{\mathrm{area}(S_{r_1})} \int_{S_{r_1}} u^2 dS \le u^2(0)$$
But I don't see how to do it.
Help would be appreciated.
We can use this to show the $n-$dimensional case.
Set $$\phi(r):=\frac1{\text{area}(S_r)}\int_{S_r}u^2(y)\,dS(y)=\frac1{\text{area}(S_1)}\int_{S_1}u^2(rz)\,dS(z).$$ Then $$\phi'(r)=\frac1{\text{area}(S_1)}\int_{S_1}2u(rz)Du(rz)\cdot z\,dS(z),$$ and consequently, using Green's formula, we compute \begin{align*} \phi'(r)&=2\frac1{\text{area}(S_r)}\int_{S_r}u(y)Du(y)\cdot\frac yr\,dS(y)\\&=2\frac1{\text{area}(S_r)}\int_{S_r}u\frac{\partial u}{\partial \nu}\,dS\\&=2\frac1{\text{area}(S_r)}\left(\int_{B_r}Du\cdot Du\,dx+\int_{B_r}u\Delta u\,dx\right)\\&=2\frac1{\text{area}(S_r)}\int_{B_r}|Du|^2\,dx\geq 0. \end{align*} Hence $\phi(r)$ is increasing.
Addendum: From above we can get $$\frac1{\text{area}(S_r)}\int_{S_r}u^2\,dS=\phi(r)\geq \lim_{r\to 0}\phi(r)=u^2(0),$$ since $\phi$ is continuous at $r=0$, which can be proven using the method in a previous post.