it might be because of the absence of pressure that changes the symmetrizers of the system. After a convenient change of variable however, the system appears under a form which is strongly reminiscent (at least for the highest order terms) of the rapidly rotating fluids. This limit has been treated in [?] .

Finally we also mention the work of Cordier & Grenier [?] , and Wang [?] , where the quasi-neutral limit ’with pressure’ is treated. The techniques used there do not apply here, and it is worth noting that the scaling obtained are not the same. The reader can also refer to [?]  and [?]  where different regimes of the Euler-Poisson system are studied.

We split the rest of the paper in two sections: the first one devoted to the study of the Euler-Poisson system, and the second devoted to the study of the Euler-Monge-Ampere system.

1.1 Notation

Hereafter $x\in {\mathbb{T}}^d={\mathbb{R}}^d/{\mathbb{Z}}^d$  and $t\in {\mathbb{R}}^{+}$  ; $v(t,x)\in {\mathbb{R}}^d$  stands for the velocity and $\rho (t,x)\in {\mathbb{R}}^{+}$  is the macroscopic density of electrons; $\phi (t,x)\in \mathbb{R}$  is the electrostatic potential; $d=2$  or $3$  .

It is always assumed that $\rho (t,\cdot )$  has total mass equal to 1.

The divergence of a vector field $v$  will be denoted by $divv$  or $\nabla \cdot v$  ; its rotational (or curl) will be denoted by $\nabla \times v$  or $curlv$  .

The components of a vector will be denoted with superindices, i.e. $v\in {\mathbb{R}}^d=(v^i{)}_{i=1..d}$  .

In all the paper, $[\cdot ]$  will denote the integer part.

We denote respectively by $\left(E\right)$  , $\left(E{P}_{\epsilon }\right)$  , $\left(EM{A}_{\epsilon }\right)$  the systems Euler incompressible, Euler-Poisson, Euler-Monge-Ampere.

2 The Euler-Poisson system

We consider the following Euler-Poisson system denoted by $\left(E{P}_{\epsilon }\right)$  :

$$\begin{array}{ccc}& & {\partial }_t\rho +\nabla \cdot \left(\rho v\right)=0,\end{array}$$

(5)

$$\begin{array}{ccc}& & {\partial }_{t}v+v\cdot \nabla v=\frac{\nabla \phi }{\epsilon },\end{array}$$

(6)

$$\begin{array}{ccc}& & \epsilon \Delta \phi =\rho -1,\end{array}$$

(7)

and consider the limit $\epsilon$  going to 0. We recall also the incompressible Euler equation $\left(E\right)$  :

$$\begin{array}{ccc}& & {\partial }_{t}v+v\cdot \nabla v=\nabla p,\end{array}$$

$$\begin{array}{ccc}& & \nabla \cdot v=0.\end{array}$$

(8)

