arxiv: v1 [math.ap] 22 Oct 2018

Transcription

1 ON GLOBAL L ESTIMATES FOR SYSTEMS WITH p-growth IN ROUGH DOMAINS MIROSLAV BULÍČEK, SUN-SIG BYUN, PETR KAPLICKÝ, JEHAN OH, AND SEBASTIAN SCHWARZACHER arxiv:8.944v [math.ap] Oct 8 Abstract. We study regularity results for nonlinear parabolic systems of p-laplacian type with inhomogeneous boundary and initial data, with p n,. We show bounds on n+ the gradient of solutions in the Lebesgue-spaces with arbitrary large integrability exponents and natural dependences on the right hand side and the boundary data. In particular, we provide a new proof of the global non-linear Calderón Zygmund theory for such systems. This extends the global result of [3] to very rough domains and more general boundary values. Our method makes use of direct estimates on the solution minus its boundary values and hence is considerably shorter than the available higher integrability results. Technically interesting is the fact that our parabolic estimates have no scaling deficit with respect to the leading order term. Moreover, in the singular case, p n,], any scaling deficit can n+ be omitted.. Introduction We are interested in higher integrability results for solutions of the following nonlinear problem: for given data, i.e., an n-dimensional domain R n with n, a length of the time interest T >, a source term f :,T R n N with N N, an initial condition u : R N, a given lateral boundary condition g :,T R N and a mapping A : R R n R n N R n N, we aim to find a function u : [,T R N satisfying t u divat,x, u = divf in,t,. u = g on,t, u, = u in. Owing to a significant number of problems originating in various applications, it is natural to reuire that A is a Carathéodory mapping, satisfying the p-coercivity, the p -growth and the strict monotonicity conditions. It means that..3 A,,η is measurable for any fixed η R n N, At,x, is continuous for almost all t,x R R n, and for some p,, there exist positive constants c and c such that for almost all t,x R R n and all η,η R n N c η p c At,x,η η At,x,η c + η p p-coercivity, p -growth, At,x,η At,x,η η η monotonicity. If for all η η the ineuality.6 is strict, then A is said to be strictly monotone. Mathematics Subject Classification. 35K5, 35K55, 35B65, 35A. Key words and phrases. asymptotically Uhlenbeck problem, nonlinear parabolic system, gradient estimate, Reifenberg domain. M. Bulíček s and S. Schwarzacher s work is supported by the project LL financed by the Ministry of Education, Youth and Sports, Czech Republic and by the University Centre for Mathematical Modelling, Applied Analysis and Computational Mathematics Math MAC. M. Bulíček is a member of the Nečas Center for Mathematical Modeling.

2 M. BULÍČEK, S. BYUN, P. KAPLICKÝ, J. OH, AND S. SCHWARZACHER The leading motive of the paper is to find assumptions on the data, i.e.,, A, u, g and f, that would guarantee not only the existence of a weak solution u L p,t;w,p ;R N but also its higher integrability in the sense that for any [, depending on the data, we have.7 u L p,t;w,p ;R N. Immediately, there appear two key difficulties when one intends to establish the validity of.7: the nonlinear and vectorial structure of the operator A and the roughness of the domain, which is supposed to be just an open bounded set. Both these difficulties will finally lead to some restrictions that we shall describe below. To introduce the novelty of our results, we give a short history on the derivation of.7 for the model elliptic case. Let g W,p ;R N and be an open bounded domain such that { div u p u = in.8 u = g on. It is classical that u = argmin v p dx v g+w,p p is the uniue solution to the system.8. After the seminal works of Ural ceva [37] and Uhlenbeck [36], solutions to systems of.8 are known to be locally in C,α for some α,. In the scalar case, this was extended up to the boundary by Lieberman [8], provided u C,α. The result by Lieberman was the starting point for many results in the scalar case. In particular, higher integrability.7 was shown for various generalizations of the non-linearity, the coefficients, the right hand sides and the boundary regularity [, 3, 9, 6,, 3, 5, 7, 3, 34, 33]. All these results relay on the so called comparison principle that is in the heart of the non-linear Calderón Zygmund theory first developed by Iwaniec [4, 5], see also [8]. Later, building on the seminal works of DiBenedetto and Friedman [, ], some higher regularity results could be transferred to the parabolic setting [, 3, 4,, 4]. However, as far as we are concerned, as of now up to the boundary, C,α regularity is not available for systems in case of inhomogeneous boundary data, as the scalar techniues are not applicable. The best regularity available for systems elliptic and parabolic are global Lipschitz bounds [4], provided the boundary of the domain is smooth. These could be used to prove global higher integrability results for uite general operators [3]. In this paper we provide a new and independent proof of the global Calderón Zygmund theory for elliptic and parabolic systems. In contrast to the available higher integrability results of the kind.7, we show estimates directly for the difference v = u g. The key observation is, that v satisfies the following PDE where the inhomogeneity of the boundary values is transferred to a right hand side: { div v p v = div v p v u p u in v = on. As was already pointed out by Uhlenbeck, local regularity for homogeneous boundary data follows by the mere fact that solutions can locally be reflected. Indeed, local estimates up to the boundary for homogeneous boundary data but inhomogeneous right hand side are in [5]. Our proof of.7 follows by combining local estimates for homogeneous boundary values with algebraic properties of the tensor A. The key estimate related to homogeneous boundary values that may well beof independentinterest is an estimate forthe gradient of theparabolic p-laplacian with no scaling deficit for the leading order terms Theorem 3.3. It can then be used to show respective estimates with no scaling deficit on the leading order terms for. see Theorem.5. Moreover, in the singular case, p <, any scaling deficit can be omitted

3 HIGHER INTEGRABILITY FOR PARABOLIC SYSTEMS 3 see Remark.6. This considerably improves the estimates in [3] uantitatively which we consider to be a major scientific contribution of the present work. Additionally we extend the theory of the global non-linear Calderón Zygmund theory also ualitatively by allowing more general assumptions on the regularity of the boundary and the boundary values in comparison to [3]. Finally we wish to remark, that the assumptions on the non-linearity A in this paper are essentially euivalent to the assumptions in [3]... Assumptions on the boundary. First, concerning the roughness of, we should be able to specify at least in which sense the boundary and the initial data, respectively, are attained. Since goes beyond the Lipschitz category, we can hardly define a trace of any Sobolev or Bochner-Sobolev function by using any kind of a trace theorem and therefore we need to proceed more carefully. Nevertheless, we can be inspired by the elliptic setting. Indeed, if one considers the problem { divax, u = in, u = g on, then a concept of weak solution attaining the trace g can benaturally expressed as follows: to look for u = g +v with v W,p ;R N such that the above euation is satisfied in the sense of distribution and g is fixed such that it attains g on in a suitable way. Notice that here we defined for any p [,, W,p ;R N := {u C ;RN },p, where,p is a standard Sobolev norm. Conseuently, we see that the natural assumption on g is that it must be extendable into the interior of such that its extension g satisfies g W,p ;R N integrability result, we need to improve also theintegrability of the extension g simultaneously. In the parabolic setting, we have to face the very similar problems on the parabolic cylinder,t, and we also reuire that the initial data are attained at least weakly. Therefore, for the nonlinear problem. having the nonlinearity A of p -growth, we are directly led to the following characterization of the boundary data, which will be done through the extension into the interior of the parabolic cylinder,t. Definition.. Let p, and [, be arbitrary given. We say that the couple u,g where u : R N and g :,T R N, is an p,-suitable couple of boundary and initial data, if there exists g W, p,t;w, p ;R N L p,t;w, p ;R N such that and lim t lim ε B εx gt,x u x ϕxdx = for all ϕ L ;R N gt,y g t,x dy = for almost all t,x,t. Before we demonstrate that the above definition is really natural for the setting, we consider in the paper, we shall briefly define all symbols used in the above definition. First, we used the notation fxdx := fxdx, Z Z where Z is a bounded subset of R n. Second, for arbitrary s,, we also defined the dual space, W,s ;R N := W,s ;R N Z

4 4 M. BULÍČEK, S. BYUN, P. KAPLICKÝ, J. OH, AND S. SCHWARZACHER where s is the dual exponent to s, i.e., s := s/s. Finally, for X being a Banach space, we introduced the Bochner or Sobolev-Bochner spaces, respectively, in a usual way. In addition, to simplify the notation, we will also use the following abbreviation in what follows: X p, := W, p,t;w, p ;R N L p,t;w, p ;R N, X p, := W, p,t;w, p ;R N L p,t;w, p ;R N. Next, let usdiscuss thechoice of the suitable coupleaccording to Definition.. Incase, we set =, then dueto the growth and coercivity assumptions on the nonlinear operator A, see.4.5, then L p -compatible data mean nothing else that, for f L p,t;l p ;R n N, we are able to construct a weak solution u L p,t;w,p ;R N which attains in the sense of Definition. boundary and initial data. Then the information that t u belongs to L p,t;w,p ;R N can be directly read from.. More precisely, we introduce the following definition of a weak solution to..6. Definition.. Let p n n+,, A satisfy..5, f Lp,T;L p ;R n N, and let u,g be an p,-suitable couple of boundary and initial data with an extension g. We say that u is a weak solution to., if there exists v such that u = v +g and.9. v X p, C[,T];L ;R N, lim vt, L t + =, and for almost all t,t and all ϕ W,p ;R N there holds. vt,ϕ + At,x, vt,x ϕxdx = ft,x ϕxdx. The duality vt,ϕ is understood in W,p,R N. We will also freuently use the euivalent reformulation of., which takes the form. vt,ϕ + At, x, vt, x+ gt, x ϕx dx = gt, ϕ + ft,x ϕxdx. Finally, to justify the setting of p,-suitable data, let us just discuss the case when we want to have u L p for some >. In view of Definition., it is natural to assume that also g and also v belong to L p and for the right hand side one should reuire f L p. Then due to.5, we also have that Ax, u L p and conseuently looking to., we see that also t u should belong to L p,t;w,p ;R N. Since u is given as a sum of g and v, it is uite obvious where the compatibility condition on g comes from. In case that the boundary of is smooth enough, one can trace the precise assumption on u,g for example by using the sharp theory for the heat euation together with the trace theorems for Sobolev-Bochner functions. Since this topic for a rough boundary goes beyond the aim of the paper and in fact is more related to the function spaces theory than to theory of PDEs, we do not investigate it in more details and rather formulate all assumptions in terms of the existence of a suitable extension g, which is included in Definition.... Assumptions on the non-linearity. Next assumptions that need to be imposed are related to the boundary and also to the structure of nonlinearity A. In case we are interested only in existence of a weak solution in the sense of Definition. for an operator A fulfilling..6, one can apply the standard Minty method see [3] to observe that for any f L p,t;l p ;R n N and any suitable L p -couple u,g there exists a function u X p,, which is a weak solution to. in the sense of Definition.. Moreover, for fixed extension g of a couple u,g the weak solution u is uniue. However, in case we want to deal with higher integrability results up to the boundary of,t, we need to impose more assumptions on A and also on the domain. Concerning the structural assumptions on A, it is well documented with counterexamples in various

5 HIGHER INTEGRABILITY FOR PARABOLIC SYSTEMS 5 dimensions [38] that in case of systems the above assumptions..6 are not enough to show higher integrability even locally in,t. The results in the higher integrability theory for systems are basically known for the model p-laplace case:.3 Ax,η = η p η, or its extension to radially symmetric case, see [36] for results in the elliptic setting and also [] for the parabolic one,.4 Ax,η = ax, η η or the most generally, to operators A s behaving asymptotically when η as.4, see the results [, 8, 9, 8] for the elliptic setting or [7, 39, 3, 7] for parabolic one. Inspired by these results we will reuire in what follows that A satisfies the following asymptotically Uhlenbeck condition. Definition.3. A nonlinearity At, x, η is said to be asymptotically Uhlenbeck if there exists a bounded nonnegative function Φ : [, [, such that.5 lim ρ Φρ = and.6 At,x,η η p η Φ η + η p for almost every t,x R R n and all η R n N. We remark that if At, x, η is asymptotically Uhlenbeck, then At,x,η η p η lim η η p = uniformly with respect to t,x R R n. Also notice that.5 and.6 directly imply that there exists a constant c > such that for almost all t,x,t and all η,η R n N we have that.7 At,x,η At,x,η c + η + η p η η. Finally, we describe our assumption on the domain. Here, we are inspired by the result for Reifenberg flat domains in scalar case, see [,, 5, 3, 33]. Hence, we will consider in the paper that the domain is δ,r -Reifenberg flat domain. More precisely: Definition.4. We say that is δ,r -Reifenberg flat if for each x and for each r,r ], there exists a coordinate system {y,,y n } such that x = in this coordinate system and that B r {y n > δr} B r B r {y n > δr}. Having collected all necessary assumptions, we can state our main theorem. Theorem.5. Let p n n+,, N N and A satisfy..6 and.5.6. Then for every [,, there exist a small constant δ > depending on n, N, p, and A such that for any [, ] and δ,δ ], we have the following properties: if R n is a bounded δ,r -Reifenberg flat domain for some r >, i.e., it satisfies Definition.4, if f L p,t;l p ;R n N and if u,g is a suitable p,-couple with an extension g X p,, then the solution u to. belongs to X p, and the following estimate holds.8 +p +p p u X p, c + g X p, + f + g p L p,t;l p ;R n N X + f, p, L p,t;l p ;R n N with a positive constant c depending only on n, N, p,, r,, T and A.

6 6 M. BULÍČEK, S. BYUN, P. KAPLICKÝ, J. OH, AND S. SCHWARZACHER Remark.6 p. Clearly, since in case p, we have + p and so.8 can be written in form without a scaling deficit: p.9 u X p, c + g X p, + f. L p,t;l p ;R n N We would like to emphasize that the nonlinearity A can go beyond Uhlenbeck category. The nonlinearity A under consideration behaves like the p-laplacian operator only near the infinity with respect to the gradient variable. In addition, it is worth pointing out that the domain under consideration is the so-called Reifenberg flat domain whose boundary is so rough that even the unit normal vector cannot be well defined. The rest of the paper is devoted to the proof of Theorem.5. First, in Section we state and prove an elliptic version of Theorem.5. In Section 3 we focus on the homogeneous case, i.e., the case when u,g =, but f. Then in Section 4 we apply the result for homogeneous case to the inhomogeneous one and prove the main theorem.. Elliptic systems Before proving the main theorem, we first prove elliptic estimates for the sake of better readability and to introduce the strategy of the proof. Let us consider the following elliptic system: { divax, u = divf in,. u = g on, whereg W,p is a suitableextension of the given boundaryvalues. Here the nonlinearity A : R n R n N R n N is a Carathéodory mapping, satisfying the p-coercivity, the p - growth and the monotonicity conditions as follows: there exist positive constants c and c such that for almost all x R n and all η,η R n N,..3.4 c η p c Ax,η η Ax,η c + η p p-coercivity, p -growth, Ax,η Ax,η η η monotonicity. We also assume that A is asymptotically Uhlenbeck, that is, there exists a bounded nonnegative function Φ : [, [, such that.5 lim ρ Φρ = and.6 Ax,η η p η Φ η + η p for almost every x R n and all η R n N. The following theorem is the main result in this section. Theorem.. Let p,, N N and A satisfy..3,..4 and.5.6. Then for every [,, there exist a small constant δ > depending on n, N, p, and A such that for any [, ] and δ,δ ], we have the following properties: if R n is a bounded δ,r -Reifenberg flat domain for some r >, i.e., it satisfies Definition.4, if f L p ;R n N and if g W,p ;R N, then the solution u to. belongs to W,p ;R N and the following estimate holds u L p ;R n N c + g L p ;R n N + f p, L p ;R n N with a positive constant c depending only on n, N, p, and on the assumptions on A and.

7 HIGHER INTEGRABILITY FOR PARABOLIC SYSTEMS 7 To prove the above theorem, we shall consider the system with homogeneous boundary data { divax, v = divf in,.7 v = on, and the p-laplacian system with homogeneous boundary data { div w p w = divf in,.8 w = on. For the p-laplacian system.8, we have the following Calderón-Zygmund estimate, see [7] and [3, Remark 4.5]. Lemma.. Let w W,p ;R N be the solution of.8. Then for every [,, there exist a small constant δ > depending on n, N, p and such that for any [, ] and δ,δ ], we have the following properties: if R n is a bounded δ,r -Reifenberg flat domain for some r > and if f L p ;R n N, then w W,p ;R N with the estimate w L p ;R n N c f p.9 L p ;R n N for some positive constant c depending only on n, N, p,, r and the assumptions on. For the asymptotically Uhlenbeck system.7 with a sufficiently smooth domain, we have the following ualitative estimate. This result is an elliptic version of the paper [3], which we repeat here to introduce our new strategy for a proof. Lemma.3. Let p,, N N and A satisfy..3,..4 and.5.6. Suppose that is a C -domain, and let v be a solution to.7. If f L p ;R n N, then we have v W,p ;R n. We can generalize this Calderón-Zygmund estimate for the weak solution to the system.7 also for Reifenberg flat domains as follows. Lemma.4. Let p,, N N and A satisfy..3,..4 and.5.6. For every [, there exist a small constant δ > depending on n, N, p, and A, such that for any [, ] and any δ,δ ] we have the following properties: if R n is a bounded δ,r -Reifenberg flat domain for some r > and if f L p ;R n N, then the solution v to.7 belongs to W,p ;R N and the following estimate holds v L p ;R n N + f c p. L p ;R n N with a positive constant c depending only on n, N, p,, r, and A. Proof. Wefirstprovetheglobalestimate.undertheaprioriassumptionv W,p ;R N. We rewrite the problem.7 as { div v p v = div f Ax, v+ v p v in, v = Since the nonlinearity A has p -growth, we have Ax, v c + v p, and hence the a priori assumption v W,p ;R N yields that Applying Lemma., we get f Ax, v+ v p v L p ;R n N. v L p ;R n N c f Ax, v+ v p v p L p ;R n N on.

8 8 M. BULÍČEK, S. BYUN, P. KAPLICKÝ, J. OH, AND S. SCHWARZACHER p c f L p ;R n N + Ax, v v p v p L p ;R n N where c is a positive constant depending on n, N, p,, r and the the assumptions on. From.5, we see that for any ε, there exists M = Mε > such that. Φρ ε, ρ M. Then.6 and. give Ax, v v p v Φ v + v p = Φ v + v p χ { v <M} +Φ v + v p χ { v M} Φ L +M p +ε+ v p Therefore, it follows that v L p ;R n N c p f +c L p ;R n N + Φ L +M p +ε v p. +c ε p v p p = c f +c L p ;R n N + Φ L +M p p p L p ;R n N +c ε p v L p ;R n N. + Φ L +M p p If we take ε, so small that c ε p, then we finally obtain the estimate 3. under the a priori assumption v W,p ;R N. Let us complete the proof of. by removing the assumption v W,p ;R N. We shall use an approximation argument as follows. From Lemma 4. in [6] see also Section 5 in [4] and a standard approximation of a Lipschitz domain by smooth domains, there exists a seuence of smooth domains { k } k= with the uniform δ,r -Reifenberg flatness property such that k k+ for all k N and. d H k, as k, where d H is the Hausdorff distance. Moreover, the approximation is achieved such that the constant in.9 is independent of k. We extend f L p ;R n N to zero outside and consider a mollifying kernel ϕ C B satisfying B ϕdx =. For ε >, we set ϕ ε x = x ε nϕ ε and define the mollifications f ε x := f ϕ ε x ϕ ε yfx ydy. R n Thenitisclearthatf ε C R n ;R n N andthatf ε f inl p ;R n N see[3, Appendix C]. Furthermore, for each k N, there exists < ε k < d H k, such that.3 f εk L p k ;R n N f L p ;R n N. We now consider the uniue weak solution v k to the problem { divax, vk = divf εk in k,.4 v k = on k.,

9 HIGHER INTEGRABILITY FOR PARABOLIC SYSTEMS 9 By Lemma.3, we know that v k L p k ;R n N ualitatively. Therefore, it follows from above and.3 that v k L p k ;R n N + f c p εk L p k ;R n N p c + f, L p ;R n N where the constant c is independent of k due to the uniform Reifenberg flatness of the domain approximation. We next let ṽ k be the zero extension of v k from k to, that is, { vk x if x ṽ k x = k, if x \ k. Then we obtain that ṽ k W,p ;R N with the estimate.5 ṽ k L p ;R n N = v k L p k ;R n N c + f p L p ;R n N Since {ṽ k } k= is uniformly bounded in W,p ;R N, there exist a subseuence, which we still denote by {ṽ k } k=, and a function ṽ W,p ;R N such that { ṽk ṽ weakly in L.6 p ;R n N ṽ k ṽ strongly in L p ;R N as k. Then we conclude from.,.5 and.6 that ṽ is a solution to.7 with the estimate ṽ L p ;R n N + f c p. L p ;R n N Finally we see from the uniueness of the weak solution to the problem.7 that ṽ = v, and the proof is complete. The proof of Theorem.. Assume that we fixed, δ and r according to Theorem.. Next, we consider arbitrary given [, ], f L p ;R n N and g W,p ;R N. Let v := u g. From the trivial ineuality u L p ;R n N v L p ;R n N + g L p ;R n N, it is enough to prove that.7 v L p ;R n N c + g L p ;R n N + f p L p ;R n N with a constant c being independent of f and g. Since u is given as u = v +g, we can rewrite the problem. as { divax, v = divax, v Ax, u+f in,.8 v = on, Since v is a solution with zero boundary data, we can apply Lemma.4 to obtain that for arbitrary [,] [, ] v p dx c + Ax, v Ax, u+f p dx.9 c + Ax, v Ax, u p dx+ f p dx whenever the right hand side is finite. The last term is bounded for any, so we shall justify that also the first integral on the right hand side of.9 is bounded and that it can.

10 M. BULÍČEK, S. BYUN, P. KAPLICKÝ, J. OH, AND S. SCHWARZACHER be estimated uniformly. Let us observe the following simple algebraic ineuality, which is a direct conseuence of the elliptic version of.7. Ax, v Ax, u c u + v p u v +c c u v + v p g +c = c v + g p g +c with a constant c possibly varying line to line but being independent of u, v and g. First, we start with the case p,]. It follows from. that Ax, v Ax, u c g p +c. Hence, substituting this ineuality into.9 and using also the fact that is bounded, we obtain v p dx c + g p dx+ f p dx and.7 then directly follows. Thus, it remains to discuss the case p,. Formally, i.e., in case we would know that v W,p, we could use the Young ineuality in. to observe that. Ax, v Ax, u p c v p p p g p +c g p +c ε v p +cε g p +c. This ineuality, used in.9 together with the fact that p and so p p, then leads to. v p dx cε v p dx+cε + g p dx+ f p dx. Hence choosing ε > sufficiently small, and assuming that the right hand side is finite, i.e., v W,p, we can use the above ineuality and absorb the term involving v to the left hand side, to obtain.3 v p dx c + g p dx+ f p dx. However, to make this procedure rigorous, we would need to know a priori that v W,p, which is not the case here. Therefore, we proceed more carefully. First, we define := p +p and set = in.9. Since +p + p =, +p we can deduce from. with the help of the Young ineuality that using also the fact that Ax, v Ax, u p c v p p p g p p +c g p +c p p = c v +p g p c+ v p + g p. +p +c g p +c Using this ineuality in.9 with := and the fact p p since p, we have.4 v p dx c + v p dx+ g p dx+ f p dx <.

11 HIGHER INTEGRABILITY FOR PARABOLIC SYSTEMS Since v p is integrable, we see that we improve the integrability of v to L p. Hence, we can now use., which has the right hand side finite and conclude the a priori estimate.3 rigorously. Next, we continue iteratively. We find for k := p k k +p that k < k+ < and k as k. Expressing now the corresponding dual exponents as we deduce that k +p + k v p p k p g p k p = v p p k k +p g p k +p =, p k k +p v p k + g p. Hence, combining this ineuality with. yields that using k again Ax, v Ax, u p k c v p p k g p k + g p k + c+ v p k + g p. Since v L p k, we can use the above ineuality in.9 to conclude improved integrability result v L p k. Conseuently, we can now use. with := k, which has now the finite right hand side and to get the uniform estimate.3 with := k, i.e.,.5 v p k dx c + g p k dx+ f p k dx. Note that the constant c > in the previous estimate may be chosen independent of k. Since k as k, we can now let k in.5, which gives the desired conclusion.7. This completes the proof. 3. Parabolic Systems: Homogeneous boundary and initial data In the following two sections we will need the following notations for parabolic intrinsic cylinders: Q λ r = Qλ r t,x = t λ p r,t B r x. If p = we have the standard parabolic cylinders Q r = Q r t,x = t r,t B r x. In this section we consider the following system: t v divat,x, v = divf in,t, 3. v = on,t, v, = in. To obtain a result for the above system, we shall consider the following parabolic p-laplacian system: t w div w p w = divf in,t, 3. w = on,t, w, = in. It can be extended to the following system t w div w p w = divfχ [,T] in,, 3.3 w = on,, w = in,].

12 M. BULÍČEK, S. BYUN, P. KAPLICKÝ, J. OH, AND S. SCHWARZACHER Clearly, we find the following natural estimate: 3.4 wt, dx+ w p dxdt c sup <t< T f p dxdt, see for instance [, ]. From [, ] one can deduce the following intrinsic local estimate: For any <, there exists R = R n,n,p,, such that for any z = t,x [,T and any parabolic cylinder Q R z T,T with < R R and, 3.5 w p dz c Q R z T w p dz Q R z T p n+p n + Q R z T for some positive constant c depending only on n, N, p and, where { p if p, d = n if n+ < p <, is the scaling deficit constant. It implies immediately the following Lemma. [ ] f p + dz Lemma 3.. Let Q R be a parabolic cylinder and < R R. Assume that for some [, ], 3.6 w p χ dz + f p χ QT dz Kλ p. Then 3.7 Q λ R z Q λ R z w p χ dz cλ p, Q λ R z T where the constant c just depends on n, N, p,, and K. Proof. The proof is by the application of 3.5 to ws,y = wλ p R s+t,ry +x, Rλ which is on its domain of definition a solution to 3.8 t w div w p w = div f, where fs,y = fλ p R s+t,ry +x λ p. Indeed, we find that in this case on the reflected system 3.9 w dz p + f p dz K, Q Q which implies an estimate that is independent of λ, and hence the result follows. Remark 3.. The above scaling implies that if we have the following solution on the whole space { t w div w p w = divf in, R n, 3. we find the homogeneous estimate 3. w X p, c f w, = in R n, p. L p, ;L p R n ;R n N d,

13 HIGHER INTEGRABILITY FOR PARABOLIC SYSTEMS 3 The same estimate holds on the half space problem. Indeed for H + = R n {x n > }, the estimate holds for solutions to t w div w p w = divf in, H +, 3. w = on, H +, w, = in H +, as we can reflect these solutions to the full space. We now state the main theorem in this section. Theorem 3.3. Let p n n+, and N N. For every [,, there exist a small constant δ > depending on n, N, p and such that for any [, ] and δ,δ ], we have the following properties: if R n is a bounded δ,r -Reifenberg flat domain for some r >, i.e., it satisfies Definition.4 and if f L p,t;l p ;R n N, then the solution w to 3.3 belongs to X p, and the following estimate holds 3.3 +p p w X p, c f +c f p L p,t;l p ;R n N L p,t;l p ;R n N with a constant c depending only on n, N, p,, r, the assumptions on and T. In addition, we have T T T + p w p dz c f p dz +c w p 3.4 dz. Remark 3.4. The famous scaling deficit that is usually necessary in order to get homogeneous estimates, as is widely known in the frame work of evolutionary p-laplace, pops up in the theorem above in a remarkable form. The above theorem shifts the scaling deficit to the lower order term. The homogeneity of the estimate with respect to the right hand side is necessary in order to include inhomogeneous boundary data. Observe that the case p = reduces to the classic parabolic estimates for the heat euation. Moreover, in the case p n n+,], one finds that + p and hence T T T 3.5 w p dz c f p dz+c w dz p +c, which possesses no scaling deficit. Proof of Theorem 3.3. First of all, we may use teh fact that we can extend w as a global in time solution to f fχ QT with related uniform estimates. We start by fixing M and M such that Let us first assume that f L p [, = M and w L p [, = M. max{m,m}. Now we introduce the parabolic maximal operator as Mgt,x := sup r> t+r and we define the upper level sets for m > as O m = t r B rx gs,y dyds, { z R n+ : M w p χ + fχ Q T p m p z > m p}.

14 4 M. BULÍČEK, S. BYUN, P. KAPLICKÝ, J. OH, AND S. SCHWARZACHER Then the Calderón-Zygmund covering [6, ] implies the following: There are parabolic cylinders Q i O m such that Q i have finite intersection, meaning that there is N N such that all Q i can be placed into N classes so that the sets Q i in each class are disjoint. w p fχ QT p χ + dz cmp m p Q i and O m Q i. i We define r i as the radius of Q i. By the weak-l estimate, we have m p Q i m p 4n+4 O m c. i It follows that ri n+ c m p, which implies that r i,r if m is chosen large enough depending only on n,r. This is the point, where the Reifenberg flatness is used. Moreover, since χ QT f p dz + w p χ dz cm p, Q i Q i we find from Lemma 3. with λ = that w p χ dz cmm. Q i Using the above we estimate, we get w p dz This implies the following: 3.6 O m c [, m p m p m p cr. w p dz + O m c [, O m c [, O m c [, If max{m,m}, then O m [, w p dz + i w p dz +cm i w p χ dz +cm i w p dz cr. w p dz Q i Q T w p dz Q i m Q i erasem p For general M,M, we use the fact that we may rescale as follows. We know that is a solution to 3.7 where Fixing λ > so that ws,y = wλ p s,y λ t w div w p w = div f, fs,y = fλ p s,y λ p. max{m p λ p p,m p λ } =,

15 HIGHER INTEGRABILITY FOR PARABOLIC SYSTEMS 5 we obtain f p dyds = w p dyds = Then we discover from 3.6 that 3.8 cr w p dyds = f p λ p dyλ p dt M p λ p p, w p λ p dyλ p dt M p λ. w p λ p p dydt. If max{m p λ p p,m p λ } = M p λ p p =, then 3.8 gives 3.9 w p dydt cm p = c f p dz. If max{m p λ p p,m p λ } = M p λ =, then 3.8 yields 3. w p p+p p dydt cm = c w p dz +p p Combining 3.9 and 3., we obtain the desired estimate 3.4. Finally, 3.3 follows by the energy ineuality and the properties of the weak time derivative. For the system 3. with a sufficiently smooth domain, we have the following ualitative estimate. Lemma 3.5. Let p n n+,, N N and A satisfy..6 and.5.6. Suppose that is a C -domain, and let v be a solution to 3.. If f L p,t;l p ;R n N, then we have v X p,. Proof. Let us consider the extended system t ṽ divat,x, ṽ = div f in T,T, 3. ṽ = on T,T, ṽ T, = in. where ṽ, f are defined by { vt,x if < t T, ṽt,x = if T t, and ft,x = { ft,x if < t T, if T t. Then it follows from [3, Theorem.3] that ṽ X p, T + ε,t for any ε,t. Therefore, we conclude from the definition of ṽ that v X p,,t. For the system 3. with a nonsmooth domain, we have the following particular result. Theorem 3.6. Let p n n+,, N N and A satisfy..6 and.5.6. Then for every [,, there exists a small constant δ > depending on n, N, p, and A such that for any [, ] and any δ,δ ], we have the following properties: If R n is a bounded δ,r -Reifenberg flat domain for some r, and if f L p,t;l p ;R n N, then the uniue weak solution v to 3. belongs to X p, and the estimate +p p 3. v X p, c + f + f p L p,t;l p ;R n N L p,t;l p ;R n N holds with a constant c depending only on n, N, p,, r, T and the assumed properties of A and..

16 6 M. BULÍČEK, S. BYUN, P. KAPLICKÝ, J. OH, AND S. SCHWARZACHER Proof. The existence of a weak solution to3. follows from the theory of monotone operators or via Galerkin approximation; we refer the reader to [, 9, 35] and the references therein. To prove the uniueness, we assume that there exist two solutions v and v to the problem 3.. Notice that v t, v t, W,p ;R N for almost all t,t. Taking this function for a test function in the weak formulation., we obtain that for almost all t,t, v i t,v v t + At,x, v i t,x v v t,xdx = ft,x v v t,xdx, i =,. Subtracting each other and integrating over,t, we have T v v t,v v t dt T + At,x, v t,x At,x, v t,x v t,x v t,x dxdt =. Since v v t,v v t, it follows that T At,x, v t,x At,x, v t,x v t,x v t,x dxdt. The strict monotonicity of A leads to v = v almost everywhere in,t. Thus v v in,t, which is due to the zero trace. We now prove the global estimate 3. under the a priori assumption v X p,. We rewrite the problem 3. as t v div v p v = div v = f At,x, v+ v p v v, = in. Notice that.5 implies At,x, v c + v p, and hence the a priori assumption v X p, yields that Applying Theorem 3.3, we get T v p dxdt c 3.3 f At,x, v+ v p v L p,t;l p ;R n N. T T +c T c +c T in,t, f At,x, v+ v p v p dxdt + p v p dxdt At,x, v v p v p dxdt on,t, T p + f p dxdt+c v p dxdt From the structural conditions.4-.6, we obtain the following energy estimate of the problem 3.: T T 3.4 sup vt, dx+ v p dxdt c f p dxdt. <t<t

17 HIGHER INTEGRABILITY FOR PARABOLIC SYSTEMS 7 On the other hand, we see from.5 that for any ε, there exists M = Mε > such that 3.5 Φρ ε, ρ M. Then.6 and 3.5 give At,x, v v p v Φ v + v p = Φ v + v p χ { v <M} 3.6 +Φ v + v p χ { v M} Φ L +M p +ε+ v p Therefore, it follows from 3.3, 3.4 and 3.6 that T T T v p dxdt c f p dxdt+c + Φ L +M p +ε v p. + p f p dxdt +c T + Φ L +M p p +cε p T v p dxdt. If we take ε, so small that cε p, then we finally obtain the estimate 3. under the a priori assumption v X p,. Let us complete the proof of 3. by removing the assumption v X p,. We shall use an approximation argument as follows. From [6, Lemma 4.] with [4, Section 5] and a standard approximation of a Lipschitz domain by smooth domains, there exists a seuence of smooth domains { k } k= with the uniform δ,r -Reifenberg flatness property such that k k+ for all k N and 3.7 d H k, as k, where d H is the Hausdorff distance. Let v k be the uniue weak solution to the problem t v k divat,x, v k = divf in,t k, 3.8 v k = on,t k, v k, = in k. By Lemma 3.5, we know that v k X p,,t k. Therefore, it follows from the a priori global estimate that +p v k X p,,t k + f c p + f p L p,t;l p k ;R n N c +p p + f + f p L p,t;l p ;R n N L p,t;l p k ;R n N L p,t;l p ;R n N where the constant c is independent of k by the uniformity of the estimates and the construction of the approximation of the domain. Next, let ṽ k be the zero extension of v k from,t k to,t, that is, { vk t,x if t,x,t ṽ k t,x = k, if t,x,t \ k. Then we obtain that ṽ k X p,,t with the estimate, 3.9 ṽ k X p,,t = v k X p,,t k +p p c + f + f p L p,t;l p ;R n N L p,t;l p ;R n N.

18 8 M. BULÍČEK, S. BYUN, P. KAPLICKÝ, J. OH, AND S. SCHWARZACHER Hence {ṽ k } k= is uniformly bounded in Xp,,T. Also it follows from 3.8,.5 and 3.9 that {ṽ k t } k= is uniformly bounded in Lp,T;W,p ;R N. In view of Aubin-Lions Lemma, there exist a subseuence, which we still denote by {ṽ k } k=, and a function ṽ X p,,t such that ṽ k ṽ weakly in L p,t;w,p ;R N 3.3 ṽ k ṽ strongly in L p,t;l p ;R N ṽ k t ṽ t weakly in L p,t;w,p ;R N as k. Then we conclude from 3.9 and 3.3 that ṽ is a solution to 3. with the estimate +p ṽ X p,,t + f c p + f p. L p,t;l p ;R n N L p,t;l p ;R n N Finally we see from the uniueness of the weak solution to the problem 3. that ṽ = v, and the proof is complete. Remark 3.7. From Theorem 3.6, we see that the uniue weak solution v to 3. belongs to X p,. Hence, it follows from the proof of Theorem 3.6 that the estimate T T T p v p dz c + f p dz + v p dz also holds with a constant c depending only on n, N, p,, r,, T and A. 4. The proof of Theorem.5 for general data Assume that we fixed, δ and r according to Theorem.5. Next, we consider any [, ], f L p,t;l p ;R n N and a suitable L p -couple u,g. For this couple, we then can find an extension g X p, and, since, we also know that there exists a uniue for a given extension weak solution u X p, to., where u = v+g with v X p, fulfilling v, =. Using the obvious ineuality u X p, v X p, + g X p,, we see that to prove Theorem.5, it is enough to show that 4. v X p, c + g X p, + f +p + g X p, p L p,t;l p ;R n N +p p + f L p,t;l p ;R n N with a constant c being independent of f and g. Thus, in what follows, we focus only on the proof of 4.. First of all, we represent t g as a divergence of a vector-valued function whose norm is controlled in an appropriate space. To do so, we look for a uniue solution to the following problem: find w : R N solving { div wt p w = t gt in, w = on. Here, we can use the standard monotone operator theory to obtain the existence of a uniue w W,p ;R N solving the above problem. Moreover, we know that it continuously depends on t g, then it is also Bochner measurable and we have an estimate wt p L p ;R n N c tgt p W,p ;R N.,

19 HIGHER INTEGRABILITY FOR PARABOLIC SYSTEMS 9 Thus, if we set G := w p w, then we have that t gt,ψ = Gt, ψ for almost every t [,T] and all ψ W,p ;R N. In addition, we have an estimate 4. T G p L p ;R n N dt c T t g p W,p ;R N dt. Next, since u is given as u = v +g, we can use the definition of G and rewrite the problem. as 4.3 t v divat,x, v = t g divat,x, v At,x, u+f = divat,x, v At,x, u+f G in,t, v = v, = in. on,t, Since v is a solution with zero initial and boundary data, we can apply Theorem 3.6 to conclude that for any [,] which is surely less or eual to, T v p dxdt { T c + At,x, v At,x, u+f G p dxdt 4.4 c { + T At,x, v At,x, u p dxdt + g p X p, + f p L p,t;l p ;R n N + T T p + } + v p dxdt + p } v p dxdt, whenever the right hand side is finite. Notice here that for the second ineuality, we used4., the Hölder ineuality and the fact that is bounded. Since v = u g, we get T T T v p dxdt u p dxdt+ g p dxdt 4.5 T T c f p dxdt+ g p dxdt, where we have used an energy estimate of the problem.. Combining 4.4 and 4.5 gives 4.6 T v p dxdt { T c + At,x, v At,x, u p dxdt + g p X p, + f p L p,t;l p ;R n N + g p+ p X p, } + f p + p p. L p,t;l p ;R n N Since the last four terms are bounded for any, it suffices to show that the first integral term in 4.6 is bounded in a universal way. We start with a simple algebraic ineuality,

20 M. BULÍČEK, S. BYUN, P. KAPLICKÝ, J. OH, AND S. SCHWARZACHER which is a direct conseuence of At,x, v At,x, u c u + v p u v +c c u v + v p g +c = c v + g p g +c, with a constant c possibly varying line to line but being independent of u, v and g. First, we consider the case that p n n+,]. Then it follows from 4.7 that At,x, v At,x, u c g p +c. Hence, substituting this ineuality into 4.6 and using also the fact that is bounded, we obtain T v p dxdt c + g p X + f p p, L p,t;l p ;R n N + g p+p + f p + p p X p,, L p,t;l p ;R n N and 4. then directly follows if for the estimate of t v, one uses 4.3 together with the growth assumption on A, see.5. We next discuss the case that p,. Formally, i.e., in case we would know that v X p,, we could use the Young ineuality in 4.7 to observe that 4.8 At,x, v At,x, u p c v p p p g p +c g p +c ε v p +cε g p +c. This ineuality, used in 4.6 together with the fact that p and so p p, then leads to T T v p dxdt cε v p dxdt 4.9 +cε + g p X p, + f p L p,t;l p ;R n N + g p+p X p, + f p + p p L p,t;l p ;R n N Hence choosing ε > sufficiently small, and assuming that the right hand side is finite, i.e., v X p,, we can use the above ineuality and absorb the term involving v to the left hand side to discover that 4. T v p dxdt c + g p X p, + f p L p,t;l p ;R n N + g p+p X p, We now remove the a priori assumption that v X p,. Write := p +p and we want to set = in 4.6. Observing that +p + p =, +p + f p + p p L p,t;l p ;R n N..

21 HIGHER INTEGRABILITY FOR PARABOLIC SYSTEMS it follows from 4.7 along with the help of the Young ineuality that using also the fact that At,x, v At,x, u p c v p p p g p p +c g p +c p p = c v +p g p c+ v p + g p. +p +c g p +c Applying this ineuality to 4.6 with := and recalling that p p since p, we find that T v p dxdt T c + v p dxdt+ g p X + f p 4. p, L p,t;l p ;R n N + g p+p + f p + p p X p, <. L p,t;l p ;R n N Because v p is integrable over T, we have improved the integrability of v to L p. Returning to 4.9, we see that the right hand side finite and deduce the a prior estimate 4.. We next continue this process iteratively for k := p k k +p such that k < k+ < and k as k. We write the corresponding dual exponents as k +p p + k k +p =, to see that v p p k p g p k p = v p p k k +p g p k k +p v p k + g p. We then combine this ineuality with 4.7 to observe that using k again At,x, v At,x, u p k c v p p k g p k + g p k + c+ v p k + g p. Since we have assumed that v L p k, we can now use the above ineuality in 4.6 to reach improved integrability result v L p k. Conseuently, we can now use 4.9 with := k that has now the finite right hand side, to further get the uniform estimate 4. with := k T v p k dxdt c + g p k X + f p k p, L p,t;l p ;R n N 4. k + g p+p X p, + f p + p pk L p,t;l p ;R n N We let k in 4. to deduce 4., as reuired. The proof is complete..

22 M. BULÍČEK, S. BYUN, P. KAPLICKÝ, J. OH, AND S. SCHWARZACHER References [] E. Acerbi and G. Mingione, Gradient estimates for the px-laplacean system, J. Reine Angew. Math , [], Gradient estimates for a class of parabolic systems, Duke Math. J. 36 7, no., [3] V. Bögelein, Global Calderón-Zygmund theory for nonlinear parabolic systems, Calc. Var. Partial Differential Euations 5 4, no. 3-4, [4] V. Bögelein, Global gradient bounds for the parabolic p-laplacian system, Proceedings of the London Mathematical Society 5, no. 3, [5] D. Breit, A. Cianchi, L. Diening, T. Kuusi, and S. Schwarzacher, The p-laplace system with right-hand side in divergence form: inner and up to the boundary pointwise estimates, Nonlinear Anal. 53 7,. MR [6] D. Breit, L. Diening, and S. Schwarzacher, Solenoidal Lipschitz truncation for parabolic PDEs, Math. Models Methods Appl. Sci. 53 3, no. 4, [7] D. Breit, B. Stroffolini, and A. Verde, Non-stationary flows of asymptotically Newtonian fluids, Commun. Contemp. Math. 8, no., 756, 6. [8] M. Bulíček, L. Diening, and S. Schwarzacher, Existence, uniueness and optimal regularity results for very weak solutions to nonlinear elliptic systems, Anal. PDE 9 6, 5 5. [9] M. Bulíček and S. Schwarzacher, Existence of very weak solutions to elliptic systems of p-laplacian type, Calc. Var. Partial Differential Euations [] S.-S. Byun and J. Oh, Global Morrey regularity for asymptotically regular elliptic euations, Appl. Math. Lett. 76 8, [] S.-S. Byun, J. Oh, and L. Wang, Global Calderón-Zygmund theory for asymptotically regular nonlinear elliptic and parabolic euations, Int. Math. Res. Not. IMRN 5, no. 7, [] S.-S. Byun and J. Ok, Nonlinear parabolic euations with variable exponent growth in nonsmooth domains, SIAM J. Math. Anal. 48 6, no. 5, [3], On W, -estimates for elliptic euations of px-laplacian type, J. Math. Pures Appl , no. 3, [4] S.-S. Byun, J. Ok, and S. Ryu, Global gradient estimates for general nonlinear parabolic euations in nonsmooth domains, J. Differential Euations 54 3, no., [5] S.-S. Byun and D.K. Palagachev, Morrey regularity of solutions to uasilinear elliptic euations over Reifenberg flat domains, Calc. Var. Partial Differential Euations 49 4, no. -, [6] S.-S. Byun and L. Wang, Parabolic euations in time dependent Reifenberg domains, Adv. Math. 7, no., [7] S.-S. Byun, L. Wang, and S. Zhou, Nonlinear elliptic euations with BMO coefficients in Reifenberg domains, J. Funct. Anal. 5 7, no., [8] L. A. Caffarelli and I. Peral, On W,p estimates for elliptic euations in divergence form, Comm. Pure Appl. Math , no.,. [9] M. Colombo and G. Mingione, Caldern-zygmund estimates and non-uniformly elliptic operators, J. Funct. Anal. 7 7, no. 4, [] E. DiBenedetto, Degenerate parabolic euations, Springer-Verlag, New York, 993. [] E. DiBenedetto and A. Friedman, Hölder estimates for nonlinear degenerate parabolic systems, J. Reine Angew. Math ,. [] L. Diening, M. Růžička, and J. Wolf, Existence of weak solutions for unsteady motions of generalized Newtonian fluids, Ann. Sc. Norm. Super. Pisa Cl. Sci. 5 9, no., 46. [3] L.C. Evans, Partial differential euations, second ed., Graduate Studies in Mathematics, vol. 9, American Mathematical Society, Providence, RI,. [4] T. Iwaniec, On L p -integrability in pde s and uasiregular mappings for large exponents, Annales Academiæ Scientiarum Fennicæ, Series A. I. Mathematica 7 98, 3 3. [5] T. Iwaniec, Projections onto gradient fields and L p -estimates for degenerated elliptic operators, Studia Math , no. 3, [6] J. Kinnunen and S. Zhou, A boundary estimate for nonlinear euations with discontinuous coefficients, Differential Integral Euations 4, no. 4, [7] T. Kuusi and G. Mingione, New perturbation methods for nonlinear parabolic problems, J. Math. Pures Appl. 9 98, no. 4, [8] G.M. Lieberman, Boundary regularity for solutions of degenerate elliptic euations, Nonlinear Anal. 988, no., 3 9. [9] J.-L. Lions, Quelues méthodes de résolution des problèmes aux limites non linéaires, Dunod; Gauthier- Villars, Paris, 969. [3] T. Mengesha and N.C. Phuc, Global estimates for uasilinear elliptic euations on Reifenberg flat domains, Arch. Ration. Mech. Anal. 3, no., 89 6.

23 HIGHER INTEGRABILITY FOR PARABOLIC SYSTEMS 3 [3] Giuseppe Mingione, The caldern-zygmund theory for elliptic problems with measure data, Ann. Sc. Norm. Super. Pisa Cl. Sci , no., [3] G. J. Minty, On a monotonicity method for the solution of non-linear euations in Banach spaces, Proc. Nat. Acad. Sci. U.S.A , [33] D.K. Palagachev and L.G. Softova, The Calderón-Zygmund property for uasilinear divergence form euations over Reifenberg flat domains, Nonlinear Anal. 74, no. 5, [34] Nguyen Cong Phuc, Nonlinear muckenhoupt-wheeden type bounds on reifenberg flat domains, with applications to uasilinear riccati type euations, Adv. Math. 5 4, [35] R.E. Showalter, Monotone operators in Banach space and nonlinear partial differential euations, Mathematical Surveys and Monographs, vol. 49, American Mathematical Society, Providence, RI, 997. [36] K. Uhlenbeck, Regularity for a class of non-linear elliptic systems, Acta Math , no. 3-4, 9 4. [37] N. N. Ural ceva, Degenerate uasilinear elliptic systems, Zap. Naučn. Sem. Leningrad. Otdel. Mat. Inst. Steklov. LOMI 7 968, 84. [38] V. Šverák and X. Yan, Non-Lipschitz minimizers of smooth uniformly convex functionals, Proc. Natl. Acad. Sci. USA 99, no. 4, [39] J. Zhang and S. Zheng, Lorentz estimate for nonlinear parabolic obstacle problems with asymptotically regular nonlinearities, Nonlinear Anal. 34 6, Mathematical Institute, Faculty of Mathematics and Physics, Charles University, Sokolovská 83, Prague, Czech Republic address: mbul86@karlin.mff.cuni.cz Department of Mathematical Sciences and Research Institute of Mathematics, Seoul National University, Seoul 886, Korea address: byun@snu.ac.kr Department of Mathematical Analysis, Faculty of Mathematics and Physics, Charles University, Sokolovská 83, Prague, Czech Republic address: kaplicky@karlin.mff.cuni.cz Fakultät für Mathematik, Universität Bielefeld, Postfach 3, D-335 Bielefeld, Germany address: joh@math.uni-bielefeld.de Institute of applied mathematics, University of Bonn, Ende nicher Allee 6, 535 Bonn, Germany & Department of Analysis, Faculty of Mathematics and Physics, Charles University, Sokolovská 83, Prague, Czech Republic address: schwarz@karlin.mff.cuni.cz