arxiv: v2 [math.ap] 24 Sep 2010

Transcription

1 LOCAL HÖLDER CONTINUITY FOR DOUBLY NONLINEAR PARABOLIC EQUATIONS TUOMO KUUSI, JUHANA SILJANDER AND JOSÉ MIGUEL URBANO arxiv:16.781v [math.ap] 4 Sep 1 Abstract. We give a proof of the Hölder continuity of weak solutions of certain degenerate doubly nonlinear parabolic equations in measure spaces. We only assume the measure to be a doubling non-trivial Borel measure which supports a Poincaré inequality. The proof discriminates between large scales, for which a Harnack inequality is used, small scales, that require intrinsic scaling methods. 1. Introduction We consider the regularity issue for nonnegative weak solutions of the doubly nonlinear parabolic equation u p 1 ) u p u) =, p <. 1.1) t This equation is a prototype of a parabolic equation of p-laplacian type. Its solutions can be scaled by nonnegative factors, but in general we cannot add a constant to a solution so that the resulting function would be a solution to the same equation. The purpose of this paper is to obtain a clear transparent proof for the local Hölder continuity of nonnegative weak solutions of 1.1). Our work is a continuation to [17], where Harnack s inequality for the same equation is proved. See also [1], [11], [1] [4]. However, since we cannot add constants to solutions, the Harnack estimates do not directly imply the local Hölder continuity. To show that our proofis based on a general principle, we consider the case where the Lebesgue measure is replaced by a more general Borel measure, which is merely assumed to be doubling to support a Poincaré inequality. In the weighted case, parabolic equations have earlier been studied in [1], [] []. See also [8]. This kind of doubly nonlinear equations have been considered by Vespri [3], Porzio Vespri [19], Ivanov [14], [15]. The known regularity proofs are based on the method of intrinsic scaling, originally introduced by DiBenedetto, they seem to depend highly on the particular form of the equation. However, the passage from one equation to another is not completely clear. For other parabolic equations, the problem has been studied at length, see [4], [3], [7] [], the references therein. The difficulty with equation 1.1) is that there is a certain kind of dichotomy in its behavior. Correspondingly, the proof has been divided in Mathematics Subject Classification. Primary 35B65. Secondary 35K65, 35D1. Key words phrases. Hölder continuity, Caccioppoli estimates, intrinsic scaling, Harnack s inequality. Research of JMU supported by CMUC/FCT project UTAustin/MAT/35/8. 1

2 KUUSI, SILJANDER AND URBANO two complementary cases: Case I : essinfu << essoscu Case II : u p u t Cu t. In large scales, i.e., in Case I, the scaling property of the equation dominates, consequently, the reduction of the oscillation follows immediately from Harnack s inequality. In small scales, on the other h, the oscillation is already very small thus the solution itself is between two constants, the infimum the supremum, whose difference is negligible. Correspondingly, the nonlinear time derivative term, which formally looks like u p u t, behaves like u t we end up with a p-parabolic type behavior. However, also in this case, westill need to modifythe known arguments. In particular, the energy estimates are not available in the usual form we need to use modified versions as in [5], [15] [5]. Our argument also applies to doubly nonlinear equations of p-laplacian type that are of the form u p 1 ) t Ax,t,u, u) =, with Ax, t,, ) satisfying the usual structure assumptions. For expository purposes, we only consider 1.1). Very recently, a direct geometric method to obtain local Hölder continuity for parabolic equations has been developed in [6] [9]. Despite the effort, the general picture remains unclear.. Preliminaries Let µ be a Borel measure be an open set in R d. The Sobolev space H 1,p ;µ) is defined to be the completion of C ) with respect to the Sobolev norm 1/p u 1,p, = u p + u )dµ) p. A function u belongs to the local Sobolev space H 1,p loc ;µ) if it belongs to H 1,p ;µ) for every. Moreover, the Sobolev space with zero boundary values H 1,p ;µ) is defined as the completion of C ) with respect to the Sobolev norm. For more properties of Sobolev spaces, see e.g. [13]. Let t 1 < t. The parabolic Sobolev space L p t 1,t ;H 1,p ;µ)) is the space of functions ux,t) such that, for almost every t, with t 1 < t < t, the function u,t) belongs to H 1,p ;µ) t u p + u p )dν <, t 1 where we denote dν = dµdt. The definition of the space L p loc t 1,t ;H 1,p loc ;µ)) is analogous.

3 LOCAL HÖLDER CONTINUITY 3 Definition.1. A function u L p loc t 1,t ;H 1,p loc ;µ)) is a weak solution of equation 1.1) in t 1,t ) if it satisfies the integral equality t ) u p u φ u p 1 φ dν =.) t for every φ C t 1,t )). t 1 Next, we recall a few definitions results from analysis on metric measure spaces. The measure µ is doubling if there is a universal constant D 1 such that µbx,r)) D µbx,r)), for every Bx,r). Here Bx,r) = {y R d : y x < r} denotes the stard open ball in R d. Let < r < R <. A simple iteration of the doubling condition implies that ) µbx,r)) R dµ µbx,r)) C, r where d µ = log D. A doubling measure in also satisfies the following annular decay property. There exist constants < α < 1 c 1 such that µbx,r)\bx,1 δ)r)) cδ α µbx,r)),.3) for all Bx,r) < δ < 1. The measure is said to support a weak 1, p)-poincaré inequality if there exist constants P > τ 1 such that 1/p u u Bx,r) dµ P r u dµ) p, Bx,r) Bx,τr) for every u H 1,p loc ;µ) Bx,τr). Here, we denote 1 u Bx,r) = udµ = udµ. µbx,r)) Bx,r) Bx,r) The word weak refers to the constant τ, that may be strictly greater than one. InR d withadoublingmeasure, theweak 1,p)-Poincaré inequality with some τ 1 implies the 1,p)-Poincaré inequality with τ = 1, see Theorem 3.4 in [1]. Hence, we may assume that τ = 1. On the other h, the weak 1, p)-poincaré inequality the doubling condition imply a weak κ, p)-sobolev-poincaré inequality with d µ p κ = d µ p, 1 < p < d µ,.4) p, p d µ, where d µ is as above. In other words, Poincaré doubling imply the Sobolev inequalities. More precisely, there are constants C > τ 1 such that ) 1/κ 1/p u u Bx,r) κ dµ Cr u dµ) p,.5) Bx,r) Bx,τ r)

