I asked a similar question before but didn't get any responses. So I will attempt again (the prior question was regarding Holder continuity).
Let $ \Omega$ denote a cube in $ R^n$ and consider solving $$- \Delta u(x) + u(x) = f(x) $$ in $ \Omega$ with $ \partial_\nu u=0$ on $ \partial \Omega$.
Suppose $ f \in L^p(\Omega)$ where $ 1 <p< \infty$. Then does one have that $ u \in W^{2,p}(\Omega)$ along with the desired estimate?
I am aware of the book by Grisvard "Elliptic_Problems_in_Nonsmooth_Domains".
Since I am on a cube it seems one can just reflect and prove the above result ??
Any comments on whether one even thinks this result is true would be greatly appreciated. Also what about if $ f $ is Holder continuous of order $ 0 < \alpha<1$. Can we get the usual $ C^{2,\alpha}$ regularity? (again i managed to "prove" this to myself but i not convinced).
thanks