ESTIMATES FOR THE MONGE-AMPERE EQUATION

Similar documents
A LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION

1. Introduction Boundary estimates for the second derivatives of the solution to the Dirichlet problem for the Monge-Ampere equation

ALEKSANDROV-TYPE ESTIMATES FOR A PARABOLIC MONGE-AMPÈRE EQUATION

Degenerate Monge-Ampère equations and the smoothness of the eigenfunction

Sobolev regularity for the Monge-Ampère equation, with application to the semigeostrophic equations

Calculus of Variations

The Monge-Ampère Equation. Connor Mooney

BOUNDARY REGULARITY FOR SOLUTIONS TO THE LINEARIZED MONGE-AMPÈRE EQUATIONS

PARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION

A VARIATIONAL METHOD FOR THE ANALYSIS OF A MONOTONE SCHEME FOR THE MONGE-AMPÈRE EQUATION 1. INTRODUCTION

MINIMAL SURFACES AND MINIMIZERS OF THE GINZBURG-LANDAU ENERGY

AFFINE MAXIMAL HYPERSURFACES. Xu-Jia Wang. Centre for Mathematics and Its Applications The Australian National University

Partial regularity for fully nonlinear PDE

REGULARITY RESULTS FOR THE EQUATION u 11 u 22 = Introduction

1 Definition of the Riemann integral

On John type ellipsoids

C 1 regularity of solutions of the Monge-Ampère equation for optimal transport in dimension two

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck

THE TWO-PHASE MEMBRANE PROBLEM REGULARITY OF THE FREE BOUNDARIES IN HIGHER DIMENSIONS. 1. Introduction

The optimal partial transport problem

Singular Solutions to the Monge-Ampère Equation. Connor Mooney

MINIMAL GRAPHS PART I: EXISTENCE OF LIPSCHITZ WEAK SOLUTIONS TO THE DIRICHLET PROBLEM WITH C 2 BOUNDARY DATA

Optimal Transportation. Nonlinear Partial Differential Equations

arxiv: v1 [math.ap] 25 Jul 2012

VISCOSITY SOLUTIONS. We follow Han and Lin, Elliptic Partial Differential Equations, 5.

REGULARITY OF MONOTONE TRANSPORT MAPS BETWEEN UNBOUNDED DOMAINS

Convex solutions to the mean curvature flow

W 2,1 REGULARITY FOR SOLUTIONS OF THE MONGE-AMPÈRE EQUATION

Homogenization and error estimates of free boundary velocities in periodic media

Continuity of convex functions in normed spaces

Everywhere differentiability of infinity harmonic functions

SHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction

NECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES

Remarks on the Gauss-Green Theorem. Michael Taylor

Regularity of flat level sets in phase transitions

PHASE TRANSITIONS: REGULARITY OF FLAT LEVEL SETS

VISCOSITY SOLUTIONS OF ELLIPTIC EQUATIONS

g 2 (x) (1/3)M 1 = (1/3)(2/3)M.

On the second differentiability of convex surfaces

Math The Laplacian. 1 Green s Identities, Fundamental Solution

On supporting hyperplanes to convex bodies

Laplace s Equation. Chapter Mean Value Formulas

A BOUNDARY VALUE PROBLEM FOR MINIMAL LAGRANGIAN GRAPHS. Simon Brendle & Micah Warren. Abstract. The associated symplectic structure is given by

Math 5052 Measure Theory and Functional Analysis II Homework Assignment 7

A few words about the MTW tensor

A GENERAL CLASS OF FREE BOUNDARY PROBLEMS FOR FULLY NONLINEAR ELLIPTIC EQUATIONS

8. Geometric problems

Perron method for the Dirichlet problem.

COMPARISON PRINCIPLES FOR CONSTRAINED SUBHARMONICS PH.D. COURSE - SPRING 2019 UNIVERSITÀ DI MILANO

THE REGULARITY OF MAPPINGS WITH A CONVEX POTENTIAL LUIS A. CAFFARELLI