4 4 KUUSI, SILJANDER AND URBANO for every Bx,τ r). The constant C depends only on p, D P. For the proof, we refer to [1]. Again, by Theorem 3.4 in [1] we may take τ = 1 in.5). For Sobolev functions with the zero boundary values, we have the following version of Sobolev s inequality. Suppose that u H 1,p Bx,r);µ). Then 1/κ 1/p u dµ) κ Cr u dµ) p..6) Bx,r) Bx,r) For the proof we refer, for example, to [18]. Moreover, by a recent result in [16], the weak 1, p)-poincaré inequality the doubling condition also imply the1, q)-poincaré inequality for some q < p, that is 1/q u u Bx,r) dµ Cr u dµ) q..7) Bx,r) Bx,r) Consequently, also.5) holds with p replaced by q. We also obtain the q,q)-poincaré inequality for some q < p. In the sequel, we shall refer to data as the set of a priori constants p, d, D, P. Our main result is the following theorem. Theorem.8. Let p < assume that the measure is doubling, supports a weak 1, p)-poincaré inequality is non-trivial in the sense that the measure of every non-empty open set is strictly positive the measure of every bounded set is finite. Moreover, let u be a weak solution of equation 1.1) in R d. Then u is locally Hölder continuous. We will use the following notation for balls cylinders, respectively: Br) = B,r) Q t s,r) = Br) t s,t). For simplicity, we will also denote Qs,r) = Q s,r) = Br) s,). Recall Harnack s inequality from [17]. Theorem.9. Let 1 < p < suppose that the measure µ is doubling supports a weak 1, p)-poincaré inequality. Moreover, let u be a weak solution to 1.1) in R d. Then there exists a constant H = H p,d,d,p,t s r p ))/r p ) such that where s > t+r p. Proof. See [17]. ess supu H essinf u, Q tr p,r) Q sr p,r) In addition, in [17] it is also proved that all solutions of equation 1.1) are locally bounded. In the sequel, we will assume this knowledge without any further comments.

5 LOCAL HÖLDER CONTINUITY 5 3. Constructing the setting Our proof is based on the known classical argument of reducing the oscillation, see [4], [7] []. However, the equation under study has some intrinsic properties which are not present, for instance, in the case of the p-parabolic equation. In large scales, the scaling property dominates the oscillation reduction follows easily from Harnack s inequality. In this case, the equation resembles the usual heat equation. In small scales, in turn, the equation changes its behavior to look more like the evolution p-laplace equation. Indeed, when we zoom in by reducing the oscillation, the infimum the supremum get closer closer to each other. Consequently, the weight u p in the time derivative term starts to behave like a constant coefficient we end up with a p-parabolic type behavior. Resembling this divide between large small scales, the proof has to be divided in two cases. We study the local) Hölder continuity in a compact set K we choose the following numbers accordingly. Let define Furthermore, choose µ µ essinf K u µ+ esssup u, K ω = µ + µ. small enough so that H +1)µ ω 3.1) holds. We will construct an increasing sequence {µ i } a decreasing sequence {µ + i } such that µ + i µ i = ω i = σ i ω for some < σ < 1. Moreover, these sequences can be chosen so that ess supu µ + i Q i ess inf Q i u µ i, for some suitable sequence {Q i } of cubes. Consequently, ess osc Q i u ω i. The cubes here will be chosen so that their size decreases in a controllable way, from which we can deduce the Hölder continuity. Observe also that if H +1)µ j ω j fails for some j, the above sequences have been chosen so that µ + j µ j < H + 3.) for all j j. We are studying the local Hölder continuity in a compact set K. Our aim is to show that the oscillation around any point in K reduces whenever we suitably decrease the size of the set where the oscillation is studied.

6 6 KUUSI, SILJANDER AND URBANO The next step is to iterate this reduction process. We end up with a sequence of cylinders Q i. For all purposes, in the sequel, it is enough to study the cylinder Q := Qηr p,r) instead of the set K. Indeed, for any point in K we can build the sequence of suitable cylinders, but since we can always translate the equation, we can, without loss of generality, restrict the study to the origin. Theequation 1.1) has its own time geometry too, that we need to respect in the arguments. This is important, especially in small scales, when the equation resembles the evolution p-laplace equation. We will use a scaling factor η = λ 1p )+1 in the time direction, where λ 1 1 is an a priori constant to be determined later. 4. Fundamental estimates We start the proof of Theorem.8 by proving the usual energy estimate in a slightly modified setting, which overcomes the problem that we cannot add constants to solutions, see [5], [15] [5]. We introduce the auxiliary function u p 1 ) Ju k) ± ) =± ξ 1/p 1) k dξ k p 1 ± =±p 1) =p 1) u k u k)± ξ k) ± ξ p dξ k±ξ) p ξdξ. Hence, we have t Ju k) ±) = ± up 1 ) u k) ±. 4.1) t In the sequel, we will also need the following estimates. Clearly, Ju k) + ) = p 1) u k)+ k +ξ) p ξdξ p 1 k +u k) +) p u k) + p 1 up u k) + u k)+ Ju k) + ) p 1)k p ξdξ p 1)k p u k) +. Observe that the assumption p is used here. On the other h, Ju k) ) = p 1) u k) p 1) u p u k). k ξ) p ξdξ 4.) 4.3) 4.4)