We recall (see [?]  for a reference on the topic) that in the periodic case, the Cauchy problem for ( 8 ) is well posed in $H^s\left({\mathbb{T}}^d\right)$  if $s>d/2+1$  . More precisely, if $d=2$  , for any divergence-free initial datum in $H^s$  , there exists a unique global solution in $L_{loc}^{\infty }(\mathbb{R},H^s({\mathbb{T}}^d\left)\right)$  ; if $d=3$  one can only prove the existence of a smooth solution in finite time, belonging to $L_{loc}^{\infty }\left(\right[0,T[,H^s({\mathbb{T}}^d\left)\right)$  for some $T>0$  . We will then prove the following:

2.1 Result

Theorem 2.1 Let $s\in \mathbb{N}$  with $s\ge [d/2]+2$  . Let ${\overline{v}}_0$  be a divergence-free vector field on ${\mathbb{T}}^d$  . Let $(\overline{v},p)$  be a smooth solution of the Euler incompressible system ( 8 ) on $[0,T]\times {\mathbb{T}}^d$  , with initial data ${\overline{v}}_0$  , satisfying $\overline{v}\in L^{\infty }\left(\right[0,T],H^{s+1}({\mathbb{T}}^d\left)\right)$  . Let $({\rho }_0^{\epsilon },v_0^{\epsilon })$  be a sequence of initial data such that ${\int }_{{\mathbb{T}}^d}{\rho }^{\epsilon }\left(x\right)dx=1$  , and such that $\left({\epsilon }^{-1}(v_0^{\epsilon }-{\overline{v}}_0),{\epsilon }^{-2}({\rho }_0^{\epsilon }-1)\right)$  is bounded in $H^s\times H^{s-1}\left({\mathbb{T}}^d\right)$  .

Then there exists a sequence $(v^{\epsilon },{\rho }^{\epsilon })$  of solutions to $\left(E{P}_{\epsilon }\right)$  with initial data $(v_0^{\epsilon },{\rho }_0^{\epsilon })$  , belonging to $L^{\infty }\left(\right[0,T_{\epsilon }],H^s\times H^{s-1}({\mathbb{T}}^d\left)\right)$  , with ${liminf}_{\epsilon \to 0}{T}_{\epsilon }\ge T$  . Moreover for any $T^{\prime }<T$  and $\epsilon$  small enough, $\left({\epsilon }^{-1}(v^{\epsilon }-\overline{v}),{\epsilon }^{-2}({\rho }^{\epsilon }-1)\right)$  is bounded in $L^{\infty }\left(\right[0,T^{\prime }],H^s\times H^{s-1}({\mathbb{T}}^d\left)\right)$  . Finally when $T=+\infty$  , $T_{\epsilon }$  goes to infinity.

Remark. The models that we consider are valid in domains without boundaries, and although stated in the space periodic case, we believe that our results hold true, with some technical adaptation, in the case of the whole space.

2.2 Proof of Theorem  2.1 

2.2.1 Heuristics.

Let us introduce $(\overline{v},p)$  the solution of the Euler incompressible system ( 8 ) and corresponding pressure. Note that by taking the divergence of ( 8 ) the pressure is given by the following:

$$\begin{array}{ccc}& & \Delta p={\sum }_{i,j=1}^d{\partial }_i{\overline{v}}^j{\partial }_j{\overline{v}}^i.\end{array}$$

We will all along the paper use the following notation: for two vector fields $u,v$  ,

$$\begin{array}{ccc}& & \nabla u:\nabla v={\sum }_{i,j=1}^d{\partial }_i{u}^j{\partial }_j{v}^i.\end{array}$$

If $v$  is solution to $\left(E{P}_{\epsilon }\right)$  , we introduce also

$$\begin{array}{ccc}& & v=\overline{v}+\epsilon v_1,\end{array}$$

$$\begin{array}{ccc}& & \rho =1+{\epsilon }^2{\rho }_1.\end{array}$$

We suppose for now that $d=2$  , and we take the curl and the divergence of the momentum equation, this yields

$$\begin{array}{ccc}& & {\partial }_{t}curl{v}_1+v\cdot \nabla curl{v}_1=R_1,\end{array}$$

$$\begin{array}{ccc}& & {\partial }_{t}div{v}_1+v\cdot \nabla div{v}_1=\frac{{\rho }_1}{\epsilon }+R_2,\end{array}$$

$$\begin{array}{ccc}& & {\partial }_t{\rho }_1+v\cdot \nabla {\rho }_1=-\frac{div{v}_1}{\epsilon }+R_3.\end{array}$$

If we assume for now that $v_1,{\rho }_1$  and their spatial derivatives remain bounded (we do not specify in what sense yet), $R_i,i=1,2,3$  are bounded terms. Hence the vorticity $curlv$  is not affected by the electric field, and the vector $\mathbf{u}=(curl{v}_1,div{v}_1,{\rho }_1)$  evolves through

$$\begin{array}{ccc}{\partial }_t\mathbf{u}+v\cdot \nabla \mathbf{u}=\frac{\mathbf{k}}{\epsilon }\times \mathbf{u}+\mathbf{R},& & \end{array}$$

with $\mathbf{R}$  bounded, and $\mathbf{k}=(1,0,0)$  . Under this form, the system looks like a rapidly rotating fluid (up to the remainder $\mathbf{R}$  , and the fact that $v\ne \mathbf{u}$  ) and the singular term $\frac{\mathbf{k}}{\epsilon }\times \mathbf{u}$  induces time oscillations of frequency ${\epsilon }^{-1}$  , but does not increase the energy of the perturbation, allowing energy estimates, as long as the remainder $\mathbf{R}$  is under control.

2.2.2 Reformulation of the system with new unknowns

For a vector $u\in {\mathbb{R}}^d$  , we denote $u^1,...,u^d$  its components. We define the new unknowns ${\omega }_1,{\beta }_1,{\rho }_1$  as follows:

$$\begin{array}{ccc}& & \nabla \cdot v=\epsilon {\beta }_1,\end{array}$$

$$\begin{array}{ccc}& & \rho =1+{\epsilon }^2{\rho }_1=1+\epsilon \Delta \phi ,\end{array}$$

$$\begin{array}{ccc}& & \text{curl}v=\omega =\overline{\omega }+\epsilon {\omega }_1,\end{array}$$

with $(\overline{v},p)$  the solution of $\left(E\right)$  , and $\overline{\omega }=\text{curl}\overline{v}$  .

We will use the following observation:

Lemma 2.2 Let $(\rho ,v)$  be a solution to $\left(E{P}_{\epsilon }\right)$  . Then, the total momentum $\int \rho (t,x)v(t,x)dx$  does not depend on time.

Proof. The momentum equation can be rewritten

$$\begin{array}{ccc}{\partial }_t\left(\rho v\right)+\nabla \cdot (\rho v\otimes v)=\rho \nabla \phi .& & \end{array}$$

We just have to show that the integral of the right hand side vanishes. For this, we use the identity $$(1+{\epsilon }^2\Delta \phi )\nabla \phi ={\epsilon }^2\left[\nabla \cdot (\nabla \phi \otimes \nabla \phi )-\frac{1}{2}\nabla |\nabla \phi {|}^2\right]+\nabla \phi ,$$  that yields 0 when integrated over ${\mathbb{T}}^d$  .

$\square$  We note also

$$\begin{array}{ccc}v=\overline{v}+\epsilon v_1,& & \end{array}$$

but $v_1$  is not really an unknown since it can be obtained from the knowledge of ${\omega }_1,{\beta }_1$  , and $\epsilon \int \rho v_1=\int {\rho }_0{v}_0-\int \rho \overline{v}$  : Indeed when $d=2$  we have

$$\begin{array}{ccc}& & {\partial }_1{\beta }_1+{\partial }_2{\omega }_1=\Delta v_1^1,\end{array}$$

(9)

$$\begin{array}{ccc}& & {\partial }_2{\beta }_1-{\partial }_1{\omega }_1=\Delta v_1^2.\end{array}$$

(10)

In the 3 dimensional case we have equations ( 9 ,  10 ) replaced by

$$\begin{array}{ccc}(\nabla \times \nabla \times v)+\nabla (\nabla \cdot v)=\Delta v,& & \end{array}$$

and thus

$$\begin{array}{ccc}(\nabla \times {\omega }_1{)}^i+{\partial }_i{\beta }_1=\Delta v_1^i.& & \end{array}$$

(11)

Note that when $d=2$  the vorticity is scalar and when $d\ne 2$  it is a vector field of ${\mathbb{T}}^d$  . We show now that $v_1$  can be estimated from ${\omega }_1$  , ${\beta }_1$  , and $\int \rho v$  .

Lemma 2.3 Let $(\rho ,v)$  be the solution at time $t$  of $\left(E{P}_{\epsilon }\right)$  with initial datum $({\rho }_0,v_0)$  , let $\overline{v}$  be a solution at time $t$  of $\left(E\right)$  with initial datum ${\overline{v}}_0$  . Then we have, for $s\in \mathbb{R}$  ,

$$\begin{array}{ccc}& & \parallel v-\overline{v}-\int [v-\overline{v}]{\parallel }_{H^{s+1}\left({\mathbb{T}}^d\right)}\le C\left(\parallel \nabla \cdot v{\parallel }_{H^s\left({\mathbb{T}}^d\right)}+\parallel \nabla \times (v-\overline{v}){\parallel }_{H^s\left({\mathbb{T}}^d\right)}\right),\end{array}$$

$$\begin{array}{ccc}& & \left|\int [v-\overline{v}]\right|\le \end{array}$$

$$\begin{array}{ccc}& & C\left(\parallel \rho -1{\parallel }_{L^2\left({\mathbb{T}}^d\right)}(\parallel \overline{v}{\parallel }_{L^2\left({\mathbb{T}}^d\right)}+\parallel \nabla \cdot v{\parallel }_{H^{-1}\left({\mathbb{T}}^d\right)}+\parallel \nabla \times (v-\overline{v}\left){\parallel }_{H^{-1}\left({\mathbb{T}}^d\right)}\right)+\left|\int [{\rho }_0{v}_0-{\overline{v}}_0]\right|\right).\end{array}$$

Proof. Let $u_1$  be the unique vector field with zero average such that

$$\begin{array}{ccc}& & \nabla \times u_1=\nabla \times (v-\overline{v}),\end{array}$$

$$\begin{array}{ccc}& & \nabla \cdot u_1=\nabla \cdot v.\end{array}$$

We have directly from identities ( 9 ,  10 ,  11 ) that $\parallel u_1{\parallel }_{H^{s+1}}\le C\left(\parallel \nabla \cdot v{\parallel }_{H^s}+\parallel \nabla \times (v-\overline{v}){\parallel }_{H^s}\right)$  , which is the first point of the lemma. The difference $v-\overline{v}-u_1=w_1$  is a constant vector field, and we have

$$\begin{array}{ccc}\int \rho v& =& \int \rho \overline{v}+\int \rho (v-\overline{v})\end{array}$$

$$\begin{array}{ccc}& =& \int \overline{v}+\int (\rho -1)\overline{v}+\int (\rho -1)u_1+w_1,\end{array}$$

which yields

$$\begin{array}{ccc}{w}_1=\int [{\rho }_0{v}_0-{\overline{v}}_0]-\int (\rho -1)\overline{v}-\int (\rho -1)u_1,& & \end{array}$$

and the result follows.

$\square$  We immediately deduce the following corollary:

Corollary 2.4 Let $(\rho ,v)$  be a solution to $\left(E{P}_{\epsilon }\right)$  , let $\overline{v}$  be a solution to $\left(E\right)$  . Then, for any $s\ge 0$  ,

$$\begin{array}{ccc}\parallel v-\overline{v}{\parallel }_{H^s}& \le & \left|\int [{\rho }_0{v}_0-{\overline{v}}_0]\right|\end{array}$$

$$\begin{array}{ccc}& +& C(1+\parallel \rho -1{\parallel }_{H^s})\left(\parallel \nabla \cdot v{\parallel }_{H^{s-1}}+\parallel \nabla \times (v-\overline{v}){\parallel }_{H^{s-1}}\right)\end{array}$$

$$\begin{array}{ccc}& +& C\parallel \rho -1{\parallel }_{H^s}\parallel \overline{v}{\parallel }_{H^s}.\end{array}$$

Taking the curl of equation ( 6 ) we recall the following identities:

$$\begin{array}{ccc}& & {\partial }_t\omega +(v\cdot \nabla )\omega +(\nabla \cdot v)\omega =0\text{when}d=2,\end{array}$$

(12)

$$\begin{array}{ccc}& & {\partial }_t\omega +(v\cdot \nabla )\omega +(\nabla \cdot v)\omega +(\omega \cdot \nabla )v=0\text{when}d=3.\end{array}$$

(13)

When $d=2$  the $\left(E{P}_{\epsilon }\right)$  system then becomes:

$$\begin{array}{ccc}& & {\partial }_t(\overline{\omega }+\epsilon {\omega }_1)+v\cdot \nabla (\overline{\omega }+\epsilon {\omega }_1)=-(\overline{\omega }+\epsilon {\omega }_1)\epsilon \beta ,\end{array}$$

(14)

$$\begin{array}{ccc}& & {\partial }_t\epsilon {\beta }_1+v\cdot \nabla \epsilon {\beta }_1+2\epsilon {\partial }_i{\overline{v}}^j{\partial }_j{v}_1^i+{\epsilon }^2{\partial }_i{v}_1^j{\partial }_j{v}_1^i=\frac{\Delta \phi }{\epsilon }-{\partial }_i{\overline{v}}^j{\partial }_j{\overline{v}}^i,\end{array}$$

(15)

$$\begin{array}{ccc}& & {\partial }_t{\epsilon }^2{\rho }_1+v\cdot \nabla {\epsilon }^2{\rho }_1=-(1+{\epsilon }^2{\rho }_1)\epsilon {\beta }_1,\end{array}$$

(16)

whereas when $d=3$  one would have to replace equation ( 14 ) by

$$\begin{array}{ccc}{\partial }_t(\overline{\omega }+\epsilon {\omega }_1)+v\cdot \nabla (\overline{\omega }+\epsilon {\omega }_1)=-(\overline{\omega }+\epsilon {\omega }_1)\epsilon \beta -\omega \cdot \nabla v.& & \end{array}$$

(17)

Noticing that $\Delta \phi =\epsilon {\rho }_1$  , if we set

$$\begin{array}{ccc}{\tilde{\rho }}_1={\rho }_1-\Delta p,& & \end{array}$$

(18)

we get the following system for $d=2$  :

$$\begin{array}{ccc}\{\begin{array}{c}{\partial }_t{\omega }_1+v\cdot \nabla {\omega }_1=-\overline{\omega }{\beta }_1-\epsilon {\omega }_1{\beta }_1-v_1\cdot \nabla \overline{\omega },\\ {\partial }_t{\beta }_1+v\cdot \nabla {\beta }_1=\frac{{\tilde{\rho }}_1}{\epsilon }-2\nabla v_1:\nabla \overline{v}-\epsilon \nabla v_1:\nabla v_1,\\ {\partial }_t{\tilde{\rho }}_1+v\cdot \nabla {\tilde{\rho }}_1=-\frac{{\beta }_1}{\epsilon }-\epsilon ({\tilde{\rho }}_1+\Delta p){\beta }_1-({\partial }_t\Delta p+v\cdot \nabla \Delta p).\end{array}& & \end{array}$$

(19)

When $d=3$  , the first equation should be replaced by

$$\begin{array}{ccc}& & {\partial }_t{\omega }_1+v\cdot \nabla {\omega }_1=\end{array}$$

$$\begin{array}{ccc}& & -\overline{\omega }{\beta }_1-v_1\cdot \nabla \overline{\omega }-\overline{\omega }\cdot \nabla v_1-\omega \cdot \nabla \overline{v}-\epsilon {\omega }_1{\beta }_1-\epsilon {\omega }_1\cdot \nabla v_1.\end{array}$$

2.2.3 Energy estimates

We handle the energy estimates when $d=2$  but the same result would hold when $d=3$  , just with more terms. The system can be written in the following way:

$$\begin{array}{ccc}& & {\partial }_t{\mathbf{u}}^{\epsilon }+{\sum }_i{v}^i{\partial }_i{\mathbf{u}}^{\epsilon }+R^{\epsilon }{\mathbf{u}}^{\epsilon }=S^{\epsilon }\left({\mathbf{u}}^{\epsilon }\right),\end{array}$$

(20)

$$\begin{array}{ccc}& & {\mathbf{u}}^{\epsilon }\left(0\right)={\mathbf{u}}_0^{\epsilon },\end{array}$$

(21)

where $v$  is still the velocity, and where

$$\begin{array}{ccc}{\mathbf{u}}^{\epsilon }=\left(\begin{array}{c}{\omega }_1\\ {\beta }_1\\ {\tilde{\rho }}_1\end{array}\right),R^{\epsilon }=\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& -\frac{1}{\epsilon }\\ 0& \frac{1}{\epsilon }& 0\end{array}\right).& & \end{array}$$

The ’source’ term $S^{\epsilon }$  is given by

$$\begin{array}{ccc}{S}^{\epsilon }=\left(\begin{array}{c}-\overline{\omega }{\beta }_1-\epsilon {\omega }_1{\beta }_1-v_1\cdot \nabla \overline{\omega }\\ -2\nabla v_1:\nabla \overline{v}-\epsilon \nabla v_1:\nabla v_1\\ -\epsilon ({\tilde{\rho }}_1+\Delta p){\beta }_1-({\partial }_t\Delta p+v\cdot \nabla \Delta p)\end{array}\right).& & \end{array}$$

We apply ${\partial }^{\gamma }$  to equation ( 20 ), where $\gamma =({\gamma }^1,...,{\gamma }^d)$  , and ${\partial }^{\gamma }$  stands for $\frac{{\partial }^{\left|\gamma \right|}}{\partial x_1^{{\gamma }^1}...\partial x_d^{{\gamma }^d}}$  , with $\left|\gamma \right|={\sum }_{i=1}^d{\gamma }^i$  .

We get

$$\begin{array}{ccc}{\partial }_t{\partial }^{\gamma }{\mathbf{u}}^{\epsilon }+v^i{\partial }_i{\partial }^{\gamma }{\mathbf{u}}^{\epsilon }+{\Sigma }^{\epsilon }+R^{\epsilon }{\partial }^{\gamma }{\mathbf{u}}^{\epsilon }={\partial }^{\gamma }{S}^{\epsilon },& & \end{array}$$

where

$$\begin{array}{ccc}{\Sigma }^{\epsilon }={\sum }_{i=1}^d{\sum }_{\left|\mu \right|\ge 1,\gamma \ge \mu }{\partial }^{\mu }{v}^i{\partial }^{\gamma -\mu }{\partial }_i{\mathbf{u}}^{\epsilon }.& & \end{array}$$

Then we have

Lemma 2.5 If $\left|\gamma \right|>d/2$  , then for ${\Sigma }^{\epsilon }$  and $S^{\epsilon }$  defined as above we have:

$$\begin{array}{ccc}\parallel {\Sigma }^{\epsilon }(t,\cdot ){\parallel }_{L^2\left({\mathbb{T}}^d\right)}\le C(1+\parallel {\mathbf{u}}^{\epsilon }(t,\cdot ){\parallel }_{H^{\left|\gamma \right|}\left({\mathbb{T}}^d\right)}+\epsilon \parallel {\mathbf{u}}^{\epsilon }(t,\cdot \left){\parallel }_{H^{\left|\gamma \right|}\left({\mathbb{T}}^d\right)}^2\right)& & \end{array}$$

and

$$\begin{array}{ccc}\parallel {\partial }^{\gamma }{S}^{\epsilon }(t,\cdot ){\parallel }_{L^2\left({\mathbb{T}}^d\right)}\le C(1+\parallel {\mathbf{u}}^{\epsilon }(t,\cdot ){\parallel }_{H^{\left|\gamma \right|}\left({\mathbb{T}}^d\right)}+\epsilon \parallel {\mathbf{u}}^{\epsilon }(t,\cdot \left){\parallel }_{H^{\left|\gamma \right|}\left({\mathbb{T}}^d\right)}^2\right).& & \end{array}$$

Proof. Point 1: the proof is a straightforward adaptation of the proof of [?]  p.151. It is based on the following estimate: [?]  Proposition 2.1.2 p. 100:

Proposition 2.6 If $u,v\in L^{\infty }\cap H^s$  $s\in \mathbb{N}$  , then for any $\delta ,\eta ,\left|\delta \right|+\left|\eta \right|=s$  , one has

$$\begin{array}{ccc}\parallel \left({\partial }^{\delta }u\right)\left({\partial }^{\eta }v\right){\parallel }_{L^2}\le C(\parallel u{\parallel }_{L^{\infty }}\parallel v{\parallel }_{H^s}+\parallel u{\parallel }_{H^s}\parallel v{\parallel }_{L^{\infty }}).& & \end{array}$$

Applying this result to ${\partial }^{\mu }{v}^i{\partial }^{\gamma -\mu }{\partial }_i{\mathbf{u}}^{\epsilon }$  , $\left|\mu \right|\ge 1$  , we obtain

$$\begin{array}{ccc}\parallel {\partial }^{\mu }{v}^i{\partial }^{\gamma -\mu }{\partial }_i{\mathbf{u}}^{\epsilon }{\parallel }_{L^2\left({\mathbb{T}}^d\right)}\le C\left(\parallel \nabla v^i{\parallel }_{L^{\infty }\left({\mathbb{T}}^d\right)}\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^{\left|\gamma \right|}\left({\mathbb{T}}^d\right)}+\parallel \nabla v^i{\parallel }_{H^{\left|\gamma \right|}\left({\mathbb{T}}^d\right)}\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{L^{\infty }\left({\mathbb{T}}^d\right)}\right).& & \end{array}$$

We know that $\parallel \cdot {\parallel }_{L^{\infty }\left({\mathbb{T}}^d\right)}\le C\parallel \cdot {\parallel }_{H^{\left|\gamma \right|}\left({\mathbb{T}}^d\right)}$  if $\left|\gamma \right|>d/2$  . Thanks to ( 9 ,  10 ) we have for any $s$ 

$$\begin{array}{ccc}\parallel \nabla v_1{\parallel }_{H^s\left({\mathbb{T}}^d\right)}\le C(\parallel {\omega }_1{\parallel }_{H^s\left({\mathbb{T}}^d\right)}+\parallel {\beta }_1{\parallel }_{H^s\left({\mathbb{T}}^d\right)}),& & \end{array}$$

(22)

thus $\parallel \nabla v_1{\parallel }_{H^{\left|\gamma \right|}\left({\mathbb{T}}^d\right)}\le C\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^{\left|\gamma \right|}\left({\mathbb{T}}^d\right)}$  . Then using that $v=\overline{v}+\epsilon v_1$  we conclude.

