Large time behavior of the Vlasov-Navier-Stokes system on the torus

TL;DR

This paper investigates the long-term behavior of solutions to the Vlasov-Navier-Stokes system on a three-dimensional torus, proving convergence of the distribution function to a velocity Dirac mass with exponential rate under small initial energy.

Contribution

It establishes the exponential convergence of the distribution function to a Dirac mass in velocity for the Vlasov-Navier-Stokes system on the torus, under small initial modulated energy.

Findings

Distribution function converges to a Dirac mass in velocity.

Convergence occurs at an exponential rate.

Global bounds on moments are achieved through bootstrap analysis.

Abstract

We study the large time behavior of Fujita-Kato type solutions to the Vlasov-Navier-Stokes system set on $\mathbb{T}^3 \times \mathbb{R}^3$. Under the assumption that the initial so-called modulated energy is small enough, we prove that the distribution function converges to a Dirac mass in velocity, with exponential rate. The proof is based on the fine structure of the system and on a bootstrap analysis allowing to get global bounds on moments.

Equations

$${\partial }_{t}f+v\cdot {\nabla }_{x}f+{\text{div}}_v[f(u-v)]$$

$${\partial }_{t}u+(u\cdot \nabla )u-\Delta u+\nabla p$$

$$\text{div}u$$

$${\rho }_f(t,x)$$

$$j_f(t,x)$$

$$\text{E}(t):=\frac{1}{2}{\int }_{{\mathbb{T}}^3}|u(t,x){|}^2\mathrm{d}x+\frac{1}{2}{\int }_{{\mathbb{T}}^3\times {\mathbb{R}}^3}f(t,x,v)|v{|}^2\mathrm{d}v\mathrm{d}x,$$

$$\text{D}(t):={\int }_{{\mathbb{T}}^3\times {\mathbb{R}}^3}f(t,x,v)|u(t,x)-v{|}^2\mathrm{d}v\mathrm{d}x+{\int }_{{\mathbb{T}}^3}|\nabla u(t,x){|}^2\mathrm{d}x.$$

$$\frac{\mathrm{d}}{\mathrm{d}t}\text{E}(t)+\text{D}(t)=0,$$

$$u_0$$

$$0\le f_0$$

$$(x,v)\mapsto f_0(x,v)|v{|}^2$$

$${\int }_{{\mathbb{T}}^3\times {\mathbb{R}}^3}{f}_0\mathrm{d}v\mathrm{d}x$$

$$u$$

$$f$$

$$j_f-{\rho }_{f}u$$

$$\text{E}(t)+{\int }_s^t\text{D}(\sigma )\mathrm{d}\sigma \le \text{E}(s),$$

$$(x,v)⟼(1+|v{|}^q)f_0(x,v)\in {\text{L}}^{\infty }({\mathbb{T}}^3\times {\mathbb{R}}^3),$$

$$N_q(f_0):=\underset{x\in {\mathbb{T}}^3,v\in {\mathbb{R}}^3}{\sup }(1+|v{|}^q)f_0(x,v).$$

$$m_{\alpha }\text{f}(t,x)$$

$$M_{\alpha }\text{f}(t)$$

$$⟨\text{h}⟩:={\int }_{{\mathbb{T}}^3}\text{h}\mathrm{d}x.$$

$${\partial }_{t}g+v\cdot {\nabla }_{x}g-{\text{div}}_v[gv]=0.$$

$$g(t,x,v)=e^{3t}{g}_0(x-(e^t-1)v,e^{t}v).$$

$$\forall \phi \in \mathcal{D}({\mathbb{R}}^3),⟨{\delta }_{\text{U}},\phi ⟩=\phi (\text{U}).$$

$$g(t,x,v){\rightharpoonup }_{t\to +\infty }\left({\int }_{{\mathbb{R}}^3}{g}_0(x-v,v)\mathrm{d}v\right)\otimes {\delta }_0.$$

$${\partial }_{t}g+v\cdot {\nabla }_{x}g+{\text{div}}_v[g(\text{U}-v)]=0,$$

$$g(t,x,v)-\left({\int }_{{\mathbb{R}}^3}{g}_0(x-v-t\text{U},v+\text{U})\mathrm{d}v\right)\otimes {\delta }_{\text{U}}{\rightharpoonup }_{t\to +\infty }0$$

$${\int }_{{\mathbb{T}}^3\times {\mathbb{R}}^3}f(t)\mathrm{d}v\mathrm{d}x=1,⟨u+j_f⟩(t)=⟨u_0+j_{f_0}⟩,$$

$$f(t,x,v)\sim {\rho }_f(t,x)\otimes {\delta }_{\frac{⟨u_0+j_{f_0}⟩}{2}}$$

$$\begin{gathered} \mathcal{E}(t):=\frac{1}{2}{\int }_{{\mathbb{T}}^3\times {\mathbb{R}}^3}f(t,x,v)|v-⟨j_f(t,x)⟩{|}^2\mathrm{d}v\mathrm{d}x \\ +\frac{1}{2}{\int }_{{\mathbb{T}}^3}|u(t,x)-⟨u(t)⟩{|}^2\mathrm{d}x+\frac{1}{4}|⟨j_f(t)⟩-⟨u(t)⟩{|}^2. \end{gathered}$$

$$\frac{\mathrm{d}}{\mathrm{d}t}\mathcal{E}(t)+\text{D}(t)=0,$$

$${\text{W}}_1\left(f(t),{\rho }_f(t)\otimes {\delta }_{\frac{⟨u_0+j_{f_0}⟩}{2}}\right)+{\Vert u(t)-\frac{⟨u_0+j_{f_0}⟩}{2}\Vert }_{{\text{L}}^2({\mathbb{T}}^3)}$$

$$\Vert u_0{\Vert }_{{\dot{\text{H}}}^{1/2}({\mathbb{T}}^3)}<\frac{1}{{\text{C}}_{\star }^2},$$

$$\begin{gathered} \phi \left(N_q(f_0)+M_{\alpha }{f}_0+\text{E}(0)+\Vert u_0{\Vert }_{{\text{H}}^{1/2}({\mathbb{T}}^3)}+1\right)\mathcal{E}(0) \\ <\min \left(1,\frac{1}{{\text{C}}_{\star }^2}-\Vert u_0{\Vert }_{{\dot{\text{H}}}^{1/2}({\mathbb{T}}^3)}^2\right), \end{gathered}$$

$$\begin{gathered} {\Vert u(t)-\frac{⟨u_0+j_{f_0}⟩}{2}\Vert }_{{\text{L}}^2({\mathbb{T}}^3)}+{\text{W}}_1\left(f(t),{\rho }^{\infty }\left(x-t\frac{⟨u_0+j_{f_0}⟩}{2}\right)\otimes {\delta }_{\frac{⟨u_0+j_{f_0}⟩}{2}}\right) \\ \le \sqrt{\mathcal{E}(0)}{\text{C}}_{\lambda }\exp \left(-\lambda t\right), \end{gathered}$$

Full text

1. 1 Introduction
2. 2 Main results
3. 3 Conservation laws, energy dissipation identities and consequences
4. 4 Changes of variables and $\textnormal{L}^{\infty}$ bounds on moments
5. 5 Regularity estimates for solutions of the Vlasov-Navier-Stokes system
6. 6 Estimates on the convection and the Brinkman force
7. 7 Exponential decay of the modulated energy
8. 8 Further description of the asymptotic state
9. 9 Appendix

1 Introduction

We consider the Vlasov-Navier-Stokes system in $\mathbb{T}^{3}\times\mathbb{R}^{3}$:

[TABLE]

where

[TABLE]

This system of nonlinear PDEs describes the transport of particles (described by their density function $f$) within a fluid (described by its velocity $u$ and its pressure $p$). It belongs to the broad family of fluid-kinetic systems, which were introduced in the pioneering works of O’Rourke [23] and Williams [25]. Among all possible couplings (we refer to the introduction of [14] for other examples), the Vlasov-Navier-Stokes has been intensively studied because of both its physical relevance (see [5] for instance) and the mathematical challenges that it offers. The Vlasov-Navier-Stokes system is fully coupled: both unknowns $f$ and $u$ depend on each other. This is due to the Brinkman force (the source term in the fluid equation) and the drag acceleration (the inertial term in the kinetic equation). We refer to [5] for the physical justification of these, and to [11, 2, 3, 18, 19] for the (partial) mathematical derivation of the former. The physical constants are all normalized in (1.1) – (1.3).

The mathematical analysis of the Vlasov-Navier-Stokes system has been for a long time focused on the existence of (weak or strong) solutions on rather academic domains [4, 10, 24] like the flat torus that we consider in this paper, or more realistic ones [16, 6]. Most of the previous results provide global existence of weak solutions in the following sense: a Leray solution for the fluid equation and a renormalized one (in the sense of DiPerna-Lions [12]) for the kinetic equation (for a more precise definition, see Definition 1.3 below). These global weak solutions are all built by an approximation-compactness argument which is based on the kinetic energy dissipation of the system. More regular solutions can also be constructed. In 2D, thanks to the uniqueness result of [17], they coincide with the weak solutions. In 3D, regular solutions are only known to exist locally (see [10] for instance). This issue is mainly due to the Navier-Stokes part of the system.

Very few articles deal with the long time behavior of this system. At the formal level, one expects a monokinetic behavior in velocity for the distribution function (in other words, concentration to a Dirac mass in velocity), due to the damping of the fluid component and the friction term acting on the kinetic phase. This behavior however has never been completely proven for the Vlasov-Navier-Stokes system. The closest attempt is the paper [10] of Choi and Kwon in which a conditional theorem is provided: the monokinetic behavior is shown to occur under a boundedness assumption that has not been established for any non-trivial global solution up to now. We intend to fill this gap by using the functional introduced by Choi and Kwon in [10] and proving that this boundedness property (in fact, a stronger one) indeed holds, for appropriate solutions of the Vlasov-Navier-Stokes system, under the assumption that the initial data are (in a sense to be made precise) sufficiently close to equilibrium.

Concerning the long-time behavior of other fluid-kinetic systems, when a Fokker-Planck dissipation is added in the kinetic equation, the situation is less involved because the equilibria are all Maxwellians, which are non-singular and (at least locally) attract all solutions. This has been investigated for instance in [15, 7]. Without this dissipation term, apart from [10], we can mention the work of Jabin [21] in which the Navier-Stokes is replaced by a stationary Stokes equation (and a different coupling term) and [14] in which a specific geometry is considered for the Vlasov-Navier-Stokes system, allowing for non-singular stationary solutions.

As far as our knowledge goes, the results that we present below constitute the first complete and rigorous proof of asymptotic monokinetic behavior for the Vlasov-Navier-Stokes system.

1.1 Weak solutions of the Vlasov-Navier-Stokes system

Let us start with a short review of the notion of weak solutions for the Vlasov-Navier-Stokes system, which will give us the opportunity to introduce some notations.

Definition 1.1**.**

The kinetic energy of the system (1.1) – (1.3) is given for $t\geq 0$ by

[TABLE]

and the dissipation is defined as

[TABLE]

The kinetic energy and dissipation stem from the seminal papers on the Vlasov-Navier-Stokes sytem [16, 4]. One can check that the identity

[TABLE]

formally holds, which paves the way for a theory of global weak solutions.

Definition 1.2**.**

We shall say that $(f_{0},u_{0})$ is an admissible initial condition if

[TABLE]

Remark 1.1**.**

The last condition does not play any role for what concerns the properties of existence, uniqueness and long time behavior that we are about to discuss. However, this normalization allows to simplify the formulas.

We shall also denote $\textnormal{H}^{1}_{\textnormal{div}}(\mathbb{T}^{3})=\textnormal{H}^{1}(\mathbb{T}^{3})\cap\textnormal{L}_{\textnormal{div}}^{2}(\mathbb{T}^{3})$.

Definition 1.3**.**

Consider an admissible initial data $(u_{0},f_{0})$ in the sense of Definition 1.2. A weak solution of the Vlasov-Navier-Stokes system with initial condition $(u_{0},f_{0})$ is a pair $(u,f)$ with the regularity

[TABLE]

with $u$ being a Leray solution of (1.2) – (1.3) (initiated by $u_{0}$) and $f$ a renormalized solution of (1.1) (initiated by $f_{0}$), and such that the following energy estimate holds for almost all $t\geq s\geq 0$ (including $s=0$),

[TABLE]

where the functionals E and D are the energy and dissipation introduced in Definition 1.1.

The existence of weak solutions $(u,f)$ (in the sense of Definition 1.3) to the Vlasov-Navier-Stokes system has been established in [4] (and even on general domains in [24, 6]).

Definition 1.4**.**

We say that an initial condition satisfies the pointwise decay assumption of order $q>0$ if

[TABLE]

and in that case we denote

[TABLE]

We finally introduce some useful notations for moments in velocity and averages on the torus.

Definition 1.5**.**

For all $\alpha\geq 0$ and any measurable non-negative function $\textnormal{f}:\,\mathbb{T}^{3}\times\mathbb{R}^{3}\to\mathbb{R}_{+}$, we set