7 LOCAL HÖLDER CONTINUITY 7 Moreover, Ju k) ) = p 1) u k) k ξ) p ξdξ u k) p 1)k p ξdξ = p 1)k p u k). Now we are ready for the fundamental energy estimate. 4.5) Lemma 4.6. Let u be a weak solution of 1.1) let k. Then there exists a constant C = Cp) > such that t ess sup Ju k) ± )ϕ p dµ+ u k) ± ϕ p dν t 1 <t<t t C t 1 u k) p ± ϕ p dν +C t 1 t for every nonnegative ϕ C t 1,t )). t 1 ) ϕ Ju k) ± )ϕ p 1 dν, t + Proof. Let t 1 < τ 1 < τ < t. We formally substitute the test function φ = ±u k) ± ϕ p in the equation obtain τ ) = u p u φ+ up 1 ) φ dν τ 1 t τ = u k) p ± u k) ± ) ±u k) ± ϕ p )dν 4.7) τ 1 τ ± τ 1 u p 1 ) u k) ± ϕ p dν. t Now the first term on the right-h side can be estimated pointwise from below as u k) p ± u k) ± ) ±u k) ± ϕ p ) u k) ± p ϕ p p u k) ± p 1 ϕ p 1 u k) ± Dϕ, the last term is estimated further by Young s inequality as p u k) ± p 1 ϕ p 1 u k) ± Dϕ 1 u k) ± p ϕ p Cu k) p ± Dϕ p. Thus we have 1 τ τ 1 τ u k) ± ϕ p dν ± τ τ 1 C u k) p ± ϕ p dν. τ 1 u p 1 ) u k) ± ϕ p dν t 4.8)

8 8 KUUSI, SILJANDER AND URBANO Using 4.1) integrating by parts, τ u p 1 ) τ ± u k) ± ϕ p dν = τ 1 t τ 1 t Ju k) ±)ϕ p dν [ ] τ = Jux,t) k) ± )ϕ p x,t)dµ t=τ 1 τ p Ju k) ± )ϕ p 1 ϕ t dν. τ 1 So we obtain τ Jux,τ ) k) ± )ϕ p x,τ )dµ+ u k) ± ϕ p dν τ 1 τ C Jux,τ 1 ) k) ± )ϕ p x,τ 1 )dµ+c u k) p ± ϕ p dν τ 1 τ ) ϕ +C Ju k) ± )ϕ p 1 dν. τ 1 t + 4.9) Now wecan dropthesecond term fromtheleft hside, let τ 1 t 1, choose τ such that Jux,τ ) k) ± )ϕ p x,τ )dµ 1 esssup Ju k) ± )ϕ p dµ 4.1) t 1 <t<t estimate the limits of integration on the right h side of 4.9). On the other h, we can also drop the first term on the left h side of 4.9) let τ 1 t 1 τ t. Summing the estimates for both terms gives the claim. Let us denote ψ ± u) := ΨH ± k,u k) ±,c) = H ± k ln c+h ± k u k) ± The following logarithmic lemma is used in forwarding information in time. Lemma Let u be a weak solution of equation 1.1). Then there exists a constant C = Cp) > such that ess sup u p ψ u)ϕp dµ k p ψ u)x,t 1)ϕ p x)dµ t 1 <t<t t +C ψ ψ ) p ϕ p dν t 1 )) ess sup k p ψ+u)ϕ p dµ u p ψ+u)x,t 1 )ϕ p x)dµ t 1 <t<t t +C ψ + ψ + ) p ϕ p dν. Above, ϕ C ) is any nonnegative time-independent test function. t 1. +

9 Proof. Choose LOCAL HÖLDER CONTINUITY 9 φ ± u) = u ψ ±u))ϕ p in the definition of weak solution observe that ψ ±) = 1+ψ ± )ψ ±). 4.1) The parabolic term will take the form t t u p 1 t 1 t up 1 φ ± u)dν = φ ± s 1/p 1) )dsdν t 1 t k [ p 1 ] u p 1 t = φ ± s 1/p 1) )dsdµ k p 1 t 1 [ u ] t = p 1) φ ± r)r p drdµ Now an integration by parts gives u k k φ ± r)r p dr = In the plus case, we have trivially u k u k t 1. ψ ±r)) r p drϕ p =ϕ p[ u ψ± r)rp ] u k p ) ψ± r)rp 3 drϕ p =ψ ± u)up ϕ p p ) φ + r)r p dr u k k ψ ± r)rp 3 drϕ p. ψ +u)u p ϕ p ψ +u)u p k p )ϕ p = p 1)ψ + u)kp ϕ p u k φ + r)r p dr ψ + u)up ϕ p, since p. Similar estimates are true also for the minus case. On the other h, by using 4.1) together with Young s inequality, we obtain u p u φ ± = u p u ψ ±u)) ϕ p ) u p 1+ψ ± )ψ ±) ϕ p p u p 1 ψ ± ψ ±ϕ p 1 ϕ 1 u p 1+ψ ± )ψ ±) ϕ p Cψ ± ψ ±) p ϕ p, almost everywhere, from which the claim follows. We will need the following notations in the next lemma, which is the most crucial part of the argument. Let r n = r + r n+1, Q± n = B n T ± n = Br n ) t γ ± r p n,t ) A ± n = { x,t) Q ± n : ±ux,t) > ±k ± n},

