Supercaloric functions for the parabolic p-Laplace equation in the fast diffusion case

We study a generalized class of supersolutions, so-called p-supercaloric functions, to the parabolic p-Laplace equation. This class of functions is defined as lower semicontinuous functions that are finite in a dense set and satisfy the parabolic comparison principle. Their properties are relatively well understood for p≥2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p\ge 2$$\end{document}, but little is known in the fast diffusion case 1<p<2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1<p<2$$\end{document}. Every bounded p-supercaloric function belongs to the natural Sobolev space and is a weak supersolution to the parabolic p-Laplace equation for the entire range 1<p<∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1<p<\infty $$\end{document}. Our main result shows that unbounded p-supercaloric functions are divided into two mutually exclusive classes with sharp local integrability estimates for the function and its weak gradient in the supercritical case 2nn+1<p<2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{2n}{n+1}<p<2$$\end{document}. The Barenblatt solution and the infinite point source solution show that both alternatives occur. Barenblatt solutions do not exist in the subcritical case 1<p≤2nn+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1<p\le \frac{2n}{n+1}$$\end{document} and the theory is not yet well understood.


Introduction
This paper studies classes of supersolutions to the parabolic p-Laplace equation ∂ t u − div |∇u| p−2 ∇u = 0.
(1.1) with respect to the solutions of (1.1), see Definition 2.6 below. The definition of supercaloric functions is the same as in classical potential theory for the heat equation when p = 2, see Watson [23]. By Juutinen et al. [13], the class of p-supercaloric functions is the same as the viscosity supersolutions to (1.1) for 1 < p < ∞. Our results can be extended to more general quasilinear equations with the p-growth, which are discussed in DiBenedetto [8], DiBenedetto et al. [9] and Wu et al. [24]. For simplicity, we discuss only the prototype case in (1.1).
A p-supercaloric function does not, in general, belong to the natural Sobolev space for (1.1). The only connection to the equation is through the comparison principle. However, Kinnunen and Lindqvist [15] proved that bounded p-supercaloric functions belong to the appropriate Sobolev space and are weak supersolutions to (1.1) for p ≥ 2. Korte et al. [17] extended the study for a more general class of parabolic equations with p-growth. In this paper we show that bounded p-supercaloric functions are weak solutions to (1.1) for the entire range 1 < p < ∞.
We are mainly interested in unbounded p-supercaloric functions. Assume that u is a p-supercaloric function in Ω T = Ω × (0, T ), where Ω is an open set in R n and T > 0. One of the main results of Kuusi et al. [22] asserts that for p > 2 there are two mutually exclusive alternatives: Either u ∈ L q loc (Ω T ) for every 0 < q < p − 1 + p n or u / ∈ L p−2 loc (Ω T ). In particular, if u ∈ L p−2 loc (Ω T ), then u ∈ L q loc (Ω T ) for every 0 < q < p − 1 + p n . For the corresponding theory for the porous medium equation, see Kinnunen et al. [19].
Examples based on the Barenblatt solution (see Barenblatt [3]) and the friendly giant show that both alternatives occur. In the first alternative the upper bound for the exponent is given by the Barenblatt  , (x, t) ∈ R n × (0, ∞), (1.2) where 2 < p < ∞, λ = n(p − 2) + p and the constant c is a positive number, which can be chosen such that is p-supercaloric in R n × R with u ∈ L q loc (R n × R) for every 0 < q < p − 1 + p n . This function solves (1.1) with a finite point source, but it fails to belong to the natural Sobolev space, since |∇u| / ∈ L p loc (R n × R). For the second alternative, we consider a bounded open set Ω in R n with a smooth boundary. By separation of variables, we obtain the friendly giant The function U is a weak solution to (1.1) in Ω × (0, ∞) and the zero extension as in We prove the corresponding result in the supercritical range 2n n+1 < p < 2 for a p-supercaloric function u in Ω T . It asserts that either u ∈ L q loc (Ω T ) for  (1.5) where λ = n(p − 2) + p and the constant c is a positive number so that (1.3) holds for every t > 0. Observe that p > 2n n+1 is equivalent with λ > 0. Barenblatt solutions do not exist in the subcritical case 1 < p ≤ 2n n+1 and, to our knowledge, the theory is not yet well understood. The Barenblatt solution is compactly supported for every t > 0 for 2 < p < ∞, but it is positive everywhere for 2n n+1 < p < 2. As above, the zero extension to the negative times is p-supercaloric in R n × R with u ∈ L q loc (R n × R) for every 0 < q < p − 1 + p n . A prime example of a p-supercaloric function for 2n n+1 < p < 2 that does not belong to the Barenblatt class is the infinite point source solution