[TABLE]

For any measurable non-negative function $\textnormal{h}:\,\mathbb{T}^{3}\to\mathbb{R}^{d}$ (for any $d\in\mathbb{N}\setminus\{0\}$), we denote its average by

[TABLE]

1.2 Heuristics for the long time behavior

In this paper, we focus on the description of the long time behavior of weak solutions to the Vlasov-Navier-Stokes system. To this end, it is enlightning to first have a look at the linear Vlasov equation with friction, around the trivial equilibrium $(0,0)$. This reads

[TABLE]

Endowed with an initial condition $g_{0}$ at $t=0$, this equation admits the explicit solution

[TABLE]

Definition 1.6**.**

For $\textnormal{U}\in\mathbb{R}^{3}$, we denote by $\delta_{\textnormal{U}}$ the Dirac measure in velocity supported at U, defined by

[TABLE]

The long time behavior of the solution to (1.16) is explicit, as we observe from (1.17) that

[TABLE]

More generally, given $\textnormal{U}\in\mathbb{R}^{3}$, for the equation

[TABLE]

the long time behavior of the solution is also explicit and described by

[TABLE]

The mechanism at stake in (1.16) and (1.18) is a competition between transport and friction. Friction always wins in the end, causing concentration to a Dirac mass in velocity. In view of this behavior, we may expect a similar concentration phenomenon in velocity for the full Vlasov-Navier-Stokes system, at least in a regime close to some equilibrium.

It is actually even possible to push the heuristics a little further. Taking for granted that the kinetic phase concentrates in velocity, with the behavior $f(t,x,v)\sim\rho_{f}(t,x)\otimes\delta_{u(t,x)}$ as $t\to+\infty$, we observe in particular that the Brinkman force in the Navier-Stokes equations vanishes as $t\to+\infty$. Since it is well known that the solution $u(t)$ of the Navier-Stokes without forcing tends to homogenize to its average in space $\langle u\rangle(t)$, we may expect that $f(t,x,v)\sim\rho_{f}(t,x)\otimes\delta_{\langle u\rangle(t)}$ as $t\to+\infty$. In particular this entails $\langle j_{f}\rangle(t)\sim\langle\rho_{f}\rangle(t)\langle u\rangle(t)$ as $t\to+\infty$. But then, by the conservation laws

[TABLE]

(see (3.2) and (3.3) in Lemma 3.1), we deduce that $\langle u\rangle(t)\sim\frac{\langle u_{0}+j_{f_{0}}\rangle}{2}$ as $t\to+\infty$. To summarize, it follows from this heuristic argument that one can expect

[TABLE]

as $t\to+\infty$, that corresponds to concentration to the constant velocity $\frac{\langle u_{0}+j_{f_{0}}\rangle}{2}$.

1.3 The modulated energy of Choi and Kwon

In [10], Choi and Kwon introduced a modulated version of the energy $\textnormal{E}(t)$ of Definition 1.1:

Definition 1.7**.**

We define the modulated energy as

[TABLE]

It is proved in [10] that the identity

[TABLE]

formally holds. Controlling the modulated energy is interesting in view of the expected long time monokinetic dynamics for the kinetic phase, because of the following statement.

Lemma 1.1**.**

With the previous notations, we have that for all $t\geq 0$,

[TABLE]

where $\textnormal{W}_{1}$ is the Wasserstein(-1) distance.

The definition and basic properties of the Wasserstein distance $\textnormal{W}_{1}$ are given in the Appendix (see Section 9.1). The proof of the previous lemma is postponed to Section 3.2.

2 Main results

Our main result provides a sharp description of the long time behavior of weak solutions to the Vlasov-Navier-Stokes system.

Theorem 2.1**.**

There exists $\textnormal{C}_{\star}>0$ and a nondecreasing onto function $\varphi:\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}$ such that the following holds. Let $(u_{0},f_{0})$ be an admissible initial condition such that $N_{q}(f_{0})<+\infty$ for some $q>4$, $M_{\alpha}f_{0}<+\infty$ for some $\alpha>3$ and $u_{0}\in\textnormal{H}^{1/2}(\mathbb{T}^{3})$. Then, if

[TABLE]

and if the initial modulated energy $\mathscr{E}(0)$ is small enough, in the sense that

[TABLE]

then for any weak solution $(u,f)$ to the Vlasov-Navier-Stokes system, there exists a profile $\rho^{\infty}\in\textnormal{L}^{\infty}(\mathbb{T}^{3})$ and $\lambda,\textnormal{C}_{\lambda}>0$ such that for all $t\geq 0$,

[TABLE]

where $\textnormal{W}_{1}$ is the Wasserstein distance.

We refer to the solutions that we consider as Fujita-Kato type, as we require small initial $\dot{\textnormal{H}}^{1/2}$ norm for the fluid velocity.

Remark 2.1**.**

The constant $\textnormal{C}_{\star}$ is the universal constant given in Proposition 9.10.

We deduce that when ${\langle u_{0}+j_{f_{0}}\rangle}=0$, the distribution function $f(t)$ weakly converges to a stationary solution, whereas when ${\langle u_{0}+j_{f_{0}}\rangle}\neq 0$, the asymptotic behavior is that of a travelling wave.

Remark 2.2**.**

As already said, existence of weak solutions follows from [4] (note by the way that both the pointwise decay assumption and the higher order Sobolev assumption are not relevant for this part).

Remark 2.3**.**

The fact that the asymptotic state for the distribution function is a Dirac mass in velocity, and thus is singular, virtually forbids the use of standard PDE techniques, such as high order Sobolev energy estimates, to prove this result.

Remark 2.4**.**

This result proves that for the Vlasov-Navier-Stokes system, the large time behavior on the torus is very different from that on a domain with partially dissipative boundary conditions (and under adequate geometric control conditions): in [14], it is indeed proved that in the latter case there exist smooth non-trivial equilibria that are locally stable.

Theorem 2.1 will be a consequence of the following result, bearing on the large time behavior of the modulated energy $\mathscr{E}(t)$.

Theorem 2.2**.**

Under the assumptions of Theorem 2.1, the following holds. There exists $\lambda,\textnormal{C}^{\prime}_{\lambda}>0$ such that for all $t\geq 0$,

[TABLE]

*Furthermore, we have the global bounds *

[TABLE]

and

[TABLE]

*where $\eta$ is a continuous nonnegative function such that $\eta(0)=0$. *

Remark 2.5**.**

The constants $\textnormal{C}_{\lambda},\textnormal{C}^{\prime}_{\lambda}$ appearing in Theorems 2.1 and 2.2 are uniform with respect to the various (semi-)norms of $u_{0}$ and $f_{0}$ that appear in the assumptions.

It is actually even possible to describe the structure of the final density $\rho^{\infty}$.

Proposition 2.3**.**

For $\delta$ small enough, under the assumptions of Theorem 2.1, if furthermore $u\in\mathscr{C}^{0}(\mathbb{R}_{+};\textnormal{W}^{1,\infty}(\mathbb{T}^{3}))$ and

[TABLE]

then there exists a vector field

[TABLE]

belonging to $\mathscr{C}^{0}(\mathbb{R}_{+};\mathscr{C}^{1}(\mathbb{T}^{3}\times\mathbb{R}^{3}))$ and such that the following holds. We have

[TABLE]

with

[TABLE]

and $s\mapsto\textnormal{Y}^{s}_{\infty,x,v}$ satisfies

[TABLE]

Remark 2.6**.**

The assumption (2.6) on $u$ is restrictive in the sense that the integral goes down to the time $t=0$. Indeed, the parabolic regularization and the estimates obtained in Section 6 prove that with the assumptions of Theorem 2.1 alone, $u\in\mathscr{C}^{0}([\varepsilon,+\infty);\textnormal{W}^{1,\infty}(\mathbb{T}^{3}))$ for all $\varepsilon>0$. The assumption (2.6) therefore requires higher regularity for the initial fluid velocity $u_{0}$. It is also possible to avoid this extra regularity assumption, replacing $f_{0}$ by the value of $f$ at time $t=1$, and all integrals starting from $s=0$ by the same starting from $s=1$. The relevant assumption replacing (2.6), namely

[TABLE]

is then obtained as a consequence of Theorem 2.2, see (2.5), when $\mathscr{E}(0)$ is taken small enough.

Remark 2.7**.**

Proposition 2.3 is proved in Section 8. For the sake of clarity the proof focuses on the case $\langle u_{0}+j_{f_{0}}\rangle=0$. The proof of the general case is similar and adds in only a few lines of computations, see Remark 8.1.

There are mainly two stabilization mechanisms at stake in the large time dynamics of solutions to the Vlasov-Navier-Stokes system. The first one is due to friction in the Vlasov equation, that forces the distribution function to concentrate in velocity. The second stabilization mechanism comes from the dissipation in the Navier-Stokes equations. There is a competition in the Navier-Stokes equations between this dissipation and the possible growth of the non-linearity and the Brinkman force $F=j_{f}-\rho_{f}u$. Loosely speaking, the smallness assumptions we make allow to tame the influence of the forcing.

As already briefly discussed in the introduction, thanks to the fine structure of the system, there happens to be a modulated energy/dissipation identity that follows from the energy identity and the conservation laws of the system, as exhibited by Choi and Kwon [10]. This identity somehow reflects the two stabilization mechanisms we have just discussed.

2.1 The case of dimension $2$

For the sake of conciseness and physical relevance, we focus in this paper on the case of dimension $3$. However, with the same method that we develop, it is possible to study the Vlasov-Navier-Stokes system on $\mathbb{T}^{2}\times\mathbb{R}^{2}$ with weaker regularity assumptions on the initial data. Namely, we can treat admissible data (in all the Definitions, statements or equations discussed in this section, one has to replace $\mathbb{T}^{3},\mathbb{R}^{3}$ by $\mathbb{T}^{2},\mathbb{R}^{2}$ when necessary), without requiring the higher $\dot{\textnormal{H}}^{1/2}$ regularity for $u_{0}$ like in dimension $3$, see Theorem 2.1 (the fact that more stringent regularity assumptions are required in dimension $3$ is due to the well-known difficulties related to the resolution of the Navier-Stokes equations). For the record, we gather in the following statement what we may obtain in dimension $2$.

Theorem 2.4**.**

There exists a nondecreasing onto function $\varphi:\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}$ such that the following holds. Let $(f_{0},u_{0})$ be an admissible initial condition such that $N_{q}(f_{0})<+\infty$ for some $q>4$. Assume that the initial modulated energy $\mathscr{E}(0)$ is small enough, in the sense that

[TABLE]

Then the weak solution $(u,f)$ to the Vlasov-Navier-Stokes system satisfies (2.2).

Note that the statement is also strengthened compared to dimension $3$ since there is uniqueness of the weak solution of the Vlasov-Navier-Stokes system: it has been indeed established in [17] that in dimension $d=2$, under the pointwise decay assumption of order $q>4$ of Definition 1.4 (and in fact an even less stringent condition is sufficient), uniqueness holds for weak solutions of the Vlasov-Navier-Stokes system.

The proof developed in dimension $3$ applies mutatis mutandis, with the following significant simplifications:

• •

the $\dot{\textnormal{H}}^{1/2}$ regularity for the fluid velocity is not required in order to get higher order energy estimates for positive times, and therefore in particular we do not need to propagate $\dot{\textnormal{H}}^{1/2}$ estimates for all times;

• •

we can rely on various estimates already proved in [17];

• •

several indices in the Sobolev embeddings are more favorable in dimension $2$.

Let us finally mention that Proposition 2.3 holds as well in dimension $2$.

2.2 Outline of the proof and organisation of the paper

To conclude this section, let us provide a (non-technical) outline of the proof of Theorem 2.2. This also gives the opportunity to describe how this paper is organized.

The purpose of the first Section 3 is to explain how Theorem 2.1 can be deduced from Theorem 2.2, and more strikingly how the proof of the latter boils down to one single uniform estimate on the local density of the kinetic phase. In Section 3.1 we gather conservation laws for the Vlasov-Navier-Stokes system. Section 3.2 emphasizes the role of the modulated energy : we prove therein Lemma 1.1 which explains how the decay of this functional leads to concentration in velocity for the particles. Sections 3.3 and Sections 3.4 detail the following key observation of [10] : up to a control of the $\textnormal{L}^{\infty}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$ norm of the local density $\rho_{f}=\int_{\mathbb{R}^{3}}f\,\mathrm{d}v$, the modulated energy is essentially controlled by its dissipation, yielding exponential decay. We also explain, following an argument of Jabin [21], how one can recover the existence of the asymptotic profile $\rho^{\infty}$ appearing in (2.2), once the exponential decay is established.

As a consequence of Section 3, the proof of Theorem 2.2 (and therefore Theorem 2.1) relies only on obtaining the following global bound for $\rho_{f}$ :

[TABLE]

and in fact our bootstrap strategy actually will prove the same estimate for $j_{f}=\int_{\mathbb{R}^{3}}fv\,\mathrm{d}v$.