10 1 KUUSI, SILJANDER AND URBANO for n =,1,,... Recall the definitions µ + i µ i = ω i = σ i ω, where µ + i esssup Q i u µ i essinf Q i u, for i 1. Observe, however, that we have to choose µ + esssup u µ essinf u, K K where the infimum the supremum are taken over K instead of Q. This is because we need the argument to be independent of the initial cylinder Q. Now we are ready to prove the fundamental lemma everything depends upon. Lemma Let < ε ± 1, k + n ) n be an increasing sequence k n ) n a decreasing sequence, both of nonnegative real numbers. Suppose u is a weak solution of equation 1.1), u k ± n) ± ε ± ω i k ± n+1 k± n ε ±ω i n+. In addition, assume further u 1 C k n 4.14) µ + i k n +, n = 1,, ) for the minus plus cases, respectively. Then there exist constants C = CD,P,C,p) > C + = CD.P,p) > such that νa ± n+1 ) νq ± n+1 ) Cn+1 ± Γ ± ) νa ± p/κ n ) νq ± 4.16) n) for every n =,1,,... Here κ is the Sobolev exponent as in.4) Γ ± = 1 ) k ± p ) ) n γ ± γ ± ε± ω p p/κ i ε ± ω i k n ± +1. Proof. Choose the cutoff functions ϕ ± n C Q± n ) so that ϕ± n 1, ϕ ± n = 1 in Q± n+1 Denote in short ϕ ± n Cn+1 r ϕ ± n t v n = u k n ) ±, k n = k ± n, ε = ε ± Cpn+1) γ ± r p. 4.17) Q n = B n T n = Q ± n, A n = A ± n, γ = γ ±, ϕ n = ϕ ± n.

11 LOCAL HÖLDER CONTINUITY 11 By Hölder s inequality, together with the Sobolev inequality.6), we obtain dν vn 1 p/κ)+p Q n+1 νq n) νq n+1 ) vn 1 p/κ)+p ϕn p1 p/κ)+p dν Q n ) 1 p/κ ) p/κ C vn ϕp n dµ v n ϕ n ) κ dµ dt T n B n B n ) 1 p/κ Cr esssup p vnϕ p ndµ v n ϕ n ) p dν. T n B n Q n Here, we applied the doubling property of the measure ν giving νq n ) νq n+1 ) C. 4.18) We continue by studying the term involving the essential supremum. By the assumption 4.4), we obtain u 1 C k n u k n) p 1 u p Ju k n) ) Ck n) p Ju k n) ). On the other h, the lower bound 4.3) gives immediately u k + n) + Ck + n) p Ju k + n) + ). Using these estimates together with the energy estimate, Lemma 4.6, yields ess sup vn ϕp n dµ Ck n) p esssup Jv n )ϕ p n dµ T n B n T n B n ) ) Ck n ) p γrn p vn ϕ p n p +Jv n )ϕ p 1 ϕn n dν. t Furthermore, the estimates 4.) 4.5) imply Q n Ju k n ) + ) Ck + n )p u k n ) +, Ju k n) ) Ck n )p u k n ). For the plus case, we used 4.15). Next, using 4.17), we arrive at ) ) rn p vn p ϕ n p +Jv n )ϕ p 1 ϕn n dν Q n t + C np v pn + k n) p ) vn dν Q n γ C np εω i ) p 1+ 1 ) ) p εωi νa n ) γ νq n ), k n ) where the last inequality follows from the fact that u k n ) ± ε ± ω i. Thus, we conclude ) p ess sup vn ϕp n dµ Cnp εω i ) εωi νa n ) γ +1) T n B n k n νq n ). 4.)

12 1 KUUSI, SILJANDER AND URBANO Furthermore, since v n ϕ n ) p dν C v n p ϕ p ndν +C vn ϕ p n p dν, Q n Q n Q n applying again the energy estimate 4.19) leads to v n ϕ n ) p dν C np εω i ) p 1+ 1 ) ) p εωi νa n ) Q n γ k n νq n ). 4.1) To finish the proof, note first that It then follows that u k ± n) ± χ {u k ± n ) ± >} u k± n) ± χ {u k ± n+1 ) ±>} vn 1 p/κ)+p Q n+1 dν k ± n+1 k± n n+) ε ± ω i. n+) εω i ) 1 p/κ)+p νa n+1 ) νq n+1 ). 4.) Inserting estimates 4.), 4.1), 4.) into 4.18) concludes the proof. Remark 4.3. If we have the extra knowledge that u kn ) =, almost everywhere in Br) at a given time level, we can choose the test functions to be independent of time the cylinder Q n so that the length in the time direction stays constant, the bottom of the cylinder stays at the given time level. In this case, by choosing a time independent test function, the right h side of the energy estimate simplifies so that we can get rid of the term +1 in the formulation ) ) Γ ± = γ ± ε± ω p 1 p/κ i k n ± in 4.16). This will get us the required extra room in the end of the first alternative of Case II. Furthermore, in the previous lemma, we chose the radii of the cylinder as r n = r + r n+1. However, the factor in the denominator can naturally be replaced by any greater number. Westarttheproofby consideringcasei. Here, weusethepreviouslemma only in the plus case. Consequently, the first case does not depend on the constant C. We recall a lemma on the fast geometric convergence of sequences from [4]. Lemma 4.4. Let Y n ) n be a sequence of positive numbers, satisfying Y n+1 Cb n Y 1+α n 4.5) where C,b > 1 α >. Then Y n ) n converges to zero as n provided Y C 1/α b 1 α. 4.6)