Weak supersolutions and supercaloric functions
We begin with notation. Let Ω ⊂ R n be an open set. For T > 0 we denote a space-time cylinders in R n+1 by Ω T = Ω× (0, T ) and Ω t1,t2 = Ω× (t 1 , t 2 ), with We denote a cube in R n by Q = (a 1 , b 1 ) × · · · × (a n , b n ) and the space-time cylinders of the form Q t1,t2 with t 1 < t 2 , are called boxes. Ω t1,t2 Ω T denotes that Ω t1,t2 is a compact subset of Ω T . Let W 1,p (Ω) denote the Sobolev space of functions u ∈ L p (Ω), whose first distributional partial derivatives ∂u ∂xi , i = 1, 2, . . . , n, exist in Ω and belong to L p (Ω). The corresponding Sobolev space with zero boundary values is denoted by W 1,p 0 (Ω). The parabolic Sobolev space L p (0, T ; W 1,p (Ω)) consists of functions u = u(x, t) such that for almost every t ∈ (0, T ), the function x → u(x, t) belongs to W 1,p (Ω) and Observe that the time derivative ∂ t u does not appear anywhere. The definition of the space Ω T . The definition for L p loc (R; W 1,p loc (R n )) is analogous. Gradient and divergence are always taken with respect to the spatial variable only.
Weak solutions to (1.1) are assumed to belong to a parabolic Sobolev space, which guarantees a priori local integrability for the function and its weak gradient. We state the definition and results in a space-time cylinder Ω T , with T > 0, but extensions to arbitrary cylinders Ω t1,t2 , with t 1 < t 2 , are obvious. Definition 2.1. Let 1 < p < ∞ and let Ω be an open set in R n . A function u ∈ L p loc (0, T ; W 1,p loc (Ω)) is called a weak solution to (1.1), if ΩT −u∂ t ϕ + |∇u| p−2 ∇u · ∇ϕ dx dt = 0 for every ϕ ∈ C ∞ 0 (Ω T ). Furthermore, we say that u is a weak supersolution if the integral above is nonnegative for all nonnegative test functions ϕ ∈ C ∞ 0 (Ω T ). If the integral is non-positive for such test functions, we call u a weak subsolution.
We consider a comparison principle for weak super-and subsolutions. The argument is similar to Kilpeläinen A pointwise limit of a sequence of uniformly bounded weak supersolutions is a weak supersolution, see Korte et al. [17,Theorem 5.3], see also Kinnunen and Lindqvist [15].
The class of weak supersolutions is closed under taking minimum of two functions, see Korte el al. [17,Lemma 3.2]. See also Kinnunen and Lindqvist [15]. The local boundedness assumption in Theorem 2.4 can be replaced with a uniform Sobolev space bound, see Kinnunen and Lindqvist [15] and Korte et al. [17,Remark 5.6]. This implies that if u ∈ L p loc (0, T ; W 1,p loc (Ω)) and min{u, k} ∈ L p loc (0, T ; W 1,p loc (Ω)) is a weak supersolution for every k = 1, 2, . . . , then u is a weak supersolution. In general, we have to consider more general class of solutions than weak supersolutions. Thus we define p-supercaloric functions as in Kilpeläinen and Lindqvist [14, Lemma 3.1], see also Kinnunen and Lindqvist [15].
Similarly, an upper semicontinuous function u is said to be a p-subcaloric function in Ω T if −u is p-supercaloric function in Ω T . Since p-supercaloric functions are lower semicontinuous, they are locally bounded from below. Thus by adding a constant, we may assume that a psupercaloric function is nonnegative, when we discuss local properties. If u is a nonnegative p-supercaloric function in Ω T , the zero extension in the past as in (1.4) is a p-supercaloric function in Ω × (−∞, T ). The following assertion is a direct consequence of the comparison principle in Definition 2.6.  Assume that u is a p-subcaloric function and v is a p-supercaloric function in We will apply an existence result for the obstacle problem. In order to guarantee continuity of the solution up to the boundary, we assume that the domain where the obstacle problem is considered has Lipschitz boundary. This condition can be relaxed, see [18,Theorem 3.1].
Theorem 2.11. Let 1 < p < ∞. Assume that Ω is an open and bounded set in R n with Lipschitz boundary and let ψ ∈ C(Ω T ). There exists u ∈ C(Ω T ) that is a weak supersolution to (1.1) in Ω T with the following properties: Next we show that every bounded p-supercaloric function u is a weak supersolution in Ω T . In this case u ∈ L p loc (0, T ; W 1,p loc (Ω)) and  Proof. Since u is lower semicontinuous, there exists a sequence of continuous functions ψ i , i = 1, 2, . . . , such that ψ 1 ≤ ψ 2 ≤ · · · ≤ u and for every (x, t) ∈ Ω T . We claim that u is a weak supersolution in Ω T . To conclude this it is sufficient to show that u is a weak supersolution in every space-time box Q t1,t2 Ω T . By Theorem 2.11, for every i = 1, 2, . . . , there exists a solution v i ∈ C(Q t1,t2 ) to the obstacle problem in Q t1,t2 with the obstacle ψ i . The function v i is a continuous weak solution in the open set ,t2 ), it follows that v i = ψ i on the boundary ∂U , except possibly when t = t 2 . By Corollary 2.10, we conclude that u ≥ v i in U for every i = 1, 2, . . . . Consequently, ψ i ≤ v i ≤ u in Q t1,t2 for every i = 1, 2, . . . and thus for every (x, t) ∈ Q t1,t2 . According to the Theorem 2.4, the function u is a weak supersolution in Q t1,t2 . Since this holds true for every Q t1,t2 Ω T we conclude that u is a weak supersolution in Ω T .

