arXiv:0810.4716v2 [math.AP] 23 Mar 2009 n o-iegnefr aaoi n litceutoswt edn c leading with equations elliptic and parabolic form non-divergence and e n h oomldrvtv rbe a salse n[]ad[] I [6]. and [5] in established was problem coefficients derivative D the consta the papers conormal Lipschitz of solvability the small the and a domain, flat with lem Reifenberg domain so-called Lipschitz a in a equations in equations elliptic form litceutostesrnetrslsu odt a efudi By in found be [17]. can Krylov date and to [6], up [5], Wang results and strongest the equations elliptic sdb .Sfnv eoeta,teslaiiyfrteDrcltan Dirichlet the we for coefficients solvability VMO for with [8] the equations in that, elliptic form com Before divergence weak lemm of the covering problems Safonov. are Vitali [6] M. the [5], by and [3], used in function, tools maximal main Hardy-Littlewood The the included. not are terms der umto ovnin vrrpae nie r nocd eas We and enforced. measurable, are and indices bounded pa are the repeated coefficients Throughout over coefficients. conventions leading rough summation with spaces Sobolev in M offiins ooe pcs ie norms. mixed spaces, Sobolev coefficients, BMO aiso rw nvriy S rn ubrDS0367fr DMS-0635607 number grant NSF DMS-0800129. University, Brown of matics hr aebe ayrsac ciiisi hsdrcin o diver For direction. this in activities research many been have There esuyteslaiiyo litcoeaosi iegneform divergence in operators elliptic of solvability the study We 1991 n[] the [3], In n[7,Kyo aeauie praho the of approach unified a gave Krylov [17], In .Dn a atal upre yasatu udn rmth from funding start-up a by supported partially was Dong H. e od n phrases. and words Key LITCEUTOSI IEGNEFR WITH FORM DIVERGENCE IN EQUATIONS ELLIPTIC ahmtc ujc Classification. Subject Mathematics isht oan o qain ntewoesaeo afs half a or space whole sp the half coefficients in a equations space, For whole domain. the Lipschitz in equations elliptic order second M einrsi h te ietos o qain nabo a in that equations assume For we directions. additionally other the in semi-norms BMO netgt litceutosi ooe pcswt mixe coefficients. with the spaces on o pro Sobolev approach assumptions derivative unified same in conormal a equations the give elliptic We and investigate domain. problem the boundary of Dirichlet boundary the of hood Abstract. C 1 , 1 W oan,adi 1 for [1] in and domains, p 1 ovblt a bandfrteDrcltpolmo divergence of problem Dirichlet the for obtained was solvability a h ovblt nSblvsae spoe o iegnef divergence for proved is spaces Sobolev in solvability The ij ATAL M COEFFICIENTS BMO PARTIALLY r sue ob esrbei n ieto n aesmall have and direction one in measurable be to assumed are a ij L OGI OGADDYO KIM DOYOON AND DONG HONGJIE eododreutos aihn enoclain part oscillation, mean vanishing equations, Second-order u r sue ohv ml M einrsadlwror- lower and semi-norms BMO small have to assumed are ( = a ij 1. a u ij x Introduction i aesalBOsm-om naneighbor- a in semi-norms BMO small have + C 51,3J5 50,35J25. 35R05, 35J15, 35K15, a 1 j 1 domains. u ) x j a + ij b L i r nfrl elliptic. uniformly are u p x i ovblt fbt divergence both of solvability + mIS n S rn number grant NSF and IAS, om cu iiino ple Mathe- Applied of Division e c,adabounded a and ace, om ne the under norms d ae h leading the pace, lm ealso We blem. ne domain, unded ohthe both f e,teusual the per, uealthe all sume rcltprob- irichlet n[] Byun [3], un Neumann d orm originally a eobtained re ec form gence oefficients pactness, al small ially those n t For nt. (1.1) 2 H. DONG AND D. KIM
VMO in the spatial variables (and measurable in the time variable in the para- bolic case). Unlike the arguments in [7], [8] and [11], which are based on certain estimates of Calder´on-Zygmund theorem and the Coifman-Rochberg-Weiss com- mutator theorem, the proofs in [17] rely mainly on pointwise estimates of sharp functions of spatial derivatives of solutions. It is worth noting that although the results in [17] are stated for equations with VMO coefficients, the proofs there only require aij to have locally small BMO semi-norms. We also remark that for di- vergence form parabolic equations a similar result was also obtained in Byun [4] by adapting the approach in [3]. Krylov’s method was later improved and gen- eralized in [9], [13]-[16], [18] and [19]. With the leading coefficients in the same class, Krylov [18] established the solvability of both divergence and non-divergence parabolic equations in mixed-norm Sobolev spaces. There are many other results in the literature regarding the Lp theory of second order parabolic and elliptic equations with discontinuous coefficients. For non- divergence form equations, we refer the reader to [2], [7], [20], [25], [27] and refer- ences therein. For divergence form equations, see also [26] and references therein. The theory of elliptic and parabolic equations with partially VMO coefficients is 2 originated in Kim and Krylov [13], where the authors proved the Wp solvability of elliptic equations in non-divergence form with leading coefficients measurable in a fixed direction and VMO in the others. Very recently, their result was generalized by Krylov [19], where the leading coefficients are assumed to be measurable in one direction and VMO in the orthogonal directions in each small ball with the direction 1,2 depending on the ball. For non-divergence parabolic equations, the Wq,p solvability was established in Kim [16], in which most leading coefficients are measurable in time variable as well as one spatial variable, and VMO in the other variables. We remark that to our best knowledge, at the time of this writing, all known results concerning Lp solvability of elliptic and parabolic equations with partially VMO/BMO coefficients are only for non-divergence form. In this paper we consider divergence form elliptic equations in the whole space, a half space and a Lipschitz domain with a small Lipschitz constant. We deal with equations with partially BMO leading coefficients with locally small BMO semi-norms (Theorem 2.2), a class of coefficients which is more general than those treated previously in [17], [3], [5] and [6]. More precisely, we assume the coefficients aij are measurable in x1 direction and BMO in the other directions with locally small BMO semi-norms (see Assumption 2.1 for a more rigorous definition). This is the same class of coefficients considered in [13], in which non-divergence form elliptic equations are studied. For equations in a Lipschitz domain, additionally we assume that aij have small BMO semi-norms in a neighborhood of the boundary 1 of the domain. Under these assumptions, we establish the unique Wp solvability of divergence form elliptic equations. We give a unified approach of both the Dirichlet boundary problem and the conormal derivative problem in a half space (Theorem 2.3, 2.4) and in a bounded Lipschitz domain (Theorem 2.7 and 2.8). We also investigate elliptic equations in Sobolev spaces with mixed norms under the same assumption on the coefficients. We point out that, as in [17] and [18], one feature of these results is that the matrix aij is not assumed to be symmetric. One of the motivations of the{ paper} is the following problem. Consider the ij i equation a ux xj = div g in B2, the ball of radius 2 centered at the origin, with zero Dirichlet boundary condition. The coefficients aij are assumed to be bounded, 3 uniformly elliptic and piecewise uniformly continuous on B1 and B2 B1. This is 1 \ a very natural problem and the W2 solvability of it follows immediately from the 1 Lax-Milgram lemma. However, for the Wp solvability when p = 2, it seems to us that none of the results above are applicable in this case. We6 will give a solution to the problem at the end of Section 2 as an application of our main results. Our approach is based on the aforementioned method from [17]. However, since aij are merely measurable in x1, we are only able to estimate the sharp function of ux′ , not the full gradient ux as in [17]. Here and throughout the paper, we denote 2 d d 1 x′ = (x , , x ) R − , so by u we mean one of u i , i =2, , d, or the whole ··· ∈ x′ x ··· 1 collection of them. Roughly speaking, the main difficulty is to bound ux by ux′ . One idea in the paper is to break the ‘symmetry’ of the coordinates so that x1 is 11 distinguished from x′. Another idea is to estimate the sharp of a ux1 instead of ux1 . This estimate together with a generalized Stein-Fefferman theorem proved in [19] enables us to bound ux1 . The main advantage of the approach is that here we can obtain the boundary estimate immediately from the estimate in the whole space since the leading coefficients are allowed to be just measurable in one direction. In a forthcoming paper, we will extend our results to systems with variably partially BMO coefficients. A brief outline of the paper: in the next section, we introduce the notation and state the main results, Theorem 2.2, 2.3, 2.4, 2.7 and 2.8. Section 3 is devoted to several auxiliary results which will be used later, in which we estimate the Lp norm 1 ux by the Lp norm of ux′ (Theorem 3.7). Then in Section 4, we give an estimate of the sharp function of ux′ . By combining this with Theorem 3.7, we are able to prove Theorem 2.2 in Section 5. Theorem 2.3 and 2.4 are proved in Section 6, while Theorem 2.7 and 2.8 are proved in Section 7, Finally, the last four sections are devoted to the mixed norm estimate.
2. Main results Before we state our assumptions and main theorems, we introduce some neces- sary notations. By Rd we mean a d-dimensional Euclidean space and a point in Rd 1 d 1 is denoted by x = (x , , x ) = (x , x′). For given two positive integers d and d ··· 1 2 such that d1 + d2 = d, we set x = (x1, , xd1 ) Rd1 , x = (xd1+1, , xd) Rd2 . 1 ··· ∈ 2 ··· ∈ d That is, for example, x1 represents the first d1 coordinates of x R . If Ω is an open subset in Rd, we define ∈
q/p 1/q p u := u x2 x1 = u(x) I (x) dx x . Lq,p(Ω) Lq Lp (Ω) Ω 1 2 k k k k Rd2 Rd1 | | Z Z ! x2 x1 Note that, in case p = q, Lp(Ω) = Lp,p(Ω) = Lp Lp (Ω). Set Rd = x Rd : x = (x1, , xd), x1 > 0 . + { ∈ ··· } 1 A function u belongs to Wq,p(Ω) if u,ux Lq,p(Ω). Unless specified otherwise, by Rd ∈ 1 Lp we mean Lp( ). Similarly, whenever we use Lq,p, Wp , Lp,loc, Wp,loc, and C0∞, we understand that Rd is omitted. For a function f in Rd, we set 1 (f) = f(x) dx = – f(x) dx, Ω Ω | | ZΩ Z Ω 4 H. DONG AND D. KIM where Ω is the d-dimensional Lebesgue measure of Ω. Throughout| | the paper we assume that the coefficients aij , ai, bi, and c are bounded by a constant K 1. Moreover, we assume the uniform ellipticity condi- tion on aij , i.e., ≥ δ ξ 2 aij (x)ξiξj | | ≤ for all x and ξ Rd, where δ (0, 1). We need a very∈ mild regularity∈ assumption on the coefficients aij . To present this assumption, let
d d 1 B (x)= y R : x y < r , B′ (x′)= y′ R − : x′ y′ < r , r { ∈ | − | } r { ∈ | − | } 1 1 Γ (x) = (x r, x + r) B′ (x′). r − × r Set Br = Br(0), Br′ = Br′ (0), and Br′ is the d 1-dimensional volume of Br′ (0). Denote | | −
x1+r ij 1 ij 1 ij 1 1 oscx′ a , Γr(x) = – a (y ,y′) – a (y ,z′) dz′ dy′ dy , 2r x1 r B (x ) − B (x ) Z − Z r′ ′ Z r′ ′ where ij 1 1 ij 1 – a (y ,z′) dz′ = a (y ,z′) dz′. B′ Z Br′ (x′) | r| ZBr′ (x′) Then we set # ij aR = sup sup sup oscx′ a , Γr(x) . x Rd r R ij ∈ ≤ The following assumption contains a parameter γ > 0, which will be specified later.
# Assumption 2.1 (γ). There is a constant R0 (0, 1] such that a γ. ∈ R0 ≤ We state the main results concerning elliptic equations in divergence form in the 1 1 usual Sobolev spaces Wp . For equations in Sobolev spaces with mixed norms Wq,p, as indicated in the introduction, our results are presented in Section 9. Theorem 2.2 (Equations in the whole space). Let p (1, ) and f, g = (g , ,g ) ∈ ∞ 1 ··· d ∈ Lp. Then there exists a constant γ = γ(d,p,δ,K) such that, under Assumption 2.1 (γ), the following hold true. (i) For any u W 1 satisfying ∈ p u λu = div g + f, (2.1) L − we have λ u + √λ u N√λ g + N f , (2.2) k kLp k xkLp ≤ k kLp k kLp provided that λ λ0, where N and λ0 0 depend only on d, p, δ, K and R0. (ii) For any λ >≥ λ , there exists a unique≥ u W 1 satisfying (2.1). 0 ∈ p (iii) If ai = bi = c =0 and aij = aij (x1), i.e., measurable functions of x1 R only ∈ with no regularity assumptions, then one can take λ0 =0. The next two theorems are about the Dirichlet problem and the conormal de- rivative problem on a half space. 5
Theorem 2.3 (Dirichlet problem on a half space). Let p (1, ) and f, g = (g , ,g ) L (Rd ). Then there exists a constant γ = γ∈(d,p,δ,K∞ ) such that, 1 ··· d ∈ p + under Assumption 2.1 (γ), for any u W 1(Rd ) satisfying ∈ p + u λu = div g + f in Rd + , (2.3) Lu =0− on ∂Rd + we have √ √ λ ux L (Rd ) + λ u L (Rd ) N λ g L (Rd ) + N f L (Rd ), (2.4) k k p + k k p + ≤ k k p + k k p + provided that λ λ0, where N and λ0 0 depend only on d, p, δ, K, and R0. Moreover, for any≥ λ > λ and g,f L ≥(Rd ), there exists a unique u W 1(Rd ) 0 ∈ p + ∈ p + satisfying (2.3). Theorem 2.4 (Conormal derivative problem on a half space). Let p (1, ) and f, g = (g , ,g ) L (Rd ). Then there exists a constant γ = γ(d,p,δ,K∈ ∞) such 1 ··· d ∈ p + that, under Assumption 2.1 (γ), for any u W 1(Rd ) satisfying ∈ p + Rd u λu = div g + f in + Li1 − 1 d , (2.5) a u i + a u = g on ∂R x 1 + we have √ √ λ ux L (Rd ) + λ u L (Rd ) N λ g L (Rd ) + N f L (Rd ), (2.6) k k p + k k p + ≤ k k p + k k p + provided that λ λ0, where N and λ0 depend only on d, p, δ, K, and R0. Moreover, for any λ > λ ≥and g,f L (Rd ), there exists a unique u W 1(Rd ) satisfying 0 ∈ p + ∈ p + (2.5). Solutions of (2.5) are understood in the weak sense. More precisely, we say u W 1(Rd ) satisfies (2.5) if we have ∈ p + ij j i a uxi φxj a uφxj + b uxi φ + (c λ)uφ dx = ( gjφxj + fφ) dx Rd − − − Rd − Z + Z + (2.7) 1 d for any φ W (R ), where p′ satisfy 1/p +1/p′ = 1. For discussions about the ∈ p′ + conormal derivative problem, we refer the reader to [22] and [23]. Next we consider the solvability of divergence form elliptic equations in domains with the homogeneous Dirichlet boundary condition: u = div g + f in Ω . (2.8) Lu =0 on ∂Ω We shall impose a little bit more regularity assumption on aij near the boundary. For any x Rd, denote ∈ dist(x, ∂Ω) = inf x y . y ∂Ω | − | ∈ d Assumption 2.5 (γ). There is a constant R1 (0, 1] such that, for any x0 R with dist(x , ∂Ω) R and any r (0, R ], we have∈ ∈ 0 ≤ 1 ∈ 1 sup – aij (x) (aij ) dx γ. | − Br (x0)| ≤ ij Z Br (x0) 6 H. DONG AND D. KIM
We also impose the same assumption on domains as in [3], i.e. the boundary ∂Ω of the domain Ω is locally the graph of a Lipschitz continuous function with a small Lipschitz constant. More precisely, we make the following assumption containing a parameter θ (0, 1], which will be specified later. ∈ Assumption 2.6 (θ). There is a constant R2 (0, 1] such that, for any x0 ∂Ω ∈ d 1 ∈ and r (0, R ], there exists a Lipschitz function φ: R − R such that ∈ 2 → 1 Ω B (x )= x B (x ): x > φ(x′) ∩ r 0 { ∈ r 0 } and φ(y′) φ(x′) sup | − | θ x ,y B (x ),x =y y′ x′ ≤ ′ ′∈ r′ 0′ ′6 ′ | − | in some coordinate system. Note that all C1 domains satisfy this assumption for any θ > 0. Theorem 2.7 (Dirichlet problem on a bounded domain). Let p (1, ) and Ω be i ∈ ∞ a bounded domain. Assume axi + c 0 in Ω in the weak sense. Then there exist γ = γ(d,p,δ,K) and θ = θ(d,p,δ,K)≤such that, under Assumption 2.1 (γ), 2.5 (γ), and Assumption 2.6 (θ), for any f, g = (g1, ,gd) Lp(Ω) there exists a unique u W 1(Ω) satisfying (2.8). Moreover, we have··· ∈ ∈ p u 1 N f + N g , (2.9) k kWp (Ω) ≤ k kLp(Ω) k kLp(Ω) where N is independent of f,g and u. Our last result is about the solvability of divergence form elliptic equations in domains with the conormal derivative boundary condition: u = div g + f in Ω Lij j j j j , (2.10) a u i n + a un = g n on ∂Ω x j where n = (n1, ,nd) is the outward normal direction of ∂Ω, which is defined almost everywhere··· on ∂Ω. Like before, solutions of (2.10) are understood in the weak sense. More precisely, we say u W 1(Ω) satisfies (2.10) if we have ∈ p ij j i a u i φ j a uφ j + b u i φ + cuφ dx = ( g φ j + fφ) dx, (2.11) − x x − x x − j x ZΩ ZΩ for any φ W 1 (Ω). ∈ p′ Theorem 2.8 (Conormal derivative problem on a bounded domain). Let p i ∈ (1, ) and Ω be a bounded domain. Assume axi + c 0 in Ω in the weak sense. Then∞ there exist γ = γ(d,p,δ,K) and θ = θ(d,p,δ,K) such≤ that, under Assumption 2.1 (γ), 2.5 (γ), and Assumption 2.6 (θ), (i) If in the weak sense ai + c 0 in Ω and aini = 0 on ∂Ω, then for any xi ≡ f, g = (g1, ,gd) Lp(Ω), the equation (2.10) has a unique up to a constant solution u ···W 1(Ω) provided∈ that bi = c =0 and f dx =0. Moreover, we have ∈ p Ω ux L (Ω) N f L (Ω) + NR g L (Ω). k k p ≤ k k p k k p (ii) Otherwise, the solution is unique and we have
u W 1(Ω) N f L (Ω) + N g L (Ω). k k p ≤ k k p k k p The constant N is independent of f,g and u. 7
Here, by ai + c 0 in Ω, we mean xi ≤ i ( a φ i + cφ) dx 0 − x ≤ ZΩ 1 i i i for any nonnegative φ C (Ω). By a i + c 0 in Ω and a n =0 on ∂Ω, we mean ∈ 0 x ≡ i ( a φ i + cφ) dx =0. (2.12) − x ZΩ 1 for any φ W2 (Ω). Restricted∈ to equations in Lipschitz domains, Theorem 2.7 and 2.8 improve the previous results in [3], [5] and [6] in two aspects: first we only assume that the lead- ing coefficients have partially small BMO semi-norms in the interior of the domain; second we also allow lower order terms. At the time of this writing it is not clear to us whether our method can be extended to deal with equations in Reifenberg flat domains. We remark that for the Poisson equation in arbitrary Lipschitz domains but with a restricted range of p, the solvability result was established by Jerison and Kenig [12] (see also [26] for a generalization to equations with VMO coefficients). We end this section by giving an example dealing with elliptic equations with piecewise continuous leading coefficients on a bounded domain. This is another nice application of Theorem 2.2, showing the possibility that the results in this paper can be applied to many different equations with not necessarily continuous coefficients. For simplicity, consider ij a u i = div g in B , u =0, (2.13) x xj 2 |∂B2 ij where each a is piecewise uniformly continuous on B1 and B2 B1. As always, aij are assumed to be uniformly elliptic. For the solvability of the\ equation (2.13) 1 ij in Wp (B2), Theorem 2.7 is not applicable because the coefficients a do not have partially small BMO semi-norms in any fixed directions. However, upon having an appropriate partition of unity and change of variables, the interior estimate is derived from the Lp-estimate of equations with piecewise continuous coefficients. Here by ‘piecewise continuous coefficients’ we mean coefficients aij continuous on Rd Rd Rd + and on +. Needless to say, this class of coefficients satisfies the assumptions of Theorem 2.2.\ The interior and boundary estimates give us u N g + N u . k xkLp(B2) ≤ k kLp(B2) k kLp(B2) Then, for p > 2, one can use the argument in the proof of Theorem 2.7 below to absorb the term N u to the left-hand side. Thus we obtain an estimate as k kLp(B2) in Theorem 2.7. The estimate when p (1, 2) follows from the duality argument. Consequently, for a given g L (B ), 1 ∈
3. Auxiliary results for equations with measurable coefficients In this section we set ij u = a u i , L0 x xj and we do not impose any regularity assumptions on the coefficients of the operator 11 11 1 0, except a . The coefficient a is assumed to be a measurable function of x onlyL or satisfying 11,# Assumption 3.1 (γ). There is a constant R0 (0, 1] such that a γ. ∈ R0 ≤ 8 H. DONG AND D. KIM
Here γ > 0 is a constant to be specified, and 11,# 11 aR = sup sup oscx′ a , Γr(x) . x Rd r R ∈ ≤ However, in Theorem 3.2 all coefficients including a11 are measurable functions of x Rd with no regularity assumptions. ∈ The first result is the classical L2-estimate for elliptic operators in divergence form with measurable coefficients. Theorem 3.2. There exists N = N(d, δ) such that, for any λ 0, ≥ √λ u + λ u N √λ g + f , k xkL2 k kL2 ≤ k kL2 k kL2 provided that u W 1, f, g = (g , ,g ) L , and ∈ 2 1 ··· d ∈ 2 u λu = div g + f. (3.1) L0 − 1 Furthermore, for any λ> 0 and f, g L2, there exists a unique solution u W2 to the equation (3.1). ∈ ∈ Proof. We present a proof for the sake of completeness. Due to the method of continuity it is enough to prove the estimate. Moreover, by the denseness of C0∞ 1 in W2 it suffices to consider u C0∞. Then from the equation and the uniform ellipticity condition it follows that∈
2 2 ij 2 δ ux dx + λ u dx a uxi uxj dx + λ u dx Rd | | Rd | | ≤ Rd Rd | | Z Z Z Z
= giuxi dx fu dx Rd − Rd Z Z 2 2 2 N 2 δ/2 uxi dx + N g dx + λ/2 u dx + f dx, ≤ Rd | | Rd | | Rd | | λ Rd | | Z Z Z Z where N = N(d, δ). This finishes the proof.
If the above operator 0 is replaced by the Laplace operator ∆, it is well known that the result as in TheoremL 3.2 holds true not only for p = 2 but also for p ∈ (1, ). More precisely, if λ> 0 and f, g Lp, then there exists a unique solution u ∞W 1 to the equation ∆u λu = div g ∈+ f. As above, we have ∈ p − 1/2 u + √λ u N g + λ− f k xkLp k kLp ≤ k kLp k kLp for all λ> 0. Using this result, we prove the following theorem. 1 Theorem 3.3. Let p (1, ), λ> 0, κ 4, and r > 0. Assume that u W loc, ∈ ∞ ≥ ∈ p, f, g = (g1, ,gd) Lp,loc, and ∆u λu = div g + f in Bκr. Then there exists a constant N ···= N(d, p∈) such that −
p p p p/2 p d p p/2 p – u (u ) dx Nκ− u + λ u +Nκ g + λ− f . x x Br x B | − | ≤ | | | | Bκr | | | | Bκr Z r Proof. We follow the idea in the proof of Theorem 7.1 in [18] taking into account the presence of λ. We can certainly assume that u, f, and g have compact supports. In addition, we assume that u, f, and g are infinitely differentiable. Indeed, if not, we take the standard mollifications and prove the estimate for the mollifications. Then we take the limit because the concerned constants are independent of the smoothness of the functions involved. 9
Take a ζ C∞ such that ∈ 0 ζ =1 on B , ζ =0 on Rd B . κr/2 \ κr Then we find a unique solution w W 1 to the equation ∈ p (∆ λ) w = div(ζg)+ ζf. − Set v := u w and observe that − (∆ λ) v = div((1 ζ)g)+(1 ζ)f. − − − The classical theory on elliptic equations in divergence form indicates that w and v are infinitely differentiable. In addition, in Bκr/2, (∆ λ) v =0. − Rd+1 Then if we view v as a function in Cloc∞ ( ) independent of t, by Lemma 7.4 in [18] p p p p/2 p vx (vx) Nκ− vx + λ v , (3.2) Br Br | − | ≤ | | | | Bκr where N depends only on d and p. On the other hand, we have 1/2 w + √λ w N ζg + λ− ζf , k xkLp k kLp ≤ k kLp k kLp which implies p d p p/2 p d p p/2 p ( w ) Nr− ζg + λ− ζf Nκ g + λ− f , x Br Lp Lp | | ≤ k k k k ≤ | | | | Bκr p p/2 p p p/2 p ( w ) + λ ( w ) N g + λ− f . x Bκr Bκr | | | | ≤ | | | | Bκr From these inequalities as well as (3.2), we see that
p p p – ux (ux) dx N ( vx (vx)B ) + N ( wx ) | − Br | ≤ | − r | Br | | Br Z Br p p p/2 p d p p/2 p Nκ− vx + λ v + Nκ g + λ− f . ≤ | | | | Bκr | | | | Bκr We also have p p/2 p p p/2 p p p/2 p vx + λ v N ux + λ u + N wx + λ w | | | | Bκr ≤ | | | | Bκr | | | | Bκr p p/2 p p p/2 p N ux + λ u + N g + λ− f . ≤ | | | | Bκr | | | | Bκr Combining the above two sets of inequalities we arrive at the desired inequality in the theorem. We frequently make use of the following change of variables to ‘break’ the sym- metry of coordinates. Let u λu = div g + f L0 − in Rd. For a number µ 1, we set ≥ ij 1 ij 1 1 1 1 1 a¯ (x , x′)= a (µ− x , x′), u¯(x , x′)= u(µ− x , x′), (3.3) 1 1 1 1 1 1 f˜(x , x′)= f(µ− x , x′), g˜(x , x′) = (µg ,g , ,g )(µ− x , x′). (3.4) 1 2 ··· d Clearlyu ¯ satisfies 2 11 1j i1 ij µ a¯ u¯ 1 + µa¯ u¯ 1 + µa¯ u¯ i + a¯ u¯ i λu¯ x x1 x xj x x1 x xj − j>1 i>1 i,j>1 X X X 10 H. DONG AND D. KIM
= divg ˜ + f.˜ ¯ 11 d If we set 0w = (¯a wx1 )x1 + ∆d 1w, where ∆d 1w = wxixi , then L − − i=2 ¯ 2 ¯ 0u¯ µ− λu¯ = divg ¯ + f, P (3.5) L − where 2 2 1 i1 f¯ = µ− f,˜ g¯ = µ− g˜ µ− a¯ u¯ i , 1 1 − x i>1 X 2 1 1j 2 ij g¯ = µ− g˜ µ− a¯ u¯ 1 µ− a¯ u¯ i +¯u j j 2. j j − x − x x ≥ i>1 X We now assume that the coefficient a11 is a measurable function of x1 R. Under this condition on a11 (no regularity assumptions on aij if ij > 1) we prove∈ 11 an estimate fora ¯ u¯x1 . 1 11 11 1 Lemma 3.4. Let λ > 0, r > 0, κ > 8Kδ− , and a = a (x ). Assume that 1 u W loc and ∈ 2, u λu = div g + f, L0 − where f, g L loc. Then there exists a constant N = N(d, δ, K) such that ∈ 2, 1/2 11 11 2 1 d/2 1 2 1/2 a¯ u¯x1 a¯ u¯x1 N κ− + κ µ− u¯x1 Br Bκr | − | Br ≤ | | 1/ 2 d/2 2 2 2 1 ˜ 2 +Nκ u¯x′ + λ u¯ + g˜ + λ− f | | | | | | | | Bκr for all µ 1, where a¯ij , u¯, f˜, and g˜ are those defined in (3.3) and (3.4). In particular,≥ if λ = f =0, i.e., u = div g, we have L0 1/2 11 11 2 1 d/2 1 2 1/2 a¯ u¯x1 a¯ u¯x1 N κ− + κ µ− u¯x1 Br Bκr | − | Br ≤ | | d/2 2 2 1/2 +Nκ u¯x + g˜ | ′ | | | Bκr for all µ 1. ≥ Proof. The second inequality in the lemma follows easily from the first. Indeed, if we write u λu = div g λu, by the first inequality L0 − − 1/2 11 11 2 1 d/2 1 2 1/2 a¯ u¯x1 a¯ u¯x1 N κ− + κ µ− u¯x1 Br Bκr | − | Br ≤ | | d/2 2 2 2 1/2 +Nκ u¯x + λ u¯ + g˜ . | ′ | | | | | Bκr Then letting λ 0 gives the result. To prove theց first inequality in the lemma, recall thatu ¯ satisfies (see (3.5)) ¯ u¯ λu¯ = divg ¯ + f¯ , L0 − λ ¯ ¯ 2 1 where λ > 0 and fλ = f + (µ− 1)λu¯. Using Theorem 3.2 we find w W2 satisfying − ∈ ¯ w λw = div (I g¯)+ I f¯ , L0 − Bκr Bκr λ where I is the indicator function of a set Ω. Then v :=u ¯ w satisfies Ω − ¯ v λv = div ((1 I )¯g)+(1 I )f¯ . L0 − − Bκr − Bκr λ In particular, ¯ v λv = 0 in B . L0 − κr 11
Now we use the following change of variables. Set 1 x 1 y1 = φ(x1) := dr, yj = xj , j 2. a¯11(r) ≥ Z0 11 1 Since δ a¯ K, we readily see that the inverse φ− exists, φ is a bi-Lipschitz function,≤ and ≤ 1 1 1 K− φ(t)/t δ− , δ φ− (t)/t K (3.6) ≤ ≤ ≤ ≤ for t = 0. We define 6 1 1 1 v¯(y ,y′)= v(φ− (y ),y′). 1 1 We also define r = √2δ− r and κ = κ/(2Kδ− ). Using the fact that ¯ v λv =0 1 1 L0 − in Bκr, κ1r1 = κr/(√2K), and (3.6), we find that, in Bκ1r1 , 11 1 11 1 v¯y1y1 +ˆa (y )∆d 1v¯ λaˆ (y )¯v =0, − − 11 1 11 1 1 wherea ˆ (y )=¯a (φ− (y )). Equivalently, in Bκ1r1 , 11 1 11 1 ∆¯v λv¯ = 1 aˆ (y ) div 0, v¯ 2 , , v¯ d λ 1 aˆ (y ) v.¯ − − y ··· y − − Then by using the change of variables as well as Theorem 3.3 (note th at κ1 4) we obtain ≥
11 2 2 – a¯ vx1 v¯y1 dx N – v¯y1 v¯y1 dy Br Br B | − 1 | ≤ B | − 1 | Z r Z r1 2 2 d 2 2 Nκ− v¯ + Nκ v¯ + λ v¯ 1 y B 1 y′ B ≤ | | κ1r1 | | | | κ1r1 2 2 d 2 2 Nκ− vx + Nκ vx + λ v , (3.7) ≤ | | Bκr | ′ | | | Bκr where N = N(d, δ, K). We also need estimates for w. By Theorem 3.2
1/2 w + √λ w N I g¯ + λ− I f¯ . k xkL2 k kL2 ≤ k Bκr kL2 k Bκr λkL2 From this and the definition of f¯ it follows that (also note that µ 1) λ ≥ 2 d 2 1 2 2 wx Nκ g¯ + λ− f¯ + λ u¯ , (3.8) | | Br ≤ | | | | | | Bκr 2 2 2 1 2 2 wx + λ w N g¯ + λ− f¯ + λu¯ . (3.9) | | | | Bκr ≤ | | | | | | Bκr This together with u ¯ = w + vyields 2 2 2 2 2 1 2 vx + λ v N u¯x + λ u¯ + g¯ + λ− f¯ . (3.10) | ′ | | | Bκr ≤ | ′ | | | | | | | Bκr To combine all the inequalities shown above, we start with 1/2 11 11 2 11 2 1/2 I := a¯ u¯x1 a¯ u¯x1 a¯ u¯x1 C , Br Br | − | Br ≤ | − | which holds true for any constant C. Upon replacing C with (¯vx1 ) and using Br1 u¯ = w + v again, we arrive at 1/2 1/2 11 2 11 2 2 1/2 I a¯ u¯x1 (¯vx1 ) N a¯ vx1 (¯vx1 ) + N wx Br1 Br1 Br ≤ | − | Br ≤ | − | Br | | =: I1 + I2. From (3.7), (3.10), and (3.9)
1 2 1/2 d/2 2 2 2 1 2 1/2 I1 Nκ− u¯x + Nκ u¯x + λ u¯ + g¯ + λ− f¯ . ≤ | | Bκr | ′ | | | | | | | Bκr 12 H. DONG AND D. KIM
Here we also usedu ¯ = w + v and κ 1. From (3.8), ≥ d/2 2 2 1 2 1/2 I2 Nκ λ u¯ + g¯ + λ− f¯ . ≤ | | | | | | Bκr Finally, notice that
2 1/2 2 2 1/2 1 2 1/2 2 1/2 g¯ Nµ− g˜ + Nµ− u¯x1 + N u¯x , | | Bκr ≤ | | Bκr | | Bκr | ′ | Bκr 2 1/2 2 2 1/2 ( f¯ ) = µ− ( f˜ ) . | | Bκr | | Bκr Therefore, 1/2 1 d/2 1 2 1/2 d/2 2 2 2 1 ˜ 2 I N(κ− + κ µ− ) u¯x1 + Nκ u¯x + λ u¯ + g˜ + λ− f Bκr ′ ≤ | | | | | | | | | | Bκr for µ 1. The lemma is proved. ≥ We recall the maximal function theorem and the Fefferman-Stein theorem. Let the maximal and sharp functions of g defined on Rd be given by
Mg(x) = sup – g(y) dy, | | r>0 Z Br(x)
# g (x) = sup – g(y) (g)Br(x) dy. r>0 | − | Z Br (x) Then g N g# , Mg N g , k kLp ≤ k kLp k kLp ≤ k kLp if g Lp, where 1