Stephan Dahlke, Lars Diening, Christoph Hartmann et al.
In this paper, we study the regularity of solutions to the $p$-Poisson equation for all $1<p<\infty$. In particular, we are interested in smoothness estimates in the adaptivity scale $ B^σ_τ(L_τ(Ω))$, $1/τ= σ/d+1/p$, of Besov spaces. The regularity in this scale determines the order of approximation that can be achieved by adaptive and other nonlinear approximation methods. It turns out that, especially for solutions to $p$-Poisson equations with homogeneous Dirichlet boundary conditions on bounded polygonal domains, the Besov regularity is significantly higher than the Sobolev regularity which justifies the use of adaptive algorithms. This type of results is obtained by combining local Hölder with global Sobolev estimates. In particular, we prove that intersections of locally weighted Hölder spaces and Sobolev spaces can be continuously embedded into the specific scale of Besov spaces we are interested in. The proof of this embedding result is based on wavelet characterizations of Besov spaces.