SOLUTION OF THE DIRICHLET PROBLEM WITH L p BOUNDARY CONDITION. Dagmar Medková

Centre for Mathematics and Its Applications The Australian National University Canberra, ACT 0200 Australia. 1. Introduction

Estimates in surfaces with positive constant Gauss curvature

for all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true

INTERIOR CURVATURE BOUNDS FOR A CLASS OF CURVATURE EQUATIONS

On the regularity of solutions of optimal transportation problems

Math 118B Solutions. Charles Martin. March 6, d i (x i, y i ) + d i (y i, z i ) = d(x, y) + d(y, z). i=1

Some lecture notes for Math 6050E: PDEs, Fall 2016

Obstacle Problems Involving The Fractional Laplacian

arxiv: v3 [math.ap] 20 Jun 2017

Constrained Optimization and Lagrangian Duality

FINITE MORSE INDEX IMPLIES FINITE ENDS

HARDY INEQUALITIES WITH BOUNDARY TERMS. x 2 dx u 2 dx. (1.2) u 2 = u 2 dx.

REGULARITY OF POTENTIAL FUNCTIONS IN OPTIMAL TRANSPORTATION. Centre for Mathematics and Its Applications The Australian National University

Analysis Finite and Infinite Sets The Real Numbers The Cantor Set

U e = E (U\E) e E e + U\E e. (1.6)

Non-Constant Stable Solutions to Reaction- Diffusion Equations in Star-Shaped Domains

Comparison Results for Semilinear Elliptic Equations in Equimeasurable Domains

NOWHERE LOCALLY UNIFORMLY CONTINUOUS FUNCTIONS

Szemerédi-Trotter theorem and applications

A GENERAL CLASS OF FREE BOUNDARY PROBLEMS FOR FULLY NONLINEAR PARABOLIC EQUATIONS

arxiv: v1 [math.mg] 26 Jun 2014

8. Geometric problems

Euler Equations: local existence

is a weak solution with the a ij,b i,c2 C 1 ( )

Variational problems with free boundaries for the fractional Laplacian

EE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 17

An introduction to Mathematical Theory of Control

Mathematical Research Letters 4, (1997) HARDY S INEQUALITIES FOR SOBOLEV FUNCTIONS. Juha Kinnunen and Olli Martio

Supplement: Universal Self-Concordant Barrier Functions

Nonlinear Programming Algorithms Handout

and BV loc R N ; R d)

The multi-layer free boundary problem for the p-laplacian in convex domains

The oblique derivative problem for general elliptic systems in Lipschitz domains

Class Meeting # 9: Poisson s Formula, Harnack s Inequality, and Liouville s Theorem

BORIS MORDUKHOVICH Wayne State University Detroit, MI 48202, USA. Talk given at the SPCOM Adelaide, Australia, February 2015

Variational problems with free boundaries for the fractional Laplacian

REGULARITY FOR INFINITY HARMONIC FUNCTIONS IN TWO DIMENSIONS

Γ -convergence of the Allen Cahn energy with an oscillating forcing term

Optimization and Optimal Control in Banach Spaces

Functional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...

Extension Properties of Pluricomplex Green Functions

ESTIMATES FOR ELLIPTIC HOMOGENIZATION PROBLEMS IN NONSMOOTH DOMAINS. Zhongwei Shen

THEOREMS, ETC., FOR MATH 516

Geometric problems. Chapter Projection on a set. The distance of a point x 0 R n to a closed set C R n, in the norm, is defined as

Regularity of Weak Solution to Parabolic Fractional p-laplacian

Pointwise estimates for marginals of convex bodies

CONSTRAINT QUALIFICATIONS, LAGRANGIAN DUALITY & SADDLE POINT OPTIMALITY CONDITIONS

LUSIN PROPERTIES AND INTERPOLATION OF SOBOLEV SPACES. Fon-Che Liu Wei-Shyan Tai. 1. Introduction and preliminaries