Barenblatt solutions
This section discusses p-supercaloric functions with a Barenblatt type behaviour in the supercritical range.  × (0, ∞). It follows that the zero extension u(x, t) as in (1.4) is a p-supercaloric function in R n × R according to Definition 2.6 and it satisfies R n ×R −u∂ t ϕ + |∇u| p−2 ∇u · ∇ϕ dx dt = Mϕ(0), with M > 0, for every ϕ ∈ C ∞ 0 (Ω T ). This means that u solves the equation ∂ t u − div(|∇u| p−2 ∇u) = Mδ in the weak sense, where δ is Dirac's delta. It follows that the weak gradient ∇u(·, t) exists, in the sense of (3.2), for almost every t and it is locally integrable to any power 0 < q < p − 1 + 1 n+1 . However, the function u is not a weak supersolution to (1.1) for every r > 0, t 1 ≤ 0 and t 2 > 0. This implies that u / ∈ L p loc (R; W 1,p loc (R n )). Observe, that the truncations min{u, k}, k = 1, 2, . . . belong to L p loc (R; W 1,p loc (R n )) and are weak solutions to (1.1) in R n+1 by Theorem 2.12.
Next we recall a Caccioppoli type inequality for nonnegative unbounded weak supersolutions, see Kuusi [20, Lemma 2.2]. Lemma 3.2. Let 1 < p < ∞, 0 < ε < 1 and let Ω be an open set in R n . Assume that u is a nonnegative weak supersolution in Ω T . There exists a constant c = c(p, ε) such that ΩT |∇u| p u −ε−1 ϕ p dx dt + ess sup 0<t<T for every nonnegative test function ϕ ∈ C ∞ 0 (Ω T ). The following version of Sobolev's inequality will be useful for us, see We prove a general local integrability result for unbounded p-supercaloric functions. The obtained integrability exponent is sharp as shown by the Barenblatt solution. We apply a Moser type iteration scheme.