Point 2: We also know thanks to Proposition  2.6 that for $s>d/2$ 

$$\begin{array}{ccc}\parallel \nabla v_1:\nabla v_1{\parallel }_{H^s}\le C\parallel \nabla v_1{\parallel }_{L^{\infty }}\parallel \nabla v_1{\parallel }_{H^s}.& & \end{array}$$

It follows that for $s>d/2$  , we have

$$\begin{array}{ccc}\parallel S^{\epsilon }{\parallel }_{H^s}\le C(1+\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^s}+\epsilon \parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^s}^2)& & \end{array}$$

where $C$  depends on the smoothness of the solution of ( 8 ).

$\square$  Having applied ${\partial }^{\gamma }$  to ( 20 ), multiplied by ${\partial }^{\gamma }{\mathbf{u}}^{\epsilon }$  , and noticed that for any $w$  , $(w,R^{\epsilon }w)=0$  , we obtain

$$\begin{array}{ccc}& & {\partial }_t\frac{1}{2}|{\partial }^{\gamma }{\mathbf{u}}^{\epsilon }{|}^2+{\sum }_{i=1}^d{\partial }_i(v^i\frac{1}{2}\left|{\partial }^{\gamma }{\mathbf{u}}^{\epsilon }{|}^2\right)\end{array}$$

$$\begin{array}{ccc}& =& \nabla \cdot v\frac{1}{2}|{\partial }^{\gamma }{\mathbf{u}}^{\epsilon }{|}^2+({\partial }_{\gamma }{S}^{\epsilon }-{\Sigma }^{\epsilon }){\partial }^{\gamma }{\mathbf{u}}^{\epsilon }.\end{array}$$

Since $\parallel \nabla \cdot v{\parallel }_{L^{\infty }}=\epsilon \parallel {\beta }_1{\parallel }_{L^{\infty }}\le C\epsilon \parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^{\gamma }}$  if $\gamma >d/2$  , using Lemma  2.5 , and integrating over ${\mathbb{T}}^d$  we have, for any $\left|\gamma \right|\le s$  , $s>d/2$  :

$$\begin{array}{ccc}\frac{d}{dt}\parallel {\partial }_{\gamma }{\mathbf{u}}^{\epsilon }(t,\cdot ){\parallel }_{L^2\left({\mathbb{T}}^d\right)}^2\le C\left(1+\parallel {\mathbf{u}}^{\epsilon }(t,\cdot ){\parallel }_{H^s\left({\mathbb{T}}^d\right)}^2+\epsilon \parallel {\mathbf{u}}^{\epsilon }(t,\cdot ){\parallel }_{H^s\left({\mathbb{T}}^d\right)}^3\right).& & \end{array}$$

By summing this over all multi-indexes $\left|\gamma \right|\le s$  , we can conclude using a standard Gronwall’s lemma, that if the solution $\overline{v},p$  of Euler is smooth (see below for the smoothness required) on the time interval $[0,T]$  , for any $T^{\prime }<T$  there exists ${\epsilon }_0$  such that the sequence $({\mathbf{u}}^{\epsilon }{)}_{\epsilon <{\epsilon }_0}$  is bounded in $L^{\infty }\left(\right[0,T^{\prime }],H^s({\mathbb{R}}^2\left)\right)$  .

Then we have

$$\begin{array}{ccc}& & {\mathbf{u}}^{\epsilon }=(curl{v}_1,div{v}_1,{\rho }_1-\Delta p),\end{array}$$

$$\begin{array}{ccc}& & v=\overline{v}+\epsilon v_1,\end{array}$$

$$\begin{array}{ccc}& & \rho =1+{\epsilon }^2{\rho }_1.\end{array}$$

We use Lemma  2.3 to get $v$  from $\mathbf{u},\overline{v}$  , and from Corollary  2.4 , the bound obtained on $\mathbf{u}$  implies a bound on ${\epsilon }^{-1}(v-\overline{v}),{\epsilon }^{-2}(\rho -1).$  (Note that from the assumption on the initial data, we have $\left|\int {\rho }_0^{\epsilon }{v}_0^{\epsilon }-{\overline{v}}_0\right|\le C\epsilon$  .) This proves Theorem  2.1 .

$\square$ 

Minimal regularity for the limiting field

In order to perform our computations, we need to have at least ${\beta }_1,{\tilde{\rho }}_1,{\omega }_1$  in $L^{\infty }\left(\right[0,T]\times {\mathbb{T}}^d)$  and thus to have an estimate on their norms in $L^{\infty }\left(\right[0,T],H^s({\mathbb{T}}^d\left)\right)$  , with $s>d/2$  . Therefore we need to control $\parallel v_1{\parallel }_{H^{s+1}},\parallel {\rho }_1{\parallel }_{H^s}$  . Then we need to apply ${\partial }^{\gamma }$  to  19 , with $\left|\gamma \right|=s>d/2$  , and we need to control $\parallel {\partial }_t\Delta p+v\cdot \nabla \Delta p{\parallel }_{H^s}$  (this is the ’worst’ term). This implies to control $\parallel \overline{v}{\parallel }_{H^{s+2}}$  . which requires $\overline{v}$  to be bounded in $L^{\infty }\left(\right[0,T],H^{s+2}({\mathbb{T}}^d\left)\right)$  with $s>d/2$  . If we take $s$  integer, we ask $s\ge [d/2]+1$  , and $\overline{v}$  must be bounded in $L^{\infty }\left(\right[0,T],H^{[d/2]+3}({\mathbb{T}}^d\left)\right)$  .

Remark. Usually, modulated energy techniques only require a bound on $\parallel \nabla \overline{v}{\parallel }_{L^{\infty }}$  . Here we need one more derivative, since the ’div-curl’ formulation of the system (performed in order to obtain energy estimates) is obtained by differentiation.

2.3 Non-prepared initial data

Here we obtain energy estimates in the general case of non-prepared initial data. What we mean by ’general case’ is the case of a generic smooth initial velocity, and smooth initial density, with finite energy, hence $\rho -1$  will be of order $\epsilon$  .

We will see that the energy estimates are the same as in the case of prepared initial data, the asymptotic $\epsilon \to 0$  is then handled similarly to [?] , although the algebra in our case is quite simple.

The solution will exhibit a good space regularity and a strongly oscillating behavior with respect to time. As explained in [?] , the motion can be decomposed along a slow and a fast manifold; the slow manifold consists of divergence-free velocities with uniform density, and the fast manifold consists of potential velocities. Due to the rapid oscillations, the potential part of $v$  will converge weakly to 0, and the divergence-free part will converge strongly to a smooth solution of the incompressible Euler equation.

We still consider $({\rho }^{\epsilon },v^{\epsilon })$  solution to $\left(E{P}_{\epsilon }\right)$  . For any vector field $v$  , we introduce its solenoidal part and potential part, which is the pair $(\Pi v,\nabla q)$  such that $v=\Pi v+\nabla q$  with $q$  periodic and $\nabla \cdot \Pi v=0$  .

A priori estimates We first rescale the density fluctuation

$$\begin{array}{ccc}& & {\rho }^{\epsilon }=1+\epsilon {\rho }_1^{\epsilon }.\end{array}$$

We use the unknown ${\mathbf{u}}^{\epsilon }$  , given by

$$\begin{array}{ccc}{\mathbf{u}}^{\epsilon }=\left(\begin{array}{c}{\omega }^{\epsilon }=\text{curl}{v}^{\epsilon }\\ {\beta }^{\epsilon }=\text{div}{v}^{\epsilon }\\ {\rho }_1^{\epsilon }\end{array}\right).& & \end{array}$$

(Note that we do not subtract $\Delta p$  to ${\rho }_1^{\epsilon }$  in this case.) We now restrict to the case $d=2$  , but one can check easily that the same results hold when $d=3$  .

We write the equation followed by ${\mathbf{u}}^{\epsilon }$  :

$$\begin{array}{ccc}& & {\partial }_t{\mathbf{u}}^{\epsilon }+{\sum }_i{v}^i{\partial }_i{\mathbf{u}}^{\epsilon }+R^{\epsilon }{\mathbf{u}}^{\epsilon }=S^{\epsilon }\left({\mathbf{u}}^{\epsilon }\right),\end{array}$$

(23)

$$\begin{array}{ccc}& & {\mathbf{u}}^{\epsilon }\left(0\right)={\mathbf{u}}_0^{\epsilon },\end{array}$$

(24)

where $R^{\epsilon }$  is a before, and the source term $S^{\epsilon }$  is now given by

$$\begin{array}{ccc}{S}^{\epsilon }=\left(\begin{array}{c}-{\beta }^{\epsilon }{\omega }^{\epsilon }\\ -\nabla v^{\epsilon }:\nabla v^{\epsilon }\\ -{\rho }_1^{\epsilon }{\beta }^{\epsilon }\end{array}\right).& & \end{array}$$

When $s>d/2$  , proceeding as in Lemma  2.5 , we have the estimate

$$\begin{array}{ccc}\parallel S^{\epsilon }{\parallel }_{H^s}\le C_s\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^s}^2.& & \end{array}$$

Arguing as in the previous case, we conclude that, given $s\ge [d/2]+2$  and a sequence of initial data $({\rho }_0^{\epsilon },v_0^{\epsilon })$  , with $\int {\rho }_0^{\epsilon }=1$  , and such that

$$\begin{array}{ccc}& & \parallel v_0^{\epsilon }{\parallel }_{H^s}\le C,\end{array}$$

$$\begin{array}{ccc}& & {\epsilon }^{-1}\parallel {\rho }_0^{\epsilon }-1{\parallel }_{H^{s-1}}\le C,\end{array}$$

there exists a sequence $({\rho }^{\epsilon },v^{\epsilon })$  of solutions to $\left(E{P}_{\epsilon }\right)$  with initial data $({\rho }_0^{\epsilon },v_0^{\epsilon })$  , and $\left({\epsilon }^{-1}\right({\rho }^{\epsilon }-1),v^{\epsilon })$  remains bounded in $L^{\infty }\left(\right[0,T],H^{s-1}\times H^s({\mathbb{T}}^d\left)\right)$  for some $T>0$  , independent of $\epsilon$  .

Convergence We use the change of variable used for rotating fluids (see [?] ), that removes the time oscillations: Considering the pair ${\tilde{\mathbf{u}}}^{\epsilon }=({\tilde{\beta }}^{\epsilon },{\tilde{\rho }}_1^{\epsilon })$  such that

$$\begin{array}{ccc}{\tilde{\beta }}^{\epsilon }+i{\tilde{\rho }}_1^{\epsilon }=e^{it/\epsilon }\left({\beta }^{\epsilon }+i{\rho }_1^{\epsilon }\right),& & \end{array}$$

we have

$$\begin{array}{ccc}{\partial }_t{\tilde{\mathbf{u}}}^{\epsilon }+v^{\epsilon }\cdot \nabla {\tilde{\mathbf{u}}}^{\epsilon }=T^{\epsilon },& & \end{array}$$

with

$$\begin{array}{ccc}& & \parallel T^{\epsilon }{\parallel }_{L^{\infty }\left(\right[0,T],H^{s-1})}\le C,\end{array}$$

$$\begin{array}{ccc}& & \parallel v^{\epsilon }{\parallel }_{L^{\infty }\left(\right[0,T],H^s)}\le C,\end{array}$$

$$\begin{array}{ccc}& & \parallel {\tilde{\mathbf{u}}}^{\epsilon }{\parallel }_{L^{\infty }\left(\right[0,T],H^{s-1})}\le C\end{array}$$

from the a priori bounds ( we still have $s\ge [d/2]+2$  ). Hence, ${\partial }_t{\tilde{\mathbf{u}}}^{\epsilon }$  is uniformly bounded at least in $L^2$  , and we deduce classically that $({\beta }^{\epsilon },{\rho }_1^{\epsilon })$  converges weakly to 0 in $[0,T]\times {\mathbb{T}}^d$  .