In Section 4, we present the main tools we used to obtain such bounds on moments. They are based on the method of characteristics, which allows, considering the characteristics curves $(\textnormal{X},\textnormal{V})$ solving the system

[TABLE]

with $(\textnormal{X}(t;t,x,v),\textnormal{V}(t;t,x,v))=(x,v)$, to write solutions to the Vlasov equation as

[TABLE]

We deduce that

[TABLE]

In order to study (2.15), we rely on a change of variables in velocity, referred to as the straightening change of variables, namely $v\mapsto\textnormal{V}(0;t,x,v)$. It is not obvious that this map is a diffeomorphism. In, Section 4, we provide a sufficient condition to ensure this: there exists a constant $\delta>0$ such that, if

[TABLE]

then indeed the straightening change of variable is admissible. Under this smallness condition, the outcome is the estimate

[TABLE]

Similar bounds for $j_{f}=\int_{\mathbb{R}^{3}}vf\,\mathrm{d}v$ can be obtained as well.

This change of variables is inspired by that used by Bardos and Degond [1] for the study of global small solutions to the Vlasov-Poisson system on $\mathbb{R}^{3}\times\mathbb{R}^{3}$.

As a consequence of Section 4, the remaining task is now to prove that for small enough initial modulated energy, the estimate (2.16) holds for all $t\geq 1$ (small times are handled by local estimates). Sections 5 to 7 are dedicated to this task. To this end, we set up a bootstrap argument. Loosely speaking, we consider

[TABLE]

and the aim is to show that $t^{\star}=+\infty$. The general strategy is as follows. Assuming $t^{\star}<+\infty$, we work on the interval of time $[1,t^{\star}]$. We shall obtain regularity estimates for $u$ using higher order energy estimates for the Navier-Stokes equations and maximal parabolic estimates for the Stokes equations. Such bounds are not relevant in terms of decay in time but on $[1,t^{\star}]$ we have thanks to the straightening change of variables

[TABLE]

Therefore, by Choi-Kwon’s key observation, $\mathscr{E}(t)$ decays exponentially fast on $[1,t^{\star})$. The idea is then to interpolate the higher regularity estimate with the pointwise $\textnormal{L}^{2}(\mathbb{T}^{3})$ bound bearing on $u-\langle u\rangle$ which is provided by the exponential decay of the modulated energy. More precisely, we use the Gagliardo-Nirenberg-Sobolev interpolation inequalities to obtain

[TABLE]

for $\alpha\in(0,1)$; we argue similarly for the control of $\nabla u$.

To apply the previous bootstrap strategy, we need enough regularity and integrability on the solutions of the Vlasov-Navier-Stokes system. We prove in Section 5 that any weak solution of the system instantaneously satisfies adequate estimates, which includes:

• •

a short time control of $\rho_{f}$ and $j_{f}$ in $\textnormal{L}^{\infty}(\mathbb{T}^{3})$, using local bounds;

• •

$\textnormal{L}^{\infty}_{t}\textnormal{H}^{1}_{x}\cap\textnormal{L}^{2}_{t}\textnormal{H}^{2}_{x}$ estimates for $u$, on time intervals away from zero, that is to say for $t\geq 1$

[TABLE]

We introduce the convenient notion of strong existence times in order to be able to propagate regularity.

In Section 6, we start to implement the interpolation strategy, relying this time on higher order maximal parabolic estimates for the Stokes equation. The outcome is a control of $\textnormal{D}^{2}u$ in $\textnormal{L}_{\textnormal{loc}}^{p}(\mathbb{R}_{+};\textnormal{L}^{q}(\mathbb{T}^{3}))$ by $(u\cdot\nabla)u$ and $j_{f}-\rho_{f}u$ in the same space.

Then Section 7 is dedicated to the proof of the global bound (2.12) : we explain therein how the previous control of $\textnormal{D}^{2}u$ can be iteratively used to produce an estimate of the form

[TABLE]

Consequently, if $\mathscr{E}(0)$ is small enough, then we must have $t^{\star}=+\infty$, which concludes the proof of Theorem 2.2.

Finally Section 8 is devoted to the proof of Proposition 2.3 which provides a sharper description of the asymptotic behavior. The analysis comes down to the study of the limit as $t\to\infty$ of characteristics (more precisely of renormalized versions of them). For the sake of clarity, the proof is written in the particular case $\langle u_{0}+j_{f_{0}}\rangle=0$ to lighten the computations (see Remark 8.1).

To conclude the paper, Section 9 is an Appendix where we provide some reminders (in particular, we shortly review some well-known basic facts about the Wasserstein distance) and justify $\textnormal{H}^{1}$ energy estimates for the Navier-Stokes equations with source.

3 Conservation laws, energy dissipation identities and consequences

3.1 Conservation laws

We discuss here some conservations laws for the Vlasov-Navier-Stokes system. We start by describing some basic ones in a first lemma: the first two ones come from the structure of the Vlasov equation alone, while the third one is a consequence of the fine structure of the complete system.

Lemma 3.1**.**

Any weak solution (in the sense of Definition (1.3)) satisfies the following conservations laws. For almost all $t\geq 0$,

[TABLE]

Proof.

Considering the results of [4], the only item to prove is (3.3). Let us assume that both $u$ and $f$ are smooth functions. Integrating the Vlasov equation against $v$, the conservation law satisfied by $j_{f}$ reads

[TABLE]

so that $\langle j_{f}\rangle$ satisfies

[TABLE]

On the other hand, from (1.2), $\langle u\rangle$ satisfies

[TABLE]

from which we deduce $\frac{\mathrm{d}}{\mathrm{d}t}\langle u+j_{f}\rangle=0$, and consequently (3.3).

In the general case, for the fluid equation we can directly use $\varphi=1$ as an admissible test function to recover a.e.

[TABLE]

For the kinetic equation we use an approximation argument relying on DiPerna-Lions theory [12] for linear transport equations : we consider a sequence of nonnegative distribution functions $(f_{n})_{n}$ solving the Vlasov equation with regularized vector fields $(u_{n})_{n}$ and regularized and truncated initial conditions $(f_{0,n})_{n}$, and such that for all $n\geq 1$ and all $t\geq 0$,

[TABLE]

By the DiPerna-Lions theory, $f$ is the (strong) limit of $(f_{n})_{n}$ in $\textnormal{L}^{\infty}(0,T;\textnormal{L}^{p}(\mathbb{T}^{3}\times\mathbb{R}^{3}))$ for all finite values of $p$ ; interpolating with the previous bound we infer that $(j_{f_{n}})_{n}\rightarrow j_{f}$ strongly in $\textnormal{L}^{\infty}(0,T;\textnormal{L}^{1}(\mathbb{T}^{3}\times\mathbb{R}^{3}))$ and $(\rho_{f_{n}}u_{n})_{n}\rightarrow\rho_{f}u$ at least in $\textnormal{L}^{1}(0,T;\textnormal{L}^{1}(\mathbb{T}^{3}))$. This is sufficient to pass to the limit in the following identity (which is justified at the regularized level)

[TABLE]

and finally, (3.3) follows for almost every $t$. ∎

A straightforward consequence of (3.3) in Lemma 3.1 is the following formula:

Lemma 3.2**.**

For almost all $t\geq 0$:

[TABLE]

Remark 3.1**.**

We shall use in this paper several times the DiPerna-Lions theory [12], in the same fashion as in the proof of Lemma 3.1. Thanks to the property of strong stability of renormalized solutions, this allows to systematically argue as if both $f$ and $u$ are smooth when looking to establish estimates for the kinetic phase. The argument, as already outlined in the proof of Lemma 3.1, is the following:

• •

consider an approximating sequence $(u_{n})_{n}$ for $u$ and $(f_{n})_{n}$ the associated solution to the Vlasov equation, with a regularized initial condition;

• •

prove the desired estimate for the solution $f_{n}$ (without explicitly using the higher regularity of $f_{n}$ or $u_{n}$);

• •

pass to the limit using the strong stability property of renormalized solutions (and Fatou’s lemma).

In the following, for brevity, we will never write down this argument explicitly but will repeatedly refer to the current remark.

3.2 The role of the modulated energy : proof of Lemma 1.1

Proof.

By the Monge-Kantorovich duality for the $\textnormal{W}_{1}$ distance (see Proposition 9.2 in the Appendix), we have

[TABLE]

We thus infer, using the Cauchy-Schwarz inequality, the normalization (3.2) and the definition of the modulated energy $\mathscr{E}(t)$

[TABLE]

Likewise,

[TABLE]

We therefore deduce, using the normalization (3.2) and the identity (3.5)

[TABLE]

so that by triangular inequality we have established

[TABLE]

On the other hand, using again (3.5), we can also estimate

[TABLE]

and the result follows. ∎

3.3 Dissipation of the modulated energy

As already said in the introduction, Choi and Kwon noticed444As a matter of fact, they consider the more general Vlasov-inhomogeneous Navier-Stokes system but we recover the system (1.1)–(1.3) as soon we stick to the case of constant fluid density. in [10] that the modulated energy (see Definition 1.7) satisfies the following formal identity

[TABLE]

At the level of weak solutions, we are only able to obtain the inequality version of (3.6), as stated in the next lemma.

Lemma 3.3**.**

For any weak solution $(u,f)$ in the sense of Definition 1.3, for almost all $t\geq 0$,

[TABLE]

In particular, we have the following modulated energy/dissipation inequality for almost all $0\leq s\leq t<+\infty$ (including $s=0$),

[TABLE]

Proof.

Let us first write

[TABLE]

that we can simplify in the following way thanks to (3.2)

[TABLE]

so that $\mathscr{E}(t)-\textnormal{E}(t)$ does not depend on $t$ thanks to (3.3). Estimate (3.7) follows then from the energy estimate (1.10). ∎

3.4 Conditional long time behavior

Definition 3.1**.**

Let $c_{P}$ be the Poincaré constant, that is the best constant such that the Poincaré-Wirtinger inequality holds:

[TABLE]

The following result relating the dissipation and the modulated energy is a variant of [10, Theorem 1.2].

Lemma 3.4**.**

There exists a continuous nonincreasing function $\psi:\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}$ such that the following holds, for any weak solution of the VNS system (in the sense of Definition 1.3) for which $\rho_{f}\in\textnormal{L}_{\textnormal{loc}}^{\infty}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$. Fix $T>0$ and define

[TABLE]

Then

[TABLE]

and we have the exponential estimate

[TABLE]

where $\lesssim$ depends only on $\lambda$.

Proof.

First we note that (3.10) $\Rightarrow$ (3.11). Indeed, combining with estimate (3.7) of Lemma 3.3, we get for almost all $0\leq s\leq t\leq T$,

[TABLE]

so that from Lemma 9.3 of the Appendix we get $\mathscr{E}(t)\lesssim\mathscr{E}(0)e^{-\lambda t}$, where $\lesssim$ depends only on $\lambda$. We therefore focus on (3.10) and try to find $\lambda>0$ of the form (3.9).

Define

[TABLE]

The Poincaré-Wirtinger inequality gives us a constant $c_{P}>0$ such that

[TABLE]

Therefore to get (3.10) for some $\lambda>0$, it is sufficient to prove that for some $\gamma,\beta>0$ we have

[TABLE]

with $\beta$ small enough (namely $\beta<c_{P}$) : in that case we have $\textnormal{D}(t)\geq\lambda\mathscr{E}(t)$ with $\lambda:=\min(\gamma,c_{P}-\beta)$.

For the sake of clarity, we omit the time variable for a few lines. We also denote $\|\rho\|_{\infty,T}:=\sup_{[0,T]}\|\rho_{f}(s)\|_{\textnormal{L}^{\infty}(\mathbb{T}^{3})}$. We start with the following identity

[TABLE]

from which we infer

[TABLE]

Now for any $\alpha\in(0,1)$, Young’s inequality entails that

[TABLE]

Combining with (3.13) we have therefore

[TABLE]

On the other hand we have

[TABLE]

from which we deduce

[TABLE]

where we used the normalization property (3.2) and

[TABLE]

In particular, we have

[TABLE]

Since $\alpha\in(0,1)$ we deduce from (3.14)

[TABLE]

which is exactly (3.12) with $\gamma:=1-\alpha$ and $\beta=(\alpha^{-1}-1)\|\rho_{f}\|_{\infty,T}$. Picking $\alpha$ close enough to $1$ (to ensure $\beta<c_{P}$), $\lambda:=\min(\gamma,c_{P}-\beta)$ satisfies (3.10). To check that $\lambda$ can indeed be chosen of the form (3.9) we have to make more explicit the choice of $\alpha$ by imposing for instance the condition $\beta=c_{P}/2$ above, that is $\alpha^{-1}=c_{P}/(2\|\rho\|_{\infty,T})+1$ which is a continuous, nonincreasing, nonvanishing function of $\|\rho\|_{\infty,T}$ : $\alpha$ is then continuous and increasing and $\lambda:=\min(1-\alpha,c_{P}/2)$ is of the form (3.9). ∎