Theorem 3.4. Let 2n
n+1 < p < 2 and let Ω be an open set in R n . Assume that u is a p-supercaloric function in Ω T . If u ∈ L s loc (Ω T ) for some s > n p (2 − p), then u ∈ L q loc (Ω T ) whenever 0 < q < p − 1 + p n . Proof. Since u is locally bounded from below, by adding a constant, we may assume that u ≥ 1. We may assume 0 < s < 1, otherwise we have u ∈ L 1 loc (Ω T ) and we can directly proceed to the last step of the proof.
for every k = 2, 3, . . . . By combining the estimates and applying Lebesgue's monotone convergence theorem we conclude Denote s 0 = s and s 1 = s 0 (1 + p n ) − (2 − p). After the first step of iteration we have u ∈ L s1 loc (Ω T ). If s 1 < 1, then in the second step of iteration we choose 1 − ε = s 1 and again combine the Caccioppoli inequality and the Sobolev inequality with m = ps1 s1− (2−p) . We continue in this way and obtain an increasing sequence of numbers s i , satisfying which can be written in terms of s 0 as After the ith step of iteration we have u ∈ L si loc (Ω T ). After a finite number of iterations we have s i ≥ 1 and thus u ∈ L 1 loc (Ω T ). In order to pass from u ∈ L 1 loc (Ω T ) to u ∈ L p−1+ p n −σ loc (Ω T ) for every σ > 0 we apply a similar argument once more. Let ε = σ By Lemma 3.2, with u k ≥ 1, this implies for every k = 2, 3, . . . . For the second term, Lemma 3.2 implies Let u be a nonnegative p-supercaloric function in Ω T . It may happen that u does not have a locally integrable weak derivative. In this case we consider the truncations u k = min{u, k} ∈ L p loc (Ω T ), k = 1, 2, . . . and define the weak gradient by which is a well defined measurable function but does not necessarily belong to L 1 loc (Ω T ). If ∇u ∈ L 1 loc (Ω T ), then ∇u(·, t) is the Sobolev gradient of u(·, t) for almost every t with 0 < t < T . Theorem 3.5. Let 2n n+1 < p < 2 and let Ω be an open set in R n . Assume that u is a p-supercaloric function in Ω T with u ∈ L s loc (Ω T ) for some s > n p (2 − p), Then ∇u ∈ L q loc (Ω T ) whenever 0 < q < p − 1 + 1 n+1 . Proof. Since u is locally bounded from below, by adding a constant, we may assume that u ≥ 1. Let 0 < t 1 < t 2 < T , Ω Ω and ε ∈ (0, 1). We consider the truncations u k = min{u, k}, k = 1, 2, . . . . By Hölder's inequality, we have where the first integral is uniformly bounded with respect to k as in (3.1) and the second integral is finite whenever q 1+ε p−q < p − 1 + p n by Theorem 3.4. From this we can conclude that the right hand side is finite for any 0 < q < p − 1 + 1 n+1 . Now we have that ∇u k is uniformly bounded in L q loc (Ω T ). Remark 3.6. If p > 2n+1 n+1 , then p − 1 + 1 n+1 > 1 and it follows that ∇u(·, t) in Lemma 3.5 is the Sobolev gradient for almost every t with 0 < t < T . However, if 2n n+1 < p ≤ 2n+1 n+1 , then p − 1 + 1 n+1 ≤ 1. In this case the Sobolev gradient may not exist and the weak gradient ∇u is interpreted as in (3.2). Remark 3.7. Let 2n n+1 < p < 2 and let Ω be an open set in R n . Assume that u is a p-supercaloric function in Ω T with u ∈ L s loc (Ω T ) for some s > n p (2 − p). By Theorem 3.5 we have ∇u ∈ L p−1 loc (Ω T ). Theorem 2.12 implies for every nonnegative ϕ ∈ C ∞ 0 (Ω T ). By the Riesz representation theorem there exists a nonnegative Radon measure μ on R n+1 such that ΩT −u∂ t ϕ + |∇u| p−2 ∇u · ∇ϕ dx dt = ΩT ϕ dμ for every ϕ ∈ C ∞ 0 (Ω T ). This means that u is a solution to the measure data problem Note that ∇u / ∈ L p loc (Ω T ), in general. The integrability of the gradient of a weak solution to a measure data problem in the fast diffusion case has been studied Baroni [4].
Remark 3.8. Let 2n n+1 < p < 2 and let Ω be an open set in R n . Assume that u is a nonnegative p-supercaloric function in Ω T with u ∈ L 2−p loc (Ω T ). Then ∇ log u ∈ L p loc (Ω T ), where ∇ stands for the usual Sobolev gradient. This property holds for u satisfying the assumptions in Theorem 3.4, since 2 − p < n p (2 − p). The logarithmic estimate can be obtained as in [16] by taking ε = p − 1 in Lemma 3.2. Remark 3.9. By the following dichotomy it follows that p-supercaloric function is either locally integrable to any power smaller than p − 1 + p n > 1, or it is not integrable to a power n p (2 − p) < 1. Observe that as p 2n n+1 we have n p (2 − p) 1 and p − 1 + p n 1. The gap between the exponents becomes smaller and shrinks to a point as p approaches the critical value from above.
We have the following characterization for Barenblatt type p-supercaloric functions.