We assume now that $\Pi v_0^{\epsilon }$  , the solenoidal part of $v_0^{\epsilon }$  , converges weakly in $L^2$  to some limit ${\overline{v}}_0$  , hence it converges strongly in $H^{s^{\prime }}$  for $s^{\prime }<s$  . We check that ${\partial }_t{\omega }^{\epsilon }$  is bounded uniformly on $L^2$  under our assumptions. Hence from the a priori bound, ${\omega }^{\epsilon }$  converges (if necessary passing to a subsequence, but see the remark below) in $C\left(\right[0,T],H^{s^{\prime }-1}({\mathbb{T}}^d\left)\right)$  for all $s^{\prime }<s$  .

Let ${\overline{v}}^{\epsilon }$  be the unique vector field with zero average such that $\nabla \times {\overline{v}}^{\epsilon }={\omega }^{\epsilon }$  . Then, $\Pi v^{\epsilon }$  , the solenoidal part of $v^{\epsilon }$  is equal to ${\overline{v}}^{\epsilon }+c^{\epsilon }$  with $c^{\epsilon }$  a constant vector field. Since $\int {\rho }^{\epsilon }{v}^{\epsilon }=\int {\rho }_0^{\epsilon }{v}_0^{\epsilon }$  , and since ${\epsilon }^{-1}({\rho }^{\epsilon }-1)$  is bounded in $H^s$  , we have $lim{c}^{\epsilon }=lim\int v_0^{\epsilon }=lim\int \Pi v_0^{\epsilon }$  . Hence $\Pi v^{\epsilon }={\overline{v}}^{\epsilon }+c^{\epsilon }$  converges in $C\left(\right[0,T],H^{s^{\prime }}({\mathbb{T}}^d\left)\right)$  for all $s^{\prime }<s$  .

Decomposing $v^{\epsilon }$  as $v^{\epsilon }=\nabla q^{\epsilon }+\Pi v^{\epsilon }$  , we have $q^{\epsilon }={\Delta }^{-1}{\beta }^{\epsilon }$  that converges weakly to 0 in $L^2\left(\right[0,T]\times {\mathbb{T}}^d)$  . Hence, using the a priori bounds, in the vorticity equation

$$\begin{array}{ccc}{\partial }_t{\omega }^{\epsilon }+v^{\epsilon }\cdot \nabla {\omega }^{\epsilon }=-{\beta }^{\epsilon }{\omega }^{\epsilon },& & \end{array}$$

we can pass to the weak limit in $L^2\left(\right[0,T]\times {\mathbb{T}}^d)$  and state that $\omega =lim{\omega }^{\epsilon }$  satisfies in ${\mathcal{D}}^{\prime }\left(\right[0,T]\times {\mathbb{T}}^d)$  ,

$$\begin{array}{ccc}& & {\partial }_t\omega +\Pi v\cdot \nabla \omega =0,\end{array}$$

$$\begin{array}{ccc}& & \omega \left(0\right)={\omega }_0,\end{array}$$

where $\Pi v$  is the limit of the solenoidal part of $v^{\epsilon }$  . Moreover we have $\omega =\nabla \times \Pi v$  , hence $\Pi v$  is a solution to the incompressible Euler equation, with initial data ${\overline{v}}_0$  the limit of the solenoidal part of $v_0^{\epsilon }$  .

Remark 1. From the regularity of $v_0$  , the solution (in the distribution’s sense) to $\left(E\right)$  with initial data $\Pi v_0$  is unique. Therefore the whole sequence $\Pi v^{\epsilon }$  is converging.

Remark 2. Here we did not introduce the solution of the limit equation (incompressible Euler in this case), and chose to argue by compactness. This method looks simpler, however, we obtain less informations concerning the ’rate of convergence’ of the sequence $({\rho }^{\epsilon },v^{\epsilon })$  with respect to $\epsilon$  .

We gather those results in the following theorem:

Theorem 2.7 Let $s\ge [d/2]+2$  , let $({\rho }_0^{\epsilon },v_0^{\epsilon })$  be a sequence of initial data, such that $\int {\rho }_0^{\epsilon }=1$  .

Assume that $v_0^{\epsilon }$  is bounded in $H^s\left({\mathbb{T}}^d\right)$  and that ${\epsilon }^{-1}({\rho }_0^{\epsilon }-1)$  is bounded in $H^{s-1}\left({\mathbb{T}}^d\right)$  . Then

• 1. There exists $T>0$  such that a sequence of solutions to $\left(E{P}_{\epsilon }\right)$  with initial data $({\rho }_0^{\epsilon },v_0^{\epsilon })$  exists on $[0,T]$  , and the sequence $\left({\epsilon }^{-1}\right({\rho }^{\epsilon }-1),v^{\epsilon })$  is bounded in $L^{\infty }\left(\right[0,T],H^{s-1}\times H^s({\mathbb{T}}^d\left)\right)$  .
• 2. The potential part of $v^{\epsilon }$  , defined above, converges weakly to 0 in $L^2\left(\right[0,T]\times {\mathbb{T}}^d)$  .
• 3. If $\Pi v_0^{\epsilon }$  , the solenoidal part of $v_0^{\epsilon }$  converges to some ${\overline{v}}_0$  weakly in $L^2$  , then, for all $s^{\prime }<s$  , $\Pi v^{\epsilon }$  , the solenoidal part of $v^{\epsilon }$  , converges in $C\left(\right[0,T],H^{s^{\prime }}({\mathbb{T}}^d\left)\right)$  to $\overline{v}$  the solution of $\left(E\right)$  with initial data ${\overline{v}}_0$  .

3 The Euler-Monge-Ampere system

We consider here the following Euler-Monge-Ampere system denoted by $\left(EM{A}_{\epsilon }\right)$  :

$$\begin{array}{ccc}& & {\partial }_t\rho +\nabla \cdot \left(\rho v\right)=0,\end{array}$$

(25)

$$\begin{array}{ccc}& & {\partial }_{t}v+v\cdot \nabla v=\frac{\nabla \psi -x}{{\epsilon }^2},\end{array}$$

(26)

$$\begin{array}{ccc}& & det{D}^2\psi =\rho .\end{array}$$

(27)

The last equation must be understood in the following weak sense: $\psi$  is the only (up to a constant) convex function with $\psi -|x{|}^2/2$  being ${\mathbb{Z}}^d$  periodic such that

$$\begin{array}{ccc}\forall f\in C^0\left({\mathbb{T}}^d\right),{\int }_{{\mathbb{T}}^d}f(\nabla \psi )d\rho ={\int }_{{\mathbb{T}}^d}f\left(x\right)dx.& & \end{array}$$

(28)

This definition will be made precise in Theorem  3.1 and Definition  3.2 .

3.1 Geometric derivation of the Euler-Monge-Ampere system

This derivation has been introduced in [?] , [?] , [?] . We reproduce it for sake of completeness, but the reader only interested in the proof of convergence may skip this section.

3.1.1 The Euler equations of an incompressible perfect fluid

The motion of an inviscid and incompressible fluid in a domain $\Omega \subset {\mathbb{R}}^d$  is described by the Euler incompressible equation $\left(E\right)$  that we recall here:

$$\begin{array}{ccc}& & {\partial }_{t}v+(v\cdot \nabla )v=\nabla p,\end{array}$$

$$\begin{array}{ccc}& & \nabla \cdot v=0.\end{array}$$

Following Arnold (see [?] ), we have a formal interpretation of the Euler incompressible equations:

introducing $G(\Omega )$  the group of all volume preserving diffeomorphisms of $\Omega$  with jacobian determinant equal to 1, the Euler equations describe the geodesics of $G(\Omega )$  with length measured in the $L^2$  sense.

3.1.2 Approximate geodesics

A general strategy to define approximate geodesics along a manifold $M$  (in our case $M=G(\Omega )$  ) embedded in a Hilbert space $H$  (here $H=L^2(\Omega ,{\mathbb{R}}^d)$  ) is to introduce a penalty parameter $\epsilon >0$  and the following $unconstrained$  dynamical system in $H$ 

$$\begin{array}{ccc}{\partial }_{tt}X+\frac{1}{2{\epsilon }^2}{\nabla }_X\left(d^2(X,M))\right)=0.& & \end{array}$$

(29)

In this equation, the unknown $t\to X\left(t\right)$  is a curve in $H$  , $d(X,M)$  is the distance (in $H$  ) of $X$  to the manifold $M$  , i.e. in our case when $M=G(\Omega )$  ,

$$\begin{array}{ccc}d(X,G(\Omega \left)\right)={inf}_{g\in G(\Omega )}\parallel X-g{\parallel }_H,& & \end{array}$$

(30)

finally, ${\nabla }_X$  denotes the gradient operator in $H$  . This penalty approach has been introduced for the Euler equations by Brenier in [?] . It is similar-but not identicalto Ebin’s slightly compressible flow theory [?] , and is a natural extension of the theory of constrained finite dimensional mechanical systems [?] . Actually if $G(\Omega )$  were a smooth manifold, the result would be exactly the one of [?] , Theorem 2.7, but this is not the case, here because the $L^2$  metric is too weak. The penalized system is formally hamiltonian in variables $(X,{\partial }_{t}X)$  with hamiltonian (or energy) given by:

$$\begin{array}{ccc}E=\frac{1}{2}\parallel {\partial }_{t}X{\parallel }_H^2+\frac{1}{2{\epsilon }^2}{d}^2(X,G(\Omega \left)\right).& & \end{array}$$

Multiplying equation ( 29 ) by ${\partial }_{t}X$  , we get immediately that the energy is formally conserved.

Therefore it is plausible that the map $X\left(t\right)$  will remain close to $G(\Omega )$  if it is close at $t=0$  . A formal computation shows that, given a point $X$  for which there is a unique closest point ${\pi }_X$  to $X$  in the $H$  closure of $G(\Omega )$  , we have:

$$\begin{array}{ccc}{\nabla }_X\left(d(X,G)\right)=\frac{1}{d(X,G)}(X-{\pi }_X).& & \end{array}$$

(31)

Thus the equation ( 29 ) formally becomes:

$$\begin{array}{ccc}{\partial }_{tt}X+\frac{1}{{\epsilon }^2}(X-{\pi }_X)=0.& & \end{array}$$

To understand why solutions to such a system may approach geodesics along $G(\Omega )$  as $\epsilon$  goes to 0, just recall that, in the simple framework of a surface $S$  embedded in the 3 dimensional Euclidean space, a geodesic $t\to s\left(t\right)$  along $S$  is characterized by the fact that for every $t$  , the plane defined by $\left\{\dot{s}\right(t),\ddot{s}(t\left)\right\}$  is orthogonal to $S$  . In our case, ${\partial }_{tt}X\left(t\right)$  is orthogonal to $G(\Omega )$  thanks to ( 30 ) and $X\left(t\right)$  remains close to $G(\Omega )$  .

Notation Since we intend to work on the flat torus ${\mathbb{T}}^d$  we might consider ${\mathbb{Z}}^d$  additive mappings, i.e. mappings such that $$\forall \vec{p}\in {\mathbb{Z}}^d,X(.+\vec{p})=X(.)+\vec{p},$$  as well as periodic mappings (i.e. mappings from ${\mathbb{T}}^d$  into itself ).

Then given $m$  an additive mapping,we denote by $\hat{m}$  the naturally associated mapping on ${\mathbb{T}}^d$  . The following polar factorization Theorem is a periodic version of [?] , it has been discovered independently by [?]  and [?] .