Transcription:

GLOBAL W 2,p ESTIMATES FOR THE MONGE-AMPERE EQUATION O. SAVIN Abstract. We use a localization property of boundary sections for solutions to the Monge-Ampere equation obtain global W 2,p estimates under natural assumptions on the domain boundary data. 1. Introduction Interior W 2,p estimates for strictly convex solutions for the Monge-Ampere equation were obtained by Caffarelli in [C] under the necessary assumption of small oscillation of the right h side. The theorem can be stated as follows. Theorem 1.1 (Interior W 2,p estimates). Let u : Ω R, (1.1) u = 0 on Ω, be a continuous convex solution to the Monge-Ampere equation (1.2) det D 2 u = f(x) in Ω, 0 < λ f Λ, for some positive constants λ, Λ. For any p, 1 < p < there exists ε(p) > 0 depending on p n such that if B ρ Ω B 1/ρ (1.3) sup log f(x) log f(y) ε(p), x y ρ for some small ρ > 0, then C depends on ρ, λ, Λ, p n. u W 2,p ({u< ρ}) C, In this short paper we answer the natural question of whether this interior estimates can be extended up to the boundary. Precisely, we obtain the following global W 2,p estimate under natural assumptions on the domain boundary data. Theorem 1.2 (Global W 2,p estimates). Let Ω be a convex bounded domain u : Ω R be a Lipschitz continuous convex solution of the Monge-Ampere equation (1.2). Assume that u Ω, Ω C 1,1, there exists ρ > 0 small such that f satisfies (1.3). If u separates quadratically on Ω from its tangent planes i.e (1.4) u(y) u(x) u(x) (y x) ρ x y 2 x, y Ω, The author was partially supported by NSF grant 0701037. 1

2 O. SAVIN then u W 2,p (Ω) C, with C depending on Ω C 1,1, u Ω C 1,1, u C 0,1, ρ, λ, Λ, p n. Remark 1.3. The gradient u(x) when x Ω is understood in the sense that y n+1 = u(x) + u(x) (y x) is a supporting hyperplane for the graph of u(y) but y n+1 = u(x) + ( u(x) δν x ) (y x) is not a supporting hyperplane for any δ > 0, where ν x denotes the exterior normal to Ω at x. In general the Lipschitz continuity of the solution can be easily obtained from the boundary data by the use of barriers. Also the quadratic separation assumption (1.4) can be checked in several situations directly from the boundary data, see Proposition 3.2 in [S2]. This is the case for example when the boundary data is more regular i.e (1.5) u Ω, Ω C 3, Ω is uniformly convex. As a consequence of Theorem 1.2 we obtain Corollary 1.4. Let u : Ω R be a continuous solution of the Monge-Ampere equation (1.2) that satisfies (1.3) (1.5). Then u W 2,p (Ω) C, with C depending on u Ω, Ω, ρ, λ, Λ, p n. In particular, if u solves the Monge-Ampere equation (1.2) with f C(Ω) (1.5) holds, then u W 2,p (Ω) for any p <. The assumptions on the boundary behavior of u Ω in Theorem 1.2 Corollary 1.4 seem to be optimal. Wang in [W] gave examples of solutions u to (1.2) with f = 1 either u C 2,1 or Ω C 2,1, that do not belong to W 2,p (Ω) for large values of p. The proof of Theorem 1.2 is based on a localization theorem for the Monge- Ampere equation at boundary points which was proved in [S1], [S2]. It states that under natural local assumptions on the domain boundary data, the sections S h (x 0 ) = {x Ω u(x) < u(x 0 ) + u(x 0 ) (x x 0 ) + h}, with x 0 Ω are equivalent to ellipsoids centered at x 0. statement below. Assume for simplicity that (1.6) B ρ (ρe n ) Ω {x n 0} B 1 ρ, We give its precise for some small ρ > 0, that is Ω (R n ) + Ω contains an interior ball tangent to Ω at 0. Let u : Ω R be continuous, convex, satisfying (1.7) det D 2 u = f, 0 < λ f Λ in Ω. After subtracting a linear function we also assume that (1.8) x n+1 = 0 is the tangent plane to u at 0, in the sense that u 0, u(0) = 0,

