Monge–Ampère equation

From HandWiki
Short description: Nonlinear second-order partial differential equation of special kind

In mathematics, a (real) Monge–Ampère equation is a nonlinear second-order partial differential equation of special kind. A second-order equation for the unknown function u of two variables x, y is of Monge–Ampère type if it is linear in the determinant of the Hessian matrix of u and in the second-order partial derivatives of u. The independent variables (x, y) vary over a given domain D of 2. The term also applies to analogous equations with n independent variables. The most complete results so far have been obtained when the equation is elliptic.

Definition

In two dimension

Given two independent variables x and y, and one dependent variable u, the general Monge–Ampère equation is of the form

L[u]=Adet(2u)+BΔu+2Cuxy+(DB)uyy+E=A(uxxuyyuxy2)+Buxx+2Cuxy+Duyy+E=0,

where A, B, C, D, and E are functions depending on the first-order variables x, y, u, ux, and uy only.

In general

Given a domain Ωn and a real-valued function u:Ω, a (real) Monge–Ampère equation is any fully nonlinear second-order equation that can be written in the form

F(x,u,Du,detD2u)=0,

for some function F. More generally, Ω can be a Riemannian manifold, since Du,D2u are well-defined on a Riemannian manifold.

If F also depends linearly on all principal minors of the Hessian matrix D2u, then it is an equation of Monge–Ampère type.

Classification

As for other second-order fully nonlinear equations, the type of a Monge–Ampère equation is defined by the linearization of the operator at a sufficiently smooth solution. Of these, the most common is the elliptic case. When people say "Monge–Ampère equation" without adjective, they usually mean the elliptic case.

Let Ωn be an open set, u:Ω be a C2 function, and consider an operator of Monge–Ampère type

L[u]=F(x,u,Du,D2u),

where F is smooth in all variables and depends on D2u only through its principal minors. The linearization of L is of form

ijaij(x)uij+lower order terms,

where

aij(x)=Fuij(x,u(x),Du(x),D2u(x)).

The quadratic form

Qx(ξ)=aij(x)ξiξj,ξn,

is the principal symbol of the linearized operator at the point x.

The equation is said to be

  • elliptic at x if Qx(ξ) if all eigenvalues are of the same sign,
  • hyperbolic at x if Qx(ξ) takes both positive and negative values (the matrix is indefinite),
  • parabolic at x if Qx(ξ) is degenerate (the matrix has vanishing determinant),
  • degenerate elliptic if it is elliptic everywhere,
  • elliptic if it is degenerate elliptic, and all eigenvalues are a bounded distance away from zero.

As follows from Jacobi's formula for the derivative of a determinant, this equation is elliptic if f is a positive function and solutions satisfy the constraint of being uniformly convex. If f is merely strictly convex, then the equation is degenerate-elliptic.({{{1}}}, {{{2}}})

Examples

The Monge–Ampère equation in its simplest form is

detD2u=f(x)

where f is a given function on a domain Ωn. This is a special case of equation (1) below with f(x,u,Du)=f(x).

The classical Liouville theorem has an analogy here. If detD2u is constant, and u is defined on all of n, then u is a quadratic function. This is the Jörgens–Calabi–Pogorelov theorem.[1]: Sec. 4.3 

If f(x) is positive and uniformly convex, and u is a solution to detD2u=f(x), then its Legendre transform u* is a solution to detD2u*=1/f(x).

Geometry

Monge–Ampère equations arise naturally in several problems in Riemannian geometry, conformal geometry, affine geometry, and CR geometry.

Given a twice-differentiable real-valued function f:Ω defined over a domain Ωn, its graph is a manifold of n dimensions. At any xΩ, the Gaussian curvature of the manifold at f(x) is detD2u(1+|Du|2)(n+2)/2. Thus, if we want to find a manifold whose Gaussian curvature is an arbitrary function we pick ourselves, then we need to solve the following Monge–Ampère equation:

detD2uK(x)(1+|Du|2)(n+2)/2=0

where K is the Gaussian curvature we want. Given such a function K, it is nontrivial to find a solution, if any. The problem of finding a solution is the Minkowski problem, or the prescribed Gaussian curvature problem.({{{1}}}, {{{2}}})

For example, the rigidity of the 2-sphere manifests as the fact that if we require n=3,K=1,u(0)=0,Du=0, then there are just two unique solutions, which is the unit 2-sphere and its reflection.

The affine spheres can be characterized by a Monge–Ampère equation.

Optimal transport

Consider the problem of optimal transport with quadratic cost (this is also called the 2-Wasserstein metric problem) on n. That is, suppose μ,ν are distributions on n with probability density functions ρμ,ρν. In this case, a map T:supp(μ)supp(ν) is a transport map iff it satisfiesh(T(x))ρμ(x)dx=h(y)ρν(y)dyfor any integrable test function hL1(supp(ν)). The problem is to find the T that minimizes the following quadratic cost function:minTxT(x)2f(x)dxBy a theorem of Brenier, the optimal transport map exists, and is the gradient of a convex function ψ:supp(μ)n, with T=Dψ. The convex function satisfies a Monge–Ampère equation:({{{1}}}, {{{2}}}): 282 [2]({{{1}}}, {{{2}}}){det(D2ψ)=ρμρνDψDψ(supp(μ))=supp(ν)The boundary condition simply states that the optimal transport maps the boundary of the source to the boundary of the target. Furthermore, the solution ψ is almost everywhere unique.

The function ψ is called the potential function of the problem in this case.

Conversely, some Monge–Ampère equations can be interpreted optimal transport. Weak-solutions of a Monge–Ampère equations obtained by optimal transport are often called Brenier solutions in the literature. Brenier solutions satisfy their corresponding Monge–Ampère equations almost everywhere.({{{1}}}, {{{2}}}): 323 

Rellich's theorem

Let Ω be a bounded domain in 3, and suppose that on Ω the coefficients A, B, C, D, and E are continuous functions of x and y only. Consider the Dirichlet problem to find u so that

L[u]=0,on Ω
u|Ω=g.

If

BDC2AE>0,

then the Dirichlet problem has at most two solutions.({{{1}}}, {{{2}}})

Ellipticity results

Suppose now that 𝐱 is a variable with values in a domain in n, and that f(𝐱,u,Du) is a positive function. Then the Monge–Ampère equation

L[u]=detD2uf(𝐱,u,Du)=0(1)

is a nonlinear elliptic partial differential equation (in the sense that its linearization is elliptic), provided one confines attention to convex solutions.


See also

References

  1. Figalli, Alessio (2017). The Monge-Ampère equation and its applications. Zurich lectures in advanced mathematics. Zürich: European Mathematical Society. ISBN 978-3-03719-170-5. 
  2. Prins, C. R.; Beltman, R.; ten Thije Boonkkamp, J. H. M.; IJzerman, W. L.; Tukker, T. W. (January 2015). "A Least-Squares Method for Optimal Transport Using the Monge--Ampère Equation" (in en). SIAM Journal on Scientific Computing 37 (6): B937–B961. doi:10.1137/140986414. ISSN 1064-8275. http://epubs.siam.org/doi/10.1137/140986414. 

Additional references