Rigidity results in generalized isothermal fluids
Abstract.
We investigate the longtime behavior of solutions to the isothermal Euler, Korteweg or quantum Navier Stokes equations, as well as generalizations of these equations where the convex pressure law is asymptotically linear near vacuum. By writing the system with a suitable timedependent scaling we prove that the densities of global solutions display universal dispersion rate and asymptotic profile. This result applies to weak solutions defined in an appropriate way. In the exactly isothermal case, we establish the compactness of bounded sets of such weak solutions, by introducing modified entropies adapted to the new unknown functions.
1. Introduction
In the isentropic case , the Euler equation on , ,
(1.1) 
enjoys the formal conservations of mass,
and entropy (or energy),
In general, smooth solutions are defined only locally in time (see [MUK86, JYC90, Xin98]). However, for some range of , if the initial velocity has a special structure and the initial density is sufficiently small, the classical solution is defined globally in time. In addition the large time behavior of the solution can be described rather precisely, as established in [Serre97]. We restate some results from [Serre97] in the following theorem:
Theorem 1.1 (From [Serre97]).
Let and . There exists such
that the following holds.
If are such that
,
then the system (1.1) with initial data and admits a unique global
solution, in the sense that , where . In addition, there exists such that
(1.2) 
Conversely, if are such that , then there exists such that the solution to (1.1) with and is global in time in the same sense as above, and (1.2) holds.
In particular, in the frame of small data (in the sense described above), the dispersion
is universal but the asymptotic profile can be arbitrary. Typically, given any function , will be allowed provided that is sufficiently small. For completeness we provide a brief proof of the above theorem in appendix.
We emphasize that the structure of the velocity is crucial: the initial velocity is a small (decaying) perturbation of a linear velocity. In a way, the above result is the Euler generalization of the global existence results for the Burgers equation with expanding data. Refinements of this result can be found in [GrSe97, Gra98, Serre2017].
The isothermal Euler equation corresponds to the value in (1.1),
(1.3) 
The mass is still formally conserved, and the energy now reads
Unlike in the isentropic case, the energy has an indefinite sign, a property which causes many technical problems. In this paper, we show that the isothermal Euler equation on , , with asymptotically vanishing density, , displays a specific large time behavior, in the sense that if the solution is global in time, then the density disperses with a rate different from the above one, and possesses a universal asymptotic Gaussian profile. This property remains when the convex pressure law satisfies , as well as for the Korteweg and quantum NavierStokes equations:
(1.4) 
with , where denotes the symmetric part of the gradient,
For this system, we still have conservation of mass and the energy
(1.5) 
where
satisfies
In the case and equation (1.4) is the precise system derived in [BruMe10], as a correction to the isothermal quantum Euler equation. We emphasize that, because of the lack of positivity of the term in the energy functional, only the barotropic variant – where with – is studied in references. Classically, a Bohm potential (corresponding to the term multiplied by in (1.4)) is also added, see [AScomp, GisVV15, Jungel, VasseurYu] for instance. In the case where the dissipation is absent (), but with capillarity (), we refer to [BDD07, AnMa09, AnMa12, AuHa17].
A loose statement of our main result reads (a more precise version is provided in the next section, see Theorem 2.9):
Theorem 1.2.
Let be a global weak solution to (1.4) with initial density/velocity satisfying
Then there exists a mapping such that
This theorem entails that, in contrast with the isentropic case, the density of solutions to (1.4) disperses as follows :
with a universal profile. This result applies to a notion of “weak solution” that is based on standard a priori estimates satisfied by smooth solutions to (1.4). We make precise the definition of such solutions in the next section, see Definition 2.1.
The main ingredient of the proof is to translate in terms of our isothermal equations a change of unknown functions introduced for the dispersive logarithmic Schrödinger equation in [CaGap]. This enables to transform (1.4) into a system with unknowns for which the associated energy is positivedefinite. A second feature of the new system is that, asymptotically in time, it reads (keeping only the dominating terms):
(1.6) 
where is the timedependent scaling mentioned in Theorem 1.2. By taking the divergence of the second equation and replace with the first one, we obtain then (keeping again only the dominating terms):
where is the FokkerPlanck operator In this last system, the first equation implies that converges to a stationary solution to the second equation. The analysis of the longtime behavior of solutions to this FokkerPlanck equation, as provided in [AMTU01], entails the expected result.
The outline of the paper is as follows. In the next section, we provide rigorous definitions of weak solutions and precise statements for our main result. Section 3 is then devoted to the longtime behavior of solutions to (1.4). In this section, we compute at first explicit solutions to (1.4) with Gaussian densities. These explicit computations motivate the introduction of the change of variable that we use afterwards. In what remains of this section we give an exhaustive proof of the precise version for Theorem 1.2. The longtime analysis mentioned here is based on the a priori existence of solutions. However, in the compressible setting, global existence of solutions is questionable. So, in the last section of the paper, we focus on the notion of weak solutions that we consider. At first, we present the a priori estimates which motivate their definition. We end the paper by proving a sequential compactness result. This sequential compactness property is a cornerstone for the proof of existence of weak solutions, see e.g. [Lio98, Fei04]. As for the large time behavior, we simply state a loose version of our result here (see Theorem 4.10 for the precise statement):
Theorem 1.3.
Assume , , with , and let . Let be a sequence of weak solutions to (1.4) on , enjoying uniformly the conservation of mass, as well as a suitable notion of energy dissipation, BDentropy dissipation, and MelletVasseur type inequality. Then up to the extraction of a subsequence, converges to a weak solution of (1.4) on .
It is for the system (2.8) in terms of , as mentioned above, that fairly natural a priori estimates are required in the above statement. Even though the notions of solution for (1.4) and (2.8) are equivalent (Lemma 2.6 below), we did not find a direct approach to express the pseudo energy, pseudo BDentropy and MelletVasseur type inequality mentioned above in a direct way in terms of , that is, without resorting to .
2. Weak solutions and large time behavior
We now state a precise definition regarding the notion of solution that we consider in this paper. Even though, in (1.4), the fluid genuine unknowns are and , the mathematical theory that we develop in Section 4 suits better to the unknowns and . Therefore we state our definition of weak solution in terms of these latter unknowns. Nevertheless, we shall keep these notations, even though no fluid velocity field underlies the computation of
Definition 2.1.
Let and Given , we call weak solution to (1.4) on any pair such that there is a collection satisfying