13 LOCAL HÖLDER CONTINUITY 13 On several occasions in the sequel, we use this lemma, together with the fundamental estimate Lemma 4.13, to conclude that a ratio of the form Y n = νa± n) νq ± n) converges to zero consequently that νa ± n) as n. This will ultimately lead to a reduction of the oscillation which is our final goal. Once a recursive inequality of type 4.16) has been established, the convergence to zero of νa ± n ) follows from the condition νa ± ) νq ± ) α±, with α ± = Γ 1/1 p/κ) ± C 1/1 p/κ)+1 1 p/κ) ±, 4.7) where the constants C ± Γ ± are the constants from the previous lemma. Note that an explicit value of α ± only follows after fixing C ± Γ ±. 5. The Case I Now we assume that 3.1) holds. Our aim is to show that the measures of certain distribution sets tend to zero that the local Hölder continuity follows from this. We start by studying the subcylinder Qr p,r) Qηr p,r). Let γ ± = 1, ε ± = 1 k n + = µ + ω 4 ω n+. Observe that, after fixing these quantities, the constant α + can be fixed as well. We will study two different alternatives which are considered in the following two lemmata, respectively. Lemma 5.1. Let λ > 1 be sufficiently large let u be a weak solution of equation 1.1). Furthermore, assume ν {x,t) Br) r p, rp ) : ux,t) µ λ + ω ) } =. 5.) Then there exists a constant σ,1) such that ess osc Q r ) p, r ) u σω. Proof. By the choices preceding the statement of this lemma, we have The assumption 3.1) implies Thus u k + n) + ε + ω. µ + = µ +ω ω. 1 k+ n ε + ω 4.

14 14 KUUSI, SILJANDER AND URBANO Plug these in Lemma 4.13 to deduce νa + n+1 ) νq + n+1 ) Cn+1 ) νa + p/κ n ) νq +. n) On the other h, by 5.) we have the trivial estimate νa + ) νq + ) 1 λ α +, choosing λ > 1 sufficiently large. By Lemma 4.4 we conclude that as n. This implies νa + n) νq + n) ess sup u µ + Q ) r p ω, ) r 4. So, if this alternative occurs, we choose These choices yield as required, with µ + 1 = µ+ ω 4 µ 1 = µ. ess osc Q r ) p, ) u 1 1 ) ω r 4 σ = 3 4. For the second possibility, we have the following lemma. Lemma 5.3. Let u be a weak solution of equation 1.1) suppose ν {x,t) Br) r p, rp ) : ux,t) µ λ + ω ) } >. 5.4) Then there exists a constant σ = σh ),1) such that ess osc ) Q r p, r λ Proof. By assumption 5.4), we have λ )u σω. ess sup u µ + ω Br) r p, rp. ) λ

15 LOCAL HÖLDER CONTINUITY 15 Now we can use Harnack s inequality Theorem.9), together with the Case I assumption 3.1), to deduce Q r ess inf λ ) p, r λ )u 1 H esssup u Br) r p, rp ) λ µ + ω H H µ + µ ω H H +1 + ω H µ + ω H H +1). Observe that the constant H depends on λ, but this does not matter since λ depends only on the data. Now, if we end up in this alternative, we choose We also obtain with as required. ess osc ) Q r p, r λ µ 1 = µ + ω H H +1) µ + 1 = µ+. λ )u ω σ = 1 ω H H +1) = σω, 1 H H +1), 6. The case II In Case II the equation looks like the evolution p-laplace equation. In this case, we need to use the scaling factor η in the time geometry of our cylinders. The difficulty is now that we cannot use the Harnack principle anymore, as the lower bound it gives might be trivial. Indeed, the infimum can be larger than the lower bound Harnack s inequality gives. On the other h, we have the following kind of elliptic Harnack s inequality. Suppose that j is the first index for which assumption 3.1) does not hold. Then we have ω j µ + j H +)µ j. 6.1) Clearly, this Harnack s inequality is valid also for every subset of the initial cylinder Q j = Qηr p,r), consequently, for every j j. Recall, that ω j = σω j 1 ω j H +) µ j µ j 1 +1 σ)ω j 1 σ)ω j 1 σ σ ω j. 6.)

16 16 KUUSI, SILJANDER AND URBANO Thus, we obtain, consequently, ) C 1 Q H +) p ηrp,r σ H +) σ) σ µ j 1 σ ω j Q C 1 µ j ω j ) p ηr p,r where C 1 = σ p / σ) p. We will consider the cylinder ) p ηr,r p Q := Q C 1 µ j ω j. Qηr p,r), By the above calculation, we have shrunk the cylinder by a factor which is controllable by the data. In the sequel, we will denote Recall the definitions θ = C 1 µ j ω j ) p. 6.3) Q ± n = Q t γ ± r p n,r n ) = B n T n = Br n ) t γ ± r p n,t ) A ± n = {x,t) Q ± n : ±u > ±k ± n}. Now the proof will follow the classical argument of DiBenedetto, see [4] [], is again divided into two alternatives. In the first one, we assume that there is a suitable cylinder for which the set where u is close to its infimum is very small. In the second alternative, we assume that this does not hold true The First Alternative. We first suppose that there exists a constant α,1) to be determined in the course of the next lemma, depending only on the data) such that ν {x,t) Q : u < µ j + ω ) j } α νq ), for a cylinder Q = Q t θrp,r) Q λ 1p )+1 θr p,r). Our aim is to use Lemma 4.13 to conclude for the reduction of the oscillation. Lemma 6.4. For every s > 3, { r ν x,t) Q θ 4 where θ is as in 6.3). ) p, r ) 4 C λ 1p ) s s 3) ν Q : ux,t) < µ j + ω }) j s r ) p, r )) θ, 4 4