GLOBAL W 2,p ESTIMATES 3 any hyperplane x n+1 = δx n, δ > 0, is not a supporting plane for u. Theorem 1.5 shows that if the boundary data has quadratic growth near {x n = 0} then, each section of u at 0 S h := S h (0) = {x Ω : is equivalent to a half-ellipsoid centered at 0. u(x) < h}, Theorem 1.5 (Localization theorem). Assume that Ω, u satisfy (1.6)-(1.8) above, ρ x 2 u(x) ρ 1 x 2 on Ω {x n ρ}. Then, for each h < c 0 there exists an ellipsoid E h of volume ω n h n/2 such that c 0 E h Ω S h c 1 0 E h Ω. Moreover, the ellipsoid E h is obtained from the ball of radius h 1/2 by a linear transformation A 1 h (sliding along the x n = 0 plane) A h E h = h 1/2 B 1 A h (x) = x νx n, ν = (ν 1, ν 2,..., ν n 1, 0), with ν c 1 0 log h. The constant c 0 > 0 above depends on ρ, λ, Λ, n. We describe the idea of the proof of Theorem 1.2 below. If the assumptions of Theorem 1.2 hold say for simplicity with f = 1 then, by rescaling, one can deduce from the Localization theorem above that the second derivatives can blow-up near Ω in a logarithmic fashion i.e D 2 u(y) C log dist(y, Ω) 2, thus D 2 u L p (Ω) for any p <. For general f we need to use the interior W 2,p estimates instead of the estimate above together with a Vitali covering lemma for sections of u near Ω. 2. Proof of Theorem 1.2 We start by remarking that under the assumptions of Theorem 1.5 above, we obtain that also the section {u < h 1/2 x n } has the shape of E h. Indeed, since S h c 1 0 E h {x n c 1 0 h1/2 } u(0) = 0, we can conclude from the convexity of u that the set satisfies for all small h F := {x Ω u < c 0 h 1/2 x n } (2.1) F S h Ω, F is tangent to Ω at 0. We show that F is equivalent to E h by bounding its volume from below. Lemma 2.1. We have F c E h for some c > 0 small depending on ρ, λ, Λ, n.

4 O. SAVIN Proof. From Theorem 1.5, there exists y S θh such that y n c 0 (θh) 1/2. We evaluate v := u c 0 h 1/2 x n, at y find v(y) θh c 0 h 1/2 c 0 (θh) 1/2 δh, for some δ > 0 provided that we choose θ small depending on c 0. Since v = 0 on F det D 2 v λ we have hence inf F v C(λ) F 2/n, ch n/2 F. Next we prove Theorem 1.2. We denote by c, C positive constants that depend on ρ, λ, Λ, p, n Ω C 1,1, u Ω C 1,1, u C 0,1. For simplicity of notation, their values may change from line to line whenever there is no possibility of confusion. Given y Ω we denote by S h (y) := {x Ω u(x) < u(y) + u(y) (x y) + h} let h(y) be the maximal value of h such that S h (y) Ω, i.e h(y) := max{h 0 S h (y) Ω}. Lemma 2.2. Let y Ω denote for simplicity h = h(y). If h c 1 then S h(y) D C h1/2 := {x Ω dist(x, Ω) C h 1/2 }, S h/2 (y) Proof. Without loss of generality assume D 2 u p dx C log h 2p hn/2. 0 S h(y) Ω. After subtracting a linear function relabeling ρ if necessary we may also assume that u satisfies the conditions of Theorem 1.5 at the origin. Since u separates quadratically from 0 on Ω we deduce that u(y) = αe n thus h > 0. Since (see [C]) we obtain from (2.1) Lemma 2.1 S h(y) = {x Ω u < αx n }, α > 0, c h n/2 S h(y) C h n/2, c h 1/2 α C h 1/2, that S h(y) is equivalent to an ellipsoid Ẽ := E C h. This implies that with y + cẽ S h(y) CẼ, Ẽ = h 1/2 A 1 B 1, Ax := x νx n,