The following regularities:
with the compatibility conditions

Euler case : The following equations in
(2.1) 
Korteweg and NavierStokes cases : The following equations in
(2.2) with the symmetric part of , and the compatibility conditions:
(2.3) (2.4)
We emphasize that the above definition is essentially the “standard” one, up to the fact that we require . The reason for this assumption will become clear in the Subsection 4.1 where we will recall the a priori estimates motivating this definition (see Lemma 2.6, as well as the definition of the pseudoenergy in (2.13)).
Several remarks are in order. When the symbol alone appears, it must be understood as , while when the symbol appears alone, it is defined by Under the compatibility condition of item this yields a welldefined vectorfield. As for the stresstensors involved in the momentum equation (2.2), we emphasize that (2.3) reads formally
An originality of the previous definition is that in the case , we do not ask for the continuity equation in terms of but in terms of However, we prove here that the usual continuity equation as written in (2.2) is a consequence to this definition thanks to the regularity of and This is the content of the following lemma:
Lemma 2.2.
Proof.
By definition, we have
Here we note that (so that ). We can then multiply this equation by We obtain:
At this point we remark that, by definition of
and, since , the products of the identity below are welldefined:
Combining these equation entails
We conclude thus that:
∎
2.1. Rewriting of (1.4) with a suitable timedependent scaling
In the case where the density is defined for all time and is dispersive (in the sense that it goes to zero pointwise), it is natural to examine the behavior of near , since it gives an “asymptotic pressure law” as time goes to infinity. A consequence of our result is that the large time behavior in (1.4) is very different according to or . Herein, we assume that with and . Typically, when , we recover the isothermal case, , and we can also consider
with no other restriction on (in any dimension), or even the exotic case . The most general class of pressure laws that we shall consider is fixed by the following assumptions:
Assumption 2.3 (Pressure law).
The pressure is convex ( for all ), and satisfies
Resuming the approach from [CaGap] (the link between Schrödinger equation and EulerKorteweg equation is formally given by the Madelung transform), we change the unknown functions as follows. Introduce solution of the ordinary differential equation
(2.5) 
The reason for considering this equation will become clear in Subsection 3.1. We find in [CaGap], for slightly more general initial data:
Lemma 2.4.
Let , . Consider the ordinary differential equation
(2.6) 
It has a unique solution , and it satisfies, as ,
where
We sketch the proof of this lemma in Appendix B, without paying attention to the quantitative estimate of the remainder term. We now introduce the Gaussian , and we set
(2.7) 
where we denote by the spatial variable for and . Denoting , (1.4) becomes, in terms of these new unknowns,
(2.8) 
The analogue of Definition 2.1 is the following:
Definition 2.5.
Let and Given , we call weak solution to (2.8) on any pair such that there exists a collection satisfying