17 LOCAL HÖLDER CONTINUITY 17 Proof. We start by using Lemma 4.13, with the choices r n = r + r n+1, k n = µ j + ω j 4 + ω j n+, ε = 1/ γ = θ. We also need the assumption 6.1) to deduce that u 1 C µ j +H +)µ j ) 1 C µ j +ω j ) 1 C k n, 6.5) with C = 3H +). This knowledge is needed in Lemma Now, after fixing ε, γ, kn C we can fix α, see 4.7). We also obtain, using 6.1) 6.3), the bounds p 1 ) k p n C 1 γ p H +) p ε ω j C 1 thus we can conclude νa n+1 ) νq n+1 ) Cn+1 νa ) p/κ n ) νq. n) By the assumption of this alternative, together with the lemma of fast geometric convergence Lemma 4.4), we have u > k almost everywhere in Q t θr/) p,r/). Thus u k) = consequently H )) k ψu) := ln c+h k u k) almost everywhere in Q t θr/) p,r/). Let r ) p t θ 6.6) 4 + = be a time level such that this is true for almost every x Br/). Now our goal is to apply Lemma 4.11 with k = µ j + ω j 4, c = ω j s H k = esssup u k), Q where Q = Qηθr p,r). Choose ϕ C Br/)) independent of time such that ϕ 1, ϕ = 1 in Br/4) In the set {u < µ j + ω j s }, we have, on the other h, ϕ C r. ψ s 3) ln, ψ s )ln ψ p ωj ) p.

18 18 KUUSI, SILJANDER AND URBANO The use of these estimates in Lemma 4.11 gives µ j ) p s 3) ln µ{x Br/4) : ux,t) < µ j + ω j s }) esssup u p ψ u)x,t)ϕ p x)dµ +C t <t< Br/) t Br/) k p ψ u)x,t )ϕ p x)dµ Br/) ψ ψ p ϕ p dµdt λ 1 ) p ω j Cs )ln θµbr/4)) λ 1 µ ) p j Cs )ln µbr/4)), for almost every t t,). Observe that, in the third inequality, we plugged in η = λ 1p )+1. The claim follows by integrating this estimate over θr/4) p,). We conclude this alternative with the following two lemmata. Lemma 6.7. Let u be a weak solution of equation 1.1) assume 6.1) holds. Then u µ j + ω j r ) p, r ) s+1 a.e. in Q θ, 8 8 where s depends only upon the data, θ is as in 6.3). Proof. Let r n = r 8 + r n+3, Q n = B n T = Br n ) t,), where t is, as in the previous lemma, such that r ) p, u k) x,t ) =, t θ 4 for a.e. x Br/). Moreover, define kn = µ j + ω j s+1 + ω j In this case, we obtain s+n+1. u k n) ε ω j, where ε = 1 s. Observe also that γ = t /r p ηθ = λ1p )+1 θ. We will substitute these in Lemma 4.13, taking into account Remark 4.3 estimate 6.5), we conclude as before that νa n+1 ) νq n+1 ) Cn+1 λ 1p ) sp ) ) 1 p/κνa ) p/κ n ) νq. n)

19 LOCAL HÖLDER CONTINUITY 19 Now choose s > λ 1. Then, by Lemma 4.4, we have νa n)/νq n) as n, provided νa )/νq ) is small enough. On the other h, by choosing s large enough, Lemma 6.4 guarantees that νa )/νq ) can be chosen to be as small as we please. This gives u µ j + ω j s+1 a.e. in Q t, r ) 8 hence the lemma is proved. Lemma 6.8. There exists < σ < 1, depending only upon the data, such that ess osc Qθ r 8) p, r 8) u σω j. Proof. By the previous lemma, ess inf Qθ r 8) p, r 8) u µ j + ω j s+1, for some s > 1, which depends only upon the data λ 1. Observe that here we used the knowledge r ) p. t θ 8 If this alternative occurs, we again choose Finally, we get as required, with µ j +1 := µ j + ω j s+1 µ + j +1 = µ+ j. ess osc Qθ 8) r p, 8) u 1 1 ) r s+1 ω j σ = 1 1 s+1. Remark 6.9. Here the choice of s is possible only after λ 1 has been determined in the second alternative. Nevertheless, both of them are a priori constants which can be assigned explicit values depending only upon the data. 6.. The Second Alternative. In the second alternative, the assumption of the first alternative is not true. In this case, for every cylinder Q t θr p,r) Qηθr p,r), we have ) ν {x,t) Q t θr p,r) : ux,t) µ j + ω j } νq t θr p < 1 α ), 6.1),r)) where α := α is the same constant as in the first alternative. This implies that, for every t η 1)θr p,), there exists a time level t with t θr p t t θα rp

20 KUUSI, SILJANDER AND URBANO for which µ { x Br) : ux,t ) > k }) 1 α 1 α µbr)). Indeed, otherwise we would have ν { x,t) Q t θr p,r) : ux,t) > k }) 6.11) t θα rp µ { x Br) : ux,t) > k }) dt t θr p > 1 α )νq t θr p,r)), which contradicts 6.1). This alternative is also based on Lemma We choose λ 1 in the definition of k n + large enough so that we can force νa+ ) to be small compared to νq + ). We start with forwarding the information of 6.11) in time. Lemma 6.1. There exists s >, depending only upon the data, such that { µ x Br) : ux,t) > µ + j ω }) 1 3α j 4 s 1 α µbr)). for almost all t t,). Proof. Let c = ω j s+n, k = µ+ j ω j s H + k = esssup u k) +, Q where s n will be chosen later Q := Qθηr p,r). Our aim is again to use Lemma 4.11 to forward the information in time. We will need some estimates for doing this. Recall the definition ψ + u) = ΨH + k,u k) +,c) = ln + H + ) k c+h + k u k). + Trivially, we have ψ + u) ln, on the other h, in the set we get ψ + u) ln The last estimate we need is 1 ψ + ) p c+h + k ωj ) s ω j = nln s+n {u > l µ + j ω j s+n}, ) ω j s ω j + ω j = n 1)ln s+n s+n ) p p ω j s ) p.