Theorem 3.1 Let $X:{\mathbb{R}}^d\to {\mathbb{R}}^d$  be ${\mathbb{Z}}^d$  additive and assume that ${\rho }_X=X_{\#}dx$  has a density in $L^1\left(\right[0,1{]}^d)$  , then there exits an a.e. unique pair $(\nabla {\phi }_X,{\pi }_X)$  satisfying $$X=\nabla {\phi }_X\circ {\pi }_X$$  with ${\phi }_X$  a convex function such that ${\phi }_X\left(x\right)-|x{|}^2/2$  is ${\mathbb{Z}}^d$  periodic, and ${\pi }_X:{\mathbb{R}}^d\to {\mathbb{R}}^d$  additive such that ${\hat{\pi }}_X$  is measure preserving in ${\mathbb{T}}^d$  . Moreover we have $$\parallel X-{\pi }_X{\parallel }_{L^2\left(\right[0,1{]}^d)}=\parallel \hat{X}-{\hat{\pi }}_X{\parallel }_{L^2\left({\mathbb{T}}^d\right)}=d(\hat{X},G({\mathbb{T}}^d\left)\right)$$  and if ${\psi }_X$  is the Legendre transform of ${\phi }_X$  then $${\pi }_X=\nabla {\psi }_X\circ X.$$ 

Remark 1. The pair $({\phi }_X,{\psi }_X)$  is uniquely defined by the density ${\rho }_X=X_{\#}dx$  .

Remark 2. Important properties of the optimal potential: The periodicity of ${\phi }_X\left(x\right)-|x{|}^2/2$  implies that $\nabla {\phi }_X$  and $\nabla {\psi }_X$  are ${\mathbb{Z}}^d$  additive, and that ${\psi }_X-|x{|}^2/2$  is also ${\mathbb{Z}}^d$  periodic. This allows the following definition:

Definition 3.2 Let $\rho$  be a probability measure on ${\mathbb{T}}^d$  , then we denote $\phi \left[\rho \right]$  (resp. $\psi \left[\rho \right]$  ) the unique up to a constant convex function such that

$$\begin{array}{ccc}& & \phi \left[\rho \right]-|\cdot {|}^2/2\text{is}{\mathbb{Z}}^d\text{periodic},\end{array}$$

$$\begin{array}{ccc}& & \forall f\in C^0\left({\mathbb{T}}^d\right),{\int }_{{\mathbb{T}}^d}f\left(\hat{\nabla }\phi \right[\rho \left]\right(x\left)\right)dx={\int }_{{\mathbb{T}}^d}f\left(x\right)d\rho \left(x\right)\end{array}$$

(resp. its Legendre transform).

Remark. We recover thus that $\psi \left[\rho \right],\phi \left[\rho \right]$  will be generalized solutions of the following Monge-Ampere equations

$$\begin{array}{ccc}& & det{D}^2\psi =\rho ,\end{array}$$

$$\begin{array}{ccc}& & \rho (\nabla \phi )det{D}^2\phi =1.\end{array}$$

3.1.3 Result

We are now ready to state the main result of this section.

Theorem 3.3 Let $s\in \mathbb{N}$  with $s\ge [d/2]+2$  . Let ${\overline{v}}_0$  be a divergence-free vector field on ${\mathbb{T}}^d$  , let $(\overline{v},p)$  be a smooth solution of the Euler incompressible system ( 8 ) on $[0,T]\times {\mathbb{T}}^d$  , with initial data ${\overline{v}}_0$  and satisfying $\overline{v}\in L^{\infty }\left(\right[0,T],H^{s+1}({\mathbb{T}}^d\left)\right)$  . Let $(v_0^{\epsilon },{\rho }_0^{\epsilon })$  be a sequence of initial data, with $\int {\rho }_0^{\epsilon }=1$  , and such that $\left({\epsilon }^{-1}(v_0^{\epsilon }-{\overline{v}}_0),{\epsilon }^{-2}({\rho }_0^{\epsilon }-1)\right)$  is bounded in $H^s\times H^{s-1}\left({\mathbb{T}}^d\right)$  . Then there exists a sequence $(v^{\epsilon },{\rho }^{\epsilon })$  of solutions to $\left(EM{A}_{\epsilon }\right)$  with initial data $(v_0^{\epsilon },{\rho }_0^{\epsilon })$  belonging to $L^{\infty }\left(\right[0,T_{\epsilon }],H^s\times H^{s-1}({\mathbb{T}}^d\left)\right)$  with ${liminf}_{\epsilon \to 0}{T}_{\epsilon }\ge T$  . Moreover for any $T^{\prime }<T$  and $\epsilon$  small enough, $\left({\epsilon }^{-1}(v^{\epsilon }-\overline{v}),{\epsilon }^{-2}({\rho }^{\epsilon }-1)\right)$  is bounded in $L^{\infty }\left(\right[0,T^{\prime }],H^s\times H^{s-1}({\mathbb{T}}^d\left)\right)$  . Finally when $T=+\infty$  , $T_{\epsilon }$  goes to infinity.

Before entering the proof of this result, we need some preliminary results concerning the linearization of the Monge-Ampere operator.

3.2 Linearization of the Monge-Ampere operator in $H^s$  norm

This section is devoted to the proof of the following Theorem:

Theorem 3.4 Let $\rho$  be a probability measure on ${\mathbb{T}}^d$  , $d\le 3$  , let $\psi$  satisfy

$$\begin{array}{ccc}det{D}^2\psi =\rho & & \end{array}$$

in the sense of Definition  3.2 . Then, there exists ${\epsilon }_0>0$  such that if $\parallel \rho -1{\parallel }_{H^2\left({\mathbb{T}}^d\right)}\le {\epsilon }_0$  , for any $s\in \mathbb{N}$  with $s>d/2$  there exists $C\left(s\right)$  that satisfies

$$\begin{array}{ccc}& & \parallel D^2\psi -I{\parallel }_{H^s\left({\mathbb{T}}^d\right)}\le C\left(s\right)\parallel \rho -1{\parallel }_{H^s\left({\mathbb{T}}^d\right)},\end{array}$$

(32)

$$\begin{array}{ccc}& & \parallel (\Delta \psi -d)-(\rho -1){\parallel }_{H^s\left({\mathbb{T}}^d\right)}\le C\left(s\right)\parallel (\rho -1){\parallel }_{H^s\left({\mathbb{T}}^d\right)}^2.\end{array}$$

(33)

3.2.1 Preliminary results

We first state the following result obtained from [?]  on the regularity of solutions to Monge-Ampere equation, adapted to the periodic case. This result will be the starting point of the proof of Theorem  3.4 .

Theorem 3.5 Let $\rho \in C^{\alpha }\left({\mathbb{T}}^d\right)$  for some $\alpha >0$  , with $0<m\le \rho \le M$  , be a probability measure on ${\mathbb{T}}^d$  , let $\psi =\psi \left[\rho \right]$  in the sense of Definition  3.2 . Then $\psi$  is a classical solution of

$$\begin{array}{ccc}det{D}^2\psi =\rho & & \end{array}$$

and satisfies for any ${\alpha }^{\prime }<\alpha$  :

$$\begin{array}{ccc}& & \parallel \nabla \psi \left(x\right)-x{\parallel }_{L^{\infty }}\le C\left(d\right)=\sqrt{d}/2,\end{array}$$

(34)

$$\begin{array}{ccc}& & \parallel D^2\psi {\parallel }_{C^{{\alpha }^{\prime }}}\le K(m,M,\parallel \rho {\parallel }_{C^{\alpha }},\alpha ,{\alpha }^{\prime }).\end{array}$$

(35)

Then we state a classical result of elliptic regularity that we will need during the course of the proof. It can be found in [?] , Theorem 9.11.

Theorem 3.6 Let $\Omega$  be an open set in ${\mathbb{R}}^d$  , $u\in W_{loc}^{2,p}(\Omega )\cap L^p(\Omega ),1<p<\infty$  , be a strong solution of the equation

$$\begin{array}{ccc}{\sum }_{i,j=1}^d{a}^{ij}{\partial }_{ij}u=f& & \end{array}$$

in $\Omega$  where the coefficients $a^{ij}$  satisfy

$$\begin{array}{ccc}& & a^{ij}\in C^0(\Omega ),f\in L^p(\Omega );\end{array}$$

$$\begin{array}{ccc}& & \lambda |\xi {|}^2\le a^{ij}{\xi }_i{\xi }_j\le \Lambda |\xi {|}^2\forall \xi \in {\mathbb{R}}^d,\end{array}$$

for $i,j=1..d$  , with $0<\lambda ,\Lambda <\infty$  . Then for any ${\Omega }^{\prime }\subset \subset \Omega$  ,

$$\begin{array}{ccc}\parallel u{\parallel }_{W^{2,p}\left({\Omega }^{\prime }\right)}\le C(\parallel u{\parallel }_{L^p\left({\Omega }^{\prime }\right)}+\parallel f{\parallel }_{L^p\left({\Omega }^{\prime }\right)}),& & \end{array}$$

where $C$  depends on $d,p,\lambda ,\Lambda ,{\Omega }^{\prime },\Omega$  and the moduli of continuity of the coefficients $a^{ij}$  on ${\Omega }^{\prime }$  .

3.2.2 Proof of Theorem  3.4 

Sketch of the proof We assume here $d=3$  . We recall that $\psi$  satisfies $det{D}^2\psi =\rho$  . We first have to prove that $\rho \in H^s$  implies $D^2\psi \in H^s$  . We will proceed by induction. We recall first that for $A,B$  two $d\times d$  matrices, we have the following expansion:

$$\begin{array}{ccc}det(A+tB)=detA+t\text{trace}\left(A_{com}^{t}B\right)+o\left(t\right),& & \end{array}$$

where $A_{com}$  is the matrix whose elements are the minors of $A$  , or co-matrix of $A$  . Hence the elements of $A_{com}$  are polynomials of degree $d-1$  in the elements of $A$  . When $A$  is invertible, we have $A_{com}=detA[A{]}^{-1}$  .

Differentiating $s$  times the Monge-Ampere equation, and denoting $M$  the co-matrix of $D^2\psi$  , we will have

$$\begin{array}{ccc}\text{trace}\left(M{D}^2{\partial }^s\psi \right)+\mathcal{T}={\partial }^s\rho ,& & \end{array}$$

where the first term contains the highest derivatives, and $\mathcal{T}$  will consist of products involving three derivatives of $\psi$  . The order of each derivative will smaller or equal to $s-1$  , and the sum of the three orders will be equal to $s$  . By a careful analysis, this product will be controlled in $L^2$  by $\parallel D^2\psi {\parallel }_{H^{s-1}}$  and $\parallel \rho {\parallel }_{H^s}$  , using Sobolev injections. From Theorem  3.5 , assuming a minimal regularity for $\rho$  (i.e.

the bound ( 37 )), $D^2\psi$  , and therefore $M$  will be continuous elliptic matrices. Hence ${\partial }^s\psi$  solves an elliptic problem, with continuous coefficients, and we will use the Theorem  3.6 to obtain $D^2{\partial }^s\psi \in L^2$  .

This intermediate step will be done in Lemma  3.7 ; using this a priori estimate and the continuity method, we will obtain the estimate ( 32 ).

Then, the expansion $det(I+D^2\phi )=1+\Delta \phi +P\left({\partial }_{ij}\phi \right)$  , where $P$  is a polynomial in ${\partial }_{ij}\phi$  whose terms are of degree two or three (when $d=3$  ), will yield ( 33 ).

Rigorous proof We recall that $\psi$  satisfies

$$\begin{array}{ccc}& & det{D}^2\psi =\rho ,\end{array}$$

(36)

$$\begin{array}{ccc}& & \parallel \rho -1{\parallel }_{H^2}\le {\epsilon }_0,\end{array}$$

(37)

for some ${\epsilon }_0$  to be chosen later. We suppose $d=3$  and the proof can be reproduced in the case $d=2$  with minor modifications. The parameter ${\epsilon }_0$  is chosen such that ( 37 ) implies

$$\begin{array}{ccc}{\lambda }_1\le \rho \le {\lambda }_2& & \end{array}$$

for some ${\lambda }_1>0,{\lambda }_2>0$  . Note also that thanks to ( 37 ), $\rho$  is in $C^{\alpha }$  for $\alpha =\frac{1}{2}$  . Then from Theorem  3.5 , $D^2\psi \in C^{{\alpha }^{\prime }}$  with ${\alpha }^{\prime }<\alpha$  . Note also that since $\rho \in [{\lambda }_1,{\lambda }_2]$  and using equation ( 36 ), $D^2\psi \in C^{{\alpha }^{\prime }}$  implies that $[D^2\psi {]}^{-1}\in C^{{\alpha }^{\prime }}$  , and thus $M^{ij}$  the co-matrix of $D^2\psi$  is uniformly elliptic and in $C^{{\alpha }^{\prime }}$  .

We first prove by induction that if $\gamma \in {\mathbb{N}}^d$  then $\rho \in H^{\left|\gamma \right|}$  implies $D^2\psi \in H^{\left|\gamma \right|}$  . It can be checked during the proof that this bound will be uniform under the condition ( 37 ) for ${\epsilon }_0$  small enough.

Lemma 3.7 Under assumption ( 37 ), for any $\gamma \in {\mathbb{N}}^d$  , $\rho \in H^{\left|\gamma \right|}$  implies that ${\partial }^{\gamma }{D}^2\psi \in L^2$  . If moreover $\rho \in W^{\left|\gamma \right|,6}$  then ${\partial }^{\gamma }{D}^2\psi \in L^6$  .

Proof. This lemma will be proved by induction. We first deal with the cases $\left|\gamma \right|=0,1,2$  .

The case $\gamma =0$  is a consequence of Theorem  3.5 .

For $\left|\gamma \right|=1$  we differentiate ( 36 ) with respect to $x_{\nu }$  , to obtain

$$\begin{array}{ccc}{M}^{ij}{\partial }_{ij}\left({\partial }_{\nu }\psi \right)={\partial }_{\nu }\rho ,& & \end{array}$$

(38)

with $M^{ij}$  the co-matrix of ${\partial }_{ij}\psi$  . Then if ${\partial }_{\nu }\rho \in L^2$  , by Theorem  3.6 , ${\partial }_{\nu }\psi \in W^{2,2}$  . If ${\partial }_{\nu }\rho \in L^6$  we also get that ${\partial }_{\nu }\psi \in W^{2,6}$  .

For $\left|\gamma \right|=2$  differentiating once more with respect to $x_{\beta }$  we obtain

$$\begin{array}{ccc}{M}^{ij}{\partial }_{ij}\left({\partial }_{\nu \beta }\psi \right)+\left({\partial }_{\beta }{M}^{ij}\right){\partial }_{ij}\left({\partial }_{\nu }\psi \right),={\partial }_{\nu \beta }\rho & & \end{array}$$

(39)

still with $M^{ij}$  the co-matrix of ${\partial }_{ij}\psi$  . Suppose that $\rho \in H^2$  , then $W^{2,2}\subset W^{1,\frac{2d}{d-2}}=W^{1,6}$  if $d=3$  , and ${\partial }_{\nu }\psi \in W^{2,6}$  . The term ${\partial }_{\beta }{M}^{ij}$  is a sum of terms of the form ${\partial }_{ij}\left({\partial }_{\beta }\psi \right){\partial }_{kl}\psi$  and the second term of the left hand side of ( 39 ) is thus bounded in $L^2$  . Then once again by Theorem  3.6 one gets that $D^2{\partial }_{\nu \beta }\psi \in L^2$  if ${\partial }_{\nu \beta }\rho \in L^2$  .

Moreover if $D^2\rho \in L^6$  then ${\partial }_{\nu }\rho \in C^{\alpha }$  for some $\alpha >0$  . Using ( 38 ) and Schauder interior estimates (see [?] , Theorem 6.2.), we obtain $D^2{\partial }_{\nu }\psi \in C^{{\alpha }^{\prime }}$  . Thus $\left({\partial }_{\beta }{M}^{ij}\right){\partial }_{ij}\left({\partial }_{\nu }\psi \right)\in C^{{\alpha }^{\prime }}$  . From ( 39 ) and Theorem  3.6 we obtain ${\partial }_{\nu \beta }{D}^2\psi \in L^6$  .

As we just saw, Lemma  3.7 is true for $\left|\gamma \right|=0,1,2$  . We assume that it holds for all $\gamma$  with $\left|\gamma \right|\le n$  for some $n\ge 2$  . Take now $\left|\gamma \right|=n+1\ge 3$  , $\rho \in H^{\left|\gamma \right|}$  , and apply ${\partial }^{\gamma }$  to ( 36 ):

$$\begin{array}{ccc}{M}^{ij}{\partial }_{ij}{\partial }^{\gamma }\psi +{\sum }_{\begin{array}{c}{\gamma }_1+{\gamma }_2+{\gamma }_3=\gamma \\ \left|\gamma \right|-1\ge \left|{\gamma }_1\right|\ge \left|{\gamma }_2\right|\ge \left|{\gamma }_3\right|\end{array}}*{\partial }_{ij}{\partial }^{{\gamma }_1}\psi {\partial }_{kl}{\partial }^{{\gamma }_2}\psi {\partial }_{mp}{\partial }^{{\gamma }_3}\psi ={\partial }^{\gamma }\rho & & \end{array}$$

(40)

with $*$  some constant coefficients. We call $\mathcal{T}$  the second term of the left hand side of ( 40 ). Since $\rho \in H^{\left|\gamma \right|}$  , $\rho \in W^{\left|\gamma \right|-1,6}$  , and we have ${\partial }^{\alpha }{D}^2\psi \in L^6\left({\mathbb{T}}^d\right)$  for any $\left|\alpha \right|\le n$  using the induction hypothesis.

Therefore $\mathcal{T}\in L^2$  and since ${\partial }^{\gamma }\rho \in L^2$  we obtain $M^{ij}{\partial }_{ij}{\partial }^{\gamma }\psi \in L^2$  . Using Theorem  3.6 it follows that ${\partial }^{\gamma }{D}^2\psi \in L^2$  .

Remember that $\left|\gamma \right|\ge 3$  thus $\left|{\gamma }_3\right|\le \frac{1}{3}\left|\gamma \right|\le \gamma -2$  , and $\left|{\gamma }_2\right|\le \frac{1}{2}\left|\gamma \right|\le \gamma -2$  . Since $d=3$  , we have $H^2\subset C^{\alpha }$  for some $\alpha >0$  and thus ${\partial }_{kl}{\partial }^{{\gamma }_2}\psi$  , ${\partial }_{mp}{\partial }^{{\gamma }_3}\psi$  are in $C^{\alpha }$  , moreover $H^1\subset L^6$  and since $\left|{\gamma }_1\right|\le \left|\gamma \right|-1$  , ${\partial }_{ij}{\partial }^{{\gamma }_1}\psi$  is in $L^6$  . Therefore $\mathcal{T}$  is in $L^6$  . If ${\partial }^{\gamma }\rho \in L^6$  we have ${\partial }^{\gamma }{D}^2\psi$  in $L^6$  . Hence the lemma holds for all $\left|\gamma \right|\le n+1$  . This achieves the proof of Lemma  3.7 .

$\square$  Now by induction on $\left|\gamma \right|$  we prove ( 32 ) and ( 33 ). From ( 37 ) we have $\parallel \rho -1{\parallel }_{L^2}\le {\epsilon }_0$  small. Take $\psi =|x{|}^2/2+\phi$  solution of $det{D}^2\psi =\rho$  with $\phi$  periodic and ${\int }_{{\mathbb{T}}^d}\phi =0$  . We begin to show that $\parallel \phi {\parallel }_{C^{2,\alpha }}$  for some $\alpha >0$  , is controlled by $\parallel \rho -1{\parallel }_{H^2}$  . Indeed, the periodic solution of

$$\begin{array}{ccc}det(I+D^2\phi )=\rho & & \end{array}$$

can be built by the continuity method (see [?] ). Starting from ${\rho }_0=1,{\phi }_0=0$  , we use the implicit function Theorem to obtain the solution ${\phi }_t$  of

$$\begin{array}{ccc}det(I+D^2{\phi }_t)=t\rho +(1-t).& & \end{array}$$

(41)

For this we differentiate ( 41 ) with respect to $t$  , to obtain

$$\begin{array}{ccc}{M}_t^{ij}{\partial }_{ij}{\partial }_t{\phi }_t=\rho -1,& & \end{array}$$

(42)

for $t\in [0,1]$  , where $M_t$  is the co-matrix of $I+D^2{\phi }_t$  . We know, from the a priori estimate of Theorem  3.5 , that for $\rho \in C^{1/2}$  , $D^2{\phi }_t$  and therefore $M_t^{ij}$  are $C^{\alpha }$  elliptic matrices for all $\alpha <1/2$  .

To see why ( 42 ) indeed admits a unique (up to a constant) periodic solution, we recall that $M$  is the comatrix of a Hessian matrix, therefore it is ’divergence-free’:

$$\begin{array}{ccc}\forall i\in [1..d],{\sum }_j{M}^{ij}=0.& & \end{array}$$

Hence equation ( 42 ) can be rewritten in divergence form

$$\begin{array}{ccc}{\sum }_{i,j}{\partial }_i\left(M_t^{ij}{\partial }_j{\partial }_t{\phi }_t\right)=\rho -1,& & \end{array}$$

and the operator $\mathcal{ℒ}=M_t^{ij}{\partial }_{ij}\cdot$  is a self adjoint operator on $H^1\left({\mathbb{T}}^d\right)$  and induces a bounded coercive bilinear form on $H_0^1\left({\mathbb{T}}^d\right)$  , where the subscript 0 means that we impose the mean value to be 0. Then the existence/uniqueness of a solution to ( 42 ) in $H_0^1$  follows by Lax-Milgram Theorem. Hence, ${\partial }_t{\phi }_t$  is the unique (up to a constant) periodic solution of the above elliptic problem, and from Schauder interior estimates we obtain

$$\begin{array}{ccc}\parallel {\partial }_t\phi {\parallel }_{C^{2,\alpha }}& \le & C\parallel \rho -1{\parallel }_{C^{1/2}}\end{array}$$

$$\begin{array}{ccc}& \le & C\parallel \rho -1{\parallel }_{H^2},\end{array}$$

uniformly in $t\in [0,1]$  , and finally for $\alpha <1/2$  ,

$$\begin{array}{ccc}\parallel \phi {\parallel }_{C^{2,\alpha }}\le C\parallel \rho -1{\parallel }_{H^2}.& & \end{array}$$

Then we have

$$\begin{array}{ccc}det(I+D^2\phi )=1+\Delta \phi +R_{ij}{\partial }_{ij}\phi & & \end{array}$$

where $R$  is a symmetric matrix whose coefficients are polynomials in ${\partial }_{ij}\phi$  of degree larger or equal to 1. The norms $\parallel R_{ij}{\parallel }_{C^{\alpha }}$  are controlled by $\parallel \rho -1{\parallel }_{C^{\alpha }}\le {\epsilon }_0$  , hence, for ${\epsilon }_0$  small enough, the matrix ${\delta }_{ij}+R_{ij}$  is uniformly bounded, elliptic, and $C^{\alpha }$  continuous. Since $\phi$  satisfies

$$\begin{array}{ccc}({\delta }_{ij}+R_{ij}){\partial }_{ij}\phi =\rho -1& & \end{array}$$

(43)

it follows from Theorem  3.6 that $\parallel {\partial }_{ij}\phi {\parallel }_{L^2}\le C\parallel \rho -1{\parallel }_{L^2}$  and this proves ( 32 ) for $\gamma =0$  .

If $\left|\gamma \right|=1$  , we have

$$\begin{array}{ccc}\left(M^{ij}\right){\partial }_{ij}{\partial }_{\nu }\phi ={\partial }_{\nu }\rho & & \end{array}$$

with $M$  uniformly bounded, elliptic and $C^{\alpha }$  continuous. For the same reasons we have $\parallel {\partial }_{ij}{\partial }_{\nu }\phi {\parallel }_{L^2}\le C\parallel {\partial }_{\nu }\rho {\parallel }_{L^2}$  .

If $\left|\gamma \right|=2$  , we do as in the Proof of Lemma  3.7 : using ( 38 ,  39 ) and keeping track of the bounds, we get

$$\begin{array}{ccc}\parallel {\partial }_{ij}{\partial }_{\nu \beta }\phi {\parallel }_{L^2}& \le & C\parallel {\partial }_{\nu \beta }\rho {\parallel }_{L^2}+C\parallel \rho -1{\parallel }_{H^2}^2\end{array}$$

$$\begin{array}{ccc}& \le & C\parallel \rho -1{\parallel }_{H^2},\end{array}$$

with $C$  uniform under the assumption of Lemma  3.7 .

If $\left|\gamma \right|\ge 3$  , we go back to equation ( 40 ): $\mathcal{T}$  is a sum of terms which contain all a product of at least two derivatives of $\psi$  of degree higher than 3. Since $D^3\psi =D^3\phi$  , we have

$$\begin{array}{ccc}{\partial }_{ij}{\partial }^{{\gamma }_1}\psi {\partial }_{kl}{\partial }^{{\gamma }_2}\psi {\partial }_{mp}{\partial }^{{\gamma }_3}\psi ={\partial }_{ij}{\partial }^{{\gamma }_1}\phi {\partial }_{kl}{\partial }^{{\gamma }_2}\phi {\partial }_{mp}{\partial }^{{\gamma }_3}\psi .& & \end{array}$$

We assume by induction that $\parallel {\partial }_{ij}\phi {\parallel }_{H^{\left|\gamma \right|-1}}\le C\parallel \rho -1{\parallel }_{H^{\left|\gamma \right|-1}}$  . Since $\rho -1\in H^{\left|\gamma \right|}$  we also have that $D^2\psi \in H^{\left|\gamma \right|}$  , with a uniform bound thanks to Lemma  3.7 . We remember that when $\left|\gamma \right|\ge 3$  , we have $\left|{\gamma }_2\right|\le \left|\gamma \right|-2,\left|{\gamma }_3\right|\le \left|\gamma \right|-2$  , therefore, using the injection of $H^2$  in $C^{1/2}$  when $d=3$  , ${\partial }_{kl}{\partial }^{{\gamma }_2}\phi$  and ${\partial }_{mp}{\partial }^{{\gamma }_3}\psi$  are uniformly bounded in $L^{\infty }$  . We obtain that

$$\begin{array}{ccc}\parallel {\partial }_{ij}{\partial }^{{\gamma }_1}\psi {\partial }_{kl}{\partial }^{{\gamma }_2}\psi {\partial }_{mp}{\partial }^{{\gamma }_3}\psi {\parallel }_{L^2}\le C\parallel \rho -1{\parallel }_{H^{\left|\gamma \right|-1}}.& & \end{array}$$

Then ${\partial }^{\gamma }\phi$  satisfies

$$\begin{array}{ccc}\left(M^{ij}\right){\partial }_{ij}{\partial }^{\gamma }\phi ={\partial }^{\gamma }\rho -\mathcal{T}& & \end{array}$$

thus, using Theorem  3.6 ,

$$\begin{array}{ccc}\parallel {\partial }^{\gamma }{D}^2\phi {\parallel }_{L^2}& \le & C(\parallel \rho -1{\parallel }_{H^{\left|\gamma \right|-1}}+\parallel {\partial }^{\gamma }\rho {\parallel }_{L^2})\end{array}$$

$$\begin{array}{ccc}& \le & C(\parallel \rho -1{\parallel }_{H^{\left|\gamma \right|}}),\end{array}$$

and we conclude that

$$\begin{array}{ccc}\parallel D^2\psi -I{\parallel }_{H^s}\le C\left(s\right)\parallel \rho -1{\parallel }_{H^s}& & \end{array}$$

for $s\in \mathbb{N},s\ge 2$  and under condition ( 37 ); thus ( 32 ) is obtained.

Using Proposition  2.6 and the fact that ${\partial }_{mp}{\partial }^{{\gamma }_3}\psi$  is uniformly bounded in $L^{\infty }$  thanks to Lemma  3.7 , we can also obtain that

$$\begin{array}{ccc}\parallel {\partial }_{ij}{\partial }^{{\gamma }_1}\phi {\partial }_{kl}{\partial }^{{\gamma }_2}\phi {\partial }_{mp}{\partial }^{{\gamma }_3}\psi {\parallel }_{L^2}& \le & C\parallel D^2\phi {\parallel }_{H^{\left|\gamma \right|}}^2\end{array}$$

$$\begin{array}{ccc}& \le & C\parallel \rho -1{\parallel }_{H^{\left|\gamma \right|}}^2\end{array}$$

for $\left|\gamma \right|\ge 3$  . (When $\left|\gamma \right|=2$  , the estimate holds also, but not using Proposition  2.6 .) Therefore, for all $\gamma$  , we have $\parallel \mathcal{T}{\parallel }_{L^2}\le C\parallel \rho -1{\parallel }_{H^{\left|\gamma \right|}}^2$  .

To conclude ( 33 ), we now write $det(I+D^2\phi )=\rho$  under the form

$$\begin{array}{ccc}\Delta \phi =\rho -1+P\left({\partial }_{ij}\phi \right),& & \end{array}$$

with $P$  consisting of products of two or three second derivatives of $\phi$  . Hence, under assumption ( 37 ), using Proposition  2.6 , we have for $s>d/2$  , $\parallel P\left({\partial }_{ij}\phi \right){\parallel }_{H^s}\le C\parallel D^2\phi {\parallel }_{H^s}^2$  . Using the bound ( 32 ), we conclude that for $s\in \mathbb{N},s>d/2$  ,

$$\begin{array}{ccc}\parallel \Delta \phi -(\rho -1){\parallel }_{H^s}\le \parallel \rho -1{\parallel }_{H^s}^2.& & \end{array}$$

and Theorem  3.4 is proved.

$\square$ 

3.3 Energy estimates and proof of the convergence

The proof of the energy estimates for Euler-Monge-Ampere is much inspired from the proof of Theorem  2.1 for the following reason: by taking the divergence of equation ( 26 ) one gets:

$$\begin{array}{ccc}{\partial }_t(\nabla \cdot v)+v\cdot \nabla (\nabla \cdot v)+{\partial }_i{v}^j{\partial }_j{v}^i=\frac{\Delta \psi -d}{{\epsilon }^2}.& & \end{array}$$

Suppose that $\rho$  is close to 1 at an order ${\epsilon }^2$  as is the case for Euler-Poisson, we guess (from Theorem  3.4 ) that we have the following:

$$\begin{array}{ccc}& & \psi =|x{|}^2/2+{\epsilon }^2\phi ,\end{array}$$

$$\begin{array}{ccc}& & \rho =det{D}^2\psi =1+{\epsilon }^2\Delta \phi +O\left({\epsilon }^4\right),\end{array}$$

and thus $\Delta \psi =d+{\epsilon }^2\Delta \phi =d+\rho -1+O\left({\epsilon }^4\right)$  . Therefore we expect that

$$\begin{array}{ccc}{\partial }_t(\nabla \cdot v)+v\cdot \nabla (\nabla \cdot v)+{\partial }_i{v}^j{\partial }_j{v}^i=\frac{\rho -1}{{\epsilon }^2}+O\left({\epsilon }^2\right),& & \end{array}$$

and that the technique of Theorem  2.1 will apply.

Before performing the div-curl decomposition for the energy estimates, we need to establish the analog of Lemma  2.2 in the present case, so that Lemma  2.3 and its Corollary  2.4 hold.

Lemma 3.8 Let $(\rho ,v)$  be a solution to $\left(EM{A}_{\epsilon }\right)$  . Then, the total momentum $\int \rho (t,x)v(t,x)dx$  does not depend on time.

Proof. We proceed similarly as in the proof of Lemma  2.2 . We need to establish that $${\int }_{{\mathbb{T}}^d}\rho (t,x)(\nabla \psi (t,x)-x)dx\equiv 0.$$  For this we use Definition  3.2 . For $f\left(x\right)=\nabla \psi \left(x\right)-x$  , we have

$$\begin{array}{ccc}{\int }_{{\mathbb{T}}^d}\rho f& =& {\int }_{{\mathbb{T}}^d}f(\nabla \hat{\phi })\end{array}$$

$$\begin{array}{ccc}& =& {\int }_{{\mathbb{T}}^d}\nabla \psi (\nabla \phi )-\nabla \phi \end{array}$$

$$\begin{array}{ccc}& =& {\int }_{{\mathbb{T}}^d}x-\nabla \phi \end{array}$$

$$\begin{array}{ccc}& =& 0,\end{array}$$

where we have used at the third line that, for $\psi ,\phi$  Legendre transform of each other, $\nabla \psi (\nabla \phi )=id$  , and at the last line that $\phi -|x{|}^2/2$  is periodic.

$\square$  Hence we have shown that one can retrieve $v$  from the initial value of $\int \rho v$  , and $\nabla \cdot v,\text{curl}v$  .

General framework We perform the same div-curl decomposition as in the Euler-Poisson case.

We then express the difference between the solution of $\left(EM{A}_{\epsilon }\right)$  and the limiting solution : either the solution of $\left(E\right)$  or the solution of $\left(EM{A}_{\epsilon }\right)$  . After having applied a proper scaling to this difference, our solution is now described by a vector $\mathbf{u}$  whose first component (that can be a vector if $d=3$  ) is the rescaled vorticity, and whose last two components are a rescaled divergence and rescaled density fluctuation. For this perturbation we will obtain

$$\begin{array}{ccc}& & {\partial }_t{\mathbf{u}}^{\epsilon }+{\sum }_i{v}^i{\partial }_i{\mathbf{u}}^{\epsilon }+R^{\epsilon }{\mathbf{u}}^{\epsilon }=Q^{\epsilon }\left({\mathbf{u}}^{\epsilon }\right),\end{array}$$

where we still use

$$\begin{array}{ccc}{R}^{\epsilon }=\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& -\frac{1}{\epsilon }\\ 0& \frac{1}{\epsilon }& 0\end{array}\right).& & \end{array}$$

For the source term $Q^{\epsilon }$  , we have $\parallel Q^{\epsilon }\parallel \le C(1+\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^s}+{\delta }_{\epsilon }\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^s})$  where ${\delta }_{\epsilon }$  goes to 0 when $\epsilon$  goes to 0, and the constant $C$  depends on the regularity of the limiting field.

