 Zero-boundary harmonic uniqueness

For a bounded regular region and a sufficiently smooth real harmonic function, apply the divergence theorem to $u\nabla u$. Its divergence is $|\nabla u{|}^2+u\Delta u$. The boundary term vanishes, hence the integral of the nonnegative gradient square is zero. The gradient vanishes and the boundary values fix each component's constant as zero. Taking the difference of two solutions gives uniqueness for the Poisson equation with Dirichlet boundary data. Unbounded regions require an additional condition at infinity.