21 LOCAL HÖLDER CONTINUITY 1 Let now ϕ C Br)) be a cutoff function which is independent of time has the properties ϕ 1, ϕ = 1 in B1 δ)r) ϕ 1 δr, where < δ < 1 is to be determined later. Apply Lemma 4.11 with these choices to conclude n 1) ln µ{x B1 δ)r) : ux,t) > l}) esssup ψ+u)x,t)ϕ p x)dµ t <t<t Br) u ) p ψ Br) k + u)x,t )ϕ p x)dµ t +Ck p ψ + ψ + ) p ϕ p dµdt t n ln ) +C nln δ p Br) µ+ j µ + j µ+ j s ω j µ + j p ) p 1 s 1 1 α 1 α µbr) {t }) ) p θµbr)) for almost every t t,t ). Observe, that in the third inequality we used 6.11). Now, by the annular decay property.3), we have µ{x Br) : ux,t) > l}) µbr)\b1 δ)r))+µ{x B1 δ)r) : ux,t) > l}) Cδ α µbr))+µ{x B1 δ)r) : ux,t) > l}). For the first term, we choose δ small enough so that α Cδ α < 16 1 α ) for the second term we use the previous estimate. Indeed, by choosing s n large enough so that ) p 1 α n 1 3α 1 α n 1) s 8 1 α ) ) Cn 1 p α ln)δ p n 1) s α ), we get the claim for almost every t t,t ) with s = s+n. Recall that t θr p t t θα rp.

22 KUUSI, SILJANDER AND URBANO Finally, since the above holds for every cylinder Q t θr p,r) Qηθr p,r), we can conclude that for almost every t ηθr p +θr p θα r p = η +1 α ) θr p, we have 1 3α µ{x Br) : ux,t) > l}) < 4 1 α µbr)). Remark Now we can choose η so large that the previous lemma holds for almost every t ηθ rp, ), i.e., η +1 α η hence η = λ1p )+1 α. But this is always guaranteed for λ 1 > 1 p. We are ready to prove the final estimate, which, together with Lemma 4.4, gives the reduction of the oscillation. Let E t) = {x Br) : ux,t) > µ + j ω j } ) ηθ E = {x,t) Q rp,r Then we have the following lemma. : ux,t) > µ + j ω j }. Lemma For every α 1,1), there exists λ 1 > such that Proof. Denote νe λ1 ) νq ηθ rp,r)) α 1. h = µ + j ω j s+1 k = µ + j ω j s, where s > will be chosen large. Let also h k, u h, v = u k, k < u < h,, u k. By the previous lemma, we can choose s large enough, namely s s, so that, for almost every t ηθ rp, ), we have µx Br) : vx,t) = }) = µ{x Br) : ux,t) k}) α 4 α µbr) {t}) α 4 µbr)).

23 LOCAL HÖLDER CONTINUITY 3 Thus, for almost every t ηθ rp, ), we obtain v Br) t) = vdµ 1 α ) h k) 4, consequently, Br) {t} h k v Br) t) α 4 h k). Using the q, q)-poincaré inequality for some q < p see.7) the remark after that), yields ) 4 q h k) q µe s+1 t)) v v α Br) t) q dµ Br) {t} Cr q v q dµ = Cr q u q dµ, Br) {t} E st)\e s+1 t) for almost every t ηθ rp, ). The constant 4/α ) q above was absorbed into the constant C. Now we integrate the above inequality over time to get h k) q νe s+1 ) Cr q E s\e s+1 u q dν. Next, we introduce a cutoff ) function ϕ C Qηθrp,r)) such that ϕ 1, ϕ = 1 in Q ηθ rp,r ϕ C ϕ r t C ηθr p. Now Hölder s inequality gives h k) q νe s+1 ) Cr q E s\e s+1 u p dν ) q/p νe s \E s+1) 1 q/p q/p Cr q u k) + p ϕ dν) p νe s \E s+1 ) 1 q/p. Qηθr p,r) By choosing λ 1 > s s in the definition of η large enough, the first factor on the right h side can be estimated by Lemma ) as u k) + p ϕ p dν Qηθr p,r) C Qηθr p,r) u k) p + ϕ p dν +C Ju k) + )ϕ p 1 ϕ Qηθr p,r) t dν ) C ωj ) p µ + j + ) p r p s u k) ηθ +dν Qηθr p,r) C ωj ) )) pν ηθ r p s Q rp,r. 6.15)

24 4 KUUSI, SILJANDER AND URBANO In the last inequality we used the doubling property of the measure ν. We obtain ωj ) qνes+1 ωj ) )) qν ηθ q/p s+1 ) C s Q rp,r νe s \E s+1) 1 q/p. Finally, summing s over s,...,λ 1 1 gives )) ηθ q/p q) )) ηθ λ 1 s )νe λ1 ) p/p q) Cν Q rp,r ν Q rp,r hence )) C ηθ νe λ1 ) λ 1 s ) p q)/pν Q rp,r. Choosing λ 1 large enough finishes the proof. Lemma Suppose that 6.1) holds. Then there exists < σ < 1, depending only upon the data, such that Proof. Let ess osc u σω j. Q ηθ rp,r) Q + n = Br n ) γ + r p n,), r n = r + r n+1 A + n as before. Substituting γ + = λ 1p ) θ, ε + = 1/ λ 1 k n + = µ+ j ω j λ 1+1 ω j λ 1+n+1 in Lemma 4.13, using 6.1) to bound 1 1 k + ) p n γ + H +) p, ε + ω j yields νa + n+1 ) νq + n+1 ) Cn+1 νa + ) p/κ n ) νq +. n) By the previous Lemma, we can choose λ 1 large enough so that νa + ) νq + ) is as small as we please. Consequently, by Lemma 4.4, we obtain ess sup u µ + Q ηθ r ) p j ω j, ) r λ 1+1, for some λ 1 > 1, which depends only upon the data. So if this alternative occurs, we choose µ + 1 = µ+ j ω j λ 1+1 µ 1 = µ j.

