Sharp global estimates for local and nonlocal porous medium-type equations in bounded domains

This paper provides a quantitative study of nonnegative solutions to nonlinear diffusion equations of porous medium-type of the form $\partial_t u + {\mathcal L}u^m=0$,$m>1$,wheretheoperator${\mathcal L}$ belongs to a general class of linear operators, and the equation is posed in a bounded domain $\mathrm{\Omega }\subset {\mathbb{R}}^N$. As possible operators we include the three most common definitions of the fractional Laplacian in a bounded domain with zero Dirichlet conditions, and also a number of other nonlocal versions. In particular, $\mathcal{L}$ can be a power of a uniformly elliptic operator with $C^1$ coefficients. Since the nonlinearity is given by $u^m$ with $m>1$, the equation is degenerate parabolic. The basic well-posedness theory for this class of equations has been recently developed in [14,15]. Here we address the regularity theory: decay and positivity, boundary behavior, Harnack inequalities, interior and boundary regularity, and asymptotic behavior. All this is done in a quantitative way, based on sharp a priori estimates. Although our focus is on the fractional models, our results cover also the local case when $\mathcal{L}$ is a uniformly elliptic operator, and provide new estimates even in this setting. A surprising aspect discovered in this paper is the possible presence of non-matching powers for the long-time boundary behavior. More precisely, when $\mathcal{L}=(-\mathrm{\Delta }{)}^s$ is a spectral power of the {Dirichlet} Laplacian inside a smooth domain, we can prove that:

• when $2s>1-1\mathrm{/}m$, for large times all solutions behave as ${\rm dist}^{1/m}$ near the boundary;
• when $2s\le 1-1\mathrm{/}m$, different solutions may exhibit different boundary behavior. This unexpected phenomenon is a completely new feature of the nonlocal nonlinear structure of this model, and it is not present in the semilinear elliptic equation $\mathcal{L}{u}^m=u$.