A fully non-linear version of the incompressible Euler equations: the semi-geostrophic system, , November 27, 2006

existence and stability of measure valued solutions, existence and uniqueness of solutions under certain continuity conditions for the density, convergence to the incompressible Euler equations. Meanwhile, a general technique to prove uniqueness of sufficiently smooth solutions to non-linearly coupled system is introduced, using optimal transportation.

Contents

1 Introduction

The semi-geostrophic equations are an approximation to the Euler equations of fluid mechanics, used in meteorology to describe atmospheric flows, in particular they are believed (see [?] ) to be an efficient model to describe frontogenesis. Different versions (incompressible [?] , shallow water [?] , compressible [?] ) of this model have been studied, and we will focus here on the incompressible 2-d and 3-d version. The 3-d model describes the behavior of an incompressible fluid in a domain

Ω \subset R^{3}

. To the evolution in

Ω

is associated a motion in a ’dual’ space, described by the following non-linear transport equation:

\begin{matrix} \partial_{t} ρ + \nabla \cdot (ρ v) = 0, \end{matrix}

\begin{matrix} v = (\nabla Ψ (x) - x)^{⊥}, \end{matrix}

\begin{matrix} det D^{2} Ψ = ρ, \end{matrix}

\begin{matrix} ρ (t = 0) = ρ^{0} . \end{matrix}

Here

ρ^{0}

is a probability measure on

R^{3}

, and for every

v = (v_{1}, v_{2}, v_{3}) \in R^{3}

v^{⊥}

stands for

(- v_{2}, v_{1}, 0)

. The velocity field is given at each time by solving a Monge-Ampere equation in the sense of the polar factorization of maps (see [?] ), i.e. in the sense that

Ψ

is convex from

R^{3}

R

and satisfies