GLOBAL W 2,p ESTIMATES 5 The inclusion above implies The rescaling ũ : Ω R of u satisfies Ω = h 1/2 A(Ω y), ν n = 0, ν C log h. S h(y) {x n C h 1/2 } D C h1/2. ũ(x) := 1 hu(y + h 1/2 A 1 x), S1 (0) = h 1/2 A(S h(y) y) Ω, B c S 1 (0) B C, where S 1 (0) represents the section of ũ at the origin at height 1. Moreover, (1.3) implies since det D 2 ũ(x) = f(x) := f(y + h 1/2 A 1 x) log f(x) log f(z) ε(p) x, z S 1 (0) h 1/2 A 1 (x z) C h 1/2 log h ρ. The interior W 2,p estimate for ũ in S 0 (1) gives D 2 ũ p dx C, S 1/2 (0) hence D 2 u p dx = A T D 2 ũ A p hn/2 dx C log h 2p hn/2. S h/2 (y) S 1/2 (0) We also need the following covering lemma. Lemma 2.3 (Vitali covering). There exists a sequence of disjoint sections such that S δ hi (y i ), hi = h(y i ), Ω S hi/2(y i ), where δ > 0 is a small constant that depends only on λ, Λ n. Proof. We choose δ such that if then i=1 S δ h(y) (y) S δ h(z) (z) 2 h(y) h(z), S δ h(z) (z) S h(y)/2 (y). The existence of δ follows from the engulfing properties of sections of solutions to Monge-Ampere equation (1.2) with bounded right h side (see [CG]).

6 O. SAVIN Now the proof is identical to the proof of Vitali s covering lemma for balls. We choose S δ h1 (y 1 ) from all the sections S δ h(y) (y), y Ω such that h(y 1 ) 1 2 sup h(y), y then choose S δ h2 (y 2 ) as above but only from the remaining sections S δ h(y) (y) that are disjoint from S δ h1 (y 1 ), then S δ h3 (y 3 ) etc. We easily obtain Ω = S δ h(y) (y) S hi/2(y i ). y Ω i=1 End of proof of Theorem 1.2 D 2 u p dx Ω i There are a finite number of sections with h i c 1, by the interior W 2,p estimate, in each such section we have D 2 u p dx C. S hi /2 (yi) S hi /2 (yi) D 2 u p dx. Next we consider the family F d of sections S hi/2(y i ) such that d/2 < h i d for some constant d c 1. By Lemma 2.2 in each such section S hi /2 (yi) D 2 u p dx C log d 2p S δ hi (y i ), since are disjoint we find i F d S δ hi (y i ) D Cd 1/2 S hi /2 (yi) D 2 u p dx C log d 2p d 1/2. We add these inequalities for the sequence d = c 1 2 k, k = 0, 1,... obtain the desired bound. References [C] Caffarelli L., Interior W 2,p estimates for solutions of Monge-Ampere equation, Ann. of Math. 131 (1990), 135-150. [CG] Caffarelli L., Gutierrez C., Properties of solutions of the linearized Monge-Ampere equations, Amer. J. Math. 119 (1997), 423-465. [S1] Savin O., A localization property at the boundary for Monge-Ampere equation, Preprint arxiv:1010.1745. [S2] Savin O, Pointwise C 2,α estimates at the boundary for Monge-Ampere equation, Preprint arxiv:1101.5436. [W] Wang X.-J., Regularity for Monge-Ampere equation near the boundary, Analysis 16 (1996) 101-107.

GLOBAL W 2,p ESTIMATES 7 Department of Mathematics, Columbia University, New York, NY 10027 E-mail address: savin@math.columbia.edu