 Solution

We prove the Interior second-derivative estimate for the Poisson equation using difference quotients, rather than differentiating the $L^2$ forcing. Fix a coordinate $e_k$ and write ${\delta }_{h}v(x)=(v(x+h{e}_k)-v(x))/h$. Choose a smooth cutoff function $\eta$ equal to one on ${\Omega }^{\prime }$ with support compactly contained in $\Omega$, leaving room for all translations with sufficiently small $|h|$. We can arrange $|D\eta |\le C(n)/d$, where $d=\mathrm{dist}({\Omega }^{\prime },\partial \Omega )$.

The standard difference quotient lemmas give discrete integration by parts, $\int a{\delta }_{-h}b=-\int ({\delta }_{h}a)b$, and $\Vert {\delta }_{h}v{\Vert }_2\le \Vert D_{k}v{\Vert }_2$ on the relevant enlarged set. In the weak formulation $\int Du\cdot D\phi =-\int f\phi$, use the admissible test function $\phi =-{\delta }_{-h}({\eta }^2{\delta }_{h}u)$. Moving a difference quotient onto $Du$ givesSet $X=\Vert \eta D{\delta }_{h}u{\Vert }_2$ and $Y=\Vert D\eta {\delta }_{h}u{\Vert }_2$. The cross term has absolute value at most $2XY$. Since $0\le \eta \le 1$, the difference quotient lemma on the right givesThe Young inequality now implies $X^2\le C(Y^2+\Vert f{\Vert }_2^2)$. Moreover $Y\le C(n)d^{-1}\Vert Du{\Vert }_{L^2(\Omega )}$, independently of $h$. Hence, for every coordinate $k$,By the Sobolev characterization by bounded difference quotients, applied to each first weak derivative on smaller interior sets, these uniform bounds give all second weak derivatives. Their bounds pass to the weak limit as $h\to 0$. Thus $u\in W_{\mathrm{loc}}^{2,2}(\Omega )$ andNo derivative of $f$ is used, and no boundary condition is needed for this interior elliptic regularity argument.