The following regularities:
with the compatibility conditions

Euler case : The following equations in
(2.9) 
Korteweg and NavierStokes cases : The following equations in
(2.10) with the symmetric part of and the compatibility conditions:
(2.11) (2.12)
Mimicking the proof of Lemma 2.2, we see that in the case , if is a weak solution to (2.10) on in the sense of Definition 2.5, then it satisfies
In view of (2.7), we check directly:
Lemma 2.6 (Equivalence of the notions of solution).
Remark 2.7.
If in Definition 2.1, we had required only , then the above equivalence would not hold. In the same spirit, the change of unknown (2.7) would make the notion of solution rather delicate in the case of the Newtonian NavierStokes equation, a case where typically . More generally, we do not consider velocities enjoying integrability properties, unless the density appear as a weight in the integral.
We define the pseudoenergy of the system (2.8) by
(2.13)  
where
which formally satisfies
(2.14) 
where the dissipation is defined by
(2.15)  
By convexity we have , and for , so . Note also the identities
(2.16) 
Recall the CsiszárKullback inequality (see e.g. [LogSob, Th. 8.2.7]): for with ,
Writing
the (formal) conservation of the mass for and the definition (2.7) imply that the pseudoenergy is nonnegative, .
As for global solutions, we have the following natural definition:
2.2. Main result: largetime behavior of weak solutions to (2.10)
With the previous definitions and remarks, a quantitative and precise statement of Theorem 1.2 reads as follows:
Theorem 2.9.
Remark 2.10.
Unlike in Theorem 1.1, no smallness assumption is made on at ( may even be linear in space), so there is no such geometrical structure on the initial velocity as in [Serre97, Gra98].
Remark 2.11.
Remark 2.12 (Wasserstein distance).
The points (b) and (c) of Theorem 2.9 imply the large time convergence of to in the Wasserstein distance , defined, for and probability measures, by
where varies among all probability measures on , and denotes the canonical projection onto the th factor. This implies, for instance, the convergence of fractional momenta (see e.g. [Vi03, Theorem 7.12])
(2.17) 
Back to the initial unknowns , Theorem 2.9 and (2.7) yield
as announced in Theorem 1.2, where the symbol means that only a weak limit is considered. However, in the special case of Gaussian initial data considered in Section 3.1, it is easy to check that all the assumptions of Theorem 2.9 are satisfied, and moreover that strongly in . Finally, another consequence of Lemma 2.4, the (proof of the) last point in Theorem 2.9, and (2.7) is
This shows that indeed, no a priori information can be directly extracted from the energy defined in (1.5).
3. From Gaussians to Theorem 2.9
This part of the paper is devoted to the large time behavior of solutions to (1.3) and its variants. We first compute explicit Gaussian solutions and then proceed to the proof of Theorem 2.9.
3.1. Explicit solution
In this section, we resume and generalize some results established in [Yuen12, CFY17]. The generalizations concern two aspects: we allow densities and velocities which are not centered at the same point (hence and below), and we consider the quantum NavierStokes equation.
3.1.1. Euler and Newtonian NavierStokes equations
We recall the compressible Euler equation for isothermal fluids on
(3.1) 
where . As noticed in [Yuen12], (3.1) has a family of explicit solutions with Gaussian densities and affine velocities centered at the same point. Allowing different initial centers for these quantities leads to considering
(3.2) 
with , . Seeking a solution of the form
and plugging this ansatz into (3.1), we obtain a set of ordinary differential equations:
(3.3)  
(3.4)  
(3.5) 
Mimicking [LiWa06], seeking and of the form
we check that the two equations in (3.3) are satisfied if and only if
(3.6) 
and we find
Remark 3.1.
Since the velocity is affine in , this computation also yields explicit solutions for the isothermal (Newtonian) NavierStokes equations, but not for its quantum counterpart, as we will see below.
3.1.2. Korteweg and quantum NavierStokes equations
As in [CFY17], we generalize (3.1) by allowing the presence of a Korteweg term (), and we extend this contribution by allowing a quantum dissipation (quantum NavierStokes equation, when ). We recall the isothermal Korteweg and quantum NavierStokes equations
(3.7) 
with , and where the Korteweg term is also equal to
which is called the Bohm’s identity. Proceeding as in the previous subsection, (3.3)–(3.5) become
(3.8)  
Again, we seek and of the form