Once exponential decay of the modulated energy is ensured, one can prove the existence of an asymptotic profile $\rho^{\infty}$ for which we have the following convergence statement.

Proposition 3.5**.**

For any weak solution $(u,f)$ to the Vlasov-Navier-Stokes system for which $\sup_{t\geq 0}\|\rho_{f}(t)\|_{\textnormal{L}^{\infty}(\mathbb{T}^{3})}<+\infty$ and $\mathscr{E}(t)\rightarrow_{t\to+\infty}0$ with exponential decay, there exists a profile $\rho^{\infty}\in\textnormal{L}^{\infty}(\mathbb{T}^{3})$ such that

[TABLE]

exponentially fast.

Proof.

We rely on an argument of Jabin [21] used in the context of the large time behavior of the Vlasov-Stokes system. The proof heavily relies on the exponential decay of the modulated energy. Recall the conservation of the mass

[TABLE]

For any smooth function $\psi\in\mathscr{C}^{\infty}(\mathbb{T}^{3})$ we have therefore for $0\leq s\leq t$

[TABLE]

Keeping in mind the definition of the Wasserstein distance (see Section 9.1), one sees that the large time convergence of $\rho_{f}$ (which would imply that the Cauchy criterion is verified for this metric) is in a way or another linked with the decay of $j_{f}(\tau)$ as $\tau\rightarrow+\infty$. In the general case, this property is not expected, as $j_{f}$ is “supposed” to converge to $\rho_{f}\langle u_{0}+j_{f_{0}}\rangle/2$. This justifies to consider the following renormalized density

[TABLE]

for which we have, denoting as well $\overline{j}_{f}:=j_{f}\left(t,x+t\frac{\langle u_{0}+j_{f_{0}}\rangle}{2}\right)$,

[TABLE]

The previous computation implies

[TABLE]

and the integrand is now expected to decay for large time. More precisely if $\|\nabla\psi\|_{\infty}\leq 1$ we have, by translation invariance of the integration over $\mathbb{T}^{3}$

[TABLE]

and we thus deduce by Cauchy-Schwarz inequality

[TABLE]

On the one hand, thanks to Lemma 3.1, the integral of $f$ over $\mathbb{T}^{3}\times\mathbb{R}^{3}$ equals $1$. On the other hand, thanks to Lemma 3.2 we have

[TABLE]

All in all, using the the Definition 1.7 of the modulated energy we have established for any $\psi\in\mathscr{C}^{\infty}(\mathbb{R}^{3})$ such that $\|\nabla\psi\|_{\infty}\leq 1$

[TABLE]

This estimate extends to Lipschitz functions $\psi$ satisfying $\|\nabla\psi\|_{\infty}\leq 1$ by a standard approximation argument and the Monge-Kantorovitch duality formula allows us to write

[TABLE]

The exponential decay of the modulated energy leads to integrability of $\mathscr{E}^{1/2}$ and therefore the Cauchy criterion for $\overline{\rho}_{f}(t)$ is verified for $t\rightarrow+\infty$ : we recover in this way the convergence of $\rho^{\infty}_{f}(t)\rightarrow\rho^{\infty}$ for some measure $\rho^{\infty}$ as $t\rightarrow+\infty$. Since $t\mapsto\rho^{\infty}(t)$ is uniformly bounded in $\textnormal{L}^{\infty}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$, we must have $\rho^{\infty}\in\textnormal{L}^{\infty}(\mathbb{T}^{3})$. Note that the convergence is indeed exponential, thanks to the exponential decay of $\mathscr{E}^{1/2}$ : this can be seen when letting $t\rightarrow+\infty$ in (3.16). Now by a change of variable we have

[TABLE]

which concludes the proof. ∎

4 Changes of variables and $\textnormal{L}^{\infty}$ bounds on moments

In this section we aim at establishing tools for obtaining bounds on the moments $\rho_{f}$ and $j_{f}$. We first obtain rough unconditional integrability results for $\rho_{f}$ and $j_{f}$ thanks to some interpolation estimates. Next, using some adequate change of variables in velocity, we get refined estimates on $\rho_{f}$ and $j_{f}$, which can be controlled along the flow in the following way. Assuming a suitable control on the quantity $\|\nabla u\|_{\textnormal{L}^{1}(0,t;\textnormal{L}^{\infty}(\mathbb{T}^{3}))}$, it is possible to prove that (cf. Lemma 4.5)

[TABLE]

which can be exploited in long time : the core of the bootstrap argument presented in Section 7 is to prove that the control on $\nabla u$ holds as long as $\mathscr{E}(0)$ is small.

Many proofs in this section rely on the representation of the solution to the Vlasov equation using characteristics, which holds at least when $u$ is a smooth vector field.

Definition 4.1**.**

Assume $u$ is smooth (say $\mathscr{C}^{1}$). We define the characteristic curves $\textnormal{X}(s;t,x,v)$ and $\textnormal{V}(s;t,x,v)$ associated with $u$ as the solution to the system of ODEs

[TABLE]

with the initial condition $(\textnormal{X}(t;t,x,v),\textnormal{V}(t;t,x,v))=(x,v)$.

By the method of characteristics, for a smooth vector field $u$, we can write the solution $f$ to the Vlasov equation as

[TABLE]

As explained in Remark 3.1, we then rely on DiPerna-Lions theory to ensure that the estimates we are able to prove with this representation formula still hold even if $u$ is not smooth enough. For instance, a rough bound on the $\textnormal{L}^{\infty}$ norm of $f$ can be directly deduced from (4.2).

Lemma 4.1**.**

For almost all $t\geq 0$,

[TABLE]

In the remaining paragraphs of this section we will systematically use the approximation procedure described in Remark 3.1, without refering to it explicitely. This is in particular the case for each of the proofs which rely on the characteristic curves.

4.1 Rough local bounds on moments

We recall the notations $M_{\alpha}$ and $m_{\alpha}$ introduced in Definition 1.5.

Lemma 4.2**.**

Consider $\alpha\geq 1$ such that $u\in\textnormal{L}_{\textnormal{loc}}^{1}(\mathbb{R}_{+};\textnormal{L}^{\alpha+3}\cap\textnormal{W}^{1,1}(\mathbb{T}^{3}))$ and $M_{\alpha}f_{0}<\infty$. Then $M_{\alpha}f(t)<\infty$ and for all $t>0$ and

[TABLE]

Proof.

Multiplying the Vlasov equation by $|v|^{\alpha}$ and integrating over $\mathbb{T}^{3}\times\mathbb{R}^{3}$, we get

[TABLE]

Recall that, for $0\leq\ell\leq k$, the following interpolation estimate

[TABLE]

holds for any non-negative $g\in L^{\infty}(\mathbb{T}^{3}\times\mathbb{R}^{3})$. In particular for $(\ell,k)=(\alpha-1,\alpha)$ we get

[TABLE]

We can control $\|g\|_{\textnormal{L}^{\infty}(\mathbb{T}^{3})}$ by Lemma 4.1, so that using Hölder’s inequality in (4.5), we infer

[TABLE]

from which we get

[TABLE]

from which (4.4) follows. ∎

Lemma 4.3**.**

Assuming $M_{3}f_{0}<+\infty$, we have the following

• (i)

$M_{3}f\in\textnormal{L}^{\infty}_{\textnormal{loc}}(\mathbb{R}_{+})$;

• (ii)

$\rho_{f}\in\textnormal{L}^{\infty}_{\textnormal{loc}}(\mathbb{R}_{+};\textnormal{L}^{2}(\mathbb{T}^{3}))$* ;*

• (iii)

$j_{f}\in\textnormal{L}^{\infty}_{\textnormal{loc}}(\mathbb{R}_{+};\textnormal{L}^{3/2}(\mathbb{T}^{3}))$.

Proof.

By Lemma 4.2, we have

[TABLE]

But, using the Sobolev embedding $\textnormal{H}^{1}(\mathbb{T}^{3})\hookrightarrow\textnormal{L}^{6}(\mathbb{T}^{3})$ and the Poincaré-Wirtinger inequality and the energy estimate (1.10), we infer

[TABLE]

This concludes the proof of $(i)$. By the interpolation estimate (4.6) for $(\ell,k)=(0,3)$ and $(\ell,k)=(1,3)$ we have

[TABLE]

We therefore obtain $(ii)$ and $(iii)$ thanks to $(i)$ and Lemma 4.1. ∎

4.2 The straightening change of variables

We discuss in this section the change of variables in velocity that will allow us, as explained at the beginning of this section, to prove long time estimates. The idea is to come down to the “free” case (that is to say to the characteristics associated with the vector field $(x,v)\to(v,-v)$ here), by using an appropriate diffeomorphism in velocity. In doing this, a smallness condition bearing on $\|\nabla u\|_{\textnormal{L}^{1}(0,t;\textnormal{L}^{\infty}(\mathbb{T}^{3}))}$ will naturally appear in our calculations.

This change of variables is close in spirit to that employed in [1] by Bardos and Degond in the study of small data solutions to the Vlasov-Poisson system on $\mathbb{R}^{3}\times\mathbb{R}^{3}$. We note however that the stabilization mechanism for Vlasov-Poisson on $\mathbb{R}^{3}\times\mathbb{R}^{3}$ is based on the dispersion properties of the free transport operator, which is significantly different from that used in our work.

We also mention that similar ideas were recently used in the context of the inertialess limit of the Vlasov-Stokes system in [20].

Lemma 4.4**.**

Fix $\delta>0$ such that $\delta e^{\delta}<1/9$. Then, for any $t\in\mathbb{R}_{+}$ satisfying

[TABLE]

and any $x\in\mathbb{R}^{3}$, the map

[TABLE]

is a $\mathscr{C}^{1}$-diffeomorphism from $\mathbb{R}^{3}$ to itself satisfying furthermore

[TABLE]

Proof.

The proof is directly inspired from the arguments outlined in [1, Proposition 1 and Corollary 1].

$(i)$ Consider a generic vector-valued flow $\textnormal{Y}^{s}_{t,z}:=\textnormal{Y}(s;t,z)$ associated with a smooth vector field $w(t,z)$ defined on $\mathbb{R}_{+}\times X$ and assume that $\|\textnormal{D}_{z}w(t)\|_{\textnormal{L}^{\infty}(X)}\leq 1+\psi(t)$, for some function $\psi\in\textnormal{L}_{\textnormal{loc}}^{1}(\mathbb{R}_{+})$. We have $\partial_{s}\textnormal{Y}_{t,z}^{s}=w(s,\textnormal{Y}_{t,z}^{s})$ which after differentiation with respect to $z$ (introducing $\Theta_{t,z}^{s}:=\textnormal{D}_{z}\textnormal{Y}_{t,z}^{s}$) leads to

[TABLE]

from which we get by Gronwall’s inequality for $s\leq t$

[TABLE]

where we used $\Theta_{t,z}^{t}=\textnormal{Id}$.

Now, let us get back to our system. Introducing the state variable $z:=(x,v)$ which belongs to $X=\mathbb{T}^{3}\times\mathbb{R}^{3}$, the vector field $w(t,z):=(v,u(t,x)-v)$ satisfies the assumption for the above abstract result, since $\|\textnormal{D}_{z}w(t)\|_{\textnormal{L}^{\infty}(X)}\leq 1+\|\nabla u\|_{\textnormal{L}^{\infty}(\mathbb{T}^{3})}$. If we denote by $(\textnormal{X}(s;t,z),\textnormal{V}(s;t,z))$ the characteristics associated with $u$, integrating the equation defining $s\mapsto\textnormal{V}(s;t,z)$ we have

[TABLE]

which leads to

[TABLE]

We thus infer from (4.9) with $\psi=\|\nabla u\|_{\textnormal{L}^{\infty}(\mathbb{T}^{3})}$ that

[TABLE]

and thus that

[TABLE]

In particular, if (4.7) holds with $\delta>0$ such that $\delta e^{\delta}\leq\frac{1}{9}$, then Lemma 9.4 applies and we can conclude. ∎

Thanks to the change of variables of Lemma 4.4, we deduce the following control on moments.

Lemma 4.5**.**

If assumption (4.7) of Lemma 4.4 is satisfied, we have for almost all $t\geq 0$,

[TABLE]

where $N_{q}(f_{0})$ is given by (1.12) and

[TABLE]

Proof.

Let $(\textnormal{X}(s,t;x,v),\textnormal{V}(s,t;x,v))$ be the characteristics (4.1) associated with $u$. We start again from the representation formula

[TABLE]

By Lemma 4.4, the mapping $v\mapsto\Gamma_{t,x}(v)=\textnormal{V}(0;t,x,v)$ defines an admissible change of variable of which we deduce

[TABLE]

which implies (the control of the jacobian is given by Lemma 4.4)

[TABLE]

For $j_{f}$ we proceed similarly and write the representation formula (valid for the same reasons)

[TABLE]

By definition of $\Gamma_{t,x}(w)$, we have the identity

[TABLE]

from which we deduce

[TABLE]

hence the claimed result. ∎

