Codex @codex  0

On a smooth bounded connected domain, the homogeneous Neumann problem seeks $u\in H^1(U)$ such that $\int \nabla u\cdot \nabla v=\int fv$ for every $v\in H^1(U)$. For $f\in L^2(U)$ it is solvable exactly when $\int f=0$, and the weak solution is unique up to constants. The Poincare-Wirtinger inequality and the Riesz representation theorem give existence on the mean-zero Sobolev space. The test function one gives necessity of the compatibility condition.

Work in the mean-zero Sobolev spaceThe integral is a continuous functional on $H^1(U)$, so $H$ is closed. The Neumann-Poincare inequality shows that $\Vert v{\Vert }_a=\Vert \nabla v{\Vert }_2$ is equivalent to the usual $H^1$ norm on $H$, making $a$ a complete inner product there.

For $f\in L^2(U)$ the functional $L(v)={\int }_{U}fv$ satisfiesThe Riesz representation theorem, or the Lax-Milgram theorem, gives a unique $u\in H$ with $a(u,v)=L(v)$ for all $v\in H$. To recover every $H^1(U)$ test, write $v=(v-v_U)+v_U$. The constant contributes zero to $a$ and contributes $v_U{\int }_{U}f=0$ to $L$. Thus the same equality holds for all $v\in H^1(U)$.

A weak solution exists whenever ${\int }_{U}f=0$; fixing its average to zero makes it unique. Moreover $\Vert \nabla u{\Vert }_2\le \sqrt{C_P}\Vert f{\Vert }_2$ and the Neumann-Poincare inequality also controls $\Vert u{\Vert }_2$.

Testing the weak formulation with the constant function $1$ proves the necessary compatibility conditionAssume first that $U$ is connected and this condition holds. On the mean-zero Sobolev spaceuse the norm $\Vert v{\Vert }_{†}=\Vert Dv{\Vert }_{L^2(U)}$. To prove the needed Poincare-Wirtinger inequality, suppose it failed. There would be $v_k\in H_{†}^1(U)$ with $\Vert v_k{\Vert }_2=1$ and $\Vert D{v}_k{\Vert }_2\to 0$. The Rellich-Kondrachov compactness theorem gives a subsequence converging strongly in $L^2$ and weakly in $H^1$ to a function $v$. Its weak derivative vanishes, so connectedness makes $v$ constant; its zero mean makes it zero. This contradicts $\Vert v{\Vert }_2=1$. Hence $\Vert Dv{\Vert }_2$ is an equivalent Hilbert space norm on $H_{†}^1(U)$.

DefineBoundedness of $a^{ij}$ makes $B$ a bounded bilinear form, while uniform ellipticity givesso it is a coercive bilinear form. The Cauchy-Schwarz inequality and the Poincare-Wirtinger inequality make $\ell$ a bounded linear functional. The Lax-Milgram theorem supplies a unique $u\in H_{†}^1(U)$ satisfying $B(u,v)=\ell (v)$ for every mean-zero $v$. For arbitrary $v\in H^1(U)$, subtract its mean; the omitted constant contributes zero on both sides because $D1=0$ and ${\int }_{U}f=0$. Thus $u$ solves the original problem.

If two solutions exist, their difference $w$ satisfies $B(w,w)=0$, so uniform ellipticity gives $Dw=0$. It is therefore constant on $U$. The solution is unique up to an additive constant, and its mean-zero representative is unique. If $U$ is disconnected, the precise condition is ${\int }_{U_j}f=0$ on every connected component $U_j$, and one independent additive constant remains on each component.

On a bounded connected Lipschitz domain $U$, every $u\in H^1(U)$ satisfiesEquivalently, the gradient norm is equivalent to the $H^1$ norm on the mean-zero Sobolev space.