Regularity of the limiting field The form of the source term will vary under the circumstances, but the general idea is that in order to bound $Q^{\epsilon }$  in $H^s$  , we will need the limiting velocity to be bounded in $H^{s+2}$  and the limiting density to be bounded in $H^{s+1}$  . Remember that the $H^s$  norm of $\mathbf{u}$  controls the norm of $(\rho ,v)$  in $H^s\times H^{s+1}$  , thus the limiting field must have one more derivative bounded than the order of the energy estimate.

A Gronwall’s lemma then yields a control on the perturbation that holds on a range of time $[0,T_{\epsilon }]$  , where $T_{\epsilon }\to T$  , $T$  being the time of existence of a smooth solution for the limiting equation.

Convergence to Euler, two dimensional case Doing the same change of variables as in the proof of Theorem  2.1 

$$\begin{array}{ccc}& & \nabla \cdot v=\epsilon {\beta }_1,\end{array}$$

$$\begin{array}{ccc}& & \rho =1+{\epsilon }^2{\rho }_1,\end{array}$$

$$\begin{array}{ccc}& & \text{curl}v=\omega =\overline{\omega }+\epsilon {\omega }_1,\end{array}$$

we obtain:

$$\begin{array}{ccc}& & {\partial }_t(\overline{\omega }+\epsilon {\omega }_1)+v\cdot \nabla (\overline{\omega }+\epsilon {\omega }_1)=-(\overline{\omega }+\epsilon {\omega }_1)\epsilon \beta ,\end{array}$$