In the next lemma, we study how the pointwise decay condition of Definition 1.4 can be locally propagated.

Lemma 4.6**.**

Let $t_{0}>0$. If $f_{0}$ satisfies (1.11) and $u\in\textnormal{L}_{\textnormal{loc}}^{1}(\mathbb{R}_{+};\textnormal{H}^{1}\cap\textnormal{L}^{\infty}(\mathbb{T}^{3}))$, then $f_{t_{0}}:=f(t_{0})$ satisfies also (1.11) and

[TABLE]

Proof.

We write

[TABLE]

Thanks to the differential equation satisfied by $s\mapsto\textnormal{V}(s;t,x,v)$ we have

[TABLE]

We deduce

[TABLE]

and therefore

[TABLE]

∎

This allows to obtain another version of Lemma 4.5 with a control like (4.7) starting only from some time $t_{0}>0$.

Lemma 4.7**.**

Let $t_{0}>0$. With the same assumptions and notations as in Lemma 4.4, except that we replace (4.7) by

[TABLE]

we have for all $t\geq t_{0}$

[TABLE]

Proof.

We can reproduce Lemma 4.4 and Lemma 4.5 replacing the initial time $t=0$ by $t=t_{0}$ and thus $f_{0}$ by $f(t_{0})$. Using Lemma 4.6, we obtain the claimed estimates. ∎

5 Regularity estimates for solutions of the Vlasov-Navier-Stokes system

This section is devoted to the following two tasks:

• •

obtaining a precise short time control for the $\textnormal{L}^{\infty}$ norm of $\rho_{f}$ and $j_{f}$ (relying on local estimates and Lemma 4.6);

• •

obtaining $\textnormal{L}^{\infty}_{t}\textnormal{H}^{1}_{x}\cap\textnormal{L}^{2}_{t}\textnormal{H}^{2}_{x}$ estimates for $u$, on time intervals away from zero, as developed in Proposition 9.10.

Such estimates will be crucial to prove Theorem 2.2, combined with the higher order estimates proved in Section 6.

We shall also introduce in this section the notion of strong existence times (see Definition 5.2). Loosely speaking, this corresponds to times $t$ for which the solution $u$ of the Navier-Stokes equation is strong on the interval of time $[0,t]$, which means in this context that it enjoys $\textnormal{H}^{1/2}(\mathbb{T}^{3})$ regularity. A smallness criterion bearing both on $u$ and on the Brinkman force $j_{f}-\rho_{f}u$ (see (5.4)) will be used.

Notation 5.1**.**

From now on, $A\lesssim_{0}B$ will mean

[TABLE]

where $\varphi:\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}$ is onto, continuous and nondecreasing, and $q>4$ and $\alpha>3$ are the exponents given in the statements of Theorem 2.1 and Theorem 2.2. Note that $\lesssim_{0}$ may depend on the integration exponents appearing in the inequality, but this will always be harmless.

Notation 5.2**.**

We will use the following notations:

[TABLE]

5.1 Local estimates

In this paragraph we establish local estimates on both the fluid and the particle densities. Namely, we prove $u\in\textnormal{L}^{1}_{\textnormal{loc}}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$ and deduce from this estimate that $\rho_{f},j_{f}\in\textnormal{L}^{\infty}_{\textnormal{loc}}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$ and then $F\in\textnormal{L}_{\textnormal{loc}}^{2}(\mathbb{R}_{+};\textnormal{L}^{2}(\mathbb{T}^{3}))$.

Proposition 5.1**.**

We have $u\in\textnormal{L}_{\textnormal{loc}}^{1}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$ and $\rho_{f},j_{f}\in\textnormal{L}_{\textnormal{loc}}^{\infty}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$. Moreover there exists a continuous nondecreasing function $\eta:\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}$ such that

[TABLE]

Proof.

In the proof we denote by $\eta$ a generic continuous function (as in the statement of the proposition), which may vary from line to line.

Since $M_{2}f_{0}<+\infty$ (see the Definition 1.2 of admissible initial data), we have also $M_{3}f_{0}\lesssim M_{2}f_{0}+M_{\alpha}f_{0}<+\infty$. We infer from the proof of Lemma 4.3 that

[TABLE]

In particular, recalling the notation $S=j_{f}-\rho_{f}u-(u\cdot\nabla)u$, we infer, using Hölder’s inequality and the Sobolev embedding $\textnormal{H}^{1}(\mathbb{T}^{3})\hookrightarrow\textnormal{L}^{6}(\mathbb{T}^{3})$ and the energy estimate (1.10),

[TABLE]

Now, if $\mathbb{P}$ stands for the Leray projector (that is the projection on divergence free vector fields), let $w$ be the unique solution of

[TABLE]

so that $u-w=e^{t\Delta}u_{0}$. Since $u_{0}\in\textnormal{H}^{\frac{1}{2}}(\mathbb{T}^{3})$, we infer from [13, Lemma 3.3] that $u-w\in\textnormal{L}^{2}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$ with the estimate

[TABLE]

Thanks to the $\textnormal{L}_{\textnormal{loc}}^{2}(\mathbb{R}_{+};\textnormal{L}^{3/2}(\mathbb{T}^{3}))$ estimate on $S$ that we obtained above, we infer from the continuity of $\mathbb{P}$ on $\textnormal{L}^{3/2}(\mathbb{T}^{3})$ and the maximal regularity of the heat operator on the torus (see Corollary 9.8) that $\Delta w\in\textnormal{L}_{\textnormal{loc}}^{2}(\mathbb{R}_{+};\textnormal{L}^{3/2}(\mathbb{T}^{3}))$. Therefore, from a standard elliptic estimate, we deduce $\textnormal{D}^{2}w\in\textnormal{L}_{\textnormal{loc}}^{2}(\mathbb{R}_{+};\textnormal{L}^{3/2}(\mathbb{T}^{3}))$ and thus $w\in\textnormal{L}_{\textnormal{loc}}^{2}(\mathbb{R}_{+};\textnormal{L}^{p}(\mathbb{T}^{3}))$ for all $p<\infty$, by Sobolev’s embedding. We have even more precisely (keeping track of the different constants)

[TABLE]

Up to now we have thus established (for any $p<\infty$) that $u\in\textnormal{L}_{\textnormal{loc}}^{2}(\mathbb{R}_{+};\textnormal{L}^{p}(\mathbb{T}^{3}))$ with

[TABLE]

In particular, we get $u\in\textnormal{L}_{\textnormal{loc}}^{1}(\mathbb{R}_{+};\textnormal{L}^{\alpha+3}(\mathbb{T}^{3}))$. Using estimate (4.4) of Lemma 4.2 we first have

[TABLE]

We use the interpolation estimate (4.6) with $k=\alpha$ and $\ell\in\{0,1\}$ to obtain this time

[TABLE]

where the integration exponents are strictly larger than $3/2$. Using (5.3) we can estimate $(u\cdot\nabla)u$ in some $\textnormal{L}_{\textnormal{loc}}^{\gamma}(\mathbb{R}_{+};\textnormal{L}^{r}(\mathbb{T}^{3}))$ for $\gamma>1$ and $r>3/2$ leading to the following estimate on the source $S$ :

[TABLE]

Since $r>3/2$, using like before the maximal regularity of the heat operator we eventually infer by the Sobolev embedding $\textnormal{W}^{2,r}(\mathbb{T}^{3})\hookrightarrow\textnormal{L}^{\infty}(\mathbb{T}^{3})$

[TABLE]

All in all, we have obtained that $u=(u-w)+w\in\textnormal{L}_{\textnormal{loc}}^{1}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$. Finally using Lemma 4.6 and the straightforward bound

[TABLE]

we infer that both $\rho_{f}$, $j_{f}$ belong to $\textnormal{L}_{\textnormal{loc}}^{\infty}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$ with the estimate

[TABLE]

∎

Lemma 5.2**.**

Recalling Notation 5.2, we have $F\in\textnormal{L}_{\textnormal{loc}}^{2}(\mathbb{R}_{+};\textnormal{L}^{2}(\mathbb{T}^{3}))$ and moreover

[TABLE]

Proof.

By Cauchy-Schwarz’s inequality, we have a.e.,

[TABLE]

from which we infer for almost all $s\geq 0$,

[TABLE]

where D is the dissipation introduced in (1.5). The estimate follows thus from the energy (1.10) and modulated energy (3.7) estimates. ∎

5.2 Parabolic regularization for the fluid

We state here a consequence of the parabolic regularization result of Proposition 9.10 of the appendix. This roughly establishes the instantaneous gain of two derivatives for the Navier-Stokes equation, if the right-hand side is square-integrable. However, such an estimate can only be obtained if a suitable smallness condition is satisfied.

Proposition 5.3**.**

Assume that for some $T>0$ there holds

[TABLE]

where $\textnormal{C}_{\star}$ is the universal constant given by Proposition 9.10. Then one has for all $1/2\leq t\leq T$ the estimate

[TABLE]

where $\lesssim$ depends only on $\textnormal{C}_{\star}$.

Proof.

If (5.4) is indeed satisfied, we can directly use the well-posedness framework given by Proposition 9.10. Thanks to Lemma 5.2 we have also (9.5) which here reduces to (5.5) because the decay of the energy (1.10) ensures $\textnormal{A}(t)\leq\textnormal{E}(0)$. ∎

5.3 Strong existence times

Thanks to Proposition 5.1, we know that $\rho_{f}$ and $j_{f}$ both belong to $\textnormal{L}_{\textnormal{loc}}^{\infty}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$. We can therefore focus on the boundedness over $[1,+\infty)$. For this purpose, the following notations will be convenient.

Definition 5.1**.**

We set for $t\geq 1$

[TABLE]

In order to use the regularization offered by Proposition 5.3, we need to ensure that the smallness condition (5.4) remains satisfied. For this reason, we introduce the following definition.

Definition 5.2** (Strong existence times).**

A real number $T\geq 0$ will be said to be a strong existence time whenever (5.4) holds.

The following lemma asserts that within our set of assumptions, we have a lower bound for strong existence times.

Lemma 5.4**.**

The smallness condition (2.1) of Theorem 2.1 suffices to ensure that $T=1$ is a strong existence time in the sense of Definition 5.2.

Proof.

Using Lemma 5.2 and estimate and (5.2), we straightforwardly have

[TABLE]

and recalling the meaning of $\lesssim_{0}$ (see Notation 5.1), one sees that the smallness condition (2.1) is indeed sufficient. ∎

6 Estimates on the convection and the Brinkman force

Our ultimate bootstrap argument requires high order estimates bearing on $u$, for which, as in the proof of Proposition 5.1, we will see the Navier-Stokes equation as

[TABLE]

where $\mathbb{P}$ is the Leray projector. We use the maximal regularity of the heat operator on the previous identity to get estimates on $\textnormal{D}^{2}u$ in terms of the Brinkman force $F$ and the the convection term $(u\cdot\nabla)u$. In this short section we explain in Proposition 6.1 the maximal regularity argument and give $\textnormal{L}^{p}_{t}\textnormal{L}^{q}_{x}$ estimates for the source terms in Lemma 6.2 and Lemma 6.3. As explained in Corollary 6.4, these estimates are already sufficient to justify the $\textnormal{L}^{1}_{t}\textnormal{W}^{1,\infty}_{x}$ regularity needed to express the condition (4.16) (and a quantitative version with the required smallness will be provided afterwards).

Proposition 6.1**.**

Fix $a,b,r\in(1,\infty)$ and $\lambda>0$. For any $t\geq 1$, and any exponent $1\leq q\leq a,b$ there holds (with a possible infinite right-hand side)

[TABLE]

where $\Phi:\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}$ is nonincreasing.

Proof.

Similarly to what we have done in the proof of Proposition 5.1, we introduce $w_{1}$ and $w_{2}$ as the unique divergence-free solutions on $[1/2,+\infty)$ of

[TABLE]

with initial conditions $w_{1}(1/2)=w_{2}(1/2)=0$ so that, denoting $u_{h}:=u-(w_{1}+w_{2})$, we have $u_{h}(t+1/2)=e^{t\Delta}u(1/2)$. Now, thanks to the maximal regularity of the heat operator (see Corollary 9.8) and the continuity of $\mathbb{P}$ on $\textnormal{L}^{r}(\mathbb{T}^{3})$, we infer for $t\geq 1/2$

[TABLE]

On the other hand, since $u_{h}(t+1/2)=e^{t\Delta}u(1/2)$, where we write

[TABLE]

we have for $t\geq 1/2$

[TABLE]

and in particular for $t\geq 1$ and $\ell\geq 1$

[TABLE]

so that for any $\ell\geq 1$ we obtain

[TABLE]

By the energy estimate (1.10), we have $\|u(1/2)\|_{\textnormal{L}^{2}(\mathbb{T}^{3})}\lesssim_{0}1$, so using (6.4) for $\ell$ large enough, we infer

[TABLE]

Using the decomposition $u=w_{1}+w_{2}+u_{h}$ and combining (6.2), (6.3) and (6.5), we infer by Hölder’s inequality the estimate (6.1). ∎