Infinite point source solutions
This section discusses p-supercaloric functions in the supercritical case that are not covered by Theorem 3.10. The leading example is the infinite point source solution in (1.6). In contrast with the friendly giant in the slow diffusion case p > 2, it is the singularity in space, not in time, that fails to be locally integrable to an appropriate power. Roughly speaking, at every time slice the singularity of the infinite point source solution is worse power type singularity than based on the fundamental solution to the elliptic p-Laplace equation. for every r > 0, t 1 > −∞ and t 2 > 0. This implies that u / ∈ L p loc (R; W 1,p loc (R n )). Observe, that the truncations min{u, k}, k = 1, 2, . . . belong to L p loc (R; W 1,p loc (R n )) and are weak solutions to (1.1) by Theorem 2.12. It is easy to see that q 2−p −q +p−n > 0 for any q ≥ p, and q 2−p ≥ −q +p+ q 2−p , so that |x| q−p− q 2−p ≤ |x| − q 2−p for every x ∈ B(0, 1). By choosing c = c(n, p, q) = (2 − p) q 2−p p−1 q 2−p − q + p − n , (4.1) we have for every (x, t) ∈ B(0, 1) × (0, ∞). Lemma 2.8 implies that u is a p-supercaloric function in B(0, 1)×(0, ∞). For any ε > 0, by choosing q > n ε (2−p), we obtain a p-supercaloric function u / ∈ L ε loc (B(0, 1) × (0, ∞)). This shows that, in general, a p-supercaloric function is not locally integrable to any positive power. The same example with the constant in (4.1) is a supersolution to (1.1) also for 1 < p ≤ 2n n+1 , with the additional requirement q > (n−p)(2−p) The function u is p-supercaloric for any c > 0 with 0 < q ≤ (n−p)(2−p) p−1 and 1 < p < 2, since then for (x, t) ∈ R n × (0, ∞).
The following intrinsic weak Harnack inequality for supersolutions can be found in [11, Proposition 3.1].

Lemma 4.3. Let 2n
n+1 < p < 2 and let Ω be an open set in R n . Assume that u is a nonnegative lower semicontinuous weak supersolution to (1.1) in Ω T . There exist constants c 1 = c 1 (n, p) ∈ (0, 1) and c 2 = c 2 (n, p) ∈ (0, 1) such that, for almost every s ∈ (0, T ), we have Next lemma will be a useful building block in the characterization of the complementary class of the Barenblatt type supercaloric functions.