Continuity of the spatial gradient of weak solutions to very singular parabolic equations involving the one-Laplacian
Shuntaro Tsubouchi
TL;DR
The paper resolves the longstanding question of spatial gradient continuity for weak solutions to a very singular parabolic equation combining a one-Laplacian and a $p$-Laplace operator. By introducing a truncation of the gradient and a parabolic approximation $E_{\varepsilon}$, the authors obtain local Hölder continuity for the truncated gradient ${\mathcal G}_{2\delta,\varepsilon}(\nabla u_{\varepsilon})$ and derive uniform estimates via De Giorgi truncation, heat-flow comparisons, and Campanato-type growth. The convergence analysis then transfers these regularity properties to the original solution, yielding $\nabla u\in C^{0}(\Omega_T)$ under the stated $p$-range and forcing conditions, with the approach avoiding intrinsic parabolic scaling. The methods extend prior stationary results to the time-dependent setting and have implications for models such as Bingham flows and crystal-surface evaporation, where facets induce nonuniform parabolicity. Overall, the work provides a robust gradient-regularity framework for highly singular parabolic equations, combining truncation techniques with classical parabolic theory to achieve sharp continuity results.
Abstract
We consider weak solutions to very singular parabolic equations involving a one-Laplace-type operator, which is singular and degenerate, and a $p$-Laplace-type operator with $\frac{2n}{n+2}<p<\infty$, where $n\ge 2$ denotes the space dimension. This type of equation is used to describe the motion of a Bingham flow. It has been a long-standing open problem of whether the spatial gradients of weak solutions are continuous in space and time. This paper aims to give an affirmative answer for a wide class of such equations. This equation becomes no longer uniformly parabolic near the facet, the place where the spatial gradient vanishes. To achieve our goal, we show local a priori Hölder continuity of gradients suitably truncated near facets. For this purpose, we consider a parabolic approximate problem and appeal to standard methods, including De Giorgi's truncation and comparisons with Dirichlet heat flows. Our method is a parabolic adjustment of our method developed to prove the corresponding statements for stationary problems.