Lemma 6.2**.**

There exists $a\in(2,4)$ and $r_{a}>2$ such that the following interpolation estimate holds for $t\geq 1$:

[TABLE]

Proof.

The proof boils down to the interpolation inequality

[TABLE]

Indeed, if the later is satisfied, since $t$ is a strong existence time, we have thanks to the regularization estimate (5.5) and the energy estimate (1.10), together with the Sobolev embedding $\textnormal{H}^{1}(\mathbb{T}^{3})\hookrightarrow\textnormal{L}^{6}(\mathbb{T}^{3})$,

[TABLE]

To justify the interpolation above, notice that for any $a>2$, we have by Hölder inequality and interpolation $[(2,6),(\infty,2)]_{\theta}$,

[TABLE]

with the following equality

[TABLE]

We deduce $\theta=2/a$. From the previous identity we also deduce the value of $r_{a}$, because $\frac{1}{r_{a}}=\frac{1}{6}(1+\frac{2}{a})+\frac{1}{2}(1-\frac{2}{a})$. In the limit case $a=2$ we get $r_{a}=3$, so that taking $|a-2|$ small enough we have indeed $r_{a}>2$ and $a\in(2,4)$. ∎

Lemma 6.3**.**

For any finite $b>4$, the following estimate holds for some $r_{b}>3$ and all strong existence times $t\geq 1$ :

[TABLE]

Proof.

Thanks to Lemma 5.2 and (5.2) we have

[TABLE]

By interpolation $[(2,2);(\infty,6)]_{\theta}$, we have

[TABLE]

where $\theta$ and $r_{b}$ are defined by the equality $(\frac{1}{b},\frac{1}{r_{b}})=\theta(\frac{1}{2},\frac{1}{2})+(1-\theta)(0,\frac{1}{6})$ from which we get $\theta=2/b$ and $\frac{1}{r_{b}}=\frac{2}{3b}+\frac{1}{6}$ ; we notice that $b>4$ implies $r_{b}>3$.

By the triangle inequality, we get

[TABLE]

Using the Sobolev embedding $\textnormal{H}^{1}(\mathbb{T}^{3})\hookrightarrow\textnormal{L}^{6}(\mathbb{T}^{3})$ together with (5.5) and the energy estimate (1.10) we have $\|u\|_{\textnormal{L}^{\infty}(1/2,t;\textnormal{L}^{6}(\mathbb{T}^{3}))}\lesssim_{0}\textnormal{M}_{\rho_{f},j_{f}}(t)^{1/2}$ which implies

[TABLE]

Combining the previous estimate with (6.8) in (6.9) we therefore get

[TABLE]

which is exactly (6.7) because $b=2/\theta$. ∎

Corollary 6.4**.**

For any strong existence time $t\geq 1$, one has $\nabla u\in\textnormal{L}^{1}(1,t;\textnormal{L}^{\infty}(\mathbb{T}^{3}))$.

Proof.

Thanks to Proposition 5.1, the right-hand sides of estimates (6.6) and 6.8 are finite. By Lemmas 6.2 and 6.3, we can therefore take $r>3/2$ in (6.1) and thus, by Sobolev’s embedding and Hölder’s inequality, we finally obtain the claimed regularity. ∎

7 Exponential decay of the modulated energy

In this section, we finish the proof of Theorem 2.2 by setting up a bootstrap procedure. Define

[TABLE]

where $\delta$ is given in Lemma 4.4. Thanks to the change of variables of Section 4, we have that $\textnormal{M}_{\rho_{f},j_{f}}(t)\lesssim_{0}1$ on for $t<t^{\star}$ (see Proposition 7.1). The main goal will be to prove that $t^{\star}=+\infty$. In order to do so, we shall combine the higher order estimates of Section 6 with the exponential decay estimates provided by Lemma 3.4.

Proposition 7.1**.**

We have $t^{\star}>1$. Moreover, for any $t<t^{\star}$, one has $\textnormal{M}_{\rho_{f},j_{f}}(t)\lesssim_{0}1$.

Proof.

By a view of the proof of Lemma 5.4 (reducing $\mathscr{E}(0)$ if necessary), we remark that for $\varepsilon>0$ small enough, $t=1+\varepsilon$ is a strong existence time, and Corollary 6.4 ensures that for $t$ close enough to $1$, the inequality $\int_{1}^{t}\|\nabla u(s)\|_{\textnormal{L}^{\infty}(\mathbb{T}^{3})}\,\mathrm{d}s<\delta$ is satisfied, ensuring $t^{\star}>1$.

For $t\in[1,t^{\star})$ we can invoke Lemma 4.7 with $t_{0}=1$ and (5.1), to obtain that $\textnormal{M}_{\rho_{f}}(t)\lesssim_{0}1$ and

[TABLE]

Thanks to Sobolev’s embedding $\textnormal{H}^{2}(\mathbb{T}^{3})\hookrightarrow\textnormal{L}^{\infty}(\mathbb{T}^{3})$ we infer

[TABLE]

and therefore (using Cauchy-Schwarz’s inequality)

[TABLE]

Thanks to (5.5) and the energy estimate (1.10) we eventually infer

[TABLE]

and we have already proved that $\textnormal{M}_{\rho_{f}}(t)\lesssim_{0}1$. We deduce that $\textnormal{M}_{j_{f}}(t)\lesssim_{0}1$ and this concludes the proof. ∎

We now combine Proposition 6.1 with Lemma 3.4.

Lemma 7.2**.**

Assume that $t_{\star}<\infty$. For any $\alpha\in[1/2,1)$, $c\in[1,\infty)$ and any finite $a,b\geq\max(1,c\alpha)$, the following estimate holds (with a possible infinite right-hand side)

[TABLE]

for $p\in[1,\infty]$ and $r\in(1,\infty)$ satisfying

[TABLE]

Proof.

Owing to Lemma 3.4, if $t^{\star}<+\infty$, there is, on $[0,t^{\star}]$, an exponential decay of the modulated energy with decay rate $\lambda^{\star}$. The Gagliardo-Nirenberg-Sobolev estimate of Theorem 9.9 for $(j,m,q)=(1,2,2)$ allows us to write for any $\alpha\in[1/2,1)$ and $s\geq 1$

[TABLE]

for $p,r$ satisfying (7.3). By definition of the modulated energy and using its exponential decay on $[1,t^{\star}]$, we have therefore

[TABLE]

for $\lambda=\lambda^{\star}(1-\alpha)$. We apply Proposition 6.1 to infer that for any exponent $c$ such that $c\alpha\leq a,b$

[TABLE]

where $\Phi$ is nonincreasing. But by Lemma 3.4, $\lambda^{\star}$ itself is a nonincreasing function of $\textnormal{M}_{\rho_{f}}(t^{\star})\lesssim_{0}1$, which yields (7.2).

∎

Lemma 7.3**.**

There exists $\gamma>0$ such that, if $t^{\star}<+\infty$, then the following estimate holds

[TABLE]

Proof.

We start by combining Lemma 7.2 with Lemma 6.2 and Lemma 6.3. Since (by Proposition 7.1) $\textnormal{M}_{\rho_{f},j_{f}}\lesssim_{0}1$ on $[1,t^{\star}]$, these results give us for some $b>4>a>2$ and $r=\min(r_{a},r_{b})>2$, the following estimate

[TABLE]

which holds for any $\alpha\in[1/2,1)$ and $p$ defined by (7.3), provided that that $\alpha c\leq\min(a,b)$.

It is important to note that $p=\infty$ is not yet reachable at this stage, due to the constraint $\alpha\in[1/2,1]$. However, we can first use Lemma 7.2 with $c=a<b$ in (7.5). In that case, going back to (7.3), we see that the limit case $\alpha=1$ leads to the equality

[TABLE]

which, since $r>2$, implies $\frac{1}{p}<\frac{1}{6}$, that is $p>6$. Taking $\alpha\in[1/2,1)$ close enough to $1$, we therefore infer the existence of $p>6$ such that,

[TABLE]

Since $p>6$, we infer from Hölder’s inequality, for some $\widetilde{r}_{a}>3$, that

[TABLE]

The point is that this last inequality can now replace Lemma 6.2 that we used earlied : we can perform the same analysis as before with the advantage that, now $\widetilde{r}_{a}>3$. This yields that $\widetilde{r}:=\min(r_{b},\widetilde{r}_{a})>3$ and hence taking

[TABLE]

we can check that $\widetilde{\alpha}\in[1/2,1)$ and satisfies

[TABLE]

So we invoke Lemma 7.2 another time with $r=\widetilde{r}>3$, $c=1$ and $\widetilde{\alpha}$ as above to infer

[TABLE]

which is an estimate of the form (7.4). ∎

We are finally in position to conclude the proof of Theorem 2.2.

Proof.

Applying Proposition 7.1, the question thus reduces to ensure $t^{\star}=+\infty$. Assuming $t^{\star}<+\infty$, we will reach a contradiction by proving (for a small enough $\mathscr{E}(0)$) the existence of $t>t^{\star}$ which is still a strong existence time and for which the inequality (7.1) is satisfied.

The first task is to exhibit strong existence times larger than $t^{\star}$. Thanks to Proposition 7.1 and Proposition 5.1, recalling the meaning of the symbol $\lesssim_{0}$ (see Notation 5.1), we have the existence of nondecreasing function $\varphi$ such for any $t\in[1,t^{\star}]$,

[TABLE]

Recall that by assumption, we have $\|u_{0}\|_{\textnormal{H}^{1/2}(\mathbb{T}^{3})}^{2}<\frac{1}{\textnormal{C}_{\star}^{2}}$. Using Lemma 5.2, we thus infer that for all strong existence times $t\leq t^{\star}$

[TABLE]

where we used the embedding $\textnormal{L}^{2}(\mathbb{T}^{3})\hookrightarrow\textnormal{H}^{-\frac{1}{2}}(\mathbb{T}^{3})$, with constant $1$. Combining this with (7.6), we get for some nondecreasing function still denoted $\varphi$

[TABLE]

Therefore, choosing $\mathscr{E}(0)$ small enough so that

[TABLE]

we deduce that

[TABLE]

hence proving by continuity the existence of strong existence times larger than $t^{\star}$.

To check that (7.1) is satisfied after $t^{\star}$ we use Lemma 7.3 to infer the existence of an universal onto nondecreasing continuous function $\varphi:\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}$ such that

[TABLE]

and we observe that a smallness condition as (2.1) ensures

[TABLE]

Therefore we can find a strong existence time $t>t^{\star}$ such that

[TABLE]

This is a contradiction with the definition of $t^{\star}$ and finally concludes the proof.

∎

8 Further description of the asymptotic state

Once the exponential decay of the modulated energy is established, Proposition 3.5 leads to the existence of a profile $\rho^{\infty}\in\textnormal{L}^{\infty}(\mathbb{T}^{3})$ which allows to describe the asymptotic behavior of $f$ in the space variable. The content of Proposition 3.5 is quite implicit as the profile is obtained by an abstract argument. It is in fact possible to describe $\rho^{\infty}$ in a finer way (but still, via implicit equations) : this is the purpose of Proposition 2.3 that we aim at proving in this last section.

Before doing so, it is interesting to compare the statement of Proposition 2.3 with the explicit asymptotic behavior of solutions to the linearized equation when $\langle u_{0}+j_{f_{0}}\rangle=0$, that is the Vlasov equation with friction

[TABLE]

for which we recall we have

[TABLE]

with

[TABLE]

From 2.9 , we therefore see that the deviation from the linearized behavior is small, as

[TABLE]

is small in $\textnormal{L}^{\infty}(\mathbb{T}^{3}\times\mathbb{R}^{3})$, as it is controlled by the initial modulated energy $\mathscr{E}(0)$ and $|\det\mathscr{A}\left(\infty,x,v\right)|-1$ is also small in $\textnormal{L}^{\infty}(\mathbb{T}^{3}\times\mathbb{R}^{3})$, as we will see in the upcoming proof.

We will detail the proof of Proposition 2.3 only in the particular case $\langle u_{0}+j_{f_{0}}\rangle=0$ for which the computations are a bit less tedious. The general case is a straightforward generalization (see Remark 8.1).

Proof of Proposition 2.3 in the case $\langle u_{0}+j_{f_{0}}\rangle=0$.

Recall the map $\Gamma_{t,x}:v\mapsto\textnormal{V}(0;t,x,v)$ that we already used in Lemma 4.4 of Section 4.2 : this very lemma ensures that, for $\delta$ small enough ($\delta e^{\delta}<1/9$ is sufficient), if (2.6) is satisfied, $\Gamma_{t,x}$ is a $\mathscr{C}^{1}(\mathbb{R}^{3})$-diffeomorphism. In order to capture the asymptotic profile of $\rho_{f}(t)$ we look at its action on a continuous function $\psi$ :

[TABLE]

Since $\rho_{f}$ does not solve a transport equation we cannot link it to the initial density $\rho_{f}(0)$, however we can write

[TABLE]

