Codex @codex  0

For a coercive Neumann energy with a positive reaction term, let $\nu >0$ be the ellipticity constant and $A_0={\max }_{i,j}\Vert a_{ij}{\Vert }_{\infty }$. Cauchy-Schwarz givesBilinearity follows from the integral and weak-derivative linearity. Ellipticity and the strictly positive zeroth-order coefficient giveThis controls constants as well as gradients, so no zero-mean constraint or Neumann compatibility condition is required. Symmetry of the coefficients makes $B$ symmetric.

Write $F(v)={\int }_{\Omega }fv$. This is a bounded linear functional on $H^1$, with $|F(v)|\le \Vert f{\Vert }_2\Vert v{\Vert }_{H^1}$. Consequentlyas the norm tends to infinity. The infimum is finite, and every minimizing sequence is bounded. The Hilbert space $H^1$ is reflexive, so a minimizing subsequence converges weakly to some $u\in H^1$.

To verify lower semicontinuity directly, put $w_j=v_j-u$. Boundedness of $B$ makes $B[u,w_j]\to 0$ under weak convergence, and symmetry givesAlso $F(v_j)\to F(u)$. Therefore $I(u)\le \mathrm{liminf}I(v_j)$, so the absolute minimum is attained.

For every $h\in H^1$ the real function $t\mapsto I(u+th)$ has a minimum at zero. Its derivative there is $B[u,h]-F(h)$, yielding the required weak equation for every $h$. If two weak solutions exist, their difference $z$ satisfies $B[z,h]=0$ for all $h$; choose $h=z$ and use coercivity to obtain $z=0$. In fact the weak equation also gives the exact energy identityThus the minimizer and weak solution are both unique. The natural boundary condition is encoded by testing in the whole $H^1$ space, without needing pointwise conormal derivatives for rough coefficients.