\nabla Ψ_{#} ρ = χ_{Ω} ℒ^{3}

, where

ℒ^{3}

is the Lebesgue measure of

R^{3}

, and

χ_{Ω}

is the indicator function of

Ω

. For compatibility

Ω

has Lebesgue-measure one. This model arises as an approximation to the primitive equations of meteorology, and we shall give a brief idea of the derivation of the model, although the reader interested in more details should refer to [?] .

In this work we will deal with various questions related to the semi-geostrophic (hereafter

S G

) system: existence and stability of measure-valued solutions, existence and uniqueness of smooth solutions, and finally convergence towards the incompressible Euler equations in 2-d. The results are stated in more details in section 1.

1.1 Derivation of the semi-geostrophic equations

We now give for sake of completeness a brief and simplified idea of the derivation of the system, inspired from [?] , and more complete arguments can be found in [?] .

Lagrangian formulation We start from the 3-d incompressible Euler equations with constant Coriolis parameter

f

in a domain

Ω

\begin{matrix} \frac{D v}{D t} + f v^{⊥} = \frac{1}{ρ} \nabla p - \nabla φ, \end{matrix}

\begin{matrix} \nabla \cdot v = 0, \frac{D ρ}{D t} = 0, \end{matrix}

\begin{matrix} v \cdot \partial Ω = 0, \end{matrix}

where

\frac{D \cdot}{D t}

stands for

\partial_{t} + v \cdot \nabla

, and we still use

v^{⊥} = (- v_{2}, v_{1}, 0)

. The term

\nabla φ

denotes the gravitational effects (here we will take

φ = g x_{3}

with constant

g

), and the term

f v^{⊥}

is the Coriolis force due to rotation of the Earth. For large scale atmospheric flows, the Coriolis force

f v^{⊥}

dominates the advection term

\frac{D v}{D t}

, and renders the flow mostly bi-dimensional. We use the hydrostatic approximation:

\partial_{x_{3}} p = - ρ g

and restrict ourselves to the case

ρ \equiv 1

Keeping only the leading order terms leads to the geostrophic balance

\begin{matrix} v_{g} = - f^{- 1} \nabla^{⊥} p, \end{matrix}

that defines

v_{g}

, the geostrophic wind. Decomposing

v = v_{g} + v_{a g}

where the second component is the ageostrophic wind, supposedly small departures from the geostrophic balance, the semi-geostrophic system reads:

\begin{matrix} \frac{D v_{g}}{D t} + f v^{⊥} = \nabla_{H} p, \end{matrix}

\begin{matrix} \nabla \cdot v = 0, \end{matrix}

where

\nabla_{H} = (\partial_{x_{1}}, \partial_{x_{2}}, 0)

. Note however that the advection operator

\partial_{t} + v \cdot \nabla

still uses the full velocity

v

. Introducing the potential

\begin{matrix} Φ = \frac{1}{2} | x_{H} |^{2} + f^{- 2} p, \end{matrix}

with

x_{H} = (x_{1}, x_{2}, 0)

, we obtain the following

\begin{matrix} \frac{D}{D t} \nabla Φ (t, x) = f (x - \nabla Φ (t, x))^{⊥} . \end{matrix}

We introduce the lagrangian map

g : Ω \times R^{+} \mapsto Ω

giving the position at time

t

of the parcel located at

x_{0}

at time 0. The previous equation means that, if for fixed

x \in Ω

we consider the trajectory in the ’dual’ space, defined by

X (t, x) = \nabla Φ (t, g (t, x))

, we have

\begin{matrix} \partial_{t} X (t, x) = f {(g (t, x) - X (t, x))}^{⊥} . \end{matrix}

By rescaling the time, we can set

f = 1

. As stated the system looks under-determined: indeed

Φ

is unknown; however we have the condition

X (t, x) = \nabla Φ (t, g (t, x))

. Moreover, the dynamic in the

x

space being incompressible and contained in

Ω

, the map

g (t, \cdot)

must be measure preserving in

Ω

for each

t

, i.e.

\begin{matrix} ℒ^{3} (g (t)^{- 1} (B)) = ℒ^{3} (B) \end{matrix}

for each

B \subset Ω

measurable (where

ℒ^{3}

denotes the Lebesgue measure of

R^{3}

We shall hereafter denote by

G (Ω)

the set of all such measure preserving maps. Then Cullen’s stability criteria asserts that the potential

Φ

should be convex for the system to be stable to small displacements of particles in the

x

space. Hence, for each

t

Φ

must be a convex function such that

\begin{matrix} X (t, \cdot) = \nabla Φ (t, g (t, \cdot)), \end{matrix}

with

g (t, \cdot) \in G (Ω)

. In the next paragraph we shall see that, under very mild assumptions on

X

, this decomposition, called polar factorization, can only happen for a unique choice of

g

and

\nabla Φ

. Now if

Φ^{*}

is the Legendre transform of

Φ

\begin{matrix} Φ^{*} (y) = {sup}_{x \in Ω} x \cdot y - Φ (x), \end{matrix}

then

\nabla Φ

and

\nabla Φ^{*}

are inverse maps of each other, and the semi-geostrophic system then reads

\begin{matrix} \frac{D X}{D t} = {(\nabla Φ^{*} (X (t)) - X (t))}^{⊥}, \end{matrix}

\begin{matrix} \nabla Φ^{*} (t) \circ X (t) \in G (Ω) . \end{matrix}

In the next paragraph, we expose the results concerning the existence and uniqueness of the gradients

\nabla Φ, \nabla Φ^{*}

1.2 Polar factorization of vector valued maps

The polar factorization of maps has been discovered by Brenier in [?] .

It has later been extended to the case of general Riemannian manifolds by McCann in [?] .

The Euclidean case Let

Ω

be a fixed bounded domain of

R^{d}

of Lebesgue measure 1 and satisfying the condition

ℒ^{d} (\partial Ω) = 0

. We consider a mapping

X \in L^{2} (Ω; R^{d})

. We will also consider the push-forward of the Lebesgue measure of

Ω

X

, that we will denote by

X_{#} χ_{Ω} ℒ^{d} = d ρ

(or, in short,

X_{#} d x

) and which is defined by

\begin{matrix} \forall f \in C_{b}^{0} (R^{d}), \int_{R^{d}} f (x) d ρ (x) = \int_{Ω} f (X (x)) d x . \end{matrix}

Let

P

be the set of probability measures

R^{d}

, and

P_{a}^{2}

the subset of

P

where the subscript

a

means absolutely continuous with respect to the Lebesgue measure (or equivalently that have a density in

L^{1} (R^{d})

), and the superscript

2

means with finite second moment. (i.e. such that

\begin{matrix} \int_{R^{d}} | x |^{2} d ρ (x) < + \infty .) \end{matrix}

Note that for

X \in L^{2} (Ω, R^{d})

, the measure

ρ = X_{#} d x

has necessarily finite second moment, and thus belongs to

P^{2}

Theorem 1.1 (Brenier, [?] ). Let

Ω

be as above,

X \in L^{2} (Ω; R^{d})

and

ρ = X_{#} d x

1. There exists a unique up to a constant convex function, that will be denoted

Φ [ρ]

, such that:

\begin{matrix} \forall f \in C_{b}^{0} (R^{d}), \int_{Ω} f (\nabla Φ [ρ] (x)) d x = \int_{R^{d}} f (x) d ρ (x) . \end{matrix}

2. Let

Ψ [ρ]

be the Legendre transform of

Φ [ρ]

, if

ρ \in P_{a}^{2}

Ψ [ρ]

is the unique up to a constant convex function satisfying

\begin{matrix} \forall f \in C_{b}^{0} (Ω), \int_{R^{d}} f (\nabla Ψ [ρ] (x)) d ρ (x) = \int_{Ω} f (x) d x . \end{matrix}

3. If $ρ \in P_{a}^{2}$ , $X$ admits the following unique polar factorization:

$\begin{matrix} X = \nabla Φ [ρ] \circ g, \end{matrix}$
with $Φ [ρ]$ convex, $g$ measure preserving in $Ω$ .

Remark:

Ψ [ρ], Φ [ρ]

depend only on

ρ

, and are solutions (in some weak sense) respectively in

R^{d}

and

Ω

, of the Monge-Ampere equations

\begin{matrix} det D^{2} Ψ = ρ, \end{matrix}

\begin{matrix} ρ (\nabla Φ) det D^{2} Φ = 1 . \end{matrix}

When

Ψ

and

Φ

are not in

C_{l o c}^{2}

these equations can be understood in the viscosity (or Alexandrov) sense or in the sense of Theorem 1.1 , which is strictly weaker. For the regularity of those solutions and the consistency of the different weak formulations the reader can refer to [?] .

The periodic case The polar factorization theorem has been extended to Riemannian manifolds in [?] (see also [?] for the case of the flat torus). In this case, we consider a mapping

X : R^{d} \mapsto R^{d}

such that for all

\vec{p} \in Z^{d}, X (\cdot + \vec{p}) = X + \vec{p}

. Then

ρ = X_{#} d x

is a probability measure on

T^{d}

. We define

Ψ [ρ], Φ [ρ]

through the following:

Theorem 1.2. Let

X : R^{d} \to R^{d}

be as above, with

ρ = X_{#} d x

1. Up to an additive constant there exists a unique convex function

Φ [ρ]

such that

Φ [ρ] (x) - x^{2} / 2

Z^{d}

-periodic (and thus

\nabla Φ [ρ] (x) - x

Z^{d}

periodic), and

\begin{matrix} \forall f \in C^{0} (T^{d}), \int_{T^{d}} f (\nabla Φ [ρ] (x)) d x = \int_{T^{d}} f (x) d ρ (x) . \end{matrix}

2. Let

Ψ [ρ]

be the Legendre transform of

Φ [ρ]

. If

ρ

is Lebesgue integrable,

Ψ [ρ]

is the unique up to a constant convex function satisfying

Ψ [ρ] (x) - x^{2} / 2

Z^{d}

-periodic (and thus

\nabla Ψ [ρ] (x) - x

Z^{d}

periodic), and

\begin{matrix} \forall f \in C^{0} (T^{d}), \int_{T^{d}} f (\nabla Ψ [ρ] (x)) d ρ (x) = \int_{T^{d}} f (x) d x . \end{matrix}

3. If $ρ$ is Lebesgue integrable, $X$ admits the following unique polar factorization:
$\begin{matrix} X = \nabla Φ [ρ] \circ g \end{matrix}$
with $g$ measure preserving from $T^{d}$ into itself, and $Φ [ρ]$ convex, $Φ [ρ] - | x |^{2} / 2$ periodic.

Remark 1: From the periodicity of

\nabla Φ [ρ] (x) - x, \nabla Ψ [ρ] (x) - x

, for every

f

Z^{d}

-periodic ,

f (\nabla Ψ [ρ]), f (\nabla Φ [ρ])

are well defined on

R^{d} / Z^{d}

Remark 2: Both in the periodic and non periodic case, the definitions of

Ψ [ρ]

and

Φ [ρ]

make sense if

ρ

is absolutely continuous with respect to the Lebesgue measure. If not, the definition and uniqueness of

Φ [ρ]

is still valid, as well as the property

\nabla Φ [ρ]_{#} ρ = χ_{Ω} ℒ^{d}

. The definition of

Ψ [ρ]

as the Legendre transform of

Φ [ρ]

is still valid also, but then the expression

\int f (\nabla Ψ [ρ] (x)) d ρ (x)

does not necessarily make sense since

\nabla Ψ

is not necessarily continuous. Moreover the polar factorization does not hold any more.

Remark 3: We have (see [?] ) the unconditional bound

∥ \nabla Ψ [ρ] (x) - x ∥_{L^{\infty} (T^{d})} \leq \sqrt{d} / 2

that will be useful later on.

1.3 Lagrangian formulation of the $S G$ system

From Theorems 1.1 , 1.2 the Lagrangian formulation of the semi-geostrophic equation then becomes

\begin{matrix} \frac{D X}{D t} = {[\nabla Ψ (X) - X]}^{⊥}, \end{matrix}

(1)

\begin{matrix} Ψ = Ψ [ρ], ρ = X_{#} d x . \end{matrix}

(2)

1.4 Eulerian formulation in dual variables

In both cases (periodic and non periodic) we thus investigate the following system that will be referred to as

S G

: we look for a time dependent probability measure

t \to ρ (t, \cdot)

satisfying

\begin{matrix} \partial_{t} ρ + \nabla \cdot (ρ v) = 0, \end{matrix}

(3)

\begin{matrix} v (t, x) = {(\nabla Ψ [ρ (t)] (x) - x)}^{⊥}, \end{matrix}

(4)

\begin{matrix} ρ (t = 0) = ρ^{0} . \end{matrix}

(5)

Weak solutions (which are defined below) of this system with

L^{p}

initial data for

p \geq 1

have been found, see [?] , [?] , [?] .

1.5 Results

In this work we deal with various mathematical problems related to this system: we first extend the notion of weak solutions that had been shown to exist for

ρ \in L^{\infty} (R_{+}, L^{q} (R^{3}))

q > 1

([?] , [?] ), and then for

ρ \in L^{\infty} (R_{+}, L^{1} (R^{3}))

([?] ), to the more general case of bounded measures. The question of existence of measure-valued solutions was raised and left unanswered in those papers, and we show here existence of global solutions to the Cauchy problem with initial data a bounded compactly supported measure, and show the weak stability/compactness of these weak measure solutions.

Then we show existence of continuous solutions, more precisely, we show local existence of solutions with Dini-continuous (see ( 8 )) density. For this solutions, the velocity field is then

C^{1}

and the Lagrangian system ( 1 , 2 ) is defined everywhere.

We also show uniqueness in the class of Holder continuous solutions (a sub-class of Dini continuous solutions). This proof uses in an original way the optimal transportation of measures by convex gradients and its regularity properties, and can be adapted to give a new proof of uniqueness for solutions of the 2-d Euler equation with bounded vorticity, but also for a broad class of non-linearly coupled system. The typical application is a density evolving through a transport equation where the velocity field depends on the gradient of a potential, the potential solving an elliptic equation with right hand side the density. Well known examples of such cases are the Vlasov-Poisson and Euler-Poisson systems.

Finally, in the 2-d case, we study the convergence of the system to the Euler incompressible equations; this convergence is expected for

ρ

close to 1, since formally expanding

Ψ = x^{2} / 2 + ε ψ

, and linearizing the determinant around the identity matrix, we get

det D^{2} Ψ = 1 + ε Δ ψ + O (ε^{2}),

and the Monge-Ampere equation turns into the Poisson equation

Δ ψ = \frac{ρ - 1}{ε} = : μ .

After a proper time scaling,

μ

satisfies

\begin{matrix} \partial_{t} μ + \nabla \cdot (μ \nabla^{⊥} ψ), \end{matrix}

\begin{matrix} Δ ψ = μ, \end{matrix}

that we recognize as the vorticity formulation of the 2-d Euler incompressible equation. The study of this ’quasi-neutral’ limit is done by two different ways:

One uses a modulated energy method similar as the one used in [?] and [?] and is valid for weak solutions. The other uses a more classical expansion of the solution, and regularity estimates, and is similar to the method used in [?] .

The second method yields also a time of existence for the smooth solution that goes to infinity, as the scaling parameter

ε

goes to 0. From a physical point of view, this asymptotic study may be seen as a justification of the consistency of the semi-geostrophic approximation.

2 Measure valued solutions

2.1 A new definition of weak solutions

We have first the following classical weak formulation of equation ( 3 ):

ρ \in C (R +, L^{1} (R^{3}) - w)

is said to be a weak solution of

S G

\begin{matrix} \forall T > 0, \forall φ \in C_{c}^{\infty} ([0, T] \times R^{2}), \end{matrix}

\begin{matrix} \int \partial_{t} φ ρ + \nabla φ \cdot (\nabla Ψ [ρ] - x)^{⊥} ρ d t d x \end{matrix}

\begin{matrix} = & \int φ (T, x) ρ (T, x) d x - \int φ (0, x) ρ (0, x) d x, \end{matrix}

where for all

t

Ψ [ρ]

is as in Theorem 1.1 . The problematic part in the case of measure valued solutions is to give sense to the product

ρ \nabla Ψ [ρ]

since at the point where

ρ

is singular

\nabla Ψ [ρ]

is unlikely to be continuous. Therefore we use the Theorem 1.1 to write for any

ρ \in P_{a}^{2} (R^{3})

\begin{matrix} \forall φ \in C_{c}^{\infty} (R^{3}), \int_{R^{3}} ρ \nabla Ψ [ρ]^{⊥} \cdot \nabla φ = \int_{Ω} x^{⊥} \cdot \nabla φ (\nabla Φ [ρ]) \end{matrix}

(the integrals would be performed over

T^{3}

in the periodic case). The property

\nabla Φ [ρ]_{#} χ_{Ω} ℒ^{3} = ρ

is still valid when

ρ

is only a measure with finite second moment (see Remark 2 after Theorem 1.2 ). Therefore, the formulation on the right hand side extends unambiguously to the case where

ρ / \in L^{1} (R^{2})

Geometric interpretation This weak formulation allows has a natural geometric interpretation: at a point where

Ψ [ρ]

is not differentiable, and thus where

\partial Ψ [ρ]

is not reduced to a single point,

\nabla Ψ [ρ]

should be replaced by

\bar{\partial} Ψ [ρ]

the center of mass of the (convex) set

\partial Ψ [ρ]

This motivates the following definition of weak measure solutions

Definition 2.1. Let, for all

t \in [0, T]

ρ (t)

be a probability measure of

R^{3}

. It is said to be a weak measure solution to

S G

with initial data

ρ^{0}

if 1The time dependent probability measure

ρ

belongs to

C ([0, T], P - w *)

, 2there exists

t \to R (t)

non-decreasing such that for all

t \in [0, T]

ρ (t, \cdot)

is supported in

B (0, R (t))

, 3for all

T > 0

and for all

φ \in C_{c}^{\infty} ([0, T] \times R^{3})

we have

\begin{matrix} \int_{[0, T] \times R^{3}} \partial_{t} φ (t, x) d ρ (d t, x) \end{matrix}

(6)

\begin{matrix} + & \int_{[0, T] \times Ω} \nabla φ (t, \nabla Φ [ρ (t)] (x)) \cdot x^{⊥} d t d x - \int_{[0, T] \times R^{3}} \nabla φ (t, x) \cdot x^{⊥} d ρ (d t, x) \end{matrix}

\begin{matrix} = & \int φ (T, x) d ρ (T, x) d x - \int φ (0, x) d ρ^{0} (x) d x . \end{matrix}

This definition is consistent with the classical definition of weak solutions if for all

t

ρ (t, \cdot)

is absolutely continuous with respect to the Lebesgue measure.

2.2 Result

Here we prove the following

Theorem 2.2.

1. Let $ρ^{0}$ be a probability measure compactly supported. There exists a global weak measure solution to the system $S G$ with initial data $ρ^{0}$ in the sense of Definition 2.1 .
2. For any $T > 0$ , if $(ρ_{n})_{n \in N}$ is a sequence of weak measure solutions on $[0, T]$ to $S G$ with initial data $(ρ_{n}^{0})_{n \in N}$ , supported in $B_{R}$ for some $R > 0$ independent of $n$ , the sequence $(ρ_{n})_{n \in N}$ is precompact in $C ([0, T], P - w *)$ and every converging subsequence converges to a weak measure solution of $S G$ .

Proof of Theorem 2.2 We first show the weak stability of the formulation of Definition ( 2.1 ), and the compactness of weak measure solutions. We then use this result to obtain global existence of solutions to the Cauchy problem with initial data a bounded measure.

Weak stability of solutions We consider a sequence

(ρ_{n})_{n \in N}

of solutions of

S G

in the sense of Definition 2.1 . The sequence is uniformly compactly supported at time 0. We first show that there exists a non-decreasing function

R (t)

such that

ρ_{n} (t)

is supported in

B (R (t))

for all

t, n

Lemma 2.3. Let

ρ \in C ([0, T], P (R^{3}) - w *)

satisfy ( 6 ), let

ρ^{0} = ρ (t = 0)

be supported in

B (0, R^{0})

, then

ρ (t)

is supported in

B (0, R^{0} + C t)

C = {sup}_{y \in Ω} {| y |}

Proof. Consider any function

ξ_{ε} (t, r) \in C^{\infty} (R)

such that

\begin{matrix} ξ_{ε} (0, r) \equiv 1 i f - \infty < r \leq R^{0}, \end{matrix}

\begin{matrix} ξ_{ε} (0, r) \equiv 0 i f r \geq R^{0} + ε, \end{matrix}

\begin{matrix} ξ_{ε} (t, r) = ξ_{ε} (r - C t), \end{matrix}

with

ξ (0, \cdot)

non increasing. Then compute

\begin{matrix} \frac{d}{d t} \int ξ_{ε} (t, | x |) d ρ (t, x) \end{matrix}

\begin{matrix} = & - \int \partial_{r} ξ_{ε} (t, | x |) C d ρ (t, x) + \int_{Ω} \partial_{r} ξ_{ε} (t, | \nabla Φ [ρ (t)] |) \frac{\nabla Φ [ρ (t)]}{| \nabla Φ [ρ (t)] |} \cdot x^{⊥} d x \end{matrix}

\begin{matrix} \geq & \int_{Ω} \partial_{r} ξ_{ε} (t, | \nabla Φ [ρ (t)] |) (- C + | x |) d x \end{matrix}

\begin{matrix} \geq & 0 \end{matrix}

since, by definition, for

x \in Ω

| x | \leq C

and

ξ

is non increasing with respect to

r

. Note also that we have used

\int \nabla_{x} [ξ (t, | x |)] \cdot x^{⊥} d ρ (t, x) d x \equiv 0

. We know on the other hand that

\begin{matrix} \int_{R^{3}} ξ_{ε} (0, | x |) d ρ (0, x) = 1, \end{matrix}

\begin{matrix} \int_{R^{3}} ξ_{ε} (t, | x |) d ρ (t, x) \leq 1, \end{matrix}

therefore we conclude that

\int_{R^{3}} ξ_{ε} (| x |, t) d ρ (t, x) \equiv 1

, which concludes the lemma by letting

ε

go to 0.

□

From this lemma, we have:

| - \int_{[0, T] \times R^{3}} \nabla φ (t, x) \cdot x^{⊥} d ρ_{n} (d t, x) + \int_{[0, T] \times Ω} \nabla φ (t, \nabla Φ [ρ_{n} (t)] (x)) \cdot x^{⊥} d t d x |

\leq C (T) ∥ φ ∥_{L^{1} ([0, T], C^{1} (B_{R (T)})} .

Thus from Definition 2.1 equation ( 6 ) we know that for any time

t \geq 0

\partial_{t} ρ_{n} (t, \cdot)

is bounded in the dual of

L^{1} ([0, T], C^{1} (R^{3}))

and thus in the dual of

L^{1} ([0, T], W^{2, p} (R^{3}))

for

p > 3

by Sobolev embeddings. Thus for some

p^{'} > 1

we have

\begin{matrix} \partial_{t} ρ_{n} \in L^{\infty} ([0, T], W^{- 2, p^{'}} (R^{2})) . \end{matrix}

With the two above results, and using classical arguments of functional analysis (see [?] ), we can obtain the following lemma:

Lemma 2.4. Let the sequence

(ρ_{n})_{n \in N}

be as above, there exists

ρ \in C ([0, T], P - w *)

and a subsequence

(ρ_{n_{k}})_{k \in N}

, such that for all

t \in [0, T]

ρ_{n_{k}} (t)

converges to

ρ (t)

in the weak-

*

topology of measures.

With this lemma, we need to show that for all

φ \in C_{c}^{\infty} ([0, T] \times R^{3})

we have

\nabla φ (t, \nabla Φ [ρ_{n} (t)])

converging to

\nabla φ (t, \nabla Φ [ρ (t)])

whenever

ρ_{n} (t)

converges weakly-

*

ρ (t)

. This last step will be a consequence of the following stability theorem:

Theorem 2.5 (Brenier, [?] ). Let

Ω

be as above. Let

(ρ_{n})_{n \in N}

be a sequence of probability measures on

R^{d}

, such that

\forall n, \int (1 + | x |^{2}) d ρ_{n} \leq C

, let

Φ_{n} = Φ [ρ_{n}]

and

Ψ_{n} = Ψ [ρ_{n}]

be as in Theorem 1.1 .

If for any

f \in C^{0} (R^{d})

such that

| f (x) | \leq C (1 + | x |^{2}),

\int f ρ_{n} \to \int ρ f,

then

Φ_{n} \to Φ [ρ]

uniformly on each compact set of

Ω

and strongly in

W^{1, 1} (Ω; R^{d}),

and

Ψ_{n} \to Ψ [ρ]

uniformly on each compact set of

R^{d}

and strongly in

W_{l o c}^{1, 1} (R^{d})

From this result, we obtain that the sequence

\nabla Φ [ρ_{n}]

converges strongly in

L^{1} (Ω)

and almost everywhere (because of the convexity of

Φ [ρ]

) to

\nabla Φ [ρ]

. Thus

\nabla φ (t, \nabla Φ [ρ_{n}])

converges to

\nabla φ (t, \nabla Φ [ρ])

L^{1} (Ω)

and one can pass to the limit in the formulation of Definition 2.1 .

This ends the proof of point 2 of Theorem 2.2 .

Existence of solutions We show briefly the existence of a solution to the Cauchy problem in the sense of Definition 2.1 . Indeed given

ρ^{0}

the initial data for the problem that we want to solve, by smoothing

ρ^{0}

, we can take a sequence

ρ_{n}^{0}

of initial data belonging to

L^{1} (R^{2})

, uniformly compactly supported and converging weakly-

*

ρ^{0}

. We know already from [?] , [?] , [?] that for every

ρ_{n}^{0}

, one can build a global weak solution of ( 3 , 4 , 5 ), that will be uniformly compactly supported on

[0, T]

for all

T \geq 0

. This sequence will also be solution in the sense of Definition 2.1 . We then use the stability result, and conclude that, up to extraction of a subsequence, the sequence

ρ_{n}

converges in

C ([0, T], P - w *)

to a weak measure solution of

S G

with initial data

ρ^{0}

. This achieves the proof of Theorem 2.2 .

□

Remark: One can prove in fact the more general result, valid for non linear functionals:

Proposition 2.6. Let

F \in C^{0} (Ω \times R^{d})

, such that

| F (x, y) | \leq C (1 + | y |^{2})

, let

(ρ_{n})_{n \in N}

be a bounded sequence of probability measures, Lebesgue integrable, with finite second moment. Let

ρ

be a probability measure with finite second moment, such that for all

f \in C^{0} (R^{d})

such that

| f (x) | \leq C (1 + | x |^{2})

\int f d ρ_{n} \to \int f d ρ

. Then as

n

goes to

\infty

, we have

\begin{matrix} \int_{R^{d}} F (\nabla Ψ [ρ_{n}] (x), x) d ρ_{n} (x) = \int_{Ω} F (y, \nabla Φ [ρ_{n}] (y)) d y \end{matrix}

\begin{matrix} \to_{n} & \int_{Ω} F (y, \nabla Φ [ρ] (y)) d y : = \int_{R^{d}} F (\bar{\partial} Ψ [ρ] (x), x) d ρ (x) . \end{matrix}

3 Continuous solutions

What initial regularity is necessary in order to guarantee that the velocity fields remains Lipschitz, or that the flow remains continuous, at least for a short time ? The celebrated Youdovich’s Theorem for the Euler incompressible equation shows that when

d = 2

, if the initial vorticity data is bounded in

L^{\infty}

, the flow is Holder continuous, with Holder index decreasing to 0 as time goes to infinity. This proof relies on the following regularity property of the Poisson equation: if

Δ φ

is bounded in

L^{\infty}

, then

\nabla φ

is Log-Lipschitz. This continuity is enough to define a Holder continuous flow for the vector field

\nabla φ^{⊥}

. Such a result is not valid for the Monge-Ampere equation. As far as we know, the optimal regularity result for Monge-Ampere equations is the following:

3.1 Regularity of solutions to Monge-Ampere equation with Dini-continuous right hand side

Theorem 3.1 (Wang, [?] ). Let

u

be a strictly convex Alexandrov solution of

\begin{matrix} det D^{2} u = ρ \end{matrix}

(7)

with

ρ

strictly positive. If

w (r)

, the modulus of continuity of

ρ

, satisfies

\begin{matrix} \int_{0}^{1} \frac{w (r)}{r} d r < \infty, \end{matrix}

(8)

then

u

is in

C_{l o c}^{2}

We will work here in the periodic case. In this case,

u

the solution of ( 7 ) will be

Ψ [ρ]

of Theorem 1.2 . The arguments of [?] , [?] , adapted to the periodic case, show that

Ψ [ρ]

is indeed a strictly convex Alexandrov solution of solution of ( 7 ). Therefore we obtain the following corollary of Theorem 3.1 :

Corollary 3.2. Let

ρ \in P (T^{d})

be such that

\begin{matrix} 0 < m \leq ρ \leq M, \end{matrix}

\begin{matrix} \int_{0}^{1} \frac{w (r)}{r} d r = C < \infty . \end{matrix}

where

m, M, C

are positive constants. Let

Ψ [ρ]

be as in Theorem 1.2 . We have, for some constant

ℋ

depending on

m, M, C

\begin{matrix} ∥ Ψ [ρ] ∥_{C^{2} (T^{d})} \leq ℋ . \end{matrix}

3.2 Result

We will now prove the following:

Theorem 3.3. Let

ρ^{0}

be a probability on

T^{3}

, such that

ρ

is strictly positive and satisfies the continuity condition ( 8 ). Then there exists

T > 0

and

C_{1}, C_{2}

depending on

ρ^{0}

, such that on

[0, T]

there exists a solution

ρ (t, x)

S G

that satisfies for all

t \in [0, T]

\begin{matrix} \int_{0}^{1} \frac{w (t, r)}{r} d r \leq C_{1}, ∥ Ψ (t, \cdot) ∥_{C^{2} (T^{3})} \leq C_{2}, \end{matrix}

where

w (t, r)

is the modulus of continuity (in space) of

ρ (t, .)

Proof of Theorem 3.3 Let us first sketch the proof: If

Ψ \in C^{2}

, then the flow

t \to X (t, x)

generated by the velocity field

[\nabla Ψ (x) - x]^{⊥}

is Lipschitz in space. Since the flow is incompressible, we have

ρ (t, x) = ρ^{0} (X^{- 1} (t, x))

Now we use the following property: If two functions

f, g

have modulus of continuity respectively

w_{f}, w_{g}

then

g \circ f

has modulus

w_{g} \circ w_{f}

Thus if

X^{- 1} (t)

is Lipschitz, we have

w_{ρ^{0} \circ X^{- 1} (t)} \leq w_{ρ^{0}} (L \cdot)

with

L

the Lipschitz constant of

X^{- 1} (t)

and condition ( 8 ) remains satisfied.

Remark 1: Note that Holder continuous functions satisfy the condition ( 8 ).

Remark 2: Note also that we do not need any integrability on

\nabla ρ

and the solution of the Eulerian system still has to be understood in the distributional sense.

A fixed point argument Let us introduce the semi-norm

\begin{matrix} ∥ μ ∥_{C} = \int_{0}^{1} \frac{w_{μ} (r)}{r} d r \end{matrix}

(9)

defined on

P (T^{3})

, where we recall that

w_{μ}

is the modulus of continuity of

μ

. We denote

P_{C}

the set

P

equipped with this semi-norm, i.e.

\begin{matrix} P_{C} = {μ \in P (T^{3}), ∥ μ ∥_{C} < \infty} . \end{matrix}

From now, we fix

ρ^{0}

a probability density in

P_{C}

, satisfying

m \leq ρ^{0} \leq M

, where

m

and

M

are strictly positive constants. Let

μ

be a time dependent probability density in

L^{\infty} ([0, T]; P_{C})

, such that

m \leq μ (t) \leq M

for all

t

, we consider the solution

ρ

of the initial value problem:

\begin{matrix} \partial_{t} ρ + (\nabla Ψ [μ] (x) - x)^{⊥} \cdot \nabla ρ = 0, \end{matrix}

(10)

\begin{matrix} ρ (t = 0) = ρ^{0} . \end{matrix}

(11)

From Theorem 3.1 and its corollary, the vector field

v [μ] = (\nabla Ψ [μ] (x) - x)^{⊥}

C^{1}

uniformly in time, therefore there exists a unique solution to this equation, by Cauchy-Lipschitz Theorem. This solution can be built by the method of characteristics as follows: Consider the flow

X (t, x)

of the vector field

v [μ]

, then

ρ (t)

ρ^{0}

pushed forward by

X (t)

, i.e.

ρ (t) = ρ^{0} \circ X^{- 1} (t)

. From the incompressibility of

v [μ]

the condition

m \leq ρ^{0} \leq M

implies that for all

t \in [0, T]

m \leq ρ (t) \leq M

The initial data

ρ^{0}

being fixed, the map

μ \mapsto ρ

will be denoted by

ℱ

The spatial derivative of

X

D_{x} X

satisfies

\begin{matrix} \partial_{t} D_{x} X = D_{x} v [μ] (X) D_{x} X, \end{matrix}

therefore we have

\begin{matrix} | D_{x} X (t) | \leq exp (t {sup}_{s \in [0, t]} | D_{x} v [μ] (s) |), \end{matrix}

and the same bound holds for

X (t)^{- 1}

. Since

w_{f \circ g} \leq w_{f} \circ w_{g}

, and writing

C_{t} = exp (t {sup}_{s \in [0, t]} | D_{x} v [μ] |)

, we obtain

w_{ρ (t)} (\cdot) \leq w_{ρ^{0}} (C_{t} \cdot)

, and

\begin{matrix} \int_{0}^{1} \frac{w_{ρ (t)} (r)}{r} d r & \leq & \int_{0}^{C_{t}} \frac{w_{ρ^{0}} (r)}{r} d r \end{matrix}

\begin{matrix} \leq & \int_{0}^{1} \frac{w_{ρ^{0}} (r)}{r} d r + (M - m) (C_{t} - 1), \end{matrix}

(using that

\forall r, w_{ρ} (r) \leq M - m

). Therefore,

∥ ρ (t) ∥_{C} \leq ∥ ρ^{0} ∥_{C} + (M - m) (C_{t} - 1) .

Now from Corollary 3.2 , and

m, M

being fixed, there exists a non-decreasing function

ℋ

such that

\begin{matrix} ∥ v [μ] ∥_{C^{1}} \leq ℋ (∥ μ ∥_{C}), \end{matrix}

and so

C_{t} \leq exp (t ℋ (∥ μ ∥_{L^{\infty} ([0, t]; P_{C}})))

. Hence we can chose

Q > 1

, and then

T

such that

\begin{matrix} ∥ ρ^{0} ∥_{C} + (M - m) (exp (T ℋ (Q ∥ ρ^{0} ∥_{C})) - 1) = Q ∥ ρ^{0} ∥_{C} . \end{matrix}

Note that for

Q > 1

, we necessarily have

T > 0

. Then the map

ℱ : μ \mapsto ρ

goes now from

A = {μ, ∥ μ ∥_{L^{\infty} ([0, T]; P_{C})} \leq Q ∥ ρ^{0} ∥_{C}, m \leq μ \leq M}

into

ℬ = {ρ, ∥ ρ (t) ∥_{C} \leq ∥ ρ^{0} ∥_{C} + (M - m) (exp (t ℋ (Q ∥ ρ^{0} ∥_{C})) - 1), \forall t \in [0, T]},

and with our choice of

T = T (Q)

, we have

ℬ \subset A

. Moreover from the unconditional bounds

\begin{matrix} ρ \leq M, \end{matrix}

\begin{matrix} ∥ v [μ] ∥_{L^{\infty} ([0, T] \times T^{3})} \leq \sqrt{3} / 2, \end{matrix}

(see the remark after Theorem 1.2 for the second bound) and using equation ( 10 ), we have also

∥ \partial_{t} ρ ∥_{L^{\infty} ([0, T]; W^{- 1, \infty})} \leq K (M)

whenever

ρ = ℱ (μ)

Call

\tilde{A}

(resp.

\tilde{ℬ}

) the set

A \cap {ρ, ∥ \partial_{t} ρ ∥_{L^{\infty} ([0, T]; W^{- 1, \infty})} \leq K (M)}

, (resp.

ℬ \cap {ρ, ∥ \partial_{t} ρ ∥_{L^{\infty} ([0, T]; W^{- 1, \infty})} \leq K (M)}

); we claim that

∙ $ℱ (\tilde{A}) \subset \tilde{ℬ} \subset \tilde{A}$ ,
∙ $\tilde{A}$ is convex and compact for the $C^{0} ([0, T] \times T^{3})$ topology,
∙ $ℱ$ is continuous for this topology,

so that we can apply the Schauder fixed point Theorem. We only check the last point, the second being a classical result of functional analysis. So let us consider a sequence

(μ_{n})_{n \in N}

converging to

μ \in A

, and the corresponding sequence

(ρ_{n} = ℱ (μ_{n}))_{n \in N}

. The sequence

ρ_{n}

is pre-compact in

C^{0} ([0, T] \times T^{3})

, from the previous point, and we see (with the stability Theorem 2.5 ) that it converges to a solution

ρ

\begin{matrix} \partial_{t} ρ + \nabla \cdot (ρ v [μ]) = 0 . \end{matrix}

But,

v [μ]

being Lipschitz, this solution is unique, and therefore

ℱ (μ_{n})

converges to

ℱ (μ)

, which proves the continuity of

ℱ

, and ends the proof of existence by the Schauder fixed point Theorem.

□

We state here without proof some consequences of the previous result:

Corollary 3.4. Let

ρ^{0} \in P (T^{3})

, such that

0 < m \leq ρ \leq M

1. If $ρ^{0} \in C^{α}, α \in] 0, 1]$ , for $T^{*} > 0$ depending on $ρ^{0}$ , a solution $ρ (t, x)$ to ( 3 , 4 , 5 ) exists in $L^{\infty} ([0, T^{*} [, C^{α} (T^{3}))$ .
2. If $ρ^{0} \in W^{1, p}, p > 3$ , for $T^{*} > 0$ depending on $ρ^{0}$ , a solution $ρ (t, x)$ to ( 3 , 4 , 5 ) exists in $L^{\infty} ([0, T [, W^{1, p} (T^{3}))$ .
3. If $ρ^{0} \in C^{k, α}, α \in] 0, 1], k \in N$ , for $T^{*} > 0$ depending on $ρ^{0}$ , a solution $ρ (t, x)$ to ( 3 , 4 , 5 ) exists in $L^{\infty} ([0, T^{*} [, C^{k, α} (T^{3}))$ .

Moreover, for these solutions, the velocity field is respectively in

C^{1, α} (T^{3})

W^{2, p} (T^{3})

, and

C^{k + 1, α} (T^{3})

[0, T^{*} [

4 Uniqueness of solutions to $S G$ with Holder continuous densities

4.1 Result

Here we prove the following theorem:

Theorem 4.1. Suppose that

ρ^{0} \in P (T^{3})

with

0 < m \leq ρ^{0} \leq M

, and belongs to

C^{α} (T^{3})

for some

α > 0

. From Theorem 3.3 , for some

T > 0

there exists a solution

\bar{ρ}

S G

L^{\infty} ([0, T], C^{α} (T^{3}))

Then every solution of

S G

L^{\infty} ([0, T^{'}], C^{β} (T^{3}))

for

T^{'} > 0, β > 0

with same initial data coincides with

\bar{ρ}

[0, inf {T, T^{'}}]

Remark 1: The uniqueness of weak solutions is still an open question.

Remark 2: Our proof of uniqueness is thus valid in a smaller class of solutions than the one found in the previous section, the reason is the following:

during the course of the proof, we will need to solve a Monge-Ampere equation, whose right-hand side is a function of the second derivatives of the solution of another Monge-Ampere equation. In Theorem 3.1 , if

u

is solution to ( 7 ) with a right hand side satisfying ( 8 ), although

u \in C^{2}

, it is not clear that the second derivatives of

u

satisfy ( 8 ). Actually, it is even known to be wrong in the case of the Laplacian (for a precise discussion on the subject, the reader may refer to [?] ). However, from Theorem 4.3 below, if

ρ \in C^{α}

then

u \in C^{2, α}

What we actually need is a continuity condition on the right hand side of ( 7 ) such that the second derivative of the solution

u

satisfies ( 8 ). This may be a weaker condition than Holder continuity, however the proof would not be affected, therefore it is enough to give it under the present form.

Proof of Theorem 4.1 Let

ρ_{1}

and

ρ_{2}

be two solutions of ( 3 , 4 , 5 ), in

L^{\infty} ([0, T], C^{β} (T^{3}))

that coincide at time 0. Let

X_{1}, X_{2}

be the two corresponding Lagrangian solutions, (i.e.

solutions of ( 1 , 2 )). The velocity field being

C^{1}

, for all

t \in [0, T]

X_{1} (t, \cdot)

and

X_{2} (t, \cdot)

are both

C^{1}

diffeomorphisms of

T^{d}

We call

v_{1}

(resp.

v_{2}

) the velocity field associated to

X_{1}

(resp.

X_{2}

v_{i} (t, x) = [\nabla Ψ_{i} (t, x) - x]^{⊥}, i = 1, 2

. We have

\begin{matrix} \partial_{t} (X_{1} - X_{2}) & = & v_{1} (X_{1}) - v_{2} (X_{2}) \end{matrix}

\begin{matrix} = & (v_{1} (X_{1}) - v_{1} (X_{2})) + (v_{1} (X_{2}) - v_{2} (X_{2})) . \end{matrix}

We want to obtain a Gronwall type inequality for

∥ X_{1} - X_{2} ∥_{L^{2}}

. Since

v_{1}

is uniformly Lipschitz in space (from Theorem 3.3 ), the first bracket is estimated in

L^{2}

norm by

C ∥ X_{1} - X_{2} ∥_{L^{2}}

We now need to estimate the second term. We first have that

\begin{matrix} \int | v_{1} (X_{2}) - v_{2} (X_{2}) |^{2} = \int ρ_{2} | \nabla Ψ_{1} - \nabla Ψ_{2} |^{2}, \end{matrix}

and since

ρ_{2}

is bounded, we need to estimate

∥ \nabla Ψ_{1} - \nabla Ψ_{2} ∥_{L^{2}}

. This will be done in the following Proposition:

Proposition 4.2. Let

X_{1}, X_{2}

be mappings from

T^{d}

into itself, such that the densities

ρ_{i} = X_{i #} d x, i = 1, 2

are in

C^{α} (T^{d})

for some

α > 0

, and satisfy

0 < m \leq ρ_{i} \leq M

. Let

Ψ_{i}, i = 1, 2

be convex such that

\begin{matrix} det D^{2} Ψ_{i} = ρ_{i} \end{matrix}

in the sense of Theorem 1.1 , i.e.

Ψ_{i} = Ψ [ρ_{i}]

. Then

∥ \nabla Ψ_{1} - \nabla Ψ_{2} ∥_{L^{2}} \leq C ∥ X_{1} - X_{2} ∥_{L^{2}},

where

C

depends on

α

(the Holder index of

ρ_{i}

∥ ρ_{i} ∥_{C^{α} (T^{d})}

m

and

M

Before giving a proof of this result, we conclude the proof of the Theorem 4.1 . The Proposition 4.2 implies immediately that

\begin{matrix} ∥ \partial_{t} (X_{1} - X_{2}) ∥_{L^{2}} \leq C ∥ X_{1} - X_{2} ∥_{L^{2}}, \end{matrix}

and we conclude the proof of the Theorem by a standard Gronwall lemma.

□

4.2 Energy estimates along Wasserstein geodesics: Proof of Proposition 4.2 .

In the proof of this result we will need the following result on optimal transportation of measures by gradient of convex functions:

Theorem 4.3 (Brenier, [?] , McCann, [?] , Cordero-Erausquin, [?] , Caffarelli,[?] ).

Let

ρ_{1}

ρ_{2}

be two probability measures on

T^{d}

, such that

ρ_{1}

is absolutely continuous with respect to the Lebesgue measure.

1. There exists a unique up to a constant convex function $φ$ such that $φ - | \cdot |^{2} / 2$ is $Z^{d}$ periodic, satisfying $\nabla φ_{#} ρ_{1} = ρ_{2}$ .

2. The map

\nabla φ

is the solution of the minimization problem

\begin{matrix} {inf}_{T_{#} ρ_{1} = ρ_{2}} \int_{T^{d}} ρ_{1} (x) | T (x) - x |_{T^{d}}^{2} d x, \end{matrix}

(12)

and for all

x \in R^{d}

| \nabla φ (x) - x |_{T^{d}} = | \nabla φ (x) - x |_{R^{d}}

3. If $ρ_{1}$ , $ρ_{2}$ are strictly positive and belong to $C^{α} (T^{d})$ for some $α > 0$ then $φ \in C^{2, α} (T^{d})$ and satisfies pointwise
$\begin{matrix} ρ_{2} (\nabla φ) det D^{2} φ = ρ_{1} . \end{matrix}$

For complete references on the optimal transportation problem ( 12 ) and its applications, the reader can refer to [?] .

Remark 1: the expression

| \cdot |_{T^{d}}

denotes the Riemannian distance on the flat torus, whereas

| \cdot |_{R^{d}}

is the Euclidian distance on

R^{d}

. The second assertion of point 2 means that, for all

x \in R^{d}

| \nabla φ (x) - x | \leq d i a m (T^{d}) = \sqrt{d} / 2

Remark 2: Here again, note that since

φ - | \cdot |^{2} / 2

is periodic, the map

x \mapsto \nabla φ (x)

is compatible with the equivalence classes of

R^{d} / Z^{d}

, and therefore is defined without ambiguity on

T^{d}

Wasserstein geodesics between probability measures In this part we use results from [?] , [?] . Using Theorem 4.3 , we consider the unique (up to a constant) convex potential

φ

such that

\begin{matrix} \nabla φ_{#} ρ_{1} = ρ_{2}, \end{matrix}

\begin{matrix} φ - | \cdot |^{2} / 2 is Z^{d} - periodic . \end{matrix}

We consider, for

θ \in [1, 2]

φ_{θ}

defined by

\begin{matrix} φ_{θ} = (2 - θ) \frac{| x |^{2}}{2} + (θ - 1) φ . \end{matrix}

We also consider, for

θ \in [1, 2]

ρ_{θ}

defined by

\begin{matrix} ρ_{θ} = \nabla φ_{θ #} ρ_{1} . \end{matrix}

Then

ρ_{θ}

interpolates between

ρ_{1}

and

ρ_{2}

. This interpolation has been introduced in [?] and [?] as the time continuous formulation of the Monge-Kantorovitch mass transfer. In this construction, a velocity field

v_{θ}

is defined

d ρ_{θ}

a.e. as follows:

\begin{matrix} \forall f \in C^{0} (T^{d}; R^{d}), \int ρ_{θ} v_{θ} \cdot f = \int ρ_{1} f (\nabla φ_{θ}) \cdot \partial_{θ} \nabla φ_{θ} . \end{matrix}

(13)

It is easily checked that the pair

ρ_{θ}, v_{θ}

satisfies

\begin{matrix} \partial_{θ} ρ_{θ} + \nabla \cdot (ρ_{θ} v_{θ}) = 0, \end{matrix}

and for any

θ \in [1, 2]

, we have (see [?] ):

\begin{matrix} \frac{1}{2} \int_{T^{d}} ρ_{θ} | v_{θ} |^{2} = \frac{1}{2} \int_{T^{d}} ρ_{1} | \nabla φ (x) - x |^{2} = W_{2}^{2} (ρ_{1}, ρ_{2}), \end{matrix}

where

W_{2} (ρ_{1}, ρ_{2})

is the Wasserstein distance between

ρ_{1}

and

ρ_{2}

, defined by

\begin{matrix} W_{2}^{2} (ρ_{1}, ρ_{2}) = {inf}_{T_{#} ρ_{1} = ρ_{2}} {\int ρ_{1} (x) | T (x) - x |_{T^{d}}^{2}} . \end{matrix}

The Wasserstein distance can also be formulated as follows:

\begin{matrix} W_{2}^{2} (ρ_{1}, ρ_{2}) = {inf}_{Y_{1}, Y_{2}} {\int_{T^{d}} | Y_{1} - Y_{2} |_{T^{d}}^{2}} \end{matrix}

where the infimum is performed over all maps

Y_{1}, Y_{2} : T^{d} \mapsto T^{d}

such that

Y_{i #} d x = ρ_{i}, i = 1, 2

. From this definition we have easily

\begin{matrix} W_{2}^{2} (ρ_{1}, ρ_{2}) \leq \int | X_{2} (t, a) - X_{1} (t, a) |^{2} d a, \end{matrix}

and it follows that, for every

θ \in [1, 2]

\begin{matrix} \int_{T^{d}} ρ_{θ} | v_{θ} |^{2} = W_{2}^{2} (ρ_{1}, ρ_{2}) \leq ∥ X_{2} - X_{1} ∥_{L^{2}} . \end{matrix}

(14)

Regularity of the interpolant measure

ρ_{θ}

From Theorem 4.3 , for

ρ_{1}, ρ_{2} \in C^{β}

and pinched between the positive postive constants

m

and

M

, we know that

φ \in C^{2, β}

and satisfies

\begin{matrix} det D^{2} φ = \frac{ρ_{1}}{ρ_{2} (\nabla φ)} . \end{matrix}

We now estimate

ρ_{θ} = ρ_{1} [det D^{2} φ_{θ}]^{- 1}

. From the concavity of

log (det (\cdot))

on symmetric positive matrices, we have

\begin{matrix} det D^{2} φ_{θ} & = & det ((2 - θ) I + (θ - 1) D^{2} φ) \end{matrix}

\begin{matrix} \geq & [det D^{2} φ]^{θ - 1} \end{matrix}

\begin{matrix} \geq & \frac{m}{M} . \end{matrix}

Moreover, since

φ \in C^{2}

det D^{2} φ_{θ}

is bounded by above. Thus

ρ_{θ}

is uniformly bounded away from 0 and infinity, and uniformly Holder continuous.

Final energy estimate If we consider, for every

θ \in [1, 2]

Ψ_{θ}

solution of

\begin{matrix} det D^{2} Ψ_{θ} = ρ_{θ}, \end{matrix}

(15)

in the sense of Theorem 1.2 , then

Ψ_{θ}

interpolates between

Ψ_{1}

and

Ψ_{2}

, and

Ψ_{θ} \in C^{2, β}

uniformly, from the regularity of

ρ_{θ}

. We will estimate

\partial_{θ} \nabla Ψ_{θ}

by differentiating ( 15 ) with respect to

θ

: for

M, N

two

d \times d

matrices,

t \in R

, we recall that

det (M + t N) = det M + t (t r a c e M_{c o}^{t} N) + o (t),

where

M_{c o}

is the co-matrix (or matrix of cofactors) of

M

. Moreover, for any

f \in C^{2} (R^{d}; R)

, if

M

is the co-matrix of

D^{2} f

, it is a common fact that

\begin{matrix} \forall j \in [1 . . d], \sum_{i = 1}^{d} \partial_{i} M_{i j} \equiv 0 . \end{matrix}

(16)

Hence, denoting

M_{θ}

the co-matrix of

D^{2} Ψ_{θ}

, we obtain that

\partial_{θ} Ψ_{θ}

satisfies

\begin{matrix} \nabla \cdot (M_{θ} \nabla \partial_{θ} Ψ_{θ}) & = & \partial_{θ} ρ_{θ} (t) \end{matrix}

\begin{matrix} = & - \nabla \cdot (ρ_{θ} v_{θ}) . \end{matrix}

(17)

From the

C^{2, β}

regularity of

Ψ_{θ}

D^{2} Ψ_{θ}

is a

C^{β}

smooth, positive definite matrix, and its co-matrix as well. Thus the problem ( 17 ) is uniformly elliptic.

If we multiply by

\partial_{θ} Ψ_{θ}

, and integrate by parts we obtain

\begin{matrix} \int \nabla^{t} \partial_{θ} Ψ_{θ} M_{θ} \nabla \partial_{θ} Ψ_{θ} = - \int \nabla \partial_{θ} Ψ_{θ} \cdot v_{θ} ρ_{θ} . \end{matrix}

Using that

M_{θ} \geq λ I

for some

λ > 0

, and combining with the inequality ( 14 ) above, we obtain

\begin{matrix} ∥ \nabla \partial_{θ} Ψ_{θ} (t) ∥_{L^{2}} & \leq & λ^{- 1} ∥ ρ_{θ} v_{θ} ∥_{L^{2}} \end{matrix}

\begin{matrix} \leq & λ^{- 1} ∥ X_{2} - X_{1} ∥_{L^{2}} {({sup}_{θ} ∥ ρ_{θ} ∥_{L^{\infty}})}^{1 / 2} . \end{matrix}

The constant

λ^{- 1}

depends on

m, M, β, {∥ ρ_{i} ∥_{C^{β}}, i = 1, 2}

, and is thus bounded under our present assumptions. We have already seen that

ρ_{θ}

is uniformly bounded, and we finally obtain that

\begin{matrix} ∥ \nabla Ψ_{1} - \nabla Ψ_{2} ∥_{L^{2}} \leq C ∥ X_{1} - X_{2} ∥_{L^{2}}, \end{matrix}

(18)

this ends the proof of Proposition 4.2 .

□

Remark 1. In [?] , the author obtains also (weaker) estimates of the type of Proposition 4.2 , for discontinuous densities

ρ_{1}, ρ_{2}

5 Uniqueness of solutions to the 2-d Euler equations with bounded vorticity: a new proof

This proof adapts easily to the case of 2-d Euler equation with bounded vorticity, giving a new proof of the uniqueness part in Youdovich’s theorem.

We start now from the following system:

\begin{matrix} \partial_{t} ρ + \nabla ψ^{⊥} \cdot \nabla ρ = 0, \end{matrix}

(19)

\begin{matrix} ρ = Δ ψ, \end{matrix}

(20)

\begin{matrix} ρ (t = 0) = ρ^{0} . \end{matrix}

(21)

For simplicity, we restrict ourselves to the periodic case, i.e.

x \in T^{d}

ρ, ψ

periodic, this implies that

ρ

has total mass equal to 0. We reprove the following classical result:

Theorem 5.1 (Youdovich, [?] ). Given an initial data

ρ^{0} \in L^{\infty} (T^{2})

satisfying

\int_{T^{2}} ρ^{0} = 0

, there exists a unique solution to ( 19 , 20 , 21 ) such that

ρ

belongs to

L_{l o c}^{\infty} (R^{+} \times T^{2})

Proof of Theorem 5.1 We consider two solutions

ρ_{1}, ψ_{1}

and

ρ_{2}, ψ_{2}

, such that

ρ_{i}, i = 1, 2

are bounded in

L^{\infty} ([0, T] \times T^{d})

. In this case the velocity fields

v_{i} = \nabla ψ_{i}^{⊥}

both satisfy

\begin{matrix} \forall (x, y) \in T^{2}, | x - y | \leq \frac{1}{2}, | v_{i} (x) - v_{i} (y) | \leq C | x - y | log \frac{1}{| x - y |} . \end{matrix}

This implies that the flows

(t, x) \mapsto X_{i} (t, x)

associated to the velocity fields

v_{i} = \nabla ψ_{i}^{⊥}

are Holder continuous, and measure preserving. Moreover, one has, for all

t \in [0, T]

ρ_{i} (t) = X_{i} (t)_{#} ρ^{0}

Applying the same technique as before, we need to estimate

∥ \nabla ψ_{1} - \nabla ψ_{2} ∥_{L^{2} (T^{2})}

in terms of

∥ X_{1} - X_{2} ∥_{L^{2} (T^{2})}

. In the present case, the energy estimate of Proposition 4.2 will hold under the weaker assumptions that the two densities are bounded.

Proposition 5.2. Let

X_{1}, X_{2}

be mappings from

T^{d}

into itself, let

ρ^{0}

be a bounded measure with a density in

L^{\infty}

with respect to the Lebsgue measure, and with

\int_{T^{d}} ρ^{0} = 0

. Let

ρ_{i} = X_{i #} ρ^{0}, i = 1, 2

. Let

ψ_{i}, i = 1, 2

be periodic solutions of

Δ ψ_{i} = ρ_{i}, i = 1, 2

, then we have

\begin{matrix} ∥ \nabla ψ_{1} - \nabla ψ_{2} ∥_{L^{2} (T^{d})} \leq {(2 ∥ ρ^{0} ∥_{L^{\infty}} max {∥ ρ_{1} ∥_{L^{\infty}}, ∥ ρ_{2} ∥_{L^{\infty}}})}^{1 / 2} ∥ X_{1} - X_{2} ∥_{L^{2} (T^{d})} . \end{matrix}

Remark: In other words, this proposition shows that for

ρ_{1}, ρ_{2}

bounded, the

H^{- 1}

norm of

ρ_{1} - ρ_{2}

is controlled by some ’generalized’ (since here we have unsigned measures) Wasserstein distance between

ρ_{1}

and

ρ_{2}

To conclude the proof of Theorem 5.1 , note first that for all

C > 0

, we can take

T

small enough so that

∥ X_{2} - X_{1} ∥_{L^{\infty} ([0, T] \times T^{2})} \leq C

. Now we have for the difference

X_{1} - X_{2}

, as long as

| X_{1} - X_{2} | \leq 1 / 2

\begin{matrix} ∥ \partial_{t} (X_{1} - X_{2}) ∥_{L^{2}} \end{matrix}

\begin{matrix} \leq & ∥ \nabla ψ_{1} (X_{1}) - \nabla ψ_{1} (X_{2}) ∥_{L^{2}} + ∥ \nabla ψ_{1} (X_{2}) - \nabla ψ_{2} (X_{2}) ∥_{L^{2}} \end{matrix}

\begin{matrix} \leq & C_{1} ∥ | X_{1} - X_{2} | log (| X_{1} - X_{2} |) ∥_{L^{2}} + C_{2} ∥ X_{1} - X_{2} ∥_{L^{2}}, \end{matrix}

where, to evaluate the second term of the second line, we have used the fact that

\begin{matrix} ∥ \nabla ψ_{1} (X_{2}) - \nabla ψ_{2} (X_{2}) ∥_{L^{2}} = ∥ \nabla ψ_{1} - \nabla ψ_{2} ∥_{L^{2}}, \end{matrix}

and then applied Proposition 5.2 .

We just need to evaluate

∥ | X_{1} - X_{2} | log (| X_{1} - X_{2} |) ∥_{L^{2}}

. We take

T

small enough so that

∥ X_{2} - X_{1} ∥_{L^{\infty} ([0, T] \times T^{2})} \leq 1 / e

and notice that

x \mapsto x {log}^{2} x

is concave for

0 \leq x \leq 1 / e

, therefore by Jensen’s inequality we have

\begin{matrix} \int_{T^{2}} | X_{2} - X_{2} |^{2} {log}^{2} (| X_{1} - X_{2} |) \end{matrix}

\begin{matrix} = & \frac{1}{4} \int_{T^{2}} | X_{2} - X_{2} |^{2} {log}^{2} (| X_{1} - X_{2} |^{2}) \end{matrix}

\begin{matrix} \leq & \frac{1}{4} \int_{T^{2}} | X_{2} - X_{1} |^{2} {log}^{2} (\int_{T^{2}} | X_{2} - X_{1} |^{2}), \end{matrix}

and some elementary computations finally yield

\begin{matrix} \partial_{t} ∥ X_{2} - X_{1} ∥_{L^{2}} \leq C ∥ X_{2} - X_{1} ∥_{L^{2}} log \frac{1}{∥ X_{2} - X_{1} ∥_{L^{2}}} . \end{matrix}

The conclusion

X_{1} \equiv X_{2}

follows then by standard arguments.

5.1 Energy estimates along Wasserstein geodesic: Proof of Proposition 5.2

The proof of this proposition is very close to the proof of Proposition 4.2 , and we will only sketch it, insisting on the specific points. Here the densities

ρ_{i}

can not be of constant sign, since their mean value is zero, hence we introduce

ρ^{0, +}

(resp.

ρ^{0, -}

) the positive (resp. negative) part of

ρ^{0}

. Then we introduce

ρ_{i}^{\pm} = X_{i #} ρ^{0, \pm}

. Note that if the mappings

X_{i}

were injective, (which is the case in our present situation) we would have

ρ_{i}^{\pm}

that coincides with the positive/negative parts of

ρ_{i}

, but this can be wrong if

X_{i}

is not injective. However what remains is that

ρ_{i} = ρ_{i}^{+} - ρ_{i}^{-}

. Now,

ρ_{i}^{\pm}, i = 1, 2

are 4 positive measures of total mass equal to say

M

, with

M < \infty

Wasserstein geodesic We interpolate between the positive parts

ρ_{i}^{+}

, and the negative part is handled in the same way. As before we introduce the density

ρ_{θ}^{+} (t)

that interpolates between

ρ_{1}^{+} (t)

and

ρ_{2}^{+} (t)

. In this interpolation, we consider

v_{θ}^{+}

such that

\begin{matrix} \partial_{θ} ρ_{θ}^{+} + \nabla \cdot (ρ_{θ}^{+} v_{θ}^{+}) = 0, \end{matrix}

and we introduce as well

ρ_{θ}^{-}, v_{θ}^{-}

. Then

ρ_{θ} = ρ_{θ}^{+} - ρ_{θ}^{-}

has mean value 0. Let the potential

ψ_{θ}

be solution to

\begin{matrix} Δ ψ_{θ} = ρ_{θ} . \end{matrix}

(22)

Note that

ρ_{θ}

has mean value zero therefore this equation is well posed on

T^{2}

, moreover

ψ_{θ}

interpolates between

ψ_{1}

and

ψ_{2}

Bound on the interpolant measure

ρ_{θ}

Instead of interpolating between two smooth densities, we interpolate between bounded densities, and use the following result from [?] :

Proposition 5.3 (McCann, [?] ). Let

ρ_{θ}^{+}

be the Wasserstein geodesic linking

ρ_{1}^{+}

ρ_{2}^{+}

defined above. Then, for all

θ \in [1, 2]

\begin{matrix} ∥ ρ_{θ}^{+} ∥_{L^{\infty}} \leq max {∥ ρ_{1}^{+} ∥_{L^{\infty}}, ∥ ρ_{2}^{+} ∥_{L^{\infty}}} . \end{matrix}

The same holds for

ρ_{i}^{-}, ρ_{θ}^{-}

Remark: This property is often referred to as ’displacement convexity’.

Energy estimates Now by differentiating ( 22 ) with respect to

θ

, we obtain

\begin{matrix} Δ \partial_{θ} ψ_{θ} = \partial_{θ} ρ_{θ} = - \nabla \cdot (ρ_{θ}^{+} v_{θ}^{+} - ρ_{θ}^{-} v_{θ}^{-}), \end{matrix}

(23)

with

v_{θ}^{\pm}

the interpolating velocity defined as in ( 13 ), and satisfying for all

θ \in [1, 2]

\begin{matrix} \int ρ_{θ}^{\pm} | v_{θ}^{\pm} |^{2} = W_{2}^{2} (ρ_{1}^{\pm} (t), ρ_{2}^{\pm} (t)) . \end{matrix}

Multiplying ( 23 ) by

\partial_{θ} ψ_{θ}

, and integrating over

θ \in [1, 2]

, we obtain

\begin{matrix} ∥ \nabla ψ_{1} - \nabla ψ_{2} ∥_{L^{2} (T^{d})} & \leq & \int_{θ = 1}^{2} ∥ ρ_{θ}^{+} v_{θ}^{+} ∥_{L^{2}} + ∥ ρ_{θ}^{-} v_{θ}^{-} ∥_{L^{2}} \end{matrix}

\begin{matrix} \leq & W_{2} (ρ_{1}^{+}, ρ_{2}^{+}) {({sup}_{θ} ∥ ρ_{θ}^{+} ∥_{L^{\infty}})}^{1 / 2} \end{matrix}

\begin{matrix} + & W_{2} (ρ_{1}^{-}, ρ_{2}^{-}) {({sup}_{θ} ∥ ρ_{θ}^{-} ∥_{L^{\infty}})}^{1 / 2} . \end{matrix}

Note that the energy estimate is easier here than in the Monge-Ampere case, since the problem is immediately uniformly elliptic.

The mappings

X_{i}

satisfy

X_{i #} ρ_{0} = ρ_{i}

, and

X_{i #} (ρ_{0}^{\pm}) = ρ_{i}^{\pm} .

Hence,

\begin{matrix} W_{2}^{2} (ρ_{1}^{\pm}, ρ_{2}^{\pm}) & \leq & \int ρ_{0}^{\pm} | X_{1} - X_{2} |^{2} . \end{matrix}

Using Proposition 5.3 , we conclude:

\begin{matrix} ∥ \nabla ψ_{1} - \nabla ψ_{2} ∥_{L^{2} (T^{d})} \end{matrix}

\begin{matrix} \leq & 2 ∥ X_{2} - X_{1} ∥_{L^{2}} {(∥ ρ^{0} ∥_{L^{\infty}} max {∥ ρ_{1} ∥_{L^{\infty}}, ∥ ρ_{2} ∥_{L^{\infty}}})}^{1 / 2} . \end{matrix}

This ends the proof of Proposition 5.2 . Note that in our specific case,

X_{i}

are Lebesgue measure preserving invertible mappings, therefore

∥ ρ_{i}^{\pm} ∥_{L^{\infty}} = ∥ ρ_{0}^{\pm} ∥_{L^{\infty}}

, and the estimate can be simplified in

\begin{matrix} ∥ \nabla ψ_{1} - \nabla ψ_{2} ∥_{L^{2} (T^{d})} \leq 2 ∥ ρ_{0} ∥_{L^{\infty}} ∥ X_{2} - X_{1} ∥_{L^{2} (T^{d})} . \end{matrix}

□

Remark: This technique can be used to conclude uniqueness for many non-linear systems, where a transport equation and an elliptic equation are coupled. The velocity field is the gradient of a potential satisfying a elliptic equation whose right hand side depends smoothly on the density.

For example, we have uniqueness of solutions to the Vlasov-Poisson system and Euler-Poisson system with bounded density and bounded velocity. The Vlasov-Monge-Ampere and Euler-Monge-Ampere systems have also been studied by the author ([?] , [?] ), and the same technique apply to yield uniqueness for solutions with

C^{α}

density and bounded velocity. Note however that to enforce uniform ellipticity, we need for the Monge-Ampere equation the density to be bounded by below which is not the case for the Poisson equation.

6 Convergence to the Euler equation

6.1 Scaling of the system

Here we present a rescaled version of the 2-d

S G

system and some formal arguments to motivate the next convergence results. Here

x \in T^{2}, t \in R^{+}

and for

v = (v_{1}, v_{2}) \in R^{2}

v^{⊥}

now means

(- v_{2}, v_{1})

. Introducing

ψ [ρ] = Ψ [ρ] - | x |^{2} / 2

, where

Ψ [ρ]

is given by Theorem 1.2 , the periodic 2-d

S G

system now reads

\begin{matrix} \partial_{t} ρ + \nabla \cdot (ρ \nabla ψ^{⊥}) = 0, \end{matrix}

\begin{matrix} det (I + D^{2} ψ) = ρ . \end{matrix}

ρ

is close to one then

ψ

should be small, and therefore one may consider the linearization

det (I + D^{2} ψ) = 1 + Δ ψ + O (| D^{2} ψ |^{2})

, that yields

Δ ψ ≃ ρ - 1

Thus for small initial data, i.e.

ρ^{0} - 1

small, one expects

ψ, μ = ρ - 1

to stay close to a solution of the Euler incompressible equation

E I

\begin{matrix} \partial_{t} \bar{ρ} + \nabla \cdot (\bar{ρ} \nabla {\bar{φ}}^{⊥}) = 0, \end{matrix}

(24)

\begin{matrix} Δ \bar{φ} = \bar{ρ} . \end{matrix}

(25)

We shall rescale the equation, in order to consider quantities of order one.

We introduce the new unknown

\begin{matrix} ρ^{ε} (t, x) = \frac{1}{ε} (ρ (\frac{t}{ε}, x) - 1), \end{matrix}

\begin{matrix} ψ^{ε} (t, x) = \frac{1}{ε} ψ (\frac{t}{ε}, x) . \end{matrix}

Then we have

\begin{matrix} ρ (t) = 1 + ε ρ^{ε} (ε t), \end{matrix}

\begin{matrix} Ψ [ρ] (t) = | x |^{2} / 2 + ε ψ^{ε} (ε t), \end{matrix}

and we define

φ^{ε}

\begin{matrix} ε φ^{ε} = | x |^{2} / 2 - Φ [ρ], \end{matrix}

so that

\begin{matrix} \nabla φ^{ε} = \nabla ψ^{ε} (\nabla Φ [ρ]) . \end{matrix}

(26)

Hence, at a point

x \in T^{2}

\nabla φ^{ε ⊥}

is the velocity of the associated dual point

\nabla Φ [ρ] (x)

. The evolution of this quantities is then governed by the system

S G_{ε}

\begin{matrix} \partial_{t} ρ^{ε} + \nabla \cdot (ρ^{ε} \nabla ψ^{ε ⊥}) = 0, \end{matrix}

(27)

\begin{matrix} det (I + ε D^{2} ψ^{ε}) = 1 + ε ρ^{ε} . \end{matrix}

(28)

Remark: Note that this system admits global weak solutions with initial data any bounded measure

ρ^{ε 0}

, as long as

\begin{matrix} \int_{T^{2}} ρ^{ε 0} = 0, \end{matrix}

(29)

\begin{matrix} ρ^{ε 0} \geq - \frac{1}{ε} . \end{matrix}

(30)

Note also that if the pair

(\bar{ρ}, \bar{φ})

is solution to the

E I

system ( 24 , 25 ), so is the pair

(\frac{1}{ε} \bar{ρ} (\frac{t}{ε}, x), \frac{1}{ε} \bar{φ} (\frac{t}{ε}, x))

We now present the convergence results. We show that solutions of

S G_{ε}

converge to solutions of

E I

in the following sense: if

ρ^{ε 0}

, the initial data of

S G_{ε}

, is close (in some sense depending on the type of convergence we wish to show) to a smooth initial data

{\bar{ρ}}^{0}

for

E I

, then

ρ^{ε}

and

\bar{ρ}

remain close for some time. This time goes to

\infty

when

ε

goes to 0.

We present two different versions of this result: the first one is for weak solutions of

S G_{ε}

, and the second one is for Lipschitz solutions.

6.2 Convergence of weak solutions

Theorem 6.1. Let

(ρ^{ε}, ψ^{ε})

be a weak solution of the

S G_{ε}

system ( 27 , 28 ). Let

(\bar{ρ}, \bar{φ})

be a smooth

C^{3} ([0, T] \times T^{2})

solution of the

E I

system ( 24 , 25 ). Let

φ^{ε}

be obtained from

ψ^{ε}

as in ( 26 ), let

H_{ε} (t)

be defined by

\begin{matrix} H_{ε} (t) = \frac{1}{2} \int_{T^{2}} {| \nabla φ^{ε} - \nabla \bar{φ} |}^{2}, \end{matrix}

then

\begin{matrix} H_{ε} (t) \leq (H_{ε} (0) + C ε^{2 / 3} (1 + t)) exp C t \end{matrix}

where

C

depends on

{sup}_{0 \leq s \leq t} {∥ D^{3} \bar{φ} (s), D^{2} \partial_{t} \bar{φ} (s) ∥_{L^{\infty} (T^{2})}}

Remark 1: Note that

\nabla φ^{ε ⊥} (t, x)

is the velocity at point

\nabla Φ [ρ] = x - ε \nabla φ^{ε}

Thus we compare the

S G_{ε}

velocity at point

x - ε \nabla φ^{ε}

(the dual point of

x

) with the

E I

velocity at point

x

. Our result allows also to compare the velocities at the same point, by noticing that

\begin{matrix} G_{ε} (t) & = & \frac{1}{2} \int_{T^{2}} ρ {| \nabla ψ^{ε} - \nabla \bar{φ} |}^{2} \end{matrix}

\begin{matrix} = & \frac{1}{2} \int_{T^{2}} {| \nabla φ^{ε} - \nabla \bar{φ} (x - ε \nabla φ^{ε}) |}^{2} \end{matrix}

\begin{matrix} \leq & C (H_{ε} (t) + ε^{2}) \end{matrix}

using the smoothness of

\bar{φ}

, and if

v_{s g_{ε}}, v_{e i}

are the respective velocities of the

S G_{ε}

and

E I

systems,

G_{ε} = \int_{T^{2}} ρ^{ε} | v_{s g_{ε}} - v_{e i} |^{2}

Remark 2: The expansion

det (I + D^{2} ψ) = 1 + Δ ψ + O (| D^{2} ψ |^{2})

, used above to justify the convergence relies a priori on the control of

D^{2} ψ

in the sup norm. But in the Theorem 6.1 , the initial data must satisfy

\nabla ψ^{ε}

close in

L^{2}

norm to a smooth divergence free velocity: this condition means that

D^{2} ψ^{ε}

is close in

H^{- 1}

norm to

D^{2} \bar{φ}

, which is smooth. This control does not allow to justify the expansion

det (I + D^{2} ψ) = 1 + Δ ψ + O (| D^{2} ψ |^{2})

, but we see that the result remains valid.

Proof of Theorem 6.1 In all the proof, we use

C

to denote any quantity that depends only on

\bar{φ}

We use the conservation of the energy of the

S G_{ε}

system, given by

\begin{matrix} E (t) = \int_{T^{2}} | \nabla φ^{ε} |^{2} . \end{matrix}

This fact, although formally easily justified, is actually not so straightforward for weak solutions, and has been proved by F. Otto in an unpublished work.

The argument is explained in [?] . Therefore

E (t) = E_{0}

. The energy of the smooth solution of

E I

is given by

\begin{matrix} \int_{T^{2}} | \nabla \bar{φ} |^{2} \end{matrix}

and also conserved. For all smooth

θ

, we will use the notation:

\begin{matrix} < D^{2} θ > (t, x) = \int_{s = 0}^{1} (1 - s) D^{2} θ (t, x - s ε \nabla φ^{ε} (t, x)) . \end{matrix}

Thus we have the identity

\begin{matrix} \int_{T^{2}} ρ^{ε} θ & = & \int_{T^{2}} θ (x - ε \nabla φ^{ε}) \end{matrix}

(31)

\begin{matrix} = & \int_{T^{2}} θ - ε \int_{T^{2}} \nabla θ \cdot \nabla φ^{ε} + ε^{2} \int_{T^{2}} < D^{2} θ > \nabla φ^{ε} \nabla φ^{ε} . \end{matrix}

(32)

Using the energy bound, the last term is bounded by

ε^{2} ∥ D^{2} θ ∥_{L^{\infty} (T^{2})} E_{0}

. Then we write

\begin{matrix} \frac{d}{d t} H_{ε} (t) & = & - \frac{d}{d t} \int_{T^{2}} \nabla \bar{φ} \cdot \nabla φ^{ε} . \end{matrix}

Using the identity ( 32 ), we have for all smooth

θ

\begin{matrix} ε \int_{T^{2}} \nabla θ \cdot \nabla φ^{ε} = - \int_{T^{2}} ρ^{ε} θ + \int_{T^{2}} θ + ε^{2} \int_{T^{2}} < D^{2} θ > \nabla φ^{ε} \nabla φ^{ε}, \end{matrix}

hence, replacing

θ

\bar{φ}

in this identity, we get

\begin{matrix} \frac{d}{d t} H_{ε} (t) & = & \frac{1}{ε} \frac{d}{d t} \int_{T^{2}} [ρ^{ε} \bar{φ} - \bar{φ} - ε^{2} < D^{2} \bar{φ} > \nabla φ^{ε} \nabla φ^{ε}] . \end{matrix}

We can suppose without loss of generality that

\int_{T^{2}} \bar{φ} (t, x) d x \equiv 0

. Then if we define

\begin{matrix} Q_{ε} (t) = \int_{T^{2}} ε < D^{2} \bar{φ} > \nabla φ^{ε} \nabla φ^{ε}, \end{matrix}

(note that

| Q_{ε} (t) | \leq C ε

), we have

\begin{matrix} \frac{d}{d t} (H_{ε} + Q_{ε}) = \frac{1}{ε} \frac{d}{d t} \int_{T^{2}} ρ^{ε} \bar{φ} . \end{matrix}

Hence we are left to compute

\begin{matrix} \frac{1}{ε} \frac{d}{d t} \int_{T^{2}} ρ^{ε} \bar{φ} & = & \frac{1}{ε} \int_{T^{2}} \partial_{t} ρ^{ε} \bar{φ} + ρ^{ε} \partial_{t} \bar{φ} \end{matrix}

\begin{matrix} = & \frac{1}{ε} \int_{T^{2}} ρ^{ε} \nabla ψ^{ε ⊥} \cdot \nabla \bar{φ} - ε \nabla φ^{ε} \cdot \nabla \partial_{t} \bar{φ} + ε^{2} < D^{2} \partial_{t} \bar{φ} > \nabla \bar{φ} \nabla \bar{φ} \end{matrix}

\begin{matrix} = & \frac{1}{ε} \int_{T^{2}} ρ^{ε} \nabla ψ^{ε ⊥} \cdot \nabla \bar{φ} - \int_{T^{2}} \nabla φ^{ε} \cdot \nabla \partial_{t} \bar{φ} + O (ε) \end{matrix}

\begin{matrix} = & T_{1} + T_{2} + O (ε), \end{matrix}

where at the second line we have used ( 27 ) for the first term and ( 32 ) with

θ = \partial_{t} \bar{φ}

for the second and third term. (Remember also that we assume

\int \partial_{t} \bar{φ} \equiv 0

.) We will now use the other formulation of the Euler equation:

v = \nabla {\bar{φ}}^{⊥}

satisfies

\begin{matrix} \partial_{t} v + v \cdot \nabla v = - \nabla p . \end{matrix}

After a rotation of

π / 2

, this equation becomes:

\begin{matrix} \partial_{t} \nabla \bar{φ} + D^{2} \bar{φ} \nabla {\bar{φ}}^{⊥} = \nabla p^{⊥}, \end{matrix}

thus for

T_{2}

we have

\begin{matrix} T_{2} & = & - \int_{T^{2}} \nabla φ^{ε} \cdot \nabla \partial_{t} \bar{φ} \end{matrix}

\begin{matrix} = & \int_{T^{2}} \nabla φ^{ε} D^{2} \bar{φ} \nabla {\bar{φ}}^{⊥} . \end{matrix}

For

T_{1}

, using ( 26 ) and ( 32 ), we have

\begin{matrix} ε T_{1} & = & \int_{T^{2}} ρ^{ε} \nabla ψ^{ε ⊥} \cdot \nabla \bar{φ} \end{matrix}

\begin{matrix} = & \int_{T^{2}} \nabla ψ^{ε ⊥} (x - ε \nabla φ^{ε}) \cdot \nabla \bar{φ} (x - ε \nabla φ^{ε}) \end{matrix}

\begin{matrix} = & \int_{T^{2}} \nabla φ^{ε ⊥} \cdot \nabla \bar{φ} - ε \nabla φ^{ε ⊥} D^{2} \bar{φ} \nabla φ^{ε} + ε Δ \end{matrix}

where

Δ

is defined by

\begin{matrix} Δ = \int_{T^{2}} \nabla φ^{ε ⊥} (D^{2} \bar{φ} - \int_{s = 0}^{1} D^{2} \bar{φ} (x - s ε \nabla φ^{ε}) d s) \nabla φ^{ε} . \end{matrix}

(33)

The term

\int_{T^{2}} \nabla φ^{ε ⊥} \cdot \nabla \bar{φ}

vanishes identically. Concerning

Δ

, we claim the following estimate:

Lemma 6.2. Let

Δ

be defined by ( 33 ), then

\begin{matrix} | Δ | \leq C (ε^{\frac{2}{3}} + H_{ε}), \end{matrix}

where

C

depends on

∥ D^{3} \bar{φ} ∥_{L^{\infty}}

We postpone the proof of this lemma after the proof of Theorem 6.1 . We now obtain

\begin{matrix} \frac{d}{d t} (H_{ε} (t) + Q_{ε} (t)) & \leq & \int_{T^{2}} (\nabla {\bar{φ}}^{⊥} - \nabla φ^{ε ⊥}) D^{2} \bar{φ} \nabla φ^{ε} + C H_{ε} + C ε^{2 / 3} . \end{matrix}

Noticing that for every

θ : T^{2} \mapsto R

we have

\begin{matrix} \int_{T^{2}} \nabla θ^{⊥} D^{2} \bar{φ} \nabla \bar{φ} = \int_{T^{2}} \nabla θ^{⊥} \cdot \nabla (\frac{1}{2} | \nabla \bar{φ} |^{2}) = 0, \end{matrix}

we find that

\begin{matrix} \int_{T^{2}} (\nabla {\bar{φ}}^{⊥} - \nabla φ^{ε ⊥}) D^{2} \bar{φ} \nabla φ^{ε} = \int_{T^{2}} (\nabla φ^{⊥} - \nabla {\bar{φ}}^{ε ⊥}) D^{2} \bar{φ} (\nabla φ^{ε} - \nabla \bar{φ}), \end{matrix}

hence

\begin{matrix} \frac{d}{d t} (H_{ε} (t) + Q_{ε} (t)) & \leq & - \int_{T^{2}} (\nabla φ^{ε ⊥} - \nabla {\bar{φ}}^{⊥}) D^{2} \bar{φ} (\nabla φ^{ε} - \nabla \bar{φ}) + C H_{ε} + C ε^{2 / 3} \end{matrix}

\begin{matrix} \leq & C (H_{ε} (t) + Q_{ε} (t) + ε^{2 / 3}) \end{matrix}

using that

Q_{ε} (t) \leq C ε

. Therefore

\begin{matrix} H_{ε} (t) + Q_{ε} (t) \leq (H_{ε} (0) + Q_{ε} (0) + C ε^{2 / 3} t) exp (C t) \end{matrix}

and finally

\begin{matrix} H_{ε} (t) \leq (H_{ε} (0) + C ε^{2 / 3} (1 + t)) exp (C t) \end{matrix}

and the result follows. Check that the constant

C

depends only on

{sup}_{0 \leq s \leq t} {∥ D^{3} \bar{φ}, D^{2} \partial_{t} \bar{φ} ∥_{L^{\infty} (T^{2})}}

. This ends the proof of Theorem 6.1

□

Proof of Lemma 6.2 First we show that if

Θ (R) = \int_{{| \nabla φ^{ε} | \geq R}} | \nabla φ^{ε} |^{2}

, then

\begin{matrix} Θ (R) \leq C \int | \nabla φ^{ε} - \nabla \bar{φ} |^{2} + \frac{C}{R^{2}} . \end{matrix}

(34)

Indeed,

\int | \nabla φ^{ε} |^{2} \leq C

, implies that

meas {| \nabla φ^{ε} | \geq R} \leq C \frac{1}{R^{2}} .

Since

| \nabla \bar{φ} (t, x) | \leq C

for

(t, x) \in [0, T^{'}] \times T^{d}

, we have

\begin{matrix} Θ (R) & \leq & \int_{{| \nabla φ^{ε} | \geq R}} | \nabla \bar{φ} |^{2} + \int_{{| \nabla φ^{ε} | \geq R}} | \nabla φ^{ε} - \nabla \bar{φ} |^{2} \end{matrix}

\begin{matrix} \leq & \frac{C}{R^{2}} + \int | \nabla φ^{ε} - \nabla \bar{φ} |^{2} . \end{matrix}

Hence ( 34 ) is proved.

Then, letting

\begin{matrix} K (x) = D^{2} \bar{φ} - \int_{s = 0}^{1} D^{2} \bar{φ} (x - s ε \nabla φ^{ε}) d s, \end{matrix}

we have

\begin{matrix} Δ \leq C Θ (R) + \int_{| \nabla φ^{ε} | \leq R} | K (x) | | \nabla φ^{ε} |^{2} \end{matrix}

with

| K (x) | \leq C ε | \nabla φ^{ε} |

thus

\begin{matrix} Δ \leq C ε \int_{| \nabla φ^{ε} | \leq R} | \nabla φ^{ε} |^{3} + C Θ (R) \end{matrix}

\begin{matrix} \leq C (ε R \int | \nabla φ^{ε} |^{2} + \frac{1}{R^{2}} + \int | \nabla φ^{ε} - \nabla \bar{φ} |^{2}) \end{matrix}

\begin{matrix} \leq C (ε R + \frac{1}{R^{2}} + \int | \nabla φ^{ε} - \nabla \bar{φ} |^{2}) \end{matrix}

for all R, so for

R = ε^{- 1 / 3}

we obtain:

\begin{matrix} Δ \leq C ε^{2 / 3} + C \int | \nabla φ^{ε} - \nabla \bar{φ} |^{2} . \end{matrix}

This proves Lemma 6.2

□

6.3 Convergence of strong solutions

We present here another proof of convergence, that holds for stronger norms. Let us consider as above the solution

(\bar{ρ}, \bar{φ})

to Euler:

\begin{matrix} \partial_{t} \bar{ρ} + \nabla \cdot (\bar{ρ} \nabla {\bar{φ}}^{⊥}) = 0, \end{matrix}

\begin{matrix} Δ \bar{φ} = \bar{ρ}, \end{matrix}

and we recall the

S G_{ε}

system

\begin{matrix} \partial_{t} ρ^{ε} + \nabla \cdot (ρ^{ε} \nabla ψ^{ε ⊥}) = 0, \end{matrix}

\begin{matrix} det (I + ε D^{2} ψ^{ε}) = 1 + ε ρ^{ε} . \end{matrix}

We have then

Theorem 6.3. Let

(\bar{ρ}, \bar{φ})

be a solution of

E I

, such that that

\bar{ρ} \in C_{l o c}^{2} (R^{+} \times T^{2})

. Let

ρ^{ε 0}

be a sequence of initial data for

S G_{ε}

satisfying ( 29 , 30 ), and such that

\frac{ρ^{ε 0} - {\bar{ρ}}^{0}}{ε}

is bounded in

W^{1, \infty} (T^{2})

. Then there exists a sequence

(ρ^{ε}, ψ^{ε})

of solutions to

S G_{ε}

that satisfies:

for all

T > 0

, there exists

ε_{T} > 0

, such that the sequence

\begin{matrix} \frac{ρ^{ε} - \bar{ρ}}{ε}, \frac{\nabla ψ^{ε} - \nabla \bar{φ}}{ε} \end{matrix}

for

0 < ε < ε_{T}

is uniformly bounded in

L^{\infty} ([0, T], W^{1, \infty} (T^{2}))

Remark: In the previous theorem, we obtained estimates in

L^{2}

norm, here we obtain estimates in Lipschitz norm. Estimates of higher derivatives follow in the same way.

Proof of Theorem 6.3 We expand the solution of

S G_{ε}

as the solution of

E I

plus a small perturbation of order

ε

and show that this perturbation remains bounded in large norms (at least Lipschitz). We first remark the the assumption on

\bar{ρ}

implies that

\forall T > 0

\bar{φ} \in L^{\infty} ([0, T]; C^{3} (T^{2}))

. Let us write

\begin{matrix} ρ^{ε} = \bar{ρ} + ε ρ_{1} \end{matrix}

\begin{matrix} ψ^{ε} = \bar{φ} + ε ψ_{1} . \end{matrix}

Rewritten in terms of

ρ_{1}, ψ_{1}

, the

S G_{ε}

system reads:

\begin{matrix} \partial_{t} ρ_{1} + (\nabla \bar{φ} + ε \nabla ψ_{1})^{⊥} \cdot \nabla ρ_{1} = - \nabla ψ_{1}^{⊥} \cdot \nabla \bar{ρ}, \end{matrix}

\begin{matrix} Δ ψ_{1} + ε trace [D^{2} ψ_{1} D^{2} \bar{φ}] + ε^{2} det D^{2} ψ_{1} = ρ_{1} - det D^{2} \bar{φ} . \end{matrix}

Differentiating the first equation with respect to space, we find the evolution equation for

\nabla ρ_{1}

\begin{matrix} \partial_{t} \nabla ρ_{1} + ((\nabla \bar{φ} + ε \nabla ψ_{1})^{⊥} \cdot \nabla) \nabla ρ_{1} \end{matrix}

\begin{matrix} = & - (D^{2} \bar{φ} + ε D^{2} ψ_{1}) \nabla ρ_{1}^{⊥} - D^{2} ψ_{1} \nabla {\bar{ρ}}^{⊥} - D^{2} \bar{ρ} \nabla ψ_{1}^{⊥} . \end{matrix}

(35)

We claim that in order to conclude the proof it is enough to have an estimate of the form

\begin{matrix} ∥ ψ_{1} (t, \cdot) ∥_{C^{1, 1} (T^{2})} \leq C (1 + ∥ ρ_{1} (t, \cdot) ∥_{C^{0, 1} (T^{2})}), \end{matrix}

(36)

where

C

depends on

\bar{φ}

. Let us admit this bound temporarily, and finish the proof of the theorem: using ( 36 ) and ( 35 ), we obtain

\begin{matrix} \frac{d}{d t} ∥ \nabla ρ_{1} ∥_{L^{\infty}} \leq C (t) (1 + ∥ \nabla ρ_{1} ∥_{L^{\infty}} + ε ∥ \nabla ρ_{1} ∥_{L^{\infty}}^{2}), \end{matrix}

where the constant

C (t)

depends on the

C^{2} (T^{2})

norm of

(\bar{ρ} (t, \cdot), \bar{φ} (t, \cdot))

. This quantity is bounded on every interval

[0, T]

Thus we conclude using Gronwall’s lemma that

∥ \nabla ρ_{1} (t, \cdot) ∥_{L^{\infty} (T^{2})}

remains bounded on

[0, T_{ε}]

with

T_{ε}

going to

T

ε

goes to 0. We then choose

T

as large as we want, since when

d = 2

the smooth solution to

E I

is global in time. From estimate ( 36 ) the

W^{1, \infty}

bound on

ρ_{1}

implies a

W^{2, \infty}

bound on

ψ_{1}

. Then, we remember that

\begin{matrix} ρ_{1} = \frac{ρ^{ε} - \bar{ρ}}{ε}, \nabla ψ_{1} = \frac{\nabla ψ^{ε} - \nabla \bar{φ}}{ε} \end{matrix}

to conclude the proof of Theorem 6.3 .

□

Proof of the estimate ( 36 ) We write the equation followed by

ψ_{1}

as follows:

\begin{matrix} Δ ψ_{1} = - trace [ε D^{2} ψ_{1} D^{2} \bar{φ}] - ε^{2} det D^{2} ψ_{1} + ρ_{1} - det D^{2} \bar{φ} . \end{matrix}

We recall that

\begin{matrix} ∥ f g ∥_{C^{2, α}} \leq ∥ f ∥_{C^{2, α}} ∥ g ∥_{C^{2, α}}, \end{matrix}

hence, using Schauder

C^{2, α}

estimates for solutions to Laplace equation (see [?] ), we have

\begin{matrix} ∥ ψ_{1} ∥_{C^{2, α}} \leq C_{1} (1 + ε ∥ ψ_{1} ∥_{C^{2, α}} + ε^{2} ∥ ψ_{1} ∥_{C^{2, α}}^{2}), \end{matrix}

(37)

where

C_{1}

depends on

∥ \bar{φ} ∥_{C^{2, α}}, ∥ ρ_{1} ∥_{C^{α}}

. The inequality ( 37 ) will be satisfied in two cases: either for

∥ ψ_{1} ∥_{C^{2, α}} \leq C_{2}

or for

∥ ψ_{1} ∥_{C^{2, α}} \geq C_{3} ε^{- 2}

where

C_{2}, C_{3}

are positive constants that depend on

C_{1}

Now we show that

ψ^{ε}

, solution of ( 28 ), is bounded in

C^{2, α}

for

ρ^{ε}

bounded in

C^{α}

norm. We consider for

t \in [0, 1]

ψ_{t}^{ε}

the unique up to a constant periodic solution of

\begin{matrix} det (I + ε D^{2} ψ_{t}^{ε}) = 1 + t ε ρ^{ε} . \end{matrix}

Diiferentiating this equation with respect to

t

, we find

\begin{matrix} M_{i j} D_{i j} \partial_{t} ψ_{t}^{ε} = ρ^{ε}, \end{matrix}

where

M

is the co-matrix of

I + ε D^{2} ψ_{t}^{ε}

. From the regularity result of Theorem 4.3 ,

M

C^{α}

and striclty elliptic. From Schauder estimates, we have then

∥ \partial_{t} ψ_{t}^{ε} ∥_{C^{2, α}} \leq C ∥ ρ^{ε} ∥_{C^{2, α}}

, and integrated over

t \in [0, 1]

, we get

\begin{matrix} ∥ ψ^{ε} ∥_{C^{2, α}} \leq C ∥ ρ^{ε} ∥_{C^{2, α}} . \end{matrix}

Hence, since

ψ^{ε} = \bar{φ} + ε ψ_{1}

, we have

ψ_{1}

bounded by

C / ε

C^{2, α}

. Hence it can not be bigger than

C_{3} / ε^{2}

, and to satsify ( 37 ), we must have

\begin{matrix} ∥ ψ_{1} ∥_{C^{2, α}} \leq C_{2}, \end{matrix}

where

C_{2}

as above depends on

∥ \bar{φ} ∥_{C^{2, α}}, ∥ ρ_{1} ∥_{C^{α}}

. This proves estimate ( 36 ).

□

Acknowledgment: The author thanks Mike Cullen for his remarks, and also Yann Brenier, since part of this work was done under his direction, during the author’s PhD thesis. He also thanks Robert McCann and the Fields Institute of Toronto for their hospitality.

G. Loeper EPFL, SB, IMA 10015 Lausanne e-mail: gregoire.loeper@epfl.ch