where $\textnormal{Y}(0;t,x,v):=\textnormal{X}(0;t,x,\Gamma_{t,x}^{-1}(v))$. Recall that

[TABLE]

hence (with the notation $\textnormal{Y}(\tau;t,x,v):=\textnormal{X}(\tau;t,x,\Gamma_{t,x}^{-1}(v))$)

[TABLE]

from which we infer

[TABLE]

All in all, introducing the variable $z:=(x,v)$ and denoting $\textnormal{Y}^{s}_{t,z}:=\textnormal{Y}(s;t,z)$, we have established

[TABLE]

where

[TABLE]

In order to understand the behavior of $\rho_{f}(t)$ as $t\rightarrow+\infty$ it is therefore natural to follow the curves $t\mapsto\textnormal{Y}^{s}_{t,z}$ as $t\rightarrow+\infty$, and this is the purpose of the following

Lemma 8.1**.**

For $\delta>0$ small enough, the following holds. For all $0\leq s\leq t$ and $z:=(x,v)\in\mathbb{T}^{3}\times\mathbb{R}^{3}$ we have

[TABLE]

Furthermore, the family of maps $(s,z)\mapsto\textnormal{Y}^{s}_{t,z}$ converges in $\mathscr{C}^{0}(\mathbb{R}_{+};\mathscr{C}^{1}(\mathbb{T}^{3}\times\mathbb{R}^{3}))$, as $t\rightarrow+\infty$, to a map $(s,z)\mapsto\textnormal{Y}^{s}_{\infty,z}$ that satisfies

[TABLE]

Proof.

We start by recalling

[TABLE]

from which, together with (8.1), we deduce the following formula for $s\leq t$

[TABLE]

From the previous expression we infer for $s\leq t$,

[TABLE]

In particular, this implies

[TABLE]

which together with the assumption (2.6) implies for $s\leq t$

[TABLE]

that implies (8.4) for $\delta\leq 1/4$. Similarly and returning to (8.6) we have for $s\leq t$

[TABLE]

and we can proceed in the same way to obtain (8.5). To establish the existence of $(s,z)\mapsto\textnormal{Y}^{s}_{\infty,z}$, we shall prove that $(s,z)\mapsto\textnormal{Y}^{s}_{t,z}$ satisfies Cauchy’s criterion as $t\rightarrow+\infty$, with respect to the local uniform metric. Since $\langle u_{0}+j_{f_{0}}\rangle=0$, we have by Lemma 3.2 and by definition of the modulated energy that $t\mapsto\langle u(t)\rangle$ is integrable over $\mathbb{R}_{+}$ (due to its exponential decay). In particular we infer the integrability over $\mathbb{R}_{+}$ of $t\mapsto\|u(t)\|_{\textnormal{L}^{\infty}(\mathbb{T}^{3})}\leq|\langle u(t)\rangle|+\|\nabla u(t)\|_{\textnormal{L}^{\infty}(\mathbb{T}^{3})}$, thanks to the assumption (2.6). In particular, by dominated convergence we infer that

[TABLE]

where the notation $\text{o}(1)$ refers a term going to [math] in $\textnormal{L}_{\textnormal{loc}}^{\infty}(\mathbb{R}_{+}\times\mathbb{T}^{3}\times\mathbb{R}^{3})$ in the limit $t\rightarrow+\infty$.

Remark 8.1**.**

In the general case $\langle u_{0}+j_{f_{0}}\rangle\neq 0$, $t\mapsto\langle u(t)\rangle$ is not integrable, as it converges to $\langle u_{0}+j_{f_{0}}\rangle/2$. One needs to replace $u(\tau,\textnormal{Y}^{\tau}_{s,t})$ by $u(\tau,\textnormal{Y}^{\tau}_{s,t})-\langle u_{0}+j_{f_{0}}\rangle/2$ in the integrand of (8.7) and, by doing so, adds a diverging drift term to the equation. In a similar fashion as the proof of Proposition 3.5, this can be counterbalanced by considering the renormalized characteristics $\textnormal{Y}(\tau;t,x+\langle u_{0}+j_{f_{0}}\rangle/2,v)$ instead of $\textnormal{Y}(\tau,t,x,v)$. The equations for these shifted trajectories are a bit different, but the convergence properties are proved in the same way, resulting in the implicit equation (2.9).

In particular, taking the difference of this identity (8.7) at times $t_{1}<t_{2}$

[TABLE]

where $\text{o}(1)$ refers here to the asymptotic $t_{1}\wedge t_{2}\rightarrow+\infty$, with the same uniformity as before. Using once more the integrability of $t\mapsto\|u(t)\|_{\textnormal{L}^{\infty}(\mathbb{T}^{3})}$, for any compact $K\subset\mathbb{R}_{+}\times\mathbb{T}^{3}\times\mathbb{R}^{3}$, if

[TABLE]

we have established

[TABLE]

with a similar (uniform) asymptotic term $\text{o}(1)$. From assumption (2.6), this proves

[TABLE]

which yields Cauchy’s criterion. We deduce the existence of $(s,z)\mapsto\textnormal{Y}^{s}_{\infty,z}$, as the (local uniform) limit of $(s,z)\mapsto\textnormal{Y}^{s}_{t,z}$ as $t\rightarrow+\infty$. By dominated convergence and continuity of $u$ for positive times, $\textnormal{Y}^{s}_{\infty,z}$ must satisfy the equation

[TABLE]

For now $\textnormal{Y}^{\tau}_{\infty,z}$ is merely continuous (as a uniform limit) in all its variables. But it turns out that the derivatives $(s;t,z)\mapsto\textnormal{D}_{z}\textnormal{Y}^{s}_{t,z}$ enjoys the same Cauchy criterion as $\textnormal{Y}^{s}_{t,z}$. Indeed, going back to (8.6), we infer, using integrability of $\tau\mapsto\|\nabla u(\tau)\|_{\textnormal{L}^{\infty}(\mathbb{T}^{3})}$ over $\mathbb{R}_{+}$ and dominated convergence

[TABLE]

where $\text{o}(1)$ refers to the asymptotic $t\rightarrow+\infty$ and is locally uniform in $s,z$, while $r_{s,z}$ is some irrelevant function which does not depend on $t$. For any $t_{1}<t_{2}$ we thus have

[TABLE]

where $\text{o}(1)$ refers to $t_{1}\wedge t_{2}\rightarrow+\infty$ and is locallly uniform in $s,z$. Owing to the integrability of $\tau\mapsto\|\nabla u(\tau)\|_{\textnormal{L}^{\infty}(\mathbb{T}^{3})}$ over $\mathbb{R}_{+}$ and the uniform bound on $(s,t,z)\mapsto\mathbf{1}_{s\leq t}\textnormal{D}_{z}\textnormal{Y}^{s}_{t,z}$ due to estimates (8.4) – (8.5), we infer

[TABLE]

Since $(\textnormal{Y}^{s}_{t,z})_{t}\rightarrow\textnormal{Y}^{s}_{\infty,z}$ pointwisely, the continuity of $\nabla u$ for positive times, its belonging to $\textnormal{L}^{1}(\mathbb{R}_{+};\textnormal{L}^{\infty}(\mathbb{T}^{3}))$ and the aforementioned uniform boundedness of $(s,t,z)\mapsto\mathbf{1}_{s\leq t}\textnormal{D}_{z}\textnormal{Y}^{s}_{t,z}$ entail, by dominated convergence,

[TABLE]

and we can then proceed as we have done for $\textnormal{Y}^{s}_{t,z}$ to establish the local uniform Cauchy criterion. ∎

If $f_{0}$ was assumed to be continuous in the space variable, we would now able to pass to the limit into formula (8.2) ; indeed we would have then by dominated convergence, using the bounds that we have established on $(s,t,z)\mapsto e^{s}\textnormal{D}_{v}\textnormal{Y}^{s}_{t,z}$ and the integrability of $v\mapsto\sup_{\mathbb{T}^{3}}f_{0}(\cdot,v)$,

[TABLE]

with

[TABLE]

Notice that here the convergence $z\mapsto\mathscr{A}(t,z)$ towards $z\mapsto\mathscr{A}(\infty,z)$ is also locally uniform in $z$. However, we are not in position to replace $f_{0}$ by a regularized version : to do so we would need a uniqueness result for the whole coupling, and such a result is only known in dimension 2 (see [17]). It turns out that the above convergence holds, but to establish it we have to use another change of variable. More precisely, in (8.2) we consider the change of variable $x\mapsto\Lambda_{t,v}(x):=\textnormal{Y}^{0}_{t,x,v}$. It is an admissible one thanks Lemma 9.4 and the estimate

[TABLE]

which itself is a consequence of (8.6), (8.4) and assumption (2.6), if $\delta$ is small enough. We have therefore

[TABLE]

The long-time behavior of $\Lambda_{t,v}^{-1}(x)$ is given by

Lemma 8.2**.**

For all $v\in\mathbb{R}^{3}$ the map $\Lambda_{\infty,v}:x\mapsto\textnormal{Y}^{0}_{\infty,x,v}$ is a $\mathscr{C}^{1}$-diffeomorphism from $\mathbb{T}^{3}$ onto itself and we have $\Lambda_{t,v}^{-1}(x)\rightarrow_{t}\Lambda_{\infty,v}^{-1}(x)$ in $\mathscr{C}^{1}(\mathbb{T}^{3}\times\mathbb{R}^{3})$, as $t\rightarrow+\infty$ and also $|\det\textnormal{D}_{x}\Lambda_{t,v}^{-1}(x)|\leq 2$ for all $x,v,t$.

Proof.

First, we infer from (8.9) the same estimate (by uniform convergence) for $\Lambda_{\infty,v}$, which is therefore also (thanks to Lemma 9.4) a $\mathscr{C}^{1}$-diffeomorphism. The same lemma gives also $\det\Lambda_{t,v}\geq 1/2$ for all $t\in[1,\infty]$. Again thanks to Lemma 9.4, we infer also uniformly in $t,x,v$, $|\textnormal{D}_{x}\Lambda_{\infty,v}^{-1}(x)|\leq 9/8$ and $\det\textnormal{D}_{x}\Lambda_{t,v}(x)\geq 1/2$. For the convergence, we write

[TABLE]

that goes to [math] locally uniformly in $x,v$ thanks to Lemma 8.1. Since the inversion map is $\mathscr{C}^{1}$ on $\text{GL}_{3}(\mathbb{R})$, using the previous lower bound on the determinants, we infer from the equatity $\textnormal{D}_{x}\Lambda_{t,v}^{-1}=(\textnormal{D}_{x}\Lambda_{t,v})^{-1}\circ\Lambda_{t,v}^{-1}$ and the previous convergence the announced convergence in $\mathscr{C}^{1}(\mathbb{T}^{3}\times\mathbb{R}^{3})$. ∎

Since $\mathscr{A}(t,z)$ is uniformly bounded and continuous and converges (locally uniformly) towards $\mathscr{A}(\infty,z)$, we infer from Lemma 8.2 and the dominated convergence theorem (using $f_{0}\in\textnormal{L}^{1}(\mathbb{T}^{3}\times\mathbb{R}^{3})$),

[TABLE]

and using back the change of variable $x\mapsfrom\Lambda_{\infty,v}(x)$ (which is admissible thanks to Lemma 8.2) we have eventually proved

[TABLE]

where

[TABLE]

which concludes the proof. ∎

9 Appendix

9.1 Wasserstein distance

To simplify the presentation, $X$ here will denote either $\mathbb{T}^{3}$ or $\mathbb{T}^{3}\times\mathbb{R}^{3}$.

Definition 9.1**.**

For $m>0$ we denote by $\mathcal{M}_{1,m}(X)$ the set of all measures $\mu$ such that

[TABLE]

Definition 9.2**.**

Fix $m>0$ and consider $\mu$ and $\nu$ in $\mathcal{M}_{1,m}(X)$. The Wassertein distance between $\mu$ and $\nu$ is

[TABLE]

where $\Gamma(\mu,\nu)$ denotes the collection of all measures on $X\times X$ with first and second marginal respectively equal to $\mu$ and $\nu$.

Proposition 9.1** ($\textnormal{W}_{1}$ metrizes the weak-$\star$ convergence).**

Fix $m>0$. Given $(\mu_{n})_{n}\in\mathcal{M}_{1,m}(X)^{\mathbb{N}}$ and $\mu\in\mathcal{M}_{1,m}(X)$, the two following facts are equivalent

• (i)

For all $f\in\mathscr{C}^{0}_{b}(X)$,

[TABLE]

• (ii)

$(\textnormal{W}_{1}(\mu_{n},\mu))_{n}\rightarrow_{n}0$.

Proposition 9.2** (Monge-Kantorovitch duality).**

Fix $m>0$ and consider $\mu$ and $\nu$ in $\mathcal{M}_{1,m}(X)$. Then

[TABLE]

9.2 Exponential decay

Lemma 9.3**.**

Consider $u:\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}$ a non-increasing integrable function satisfying for some $\lambda>0$ and almost all $t\geq 0$

[TABLE]

Then for $t\geq 0$ there holds

[TABLE]

