My PDE professor showed the following result:
Let $f \in L^{p}(\Omega)$, for $1 < p < \infty$. Also consider $u \in L^{1}_{loc}(\Omega)$ a solution of
\begin{align} -\Delta u + a(x)u&= f, \quad x\in \Omega\\ \hspace{2.7cm}u &= 0, \quad x\in \partial \Omega . \end{align} Then $u \in W^{2,p}(\Omega).$
In Evans, page 332, Theorem 2, is a result similar result, however for $f \in H^{m}(\Omega)$. But I really need this for $f \in L^{p}(\Omega)$. I've looked in several books for this result but haven't found it. Anyone know where I can find this result, or something similar?
I'm no expert on elliptic PDEs, but I think I found two books:
Wu, Zhuoqun; Yin, Jingxue; Wang, Chunpeng, Elliptic and parabolic equations, Hackensack, NJ: World Scientific (ISBN 981-270-025-0/hbk; 981-270-026-9/pbk). xv, 408 p. (2006). ZBL1108.35001.
It seems Gilbarg and Trudinger also covers this in chapter 9 (Strong solutions).
For the first book, the introduction for Chapter 9 "$L^{p}$ Estimates for Linear Equations and Existence of Strong Solutions" reads:
Here, for $$-a_{i j}(x) D_{i j} u+b_{i}(x) D_{i} u+c(x) u=f(x), \quad x \in \Omega\tag{9.1.7}$$ $$a_{i j}(x) \xi_{i} \xi_{j} \geq \lambda|\xi|^{2}, \quad \forall \xi \in \mathbb{R}^{n}, x \in \Omega\tag{9.1.8}$$ $$\sum_{i, j=1}^{n}\left\|a_{i j}\right\|_{L^{\infty}(\Omega)}+\sum_{i=1}^{n}\left\|b_{i}\right\|_{L^{\infty}(\Omega)}+\|c\|_{L^{\infty}(\Omega)} \leq M .\tag{9.1.9}$$ they prove
This result is used to prove for your equation ((9.1.15) in the book) the following:
Since the equation is in divergence form they also proceed with energy estimates like in the $L^2$ case (testing with $|u|^{p-2}u$).
For $p<2$ they have
But they don't elaborate. (Perhaps it is clear why, but I don't know and I am only glancing through the book) Hope this helps.