25 LOCAL HÖLDER CONTINUITY 5 We also obtain as required, with ess osc Q ηθ ) r p, ) u 1 1 ) r λ ω 1+1 j, σ = 1 1 λ 1+1. Now the Hölder continuity follows by stard iterative real analysis methods, see chapter III of [4], or []. Acknowledgments. The authors would like to thank Juha Kinnunen for interesting discussions on the subject of this paper. References [1] Filippo Chiarenza Raul Serapioni. Degenerate parabolic equations Harnack inequality. Ann. Mat. Pura Appl. 4), 137:139 16, [] Filippo M. Chiarenza Raul P. Serapioni. A Harnack inequality for degenerate parabolic equations. Comm. Partial Differential Equations, 98): , [3] Emmanuele DiBenedetto. On the local behaviour of solutions of degenerate parabolic equations with measurable coefficients. Ann. Scuola Norm. Sup. Pisa Cl. Sci. 4), 133): , [4] Emmanuele DiBenedetto. Degenerate parabolic equations. Universitext. Springer- Verlag, New York, [5] Emmanuele DiBenedetto Avner Friedman. Hölder estimates for nonlinear degenerate parabolic systems. J. Reine Angew. Math., 357:1, [6] Emmanuele DiBenedetto, Ugo Gianazza, Vincenzo Vespri. A geometric approach to the hölder continuity of solutions to certain singular parabolic partial differential equations. Preprint, 1. [7] Emmanuele DiBenedetto, José Miguel Urbano, Vincenzo Vespri. Current issues on singular degenerate evolution equations. In Evolutionary equations. Vol. I, Hb. Differ. Equ., pages North-Holl, Amsterdam, 4. [8] Eugene B. Fabes, Carlos E. Kenig, Raul P. Serapioni. The local regularity of solutions of degenerate elliptic equations. Comm. Partial Differential Equations, 71):77 116, 198. [9] Ugo Gianazza, Mikhail Surnachev, Vincenzo Vespri. A new proof of the hölder continuity of solutions to p-laplace type parabolic equations. Preprint, 1. [1] Ugo Gianazza Vincenzo Vespri. A Harnack inequality for solutions of doubly nonlinear parabolic equations. J. Appl. Funct. Anal., 13):71 84, 6. [11] Ugo Gianazza Vincenzo Vespri. Parabolic De Giorgi classes of order p the Harnack inequality. Calc. Var. Partial Differential Equations, 63): , 6. [1] Piotr Hajlasz Pekka Koskela. Sobolev met Poincaré. Mem. Amer. Math. Soc., ):x+11,. [13] Juha Heinonen, Tero Kilpeläinen, Olli Martio. Nonlinear potential theory of degenerate elliptic equations. Oxford Mathematical Monographs. The Clarendon Press Oxford University Press, New York, 1993., Oxford Science Publications. [14] Alexer V. Ivanov. Uniform hölder estimates for generalized solutions of quasilinear parabolic equations that admit double degeneration. Algebra i Analiz, 3): , Translation in St. Petersburg Math. J ), no., [15] Alexer V. Ivanov. Hölder estimates for a natural class of equations of fast diffusion type hölder estimates for equations of fast diffusion type. Zap. Nauchn. Sem. S.- Peterburg. Otdel. Mat. Inst. Steklov. POMI), 911):9 3, Translation in J. Math. Sci. New York) ), no. 6,

26 6 KUUSI, SILJANDER AND URBANO [16] Stephen Keith Xiao Zhong. The Poincaré inequality is an open ended condition. Ann. of Math. ), 167): , 8. [17] Juha Kinnunen Tuomo Kuusi. Local behaviour of solutions to doubly nonlinear parabolic equations. Math. Ann., 3373):75 78, 7. [18] Juha Kinnunen Nageswari Shanmugalingam. Regularity of quasi-minimizers on metric spaces. Manuscripta Math., 153):41 43, 1. [19] Maria M. Porzio Vincenzo Vespri. Hölder estimates for local solutions of some doubly nonlinear degenerate parabolic equations. J. Differential Equations, 131): , [] Mikhail Surnachev. A harnack inequality for weighted degenerate parabolic equations. To appear in J. Differential Equations, 1. [1] Neil S. Trudinger. Pointwise estimates quasilinear parabolic equations. Comm. Pure Appl. Math., 1:5 6, [] José Miguel Urbano. The method of intrinsic scaling, volume 193 of Lecture Notes in Mathematics. Springer-Verlag, Berlin, 8. A systematic approach to regularity for degenerate singular PDEs. [3] Vincenzo Vespri. On the local behaviour of solutions of a certain class of doubly nonlinear parabolic equations. Manuscripta Math., 751):65 8, 199. [4] Vincenzo Vespri. Harnack type inequalities for solutions of certain doubly nonlinear parabolic equations. J. Math. Anal. Appl., 1811):14 131, [5] Shulin Zhou. Parabolic Q-minima their application. J. Partial Differential Equations, 74):89 3, Addresses: T.K.: Aalto University, Institute of Mathematics, P.O. Box 111, FI-76 Aalto, Finl. tuomo.kuusi@tkk.fi J.S.: Aalto University, Institute of Mathematics, P.O. Box 111, FI-76 Aalto, Finl. juhana.siljer@tkk.fi J.M.U.: CMUC, Department of Mathematics, University of Coimbra, Coimbra, Portugal. jmurb@mat.uc.pt