Proof.

If $v(t)$ denotes the integral in the estimate, $v\in\textnormal{W}^{1,\infty}(\mathbb{R}_{+})$ satisfies $v^{\prime}\leq-\lambda v$, so the standard version of the Gronwall Lemma implies $v(t)\leq v(0)e^{-\lambda t}$. Since $u\leq u(0)$ w.l.o.g. we can assume $t\geq 1$ and since $u$ is non-increasing, we have

[TABLE]

9.3 Perturbation of the identity map

We use in this work the following version of the inverse function theorem.

Lemma 9.4**.**

For $\Omega=\mathbb{T}^{3}$ or $\Omega=\mathbb{R}^{3}$, if $\phi:\Omega\rightarrow\Omega$ is $\mathscr{C}^{1}$ and satisfies $\|\nabla\phi\|_{\infty}<1$, then $f:=\textnormal{Id}+\phi$ is a $\mathscr{C}^{1}$-diffeomorphism of $\Omega$ onto itself satisfying $\|\nabla f\|_{\infty}\leq(1-\|\nabla\phi\|_{\infty})^{-1}$. If furthermore $\|\nabla\phi\|_{\infty}\leq 1/9$, then $\det\nabla f\geq 1/2$.

9.4 Maximal regularity

The maximal regularity estimate for the heat equation, on the whole space, can be stated in the following way.

Theorem 9.5**.**

For $p,q\in(1,\infty)$, and $\varphi\in\mathcal{S}(\mathbb{R}\times\mathbb{R}^{3})$ such that $\varphi(0,\cdot)=0$, there holds

[TABLE]

This estimate is for instance a consequence of [22]. Naturally, one expects an analogous estimate on the torus, but we did not manage to exhibit a precise reference in the literature. For the sake of completeness we give therefore a proof of the following corollary.

Corollary 9.6**.**

For $p,q\in(1,\infty)$ and $\psi\in\mathcal{S}(\mathbb{R}\times\mathbb{R}^{3})$ which is $\mathbb{Z}^{3}$-periodic in the space variable and such that $\psi(0,\cdot)=0$, there holds

[TABLE]

Proof.

Let’s use the Dirac comb $\Sh:=\sum_{n\in\mathbb{Z}^{3}}\delta_{n}$ as a getway between functions defined on $\mathbb{R}^{3}$ and $\mathbb{Z}^{3}$-periodic functions (identified as functions defined on the torus $\mathbb{T}^{3}$). In the sequel $C$ denotes the open unit cube $(0,1)^{3}$.

Lemma 9.7**.**

For any $g\in\mathcal{S}(\mathbb{R}^{3})$ which is $\mathbb{Z}^{3}$-periodic there exists $h\in\mathscr{D}(C)$ such that $g=\Sh\star h$. Furthermore, for any such function $h$, and for any $\in[1,\infty]$ there holds $\|g\|_{\textnormal{L}^{q}(\mathbb{R}^{3})}=\|h\|_{\textnormal{L}^{q}(\mathbb{T}^{3})}$.

Proof.

Fix a non-zero $\theta\in\mathscr{D}(C)$, then $h:=g\theta/(\Sh\star\theta)$ is a well-defined element of $\mathscr{D}(C)$ satisfying $g=\Sh\star h=\sum_{n\in\mathbb{Z}^{3}}\tau_{n}h$. Since $h\in\mathscr{D}(C)$, the functions $\tau_{n}h$ have disjoint supports which justifies the equality of the $\textnormal{L}^{q}$-norms. ∎

Obviously the previous lemma holds also when adding a time variable. In particular we have the existence of $\varphi\in\mathcal{S}(\mathbb{R}\times\mathbb{R}^{3})$, such that $\varphi(0,\cdot)=0$, $\varphi(t,\cdot)\in\mathscr{D}(C)$ for all $t$ and $\psi=\Sh\star\varphi$ (here the convolution is to be understood in the space variable only). From this representation formula we also deduce $\Delta\psi=\Sh\star\Delta\varphi$ and $\partial_{t}\psi-\Delta\psi=\Sh\star(\partial_{t}\varphi-\Delta\varphi)$, where the spatial support of $\Delta\varphi$ and $\partial_{t}\varphi-\Delta\varphi$ are still included in $C$ : the previous lemma applies therefore to write

[TABLE]

where the inequality is obtained from (9.5). ∎

In the current article we will use the following consequence of Corollary 9.6, which is obtained by a standard approximation argument.

Corollary 9.8**.**

For $p,q\in(1,\infty)$ and $T>0$ if $S\in\textnormal{L}^{p}(0,T;\textnormal{L}^{q}(\mathbb{T}^{3}))$, the unique tempered solution $u$ of

[TABLE]

satisfies

[TABLE]

9.5 Interpolation

The following classical interpolation estimate can be for instance found in [9, Thm 1.5.2].

Theorem 9.9** (Gagliardo-Nirenberg-Sobolev).**

Consider $1\leq p,q,r\leq\infty$ and $m\in\mathbb{N}$. Assume that $j\in\mathbb{N}$ and $\alpha\in\mathbb{R}$ satisfy

[TABLE]

with the exception $\alpha<1$ if $m-j-3/r\in\mathbb{N}$. Then, the following holds. For any $g\in\textnormal{L}^{q}(\mathbb{T}^{3})$, if $\textnormal{D}^{m}g\in\textnormal{L}^{r}(\mathbb{T}^{3})$, then $\textnormal{D}^{j}g\in\textnormal{L}^{p}(\mathbb{T}^{3})$ and we have the following estimate for $g$

[TABLE]

where the constant behind $\lesssim$ does not depend on $g$. If $\langle\textnormal{D}^{j}g\rangle=0$, then the term $\|g\|_{\textnormal{L}^{q}(\mathbb{T}^{3})}$ in the right-hand side can be dispensed with.

9.6 Parabolic regularization for the Navier-Stokes equations with a source term

The main result of this section is Proposition 9.10, which gives higher order energy estimates for the Navier-Stokes system together with a form of regularization along time. These estimates seem to be folklore but we give here the proof for the sake of completeness.

Proposition 9.10**.**

There exists a universal constant $\textnormal{C}_{\star}>0$ such that the following holds. Consider $u_{0}\in\textnormal{H}_{\textnormal{div}}^{1/2}(\mathbb{T}^{3})$, $F\in\textnormal{L}_{\textnormal{loc}}^{2}(\mathbb{R}_{+};\textnormal{H}^{-1/2}(\mathbb{T}^{3}))$ and $T>0$ such that

[TABLE]

Then, there exists on $[0,T]$ a unique Leray solution of the Navier-Stokes system with source $F$ and with initial data $u_{0}$. This solution $u$ belongs to $\textnormal{L}^{\infty}([0,T];\textnormal{H}^{1/2}(\mathbb{T}^{3}))\cap\textnormal{L}^{2}(0,T;\textnormal{H}^{3/2}(\mathbb{T}^{3}))$ and satisfies for a.e. $0\leq t\leq T$

[TABLE]

Furthermore, if $F\in\textnormal{L}_{\textnormal{loc}}^{2}(\mathbb{R}_{+};\textnormal{L}^{2}(\mathbb{T}^{3}))$, we have for a.e. $1/2\leq t\leq T$

[TABLE]

where $\lesssim$ depends only on $\textnormal{C}_{\star}$, and A is defined by

[TABLE]

Proof.

The proof proceeds in two different steps. First, if such a Leray solution of the Navier-Stokes exists, because of the interpolation estimate

[TABLE]

we have in particular $u\in\textnormal{L}^{4}(0,T;\textnormal{H}^{1}(\mathbb{T}^{3}))$. This is a known case of weak-strong uniqueness, see for instance the stability result [8, Theorem 3.3]. The second step is to prove that such a solution indeed exists. This follows by a simple compactness argument, using Proposition 9.11 and Proposition 9.12 below, choosing for $\gamma$ an appropriate regularization of $t\mapsto 2t\mathbf{1}_{0\leq t\leq 1/2}+\mathbf{1}_{t>1/2}$. ∎

In order to prove the existence of a solution as in Proposition 9.10, we rely on the following standard approximation procedure: we consider, for $\chi\in\mathscr{C}^{\infty}(\mathbb{T}^{3})$, the regularized system:

[TABLE]

where $\widetilde{u}_{\chi}:=u\star\chi$. When $u_{0}$ and $F$ are smooth, the existence of a unique smooth solution to system (9.7) – (9.9) is standard.

Proposition 9.11**.**

Consider a nondecreasing function $\gamma\in\mathscr{C}^{1}_{b}(\mathbb{R})$ vanishing at [math] and such that $\|\gamma\|_{\textnormal{W}^{1,\infty}(\mathbb{R})}\leq 1$. There exists $\textnormal{C}>0$ and an onto nondecreasing continuous function $\varphi:\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}$, such that for any $u_{0}\in\mathscr{C}^{\infty}_{\textnormal{div}}(\mathbb{T}^{3})$, $F\in\mathscr{C}^{\infty}(\mathbb{R}_{+}\times\mathbb{T}^{3})$ and any $\chi\in\mathscr{C}^{\infty}(\mathbb{T}^{3})$ such that $\|\chi\|_{1}=1$, the unique solution $u$ of (9.7) – (9.9) satisfies for $t\geq 0$,

[TABLE]

where the constant behind $\lesssim$ is universal, $\Phi(z):=(1+z)e^{z}$, A is given by (9.6) and

[TABLE]

Proof.

We multiply the equation by $-\gamma(t)\Delta u$, and use adequate integrations by parts together with Young’s and Hölder’s inequality, to get

[TABLE]

We use then another time Young’s inequality and the Sobolev embedding $\textnormal{H}^{1}(\mathbb{T}^{3})\hookrightarrow\textnormal{L}^{6}(\mathbb{T}^{3})$ to write

[TABLE]

Using the definition (9.6) of $\textnormal{A}(t)$ and the fact that $\|\gamma\|_{\textnormal{W}^{1,\infty}(\mathbb{R})}\leq 1$, we infer, introducing $\ell(t):=\gamma(t)\|\nabla u(t)\|_{\textnormal{L}^{2}(\mathbb{T}^{3})}^{2}$

[TABLE]

which implies by Gronwall’s inequality (since $\ell(0)=0$), using once again the definition of $\textnormal{A}(t)$,

[TABLE]

where $h$ is given by (9.11), for some univeral constant $\textnormal{C}>0$ ; this last estimate can be recasted into (9.10). ∎

Recall the notation $\|\cdot\|_{\dot{\textnormal{H}}^{s}(\mathbb{T}^{3})}$ for the $\textnormal{L}^{2}$ norm associated with the multiplier $|\xi|^{s}$.

Proposition 9.12**.**

There exists a universal constant $\textnormal{C}_{\star}$ such that the following holds. For any $u_{0}\in\mathscr{C}^{\infty}_{\textnormal{div}}(\mathbb{T}^{3})$, $F\in\mathscr{C}^{\infty}(\mathbb{R}_{+}\times\mathbb{T}^{3})$ and any $\chi\in\mathscr{C}^{\infty}(\mathbb{T}^{3})$ such that $\|\chi\|_{1}=1$, if for some $T>0$ one has

[TABLE]

then the unique solution $u$ of (9.7) – (9.9) satisfies for $t\in[0,T]$,

[TABLE]

In particular, recalling the definition (9.11), on $[0,T]$ we have $h\leq\textnormal{C}/\textnormal{C}_{\star}$ where C is the universal constant given in Proposition 9.11.

Proof.

Let us first recall the fundamental energy estimate

[TABLE]

Consider $\Lambda$ the Fourier multiplier associated with $|\xi|$. After taking the scalar product with $\Lambda u$, thanks to Plancherel’s formula, Hölder’s inequality, to the continuity of the Leray projector $\mathbb{P}$ on $\textnormal{L}^{3/2}(\mathbb{T}^{3})$, we can obtain as well

[TABLE]

Using Young’s inequality and combining with (9.14) we infer

[TABLE]

we therefore have

[TABLE]

We have by Sobolev embedding

[TABLE]

Since $\Lambda u$ and $\nabla u$ have a vanishing mean, we therefore have

[TABLE]

This is a differential inequality of the form

[TABLE]

where C is some universal constant and

[TABLE]

After integration, we hence have

[TABLE]

In particular, if for some $T>0$ one has (this precisely corresponds to the assumption (9.3))

[TABLE]

then by a standard continuity argument we can show that the inequality $x(t)^{1/2}\leq 1/C$ which is true for $t=0$ remains valid up to $t=T$, entailing on $[0,T]$,

[TABLE]

which corresponds to the desired inequality, recalling (9.15). ∎

Acknowledgements. We thank Young-Pil Choi for pointing out a mistake in a previous version of the paper. DHK and AM were partially supported by the grant ANR-19-CE40-0004. IM was partially supported by the European Research Council (ERC) MAFRAN grant under the European’s Union Horizon 2020 research and innovation programme (grant agreement No 726386) while he was a research associate in PDEs at DPMMS, University of Cambridge.