(44)

$$\begin{array}{ccc}& & {\partial }_t\epsilon {\beta }_1+v\cdot \nabla \epsilon {\beta }_1+2\epsilon {\partial }_i{\overline{v}}^j{\partial }_j{v}_1^i+\epsilon {\partial }_i{v}_1^j{\partial }_j{v}_1^i=\frac{\Delta \psi -d}{{\epsilon }^2}-{\partial }_i{\overline{v}}^j{\partial }_j{\overline{v}}^i,\end{array}$$

(45)

$$\begin{array}{ccc}& & {\partial }_t{\epsilon }^2{\rho }_1+v\cdot \nabla {\epsilon }^2{\rho }_1=-(1+{\epsilon }^2{\rho }_1)\epsilon {\beta }_1.\end{array}$$

(46)

Now we define $\Xi$  by

$$\begin{array}{ccc}\Delta \psi -d={\epsilon }^2{\rho }_1+{\epsilon }^4\Xi ,& & \end{array}$$

and from Theorem  3.4 inequality ( 33 ), we have, if $s\ge 2$  , $\parallel \Xi {\parallel }_{H^s}\le C\parallel {\rho }_1{\parallel }_{H^s}^2$  . The system can here be written in the following way:

$$\begin{array}{ccc}& & {\partial }_t{\mathbf{u}}^{\epsilon }+{\sum }_i{v}^i{\partial }_i{\mathbf{u}}^{\epsilon }+R^{\epsilon }{\mathbf{u}}^{\epsilon }=S^{\epsilon }\left({\mathbf{u}}^{\epsilon }\right)+V^{\epsilon },\end{array}$$

$$\begin{array}{ccc}& & {\mathbf{u}}^{\epsilon }\left(0\right)={\mathbf{u}}_0^{\epsilon },\end{array}$$

still with

$$\begin{array}{ccc}{\mathbf{u}}^{\epsilon }=\left(\begin{array}{c}{\omega }_1\\ {\beta }_1\\ {\tilde{\rho }}_1\end{array}\right),R^{\epsilon }=\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& -\frac{1}{\epsilon }\\ 0& \frac{1}{\epsilon }& 0\end{array}\right),& & \end{array}$$

with the same $S^{\epsilon }$  as in the Euler-Poisson case, ${\tilde{\rho }}_1={\rho }_1-\Delta p$  , and with

$$\begin{array}{ccc}{V}^{\epsilon }=\left(\begin{array}{c}0\\ \epsilon \Xi \\ 0\end{array}\right).& & \end{array}$$

We have $\parallel V^{\epsilon }{\parallel }_{H^s}\le C\epsilon (1+\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^s}^2)$  , for $s$  large enough. Then the energy estimates are the same as in the first proof, the solution ${\mathbf{u}}^{\epsilon }$  satisfying a control of the form:

$$\begin{array}{ccc}\frac{d}{dt}\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^s}^2\le C\left(1+\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^s}^2+\epsilon \parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^s}^3\right)& & \end{array}$$

and the same conclusion holds true. Then from Corollary  2.4 , $v-\overline{v},\rho -1$  can be retrieved from ${\mathbf{u}}^{\epsilon }$  , and we obtain the expected conclusion.

$\square$  Convergence to Euler, three dimensional case In the 3-d case, equation ( 44 ) should be replaced by

$$\begin{array}{ccc}{\partial }_t(\overline{\omega }+\epsilon {\omega }_1)+v\cdot \nabla (\overline{\omega }+\epsilon {\omega }_1)=-(\overline{\omega }+\epsilon {\omega }_1)\epsilon \beta -\omega \cdot \nabla v.& & \end{array}$$

Note that the vorticity equation is the same as in the Euler-Poisson case. This change would not affect the energy estimates.

Higher order approximation Here we prove that the the Euler-Poisson system and the Euler-Monge-Ampere system are closer as $\epsilon$  goes to 0 than Euler-Poisson and Euler. We fix $s\ge s_0=[d/2]+2$  . For ${\overline{v}}_0$  a $H^{s+2}$  smooth divergence-free vector field, we consider $(\overline{v},p)$  a solution of the Euler incompressible system ( 8 ) such that $\overline{v}\in L^{\infty }\left(\right[0,T],H^{s+2}({\mathbb{T}}^d\left)\right)$  for some $T>0$  . We consider also a sequence $(v_{ep}^{\epsilon },{\rho }_{ep}^{\epsilon })$  of solutions of the $\left(E{P}_{\epsilon }\right)$  system with initial data $(v_{ep,0}^{\epsilon },{\rho }_{ep,0}^{\epsilon })$  (with $\int {\rho }_0^{\epsilon }=1$  ) such that ${\epsilon }^{-1}(v_{ep}^{\epsilon }-\overline{v}),{\epsilon }^{-2}({\rho }_{ep}^{\epsilon }-1)$  is bounded in $L^{\infty }\left(\right[0,T^{\prime }],H^{s+1}\times H^s({\mathbb{T}}^d\left)\right)$  , for any $0<T^{\prime }<T$  , if $\epsilon$  is small enough. Thanks to Theorem  2.1 , and from the regularity assumption made on $\overline{v}$  , such a sequence exists for any sequence of well prepared initial data.

Theorem 3.9 Let $s\in \mathbb{N}$  with $s\ge [d/2]+2$  . Let $\overline{v},v_{ep}^{\epsilon },{\rho }_{ep}^{\epsilon }$  be as above. Let $(v_0^{\epsilon },{\rho }_0^{\epsilon })$  be a sequence of initial data such that $\left({\epsilon }^{-2}(v_0^{\epsilon }-v_{ep,0}^{\epsilon }),{\epsilon }^{-3}({\rho }_0^{\epsilon }-{\rho }_{ep,0}^{\epsilon })\right)$  is bounded in $H^s\times H^{s-1}\left({\mathbb{T}}^d\right)$  .

Then there exists a sequence $(v^{\epsilon },{\rho }^{\epsilon })$  of solutions to $\left(EM{A}_{\epsilon }\right)$  with initial data $(v_0^{\epsilon },{\rho }_0^{\epsilon })$  , belonging to $L^{\infty }\left(\right[0,T_{\epsilon }],H^s\times H^{s-1}({\mathbb{T}}^d\left)\right)$  , with ${liminf}_{\epsilon \to 0}{T}_{\epsilon }\ge T$  . Moreover for any $T^{\prime }<T$  and $\epsilon$  small enough, the sequence $\left({\epsilon }^{-2}(v^{\epsilon }-v_{ep}^{\epsilon }),{\epsilon }^{-3}({\rho }^{\epsilon }-{\rho }_{ep}^{\epsilon })\right)$  is bounded in $L^{\infty }\left(\right[0,T^{\prime }],H^s\times H^{s-1}({\mathbb{T}}^d\left)\right)$  .

Remark. We see here that the difference between solutions of $\left(E{P}_{\epsilon }\right)$  and $\left(EM{A}_{\epsilon }\right)$  is of order ${\epsilon }^3$  for the density and of order ${\epsilon }^2$  for the velocity whereas the difference between solutions of $\left(E{P}_{\epsilon }\right)$  (or $\left(EM{A}_{\epsilon }\right)$  ) and Euler was of order ${\epsilon }^2$  for the density and of order $\epsilon$  for the velocity.

Proof. We give the proof when $d=2$  , the proof would be the same when $d=3$  , just with more terms. We introduce $(v_{ep},{\rho }_{ep}=1+{\epsilon }^2{\rho }_1)$  solution to $\left(E{P}_{\epsilon }\right)$  , and $({\beta }_{ep},{\omega }_{ep})=(\nabla \cdot v_{ep},\nabla \times v_{ep})$  . Then we set

$$\begin{array}{ccc}& & v=v_{ep}+{\epsilon }^2{v}_2,\end{array}$$

$$\begin{array}{ccc}& & \nabla \cdot v={\beta }_{ep}+{\epsilon }^2{\beta }_2,\end{array}$$

$$\begin{array}{ccc}& & \rho ={\rho }_{ep}+{\epsilon }^3{\rho }_2,\end{array}$$

$$\begin{array}{ccc}& & \text{curl}v=\omega ={\omega }_{ep}+{\epsilon }^2{\omega }_2.\end{array}$$

The system $\left(EM{A}_{\epsilon }\right)$  now reads:

$$\begin{array}{ccc}& & {\partial }_t({\omega }_{ep}+{\epsilon }^2{\omega }_2)+v\cdot \nabla ({\omega }_{ep}+{\epsilon }^2{\omega }_2)=-({\omega }_{ep}+{\epsilon }^2{\omega }_2)({\beta }_{ep}+{\epsilon }^2{\beta }_2),\end{array}$$

(47)

$$\begin{array}{ccc}& & {\partial }_t({\beta }_{ep}+{\epsilon }^2{\beta }_2)+v\cdot \nabla ({\beta }_{ep}+{\epsilon }^2{\beta }_2)+\nabla (v_{ep}+{\epsilon }^2{v}_2):\nabla (v_{ep}+{\epsilon }^2{v}_2)=\frac{\Delta \psi -d}{{\epsilon }^2},\end{array}$$

(48)

$$\begin{array}{ccc}& & {\partial }_t({\rho }_{ep}+{\epsilon }^3{\rho }_2)+v\cdot \nabla ({\rho }_{ep}+{\epsilon }^3{\rho }_2)=-({\rho }_{ep}+{\epsilon }^3{\rho }_2)({\beta }_{ep}+{\epsilon }^2{\beta }_2).\end{array}$$

(49)

We still define $\Xi$  by

$$\begin{array}{ccc}\Delta \psi -d=\rho -1+{\epsilon }^4\Xi ,& & \end{array}$$

and from Theorem  3.4 we will have

$$\begin{array}{ccc}\parallel \Xi {\parallel }_{H^s\left({\mathbb{T}}^d\right)}\le C{\epsilon }^{-4}\parallel \rho -1{\parallel }_{H^s\left({\mathbb{T}}^d\right)}^2\le C(\parallel {\rho }_1{\parallel }_{H^s\left({\mathbb{T}}^d\right)}^2+{\epsilon }^2\parallel {\rho }_2{\parallel }_{H^s\left({\mathbb{T}}^d\right)}^2),& & \end{array}$$

(we use the notation ${\rho }_{ep}=1+{\epsilon }^2{\rho }_1$  ).

Setting

$$\begin{array}{ccc}{\mathbf{u}}^{\epsilon }=\left(\begin{array}{c}{\omega }_2\\ {\beta }_2\\ {\rho }_2\end{array}\right),& & \end{array}$$

we obtain that

$$\begin{array}{ccc}{\partial }_t{\mathbf{u}}^{\epsilon }+v\cdot \nabla {\mathbf{u}}^{\epsilon }+R^{\epsilon }{\mathbf{u}}^{\epsilon }=T_{\epsilon },& & \end{array}$$

with $R^{\epsilon }$  as before and $T_{\epsilon }$  defined by

$$\begin{array}{ccc}{T}_{\epsilon }=\left(\begin{array}{c}-v_2\cdot \nabla w_{ep}-{\beta }_{ep}{\omega }_2-{\beta }_2{\omega }_{ep}-{\epsilon }^2{\omega }_2{\beta }_2\\ -v_2\cdot \nabla {\beta }_{ep}-2\nabla v_{ep}:\nabla v_2-{\epsilon }^2\nabla v_2:\nabla v_2+\Xi \\ -\epsilon v_2\cdot \nabla {\rho }_1-{\beta }_{ep}{\rho }_2-\epsilon {\beta }_2{\rho }_1-{\epsilon }^2{\beta }_2{\rho }_2\end{array}\right).& & \end{array}$$

Using again Proposition  2.6 as in Lemma  2.5 we obtain that, for $s>d/2$  ,

$$\begin{array}{ccc}\parallel T_{\epsilon }{\parallel }_{H^s\left({\mathbb{T}}^d\right)}\le C_s(1+\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^s\left({\mathbb{T}}^d\right)}+\epsilon \parallel {\mathbf{u}}^{\epsilon }{\parallel }_{H^s\left({\mathbb{T}}^d\right)}^2),& & \end{array}$$

where the constant $C_s$  is controlled by $\parallel v_{ep}{\parallel }_{H^{s+2}},\parallel {\rho }_1{\parallel }_{H^{s+1}}$  (still with ${\rho }_1={\epsilon }^{-2}({\rho }_{ep}-1)$  ). From Theorem  2.1 these quantities are controlled for $0\le t\le T^{\prime }<T$  , $T$  being the time on which the solution of $\left(E\right)$  is smooth. Hence we have by Gronwall’s lemma a bound on $\parallel {\mathbf{u}}^{\epsilon }{\parallel }_{L^{\infty }\left(\right[0,T^{\prime }],H^s)}$  .

Arguing as in the previous proofs, and using Corollary  2.4 , we obtain that $(v_2,{\rho }_2)$  remains bounded in $L^{\infty }\left(\right[0,T^{\prime \prime }],H^{s+1}\times H^s({\mathbb{T}}^d\left)\right)$  for any $T^{\prime \prime }<T^{\prime }$  and for $\epsilon <{\epsilon }_0$  small enough. It follows that $\left({\epsilon }^{-2}\right(v_{ep}-v_{ema}),{\epsilon }^{-3}({\rho }_{ep}-{\rho }_{ema}\left)\right)$  remains bounded in $L^{\infty }\left(\right[0,T^{\prime \prime }],H^s\times H^{s-1}({\mathbb{T}}^d\left)\right)$  . This achieves the proof of Theorem  3.9 .

$\square$ 

3.4 Non-prepared initial data

In this case, we obtain exactly the same result as for Euler-Poisson, using the same techniques.

We follow closely section 2.3, and we only have to estimate the additional source term that will appear in the equation followed by ${\beta }_1$  , due to the Monge-Ampere coupling. We recall that $det{D}^2\psi =\rho$  , and we will have to estimate the difference

$$\begin{array}{ccc}\frac{1}{{\epsilon }^2}\left[(\Delta \psi -d)-(\rho -1)\right]& & \end{array}$$

in $H^s$  when we know that ${\epsilon }^{-1}(\rho -1)$  is bounded in $H^s$  . Thanks to Theorem  3.4 , we conclude that for $s\ge 2$  this term is controlled by $\parallel {\epsilon }^{-1}(\rho -1){\parallel }_{H^s}^2=\parallel {\rho }_1{\parallel }_{H^s}$  . Hence the energy estimate can be handled similarly just with an additional term, and the conclusion remains true.

Theorem 3.10 The Theorem  2.7 holds also when replacing the $\left(E{P}_{\epsilon }\right)$  system by the $\left(EM{A}_{\epsilon }\right)$  system.

$\square$ 

Gregoire Loeper EPFL, SB-IMA 1015 Lausanne, Switzerland gregoire.loeper@epfl.ch