Dissipative solutions and semiflow selection for the complete Euler system

TL;DR

This paper introduces a dynamical systems approach to select physically relevant solutions for the complete Euler system, ensuring consistency with admissibility criteria and stability of stationary states.

Contribution

It proposes a novel selection method based on dissipative solutions that satisfies the semiflow property and respects entropy maximization and weak-strong uniqueness.

Findings

Selection of a unique solution system for the Euler equations.

Ensures strong solutions are always selected when they exist.

Includes stationary states in the solution set, maintaining stability.

Abstract

To circumvent the ill-posedness issues present in various models of continuum fluid mechanics, we present a dynamical systems approach aiming at selection of physically relevant solutions. Even under the presence of infinitely many solutions to the full Euler system describing the motion of a compressible inviscid fluid, our approach permits to select a system of solutions (one trajectory for every initial condition) satisfying the classical semiflow property. Moreover, the selection respects the well accepted admissibility criteria for physical solutions, namely, maximization of the entropy production rate and the weak--strong uniqueness principle. Consequently, strong solutions are always selected whenever they exist and stationary states are stable and included in the selection as well. To this end, we introduce a notion of dissipative solution, which is given by a triple of density, momentum and total entropy defined as expectations of a suitable measure--valued solution.

Equations

$$\begin{array}{rl}{\partial }_t\varrho +{\mathrm{div}}_x\mathbf{m}& =0,\\ {\partial }_t\mathbf{m}+{\mathrm{div}}_x\left(\frac{\mathbf{m}\otimes \mathbf{m}}{\varrho }\right)+{\nabla }_{x}p& =0,\\ {\partial }_t\mathcal{E}+{\mathrm{div}}_x\left[\left(\mathcal{E}+p\right)\frac{\mathbf{m}}{\varrho }\right]& =0.\end{array}$$

$$\mathcal{E}=\frac{1}{2}\frac{|\mathbf{m}{|}^2}{\varrho }+\varrho e,$$

$$(\gamma -1)\varrho e=p,\text{mbox}where\gamma >1\text{mbox}istheadiabaticconstant.$$

$$p=\varrho \vartheta \text{mbox}yieldinge=c_v\vartheta ,c_v=\frac{1}{\gamma -1}.$$

$$\mathbf{m}\cdot \mathbf{n}{|}_{\partial \Omega }=0.$$

$$\varrho (0,\cdot )={\varrho }_0,\mathbf{m}(0,\cdot )={\mathbf{m}}_0,\mathcal{E}(0,\cdot )={\mathcal{E}}_0.$$

$${\partial }_t(\varrho s)+{\mathrm{div}}_x(s\mathbf{m})=0\text{mbox}or{\partial }_{t}s+\left(\frac{\mathbf{m}}{\varrho }\right)\cdot {\nabla }_{x}s=0,$$

$$s(\varrho ,\vartheta )=\log ({\vartheta }^{c_v})-\log (\varrho ).$$

$${\partial }_t(\varrho s)+{\mathrm{div}}_x(s\mathbf{m})\ge 0$$

$$S_i={\varrho }_{i}s({\varrho }_i,{\mathbf{m}}_i,{\mathcal{E}}_i),i=1,2.$$

$$[{\varrho }_1,{\mathbf{m}}_1,{\mathcal{E}}_1]{\succ }_{\mathcal{D}}[{\varrho }_2,{\mathbf{m}}_2,{\mathcal{E}}_2]$$

$$[{\varrho }_1(t,\cdot ),{\mathbf{m}}_1(t,\cdot ),{\mathcal{E}}_1(t,\cdot )]=[{\varrho }_2(t,\cdot ),{\mathbf{m}}_2(t,\cdot ),{\mathcal{E}}_2(t,\cdot )]\text{mbox}foranyt\in [0,\tau ];$$

$${\int }_{\Omega }{S}_1(t+,x)\mathrm{d}x\ge {\int }_{\Omega }{S}_2(t+,x)\mathrm{d}x\text{mbox}forallt\in (\tau ,\tau +\delta ).$$

$${\int }_{\Omega }{S}_1({\tau }_n+,x)\mathrm{d}x\ge {\int }_{\Omega }{S}_2({\tau }_n+,x)\mathrm{d}x\text{mbox}foralln\in \mathbb{N}.$$

$$\partial_{t}S+{\rm div}_{x}\left(S\frac{{\bf m}}{\varrho}\right)={\color[rgb]{1,0,0}\Sigma},\ S=\varrho s,$$

$$[{\varrho }_1,{\mathbf{m}}_1,{\mathcal{E}}_1]\succ [{\varrho }_2,{\mathbf{m}}_2,{\mathcal{E}}_2]\iff {\Sigma }_1\ge {\Sigma }_2\text{mbox}onanycompactsubsetof[0,\infty )\times \overline{\Omega }.$$

$$\text{mbox}thedensity\varrho (t,\cdot ),\text{mbox}themomentum\mathbf{m}(t,\cdot ),\text{mbox}thetotalentropyS(t,\cdot )=\varrho s(t,\cdot ),$$

$$t\mapsto {\int }_{\Omega }\varrho (t,\cdot )\phi \mathrm{d}x,t\mapsto {\int }_{\Omega }\mathbf{m}(t,\cdot )\cdot \mathbf{\phi }\mathrm{d}x$$

$$t\mapsto {\int }_{\Omega }S(t,\cdot )\phi \mathrm{d}x\in {\mathrm{BV}}_{\mathrm{loc}}[0,\infty );$$

$${\int }_{\Omega }S(t+,\cdot )\phi \mathrm{d}x,{\int }_{\Omega }S(t-,\cdot )\phi \mathrm{d}x\text{mbox}arewelldefinedfort\in [0,\infty )$$

$$S(0-)=S_0.$$

$${\mathbf{U}}_0=[{\varrho }_0,{\mathbf{m}}_0,S_0],\text{mbox}togetherwiththetotalenergy{E}_0\ge {\int }_{\Omega }\left[\frac{|{\mathbf{m}}_0{|}^2}{{\varrho }_0}+{\varrho }_{0}e({\varrho }_0,S_0)\right]\mathrm{d}x,$$

$$\displaystyle\mathbb{U}[t;({\bf U}_{0},E_{0})]=\Big{(}[\varrho(t,\cdot),{\bf m}(t,\cdot),S(t-,\cdot)],E_{0}\Big{)},$$

$$\displaystyle E_{0}\geq\int_{\Omega}\left[\frac{|{\bf m}(t,\cdot)|^{2}}{\varrho(t,\cdot)}+\varrho(t,\cdot)e\Big{(}\varrho(t,\cdot),S(t\pm,\cdot)\Big{)}\right]\,\,{\rm d}{x}.$$

$$\mathbb{U}[0;({\mathbf{U}}_0,E_0)]=({\mathbf{U}}_0,E_0)$$

$$\mathbb{U}\Big{[}t_{1}+t_{2};({\bf U}_{0},E_{0})\Big{]}=\mathbb{U}\Big{[}t_{1};\mathbb{U}[t_{2},({\bf U}_{0},E_{0})]\Big{]}\ \mbox{for any}\ t_{1},\ t_{2}\geq 0.$$

$$t\mapsto\mathbb{U}\Big{[}t;({\bf U}_{0},E_{0})\Big{]}=\Big{(}[\varrho(t,\cdot),{\bf m}(t,\cdot),S(t-,\cdot)],E_{0}\Big{)}$$

$$E_0={\int }_{\Omega }{\mathcal{E}}_0\mathrm{d}x.$$

$$\mathbb{U}\Big{[}t;({\bf U}_{0},E_{0})\Big{]}=\left(\Big{[}\varrho(t,\cdot),{\bf m}(t,\cdot),S(\varrho(t,\cdot),{\bf m}(t,\cdot),\mathcal{E}(t,\cdot))\Big{]},E_{0}\right)$$

$$\left([\varrho_{1},{\bf m}_{1},S_{1}],E_{0}\right)=\mathbb{U}\Big{[}\cdot\,;({\bf U}_{0},E_{0})\Big{]}$$

$${\sigma }_1(t):={\int }_{\Omega }(S_1(t)-S_0)\mathrm{d}x\le {\sigma }_2(t):={\int }_{\Omega }(S_2(t)-S_0)\mathrm{d}x\text{on}[0,\infty ),$$

Full text

The research of E.F. leading to these results has received funding from the Czech Sciences Foundation (GAČR), Grant Agreement 18–05974S. The Institute of Mathematics of the Academy of Sciences of the Czech Republic is supported by RVO:67985840.

M.H. gratefully acknowledges the financial support by the German Science Foundation DFG via the Collaborative Research Center SFB1283.

1. 1 Introduction

2. 1.1 Admissibility criteria

3. 1.2 Semiflow solutions

4. 1.3 Selection procedure

5. 2 Dynamical systems approach

6. 2.1 Trajectory space

7. 2.2 Solution space

8. 3 Dissipative solutions

9. 3.1 Phase space

10. 3.2 Dissipative measure–valued solutions

11. 3.3 Dissipative solution

12. 3.4 Sequential stability

13. 3.5 Existence of dissipative solutions

14. 4 Semiflow selection

15. 5 Concluding discussion

16. 5.1 Stability of equilibrium states

17. 5.2 Dafermos’ admissibility criteria

1. Introduction

The Euler system describing the motion of a general compressible inviscid fluid represents one of the basic models in the framework of continuum fluid mechanics. The unknown fields are the fluid density $\varrho=\varrho(t,x)$, the momentum ${\bf m}={\bf m}(t,x)$, and the energy $\mathcal{E}=\mathcal{E}(t,x)$ satisfying the system of partial differential equations:

[TABLE]

Writing the energy as a sum of its kinetic and internal components,

[TABLE]

we suppose that the pressure $p$ and the internal energy $e$ satisfy the caloric equation of state in the form

[TABLE]

In addition, we introduce the absolute temperature $\vartheta$ through the Boyle–Mariotte thermal equation of state:

[TABLE]

Supposing that the fluid occupies a bounded spatial domain $\Omega\subset R^{N}$, $N=1,2,3$ we impose the impermeability boundary condition

[TABLE]

Finally, the initial state of the fluid is given through the initial conditions

[TABLE]

The Second law of thermodynamics is enforced through the entropy balance equation

[TABLE]

where the entropy $s$ is given as

[TABLE]

There is a vast amount of literature dedicated to the mathematical theory of the Euler system. In particular, it is known that the initial–value problem is well posed locally in time in the class of smooth solutions, see e.g. the monograph by Majda [17] or the more recent treatment by Benzoni–Gavage and Serre [1]. Smooth solutions, however, develop singularities in a finite time for a fairly general class of initial data, see e.g. Smoller [19]. Thus if the Euler system is accepted as an adequate description of the fluid motion in a long run, a concept of generalized solution is needed.

The modern theory of partial differential equations is based on the concept of weak solution, where the derivatives are understood in the sense of distributions. This gives rise to a large class of objects in which uniqueness might be lost. Several admissibility criteria have been proposed to select the physically relevant weak solution, among which the entropy inequality

[TABLE]

reflecting the Second law of thermodynamics. The recent adaptation of the method of convex integration, developed in the context of incompressible fluids by De Lellis and Székelyhidi [10], gave rise to numerous examples of ill–posedness also in the class of compressible fluids, see Chiodaroli, De Lellis, Kreml [6], Chiodaroli and Kreml [7] , Chiodaroli et al. [8], among others. In particular, it was shown in [12] that the Euler system (1.1)–(1.5) is ill–posed, specifically it admits infinitely many weak solutions for a large class of initial data. Moreover, these solutions satisfy the entropy inequality (1.7). In addition, examples of regular initial data producing infinitely many weak solutions in the long run have been also obtained in [12].

1.1. Admissibility criteria

In view of these facts, more refined admissibility criteria are needed in order to select the physically relevant solutions. Dafermos [9] proposed a selection criterion based on maximality of the total entropy. We formulate it in a slightly stronger form, avoiding the issue of existence of right-sided time derivatives. Specifically, suppose that $[\varrho_{i},{\bf m}_{i},\mathcal{E}_{i}]$, $i=1,2$ are two solutions of the Euler system (1.1)–(1.5) with the associated total entropies

[TABLE]

We say that

[TABLE]

if:

• •

there exists $\tau\geq 0$ such that

[TABLE]

• •

there exists $\delta>0$ such that

[TABLE]

A weak solution $[\varrho,{\bf m},E]$ is called maximal (admissible) if it is maximal with respect to the relation ${\succ}_{\mathcal{D}}$.

A modification of this criterion was further investigated in the context of the barotropic Euler system in [11]. However, in this case it is rather the total energy that shall be minimized, according to the principle of maximal energy dissipation. The criterium of [11] translated to our setting of the full Euler system, leads to a weaker version of the condition (1.8). Namely, we require that there exists a sequence $(\tau_{n})_{n}$, $\tau_{n}>\tau$, $\tau_{n}\to\tau$ such that

[TABLE]

It was shown in [11] that the solutions constructed there by the method of convex integration do not fulfill the corresponding criterium for maximal energy dissipation, suggesting that such a criterium shall be retained in order to exclude nonphysical solutions.

To compare these two criteria, let us denote by ${\succ}_{\mathcal{F}}$ the partial ordering induced by (1.9). The following holds true: If a solution is maximal with respect to ${\succ}_{\mathcal{F}}$ then it is also maximal with respect to ${\succ}_{\mathcal{D}}$. This may seem surprising at first sight since (1.8) is obviously a stronger condition than (1.9). But it is exactly for this reason why the implication of maximality is valid. More precisely, the weaker condition (1.9) allows to compare more solutions, for instance also those that oscillate around each other and that are therefore not comparable by the condition (1.8).

An alternative criterion enforcing the Second law of thermodynamics is maximility of the global entropy production proposed in [3]. In accordance with the Schwartz representation principle, the inequality in the entropy balance (1.7) can be interpreted as

[TABLE]

where $\Sigma$ is a non–negative Borel measure on $[0,\infty)\times\overline{\Omega}$. Similarly to the above, we say that

[TABLE]

It was shown that maximal solutions with respect to the relation $\succ$ exist for the full Euler system under rather general hypotheses on the data as well as the constitutive relations.

Despite these efforts, however, none of the above selection criteria proved sufficient to guarantee the desired well-posedness result. The major open question therefore remains: how can physically relevant solutions be distinguished from the nonphysical ones? The aim of the present paper is to take into account another physical property of an evolution system, namely, the so-called semiflow property: starting the system at time [math], letting it run to time $s>0$ and then restarting and continue for the time $t>0$, the state of the system at the final time $s+t$ should be the same as if the system ran directly from [math] to $s+t$. If uniqueness holds, the semiflow property immediately follows. However, for systems where uniqueness is unknown or not valid, it is generally not clear whether such a semiflow even exists.

1.2. Semiflow solutions

On the following pages we show how to approach this problem and we construct a solution semiflow to the complete Euler system (1.1)–(1.4). In particular, this leads to a well-defined dynamical system associated to (1.1)–(1.4), which depends in a measurable way on the initial condition. To this end, it is in the first place necessary to identify the correct phase variables of the system together with a suitable notion of solution. Even though the system (1.1) describes the evolution of the density, momentum and energy, it turns out to be beneficial to replace the energy by the total entropy. In other words, the state of the fluid at a given time $t\geq 0$ will be determined by the values of three phase variables,

[TABLE]

interpreted through their spatial integral means as quantities in suitable abstract function spaces.

The reason why we prefer the entropy $S$ instead of the energy $\mathcal{E}$ is the lack of suitable a priori bounds for the latter. The integral means

[TABLE]

will be continuous for $t\in[0,\infty)$ for any smooth $\varphi$, $\boldsymbol{\varphi}$, while

[TABLE]

whence the limits

[TABLE]

with the convention

[TABLE]

Hence, another difficulty regarding the classical theory of dynamical systems stems from the fact that solutions are not continuous in time, which in particular holds here for the total entropy. We overcome this issue by replacing continuity by existence of one–sided limits at every time $t\geq 0$.

As the next step, we shall determine what input information on the initial state of the system is necessary. Apart from the initial state for the density, momentum and total entropy, we shall be also given the total energy $E_{0}$ which is a constant of the motion and provides various bounds for all the corresponding quantitities. Given the initial state of the system $\mathbb{U}_{0}=({\bf U}_{0},E_{0})$, where

[TABLE]

we then identify in a unique way the state of the system at the time $t\geq 0$,

[TABLE]

In addition, the mapping $t\mapsto\mathbb{U}[t;({\bf U}_{0},E_{0})]$ will enjoy the semiflow property:

• •

[TABLE]

• •

[TABLE]

In particular, the state of the system at a time $t\geq 0$ is uniquely determined in terms of the initial data $[\varrho_{0},{\bf m}_{0},S_{0}]$ and the initial energy $E_{0}$.

The trajectories

[TABLE]

represent a generalized - dissipative - solution to the Euler system (1.1)–(1.5), which we introduce in Section 3. It will be shown that they comply with the following stipulations:

• •

Weak–strong uniqueness. Suppose that the Euler system (1.1)–(1.5) admits a classical solution $[\varrho,{\bf m},\mathcal{E}]$ on a time interval $[0,T_{\rm max})$,

[TABLE]

Then for ${\bf U}_{0}=[\varrho_{0},{\bf m}_{0},S(\varrho_{0},{\bf m}_{0},\mathcal{E}_{0})]$ we have

[TABLE]

for all $t\in[0,T_{\rm max})$.

• •

Maximal entropy production. Suppose that

[TABLE]

and that $\left([\varrho_{2},{\bf m}_{2},S_{2}],E_{0}\right)$ is another dissipative solution starting from the same initial data $({\bf U}_{0},E_{0})$ and such that

[TABLE]

where $\sigma_{i}$ is the entropy production rate associated to $[\varrho_{i},{\bf m}_{i},S_{i}]$, $i=1,2$, respectively.

Then

[TABLE]

• •

Stability of stationary states. Let $\varrho=\overline{\varrho}>0$, ${\bf m}=0$, $\overline{\mathcal{E}}=\frac{1}{|\Omega|}E_{0}$ be a stationary solution of the Euler system (1.1). Suppose that the system reaches the equilibrium state at some time $T\geq 0$, i.e.

[TABLE]

Then

[TABLE]

for all $t\geq T$.

1.3. Selection procedure

The semiflow $\mathbb{U}$ will be constructed by means of the selection procedure originally proposed in the context of stochastic Markov processes by Krylov [16]. The method was later adapted to deterministic evolutionary problems with time continuous solutions by Cardona and Kapitanski [5]. We have developed a similar approach to the isentropic Euler system in [2], where we also relaxed the continuity assumption on the trajectories. We refer the reader to Section 2 for details of the abstract formulation.

In Section 3, we specify a class of generalized solutions to the Euler system (1.1)–(1.5) termed dissipative solutions. They are, loosely speaking, the expected values of suitable measure–valued solutions as introduced in [4] but with an additional refinement regarding the associated concentration defect measures.

Applying carefully the selection procedure of [16], [5], [2], we obtain the desired semiflow solution in Section 4. The paper is concluded by a discussion of further extensions in Section 5. In particular, in Section 5.2 we compare our construction to the two versions of Dafermos’ criterium introduced in Section 1.1.

2. Dynamical systems approach

We start by adapting the general approach of [16], [5] to problems with discontinuous solutions paths, see also [2]. Suppose that the state of a physical system at each instant $t\geq 0$ is characterized by an abstract vector $\mathbb{U}(t)\in\mathcal{S}$ ranging in a suitable phase space $\mathcal{S}$. In our setting, the state space will be a separable Hilbert space or a Polish topological space.

Consider a mapping

[TABLE]

meaning $\mathbb{U}[t;\mathbb{U}_{0}]$ is the state of the system emanating from the initial state $\mathbb{U}_{0}$ at a time $t\geq 0$.

Definition 2.1** (Semiflow).**

We say that a mapping $\mathbb{U}:[0,\infty)\times\mathcal{S}\to\mathcal{S}$ is a semiflow if:

1. (a)

The flow starts at the initial datum, that is

[TABLE] 2. (b)

It has the semigroup property, that is

[TABLE]

In the classical theory of dynamical systems, continuity of the semiflow $\mathbb{U}$ is required both in $t$ and in $\mathbb{U}_{0}$. In view of the issues discussed in the previous section, we develop a generalized theory, where continuity in time is relaxed to the existence of the one–sided limits at $t\pm$, while continuity in $\mathbb{U}_{0}$ is replaced by measurability with respect to Borel sets generated by suitable topologies.

2.1. Trajectory space

The trajectories space can be defined as $\mathcal{T}=\mathrm{BV}_{\rm loc}([0,\infty);\mathcal{S}).$

For trajectories we can define the following operations:

• •

Time shift. For $\xi\in\mathcal{T}$ and $T>0$ we set

[TABLE]

• •

Continuation. For $\xi^{1},\ \xi^{2}\in\mathcal{T}$ and $T>0$ we set

[TABLE]

2.2. Solution space

The solution space can be loosely described as the family of all solutions of a given system of equations emanating from a fixed initial datum in $\mathcal{S}$. As such the solution operator $\mathcal{U}$ can be understood as a mapping

[TABLE]

Definition 2.2** (Solution operator).**

We say that a mapping $\mathcal{U}:\mathcal{S}\to 2^{\mathcal{T}}$ is solution operator if the following properties are satisfied:

1. [A1]

Existence, compactness. For each $\mathbb{U}_{0}\in\mathcal{S}$ the set $\mathcal{U}[\mathbb{U}_{0}]$ is a non–empty bounded subset of $\mathrm{BV}_{\rm loc}([0,\infty);\mathcal{S})$, compact in the topology of $L^{1}_{{\rm loc}}([0,\infty);\mathcal{S})$,

[TABLE] 2. [A2]

Measurability. The set valued mapping

[TABLE]

is Borel measurable; where $2^{\mathcal{T}}$ is endowed with the Hausdorff topology defined on compact subsets of the metric space $L^{1}_{{\rm loc}}([0,\infty);\mathcal{S})$. 3. [A3]

Shift invariance. For any $\mathbb{U}\in\mathcal{U}[\mathbb{U}_{0}]$ and any $T>0$, we have

[TABLE] 4. [A4]

Continuation. For any $\mathbb{U}^{1}\in\mathcal{U}[\mathbb{U}_{0}]$, $T>0$, and $\mathbb{U}^{2}\in\mathcal{U}[\mathbb{U}^{1}(T-)]$, we have

[TABLE]

A version of the following result for the case of $\mathcal{T}$ being a set of continuous trajectories with values in a Polish space was proved by Cardona and Kapitanski [5]. However, for applications in fluid dynamics, where time continuity of the energy or entropy is not valid, a suitable modification is necessary which was proved in [2, Section 5].

Proposition 2.3**.**

Let $\mathcal{S}$ be a Polish space and let $\mathcal{T}=\mathrm{BV}_{\rm loc}([0,\infty);\mathcal{S})$. Let $\mathcal{U}:\mathcal{S}\to 2^{\mathcal{T}}$, be a solution operator enjoying the properties [A1]–[A4] from Definition 2.2.

Then $\mathcal{U}$ admits a measurable semiflow selection, specifically, for any $\mathbb{U}_{0}\in\mathcal{S}$ there exists a single trajectory

[TABLE]

such that

[TABLE]

is a semiflow in the sense of Definition 2.1.

In addition, if $\beta:\mathcal{S}\to R$ is a bounded continuous function and $\lambda>0$, the selection can be chosen to satisfy

[TABLE]

We recall that the selection procedure of Proposition 2.3 relies on a subsequent minimization (or alternatively maximization) of a sequence of suitable continuous functionals, as e.g. the functional in (2.1), over the set of all solutions. The functionals are chosen in a way to separate points of the trajectory space $\mathcal{T}$. For this purpose the Laplace transform with respect to the time variable proved to be beneficial. In general, the selection depends on the particular choice of a sequence $(\lambda_{n})_{n\in\mathbb{N}}$ representing the points where the Laplace transform is evaluated. The functional in (2.1) shall then be chosen as the first one to minimize/maximize. In our application to the complete Euler system in the next section we chose (2.1) in order to maximize the entropy production rate.

In the remaining part of the paper, we identify a suitable solution operator $\mathcal{U}$ associated to the Euler system that complies with the hypotheses of Proposition 2.3. In Section 5, we present further discussion on the question of optimality of the choice of the functionals and in particular the relation to the Dafermos criterion discussed in the Introduction.

3. Dissipative solutions

Our first goal in this section is to identify a suitable class of generalized solutions to the Euler problem (1.1)–(1.5). We start by introducing the basic state variables, then we recall the (slightly modified) notion of dissipative measure–valued solutions from [4]. This leads us to our notion of dissipative solution which we introduce in Section 3.3. Afterwards we study stability of dissipative solutions as well as their existence. The semiflow is then constructed in Section 4.

3.1. Phase space

In what follows, we plan to work with the state variables $\varrho$, ${\bf m}$, and the total entropy $S=\varrho s$. In accordance with hypotheses (1.2), (1.3), the pressure $p$ and the internal energy $e$ can be written in the form

[TABLE]

Lemma 3.1**.**

The mapping

[TABLE]

is strictly convex.

Proof.

This is a matter of direct computation of the Hessian matrix. We have

[TABLE]

and

[TABLE]

Obviously the Hessian matrix has positive trace, while its determinant reads

[TABLE]

which completes the proof. ∎

Consequently, we may define

[TABLE]

which is a convex lower semi–continuous function on $[0,\infty)\times R$.

Similarly, we define the kinetic energy

[TABLE]

which is a convex lower semi–continuous function on $[0,\infty)\times R^{N}$. We conclude that the total energy

[TABLE]

may be viewed as a convex lower semi–continuous function of $\varrho\geq 0$, ${\bf m}\in R^{N}$, and $S\in R$.

3.2. Dissipative measure–valued solutions

Following [4] we introduce the concept of dissipative measure–valued solution to the Euler problem (1.1)–(1.5). In addition, similarly to [2], we refine the definition of the measures describing the concentration defect. The reader may consult [4] for the physical background and mathematical objects like Young measures used in what follows.

We start by introducing the state space of “dummy variables”:

[TABLE]

By $\mathcal{P}(\mathcal{Q})$ we denote the set of probability measures on $\mathcal{Q}$, whereas $\mathcal{M}^{+}(\overline{\Omega})$ and $\mathcal{M}^{+}(S^{N-1}\times\overline{\Omega})$ denotes the set of positive bounded Radon measures on $\overline{\Omega}$ and $S^{N-1}\times\overline{\Omega}$, respectively, where $S^{N-1}\subset R^{N}$ is the unit sphere.

A dissipative measure–valued solution of the Euler system (1.1)–(1.4) with the initial data

[TABLE]

consists of the following objects:

• •

a parametrized family of probability measures

[TABLE]

• •

kinetic and internal energy concentration defect measures

[TABLE]

• •

convective concentration defect measure

[TABLE]

The Euler equations are satisfied in the following sense:

• •

Total energy is a constant of motion:

[TABLE]

• •

Mass conservation or continuity equation $\mbox{\eqref{I1}}_{1}$ reads

[TABLE]

for any $\varphi\in C^{1}_{c}([0,\infty)\times\overline{\Omega})$.

• •

Momentum balance (equation $\mbox{\eqref{I1}}_{2}$) reads

[TABLE]

for any $\boldsymbol{\varphi}\in C^{1}_{c}([0,\infty)\times\overline{\Omega};R^{N})$, $\boldsymbol{\varphi}\cdot{\bf n}|_{\partial\Omega}=0$.

• •

Entropy balance (inequality (1.7)) is rewritten in the renormalized form:

[TABLE]

for any $\varphi\in C^{1}_{c}([0,\infty)\times\overline{\Omega})$, $\varphi\geq 0$, and any $Z$,

[TABLE]

It follows from (3.2), (3.3) that

[TABLE]

are continuous functions of time. Accordingly, equations (3.2), (3.3) can be written as

[TABLE]

for any $0\leq\tau_{1}<\tau_{2}<\infty$, and any $\varphi\in C^{1}_{c}([0,\infty)\times\overline{\Omega})$, where $\left<\mathcal{V}_{0,x};\tilde{\varrho}\right>=\varrho_{0}(x)$;

[TABLE]

for any $0\leq\tau_{1}<\tau_{2}<\infty$, and any $\boldsymbol{\varphi}\in C^{1}_{c}([0,\infty)\times\overline{\Omega};R^{N})$, $\boldsymbol{\varphi}\cdot{\bf n}|_{\partial\Omega}=0$.

As shown in [4, Section 2.1.1], the renormalized entropy inequality (3.4) implies that

[TABLE]

as soon as

[TABLE]

This is the minimum principle for the entropy $s=\frac{S}{\varrho}\geq s_{0}$. From now on, we fix $s_{0}$ and consider only solutions satisfying (3.7). This corresponds to having a new entropy $s-s_{0}\geq 0$. Then one can perform the limit passage $Z(s)\nearrow s$ in (3.4) obtaining the entropy balance

[TABLE]

cf. [3]. In particular,

[TABLE]

where $\chi^{1}$ is continuous and $\chi^{2}$ non–decreasing. Thus (3.8) can be rewritten in the form

[TABLE]

for any $0\leq\tau_{1}<\tau_{2}<\infty$, and any $\varphi\in C^{1}([0,\infty)\times\overline{\Omega})$, $\varphi\geq 0$, where $\left<\mathcal{V}_{0-,x};{\tilde{S}}\right>=S_{0}(x)$.

Now we have all in hand to formulate the definition of dissipative measure–valued solution.

Definition 3.2** (Dissipative measure–valued solution).**

A dissipative measure–valued solution of the Euler system (1.1)–(1.4) with the initial data

[TABLE]

is a parameterized family of probability measures

[TABLE]

together with the energy concentration defect measures

[TABLE]

and the convection concentration defect measure

[TABLE]

satisfying the integral identities (3.1), (3.5), (3.6), and (3.9).

Next, we list certain bounds that can be derived from hypothesis (3.7) and the energy equality (3.1). As the entropy is bounded below by $s_{0}$, we deduce from (3.1) that

[TABLE]

Similarly, as

[TABLE]

we conclude

[TABLE]

Finally, we derive some bounds on the total entropy $S$. Recall that $\tilde{S}\geq s_{0}\tilde{\varrho}$; whence

[TABLE]

If $\tilde{S}>0$, we compute

[TABLE]

Consequently,

[TABLE]

Next, we estimate the quantity $\tilde{S}/\sqrt{\tilde{\varrho}}$. If $\tilde{S}\leq 0$, we get, repeating the above argument,

[TABLE]

If $\tilde{S}>0$, we write

[TABLE]

We may therefore infer that

[TABLE]

Finally, we recall the weak–strong uniqueness principle proved in [4, Theorem 3.3].

Proposition 3.3**.**

Let $\Omega\subset R^{N}$, $N=1,2,3$ be a bounded domain with smooth boundary. Suppose that the Euler system (1.1)–(1.5) admits a classical solution $[\varrho,{\bf m},\mathcal{E}]$ in the class

[TABLE]

with the initial data

[TABLE]

Let $\mathcal{V}_{t,x}$ be a dissipative measure valued solution specified in Definition 3.2 starting from the data

[TABLE]

Then

[TABLE]

and

[TABLE]

for a.a. $(t,x)\in[0,T)\times\Omega$.

Note that existence of a local–in–time classical solution in the class (3.14) was established by Schochet [18].

3.3. Dissipative solution

Having collected the necessary preliminary material, we are ready to introduce the central object of the present paper - the dissipative solutions to the Euler system.

Definition 3.4** (Dissipative solution).**

The quantity $([\varrho,{\bf m},S],E_{0})$, where

[TABLE]

is a dissipative solution of the Euler system (1.1) with the initial data

[TABLE]

if there exists a dissipative measure–valued solution $\mathcal{V}_{t,x}$ as specified in Definition 3.2 such that

[TABLE]

Remark 3.5**.**

In accordance with the bounds (3.1), (3.10), (3.11), and (3.12), any dissipative solution belongs to the class

[TABLE]

whence to bounded balls in the afore mentioned spaces. These are compact metric (Polish) spaces with respect to the weak topology. In particular, condition (3.15) reduces to

[TABLE]

Finally, we introduce a subclass of dissipative solutions that reflect the physical principle of maximal dissipation defined via the entropy production rate

[TABLE]

discussed in Section 1.2. Let $([\varrho^{i},{\bf m}^{i},S^{i}],E_{0})$, $i=1,2$, be two dissipative solutions starting from the same initial data $([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})$ with entropy production rates $\sigma^{1}$ and $\sigma^{2}$. Similarly to (1.10), we introduce the relation

[TABLE]

Definition 3.6** (Maximal dissipative solution).**

A dissipative solution $([\varrho,{\bf m},S],E_{0})$ starting from the initial data $([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})$ is a maximal dissipative solution if it is maximal with respect to the relation $\succ$. Specifically, if

[TABLE]

where $([\tilde{\varrho},\tilde{{\bf m}},\tilde{S}],{E}_{0})$ is another dissipative solution starting from $([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})$, then

[TABLE]

Here $\sigma$ and $\tilde{\sigma}$ are the entropy production rates of $([\tilde{\varrho},\tilde{{\bf m}},\tilde{S}],{E}_{0})$ and $([\varrho,{\bf m},S],E)$ respectively.

3.4. Sequential stability

We start by introducing suitable topologies on the space of the initial data and the space of dissipative solutions. Fix $\ell>N/2+1$ and consider the Hilbert space

[TABLE]

together with its subset containing the initial data

[TABLE]

Note that $\mathcal{D}_{\rm Euler}$ is a closed convex subset of $\mathcal{S}_{\rm Euler}$. We also define the trajectory space

[TABLE]

In the following we are going to show sequential stability (compactness) of the set of dissipative solutions. This will be subsequently used in the proof of existence of dissipative solutions (using a suitable approximation) as well as measurability of the mapping

[TABLE]

where $\mathcal{U}([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})$ denotes the solution set

[TABLE]

for the initial data $([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})\in\mathcal{D}_{\rm Euler}$. We have the following result.

Proposition 3.7**.**

Suppose that $\{([\varrho_{0,\varepsilon},{\bf m}_{0,\varepsilon},S_{0,\varepsilon}],E_{0,\varepsilon})\}_{\varepsilon>0}\subset\mathcal{D}_{\rm Euler}$ is a sequence of data giving rise to a family of dissipative solutions $\{([\varrho_{\varepsilon},{\bf m}_{\varepsilon},S_{\varepsilon}],E_{0,\varepsilon})\}_{\varepsilon>0}$, that is,

[TABLE]

Moreover, we assume that there exists $\overline{E}>0$ such that $E_{0,\varepsilon}\leq\overline{E}$ for all $\varepsilon>0$.

Then, at least for suitable subsequences,

[TABLE]

and

[TABLE]

where

[TABLE]

Proof.

We mainly follow the ideas from [2, Proposition 3.1], which we refer to for further details. Some modifications are needed to accommodate the additional variable $S$ (in the internal energy and the entropy balance).

First, we claim that the convergence in (3.17) follows immediately from the boundedness of the initial energy $E_{0,\varepsilon}$ and the uniform bounds (3.10)–(3.12). Indeed, using Jensen’s inequality and convexity of the energy established in Lemma 3.1, we deduce that

[TABLE]

where $\mathcal{V}^{\varepsilon}$ is the (Young) measure associated with the solution $([\varrho_{\varepsilon},{\bf m}_{\varepsilon},S_{\varepsilon}],E_{0,\varepsilon})$. Integrating over $\Omega$ we can see that the right hand side is bounded by $E_{0,\varepsilon}\leq\overline{E}$ using (3.1).

On the other, we can use (3.10)–(3.13) to obtain the estimates

[TABLE]

uniformly in $\varepsilon$. Consequently, we deduce from equations (3.2), (3.3), and the energy balance (3.1), that

[TABLE]

where

[TABLE]

Similarly, we get from (3.19) that

[TABLE]

Moreover, we deduce from the entropy balance (3.9) and Helly’s selection theorem that

[TABLE]

modifying $S$ on the set of times of measure zero as the case may be (we split $\varphi=\varphi^{+}-\varphi^{-}$ and argue for both separately). Note that it is enough to show the former one for a dense subset of $W^{\ell,2}(\Omega)$ which follows from a diagonal argument. Seeing that $L^{\gamma}(\Omega)$ endowed with the weak topology is compactly embedded in $W^{-\ell,2}(\Omega)$ we obtain the desired conclusion

[TABLE]

Obviously,

[TABLE]

passing to a subsequence as the case may be. Next, it is easy to perform the limit in the equation of continuity (3.2) to obtain (3.5).

Next, as a consequence of the energy balance (3.1), the Young measures $\mathcal{V}^{\varepsilon}_{t,x}$ have uniformly bounded first moments; whence

[TABLE]

at least for a subsequence. In addition, using estimates (3.11)–(3.13) we can pass to the limit in the entropy balance (3.9).

The limit in the total energy balance (3.1) and the momentum balance (3.3) involves the concentration measures and is more technical. However, the procedure is exactly the same as in [2, Section 3] therefore we omit the proof here. ∎

3.5. Existence of dissipative solutions

The sequential stability from the previous part combined with a suitable approximation implies the existence of a dissipative solution. The precise statement is the content of the following assertion.

Proposition 3.8**.**

Let $([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})\in\mathcal{D}_{\rm Euler}$ be given. Then the Euler system (1.1)–(1.4) admits a dissipative solution in the sense of Definition 3.4 with the initial data $([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})$.

Proof.

Given initial data $([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})\in\mathcal{D}_{\rm Euler}$ implies

[TABLE]

with the respective bounds in terms of $E_{0}$, recall the lower bounds for the energy obtained in (3.10)–(3.12). It is standard to construct smooth functions $\varrho_{0,\varepsilon}$, ${\bf m}_{0,\varepsilon}$ and $S_{0,\varepsilon}$, where $\varrho_{0,\varepsilon}$ is strictly positive and ${\bf m}_{0,\varepsilon}$ compactly supported in $\Omega$, $S_{0,\varepsilon}\geq\varrho_{0,\varepsilon}s_{0}$ such that

[TABLE]

and

[TABLE]

as $\varepsilon\to 0$. Finally, we set $s_{0,\varepsilon}=S_{0,\varepsilon}/\varrho_{0,\varepsilon}$ and ${\bf u}_{0,\varepsilon}={\bf m}_{0,\varepsilon}/\varrho_{0,\varepsilon}$.

Now, similarly to Kröner and Zajaczkowski [15], we consider the approximate system

[TABLE]

with the initial data

[TABLE]

Here, the symbol $\mathcal{L}$ denotes a suitable “viscosity” operator, with the associated set of boundary conditions to be imposed on the velocity field ${\bf u}$. In [15], the authors consider $\mathcal{L}=\Delta_{x}^{3}$ with the Dirichlet boundary conditions. This is not convenient here, as the resulting admissible test functions $\boldsymbol{\varphi}$ in the momentum equation (3.3) would be compactly supported in $\Omega$. In order to preserve the weak–strong uniqueness principle (Proposition 3.3), however, we need a larger class of test functions satisfying merely the impermeability condition $\boldsymbol{\varphi}\cdot{\bf n}|_{\partial\Omega}=0$. To this end, we use a different ansatz inspired by Kato and Lai [14].

Let $W^{3,2}_{n}(\Omega;R^{N})$ be the Hilbert space,

[TABLE]

We suppose that $\Omega\subset R^{N}$, $N=1,2,3$ is a bounded regular domain so that $W^{3,2}_{n}(\Omega;R^{N})$ is compactly embedded in $C^{1}(\overline{\Omega};R^{N})$. Moreover, $W^{3,2}_{n}$ being a closed subset of the separable Hilbert space $W^{3,2}$ is separable. Let $((\ ;\ ))$ be the scalar product on $W^{3,2}_{n}$, that is,

[TABLE]

Similarly to [14], we consider $\mathcal{L}$ to be the self–adjoint operator on $W^{3,2}_{n}$ associated to the bilinear form $((\ ;\ ))$, namely

[TABLE]

The associated variational formulation of the momentum equation in (3.22) then reads

[TABLE]

for any $\tau>0$, and any $\boldsymbol{\varphi}\in W^{3,2}_{n}(\Omega;R^{N})$.

As $W^{3,2}_{n}$ is a separable Hilbert space, the existence proof used in [15, Theorem 4.1] applies without essential modifications (see also Hoff [13] for how to handle the transport equations for $\varrho$ and $s$). Given $\varepsilon>0$, there exists a trio $(\varrho_{\varepsilon},{\bf u}_{\varepsilon},s_{\varepsilon})$ of the following class which is a classical solution to (3.22) in the following sense:

• •

The balance of mass (3.22)1 and the entropy balance (3.22)3 hold pointwise in $(0,\infty)\times\Omega$.

• •

The balance of momentum (3.22)2 holds as (3.24).

• •

The solution possesses enough regularity to rigorously justify the standard energy estimate

[TABLE]

Finally, we set $S_{\varepsilon}=\varrho_{\varepsilon}s_{\varepsilon}$ and ${\bf m}_{\varepsilon}=\varrho_{\varepsilon}{\bf u}_{\varepsilon}$. Using the arguments of the preceding section, it is easy to perform the limit $\varepsilon\to 0$ in the sequence of approximate solutions

[TABLE]

to obtain the desired dissipative solution, with the exception of the energy equality (3.1) that now reads

[TABLE]

Note that, on account of (3.25), it is easy to see that viscous term in the momentum equation vanishes as $\varepsilon\to 0$. In contrast with (3.26), the entropy balance (3.9) holds as equality. To convert (3.26) to equality, it is enough to augment $\mathfrak{C}_{\rm int}(t)$ by $h(t)\,{\rm d}{x}$ with a suitable spatially homogeneous $h\geq 0$. Note that the momentum equation (3.3) remains valid as $\mathfrak{C}_{\rm int}(t)$ acts on ${\rm div}_{x}\boldsymbol{\varphi}$, where $\boldsymbol{\varphi}\cdot{\bf n}|_{\partial\Omega}=0$; whence

[TABLE]

∎

4. Semiflow selection

The goal of this section is to show that there is a semiflow selection to the Euler system (1.1). We recall that the precise definition of a semiflow - in the abstract framework - is given in Definition 2.1. The following is the main result of the present paper.

Theorem 4.1** (Semiflow solution).**

The Euler system (1.1)–(1.4) admits a semiflow solution in the class of dissipative solutions in the sense of Definition 3.4. More specifically, for any initial data

[TABLE]

there exists a dissipative solution

[TABLE]

enjoying the following properties:

• •

for each $\mathbb{U}_{0}\in\mathcal{D}_{\rm Euler}$ the solution $\mathbb{U}[\cdot,\mathbb{U}_{0}]$ is maximal in the sense of Definition 3.6;

• •

$\mathcal{D}_{\rm Euler}$* is an invariant set, meaning*

[TABLE]

for any $t\geq 0$;

• •

the mapping

[TABLE]

is Borel measurable;

• •

the mapping

[TABLE]

is a semiflow, specifically,

[TABLE]

for any $t_{1},t_{2}\geq 0$, and any $\mathbb{U}_{0}\in\mathcal{D}_{\rm Euler}$.

The claim of Theorem 4.1 will follow from Proposition 2.3 as soon as we verify [A1]–[A4] from Definition 2.2 for the solution set $\mathcal{U}([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})$ with $([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})\in\mathcal{D}_{\rm Euler}$. This will be done in the following lemmas.

Lemma 4.2**.**

For each $([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})\in\mathcal{D}_{\rm Euler}$ the set $\mathcal{U}[\varrho_{0},{\bf m}_{0},S_{0},E_{0}]$ is a non–empty bounded subset of $\mathrm{BV}_{\rm loc}([0,\infty);\mathcal{S}_{\rm Euler})$ and it is compact in the topology of $L^{1}_{{\rm loc}}([0,\infty);\mathcal{S}_{\rm Euler})$.

Proof.

The claim follows from Propositions 3.7 and 3.8. ∎

Lemma 4.3**.**

Let $([\varrho,{\bf m},S],E_{0})$ be a dissipative solution to the Euler system in the sense of Definition 3.4. Then $\mathfrak{S}_{T}\circ([\varrho,{\bf m},S],E_{0})$ is a dissipative solution corresponding to the data $([\varrho,{\bf m},S](T-),E_{0})$.

Proof.

We recall that the operation of time shift is given by

[TABLE]

Since a shift of a test function in (3.1)–(3.4) produces a test function in the same class, the claim easily follows. ∎

Lemma 4.4**.**

Let $([\varrho^{1},{\bf m}^{1},S^{1}],E_{0})$ be a dissipative solution of the Euler system, $T>0$ and let $([\varrho^{2},{\bf m}^{2},S^{2}],E_{0})$ be another dissipative solution with the initial data $([\varrho^{1},{\bf m}^{1},S^{1}](T-),E_{0})$.

Then

[TABLE]

is a dissipative solution of the Euler system.

Proof.

The concatenation property for (3.1) and (3.2) is obvious. As far as (3.3) and (3.4) are concerned we observe the following. The energy is convex lower semicontinuous, in as much we obtain

[TABLE]

The function

[TABLE]

is concave such that

[TABLE]

So, (3.3) and (3.4) have the concatenation property as well. ∎

Lemma 4.5**.**

The mapping

[TABLE]

is Borel measurable; where $2^{\mathcal{T}_{\rm Euler}}$ is endowed with the Hausdorff topology defined on compact subsets of the metric space $L^{1}_{{\rm loc}}([0,\infty);\mathcal{S}_{\rm Euler})$.

Proof.

As a consequence of the sequential stability from Propositions 3.7 the claim follows from [20, Thm. 12.1.8]. ∎

Proof of Theorem 4.1.

In view of Proposition 2.3, the first claim of Theorem 4.1 - the existence of the semiflow $\mathbb{U}$ - follows now from Lemmas 4.2–4.5. Let us finally explain why $\mathbb{U}([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})$ can be selected to be maximal in the sense of Definition 3.6. As stated in Proposition 2.3 the selection can be chosen to satisfy

[TABLE]

for any $([\tilde{\varrho},\tilde{{\bf m}},\tilde{S}],E_{0})\in\mathcal{U}([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})$. Here $\beta:\mathcal{S}_{\rm Euler}\to R$ is a bounded and continuous function. We suppose that

[TABLE]

where $\alpha:R\rightarrow R$ is smooth, bounded and strictly decreasing. We proceed by contradiction. Let $([\tilde{\varrho},\tilde{{\bf m}},\tilde{S}],{E}_{0})\in\mathcal{U}([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})$ be such that $([\tilde{\varrho},\tilde{{\bf m}},\tilde{S}],{E}_{0})\succ([\varrho,{\bf m},S],E_{0}),$ that is, $\int_{\Omega}\tilde{S}\,{\rm d}{x}\geq\int_{\Omega}S\,{\rm d}{x}$ in $(0,\infty),$ where we denote $([\varrho,{\bf m},S],E_{0})=\mathbb{U}([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})$. Then we get

[TABLE]

whence $\int_{\Omega}S\,{\rm d}{x}=\int_{\Omega}\tilde{S}\,{\rm d}{x}$ a.a. in $(0,\infty)$ since $\alpha$ is strictly decreasing. Hence $\mathbb{U}([\varrho_{0},{\bf m}_{0},S_{0}],E_{0})$ is maximal with respect to $\succ$. The proof of Theorem 4.1 is hereby complete. ∎

5. Concluding discussion

We have shown the existence of a semiflow solution to the complete Euler system in the class of dissipative solutions. The semiflow solution has been selected among other solutions starting from the same initial data. The major issue is, of course, uniqueness of such a selection. In accordance with Proposition 3.3 the selected solution is unique and coincides with the strong solutions emanating from the same initial state as long as the latter exists. There is a special class of strong solutions - the equilibrium states. As we shall see below, they enjoy certain stability properties thanks to the fact that the selected solutions produce the maximal amount of entropy.

5.1. Stability of equilibrium states

The equilibrium states are stationary solutions to the Euler system (1.1). Specifically, the momentum vanishes identically ${\bf m}\equiv 0$, while the density $\varrho=\overline{\varrho}>0$ and the entropy $S=\overline{S}$ are constant. Note that $\overline{\varrho}$ is uniquely determined by the total mass

[TABLE]

Given $\overline{\varrho}$ and the total energy $E_{0}$, the (constant) equilibrium entropy $\overline{S}$ is the unique maximizer of the total entropy among all states with given mass and energy:

[TABLE]

Indeed, suppose that $\overline{\varrho}$, $\overline{S}$ are constant, and

[TABLE]

Normalizing the above relation by a factor $\frac{1}{|\Omega|}$ we may suppose $|\Omega|=1$. Using convexity of the function $\Phi(\varrho,S)\equiv\varrho e(\varrho,S)$ (cf. Lemma 3.1) and Jensen’s inequality, we get

[TABLE]

As $\Phi$ is strictly increasing in $S$ and strictly convex in $(\varrho,S)$, we get successively $\int_{\Omega}S\,\,{\rm d}{x}=\int_{\Omega}\overline{S}\,\,{\rm d}{x}$, and $\varrho=\overline{\varrho}$, $S=\overline{S}$.

We claim that a maximal solution can only see the equilibria with the energy $E_{0}$.

Proposition 5.1** (Stability of equilibria).**

Let $([\varrho,{\bf m},S],E_{0})$ be a dissipative solution to the Euler system emanating from the initial data

[TABLE]

Suppose that $([\varrho,{\bf m},S],E_{0})$ is maximal in the sense of Definition 3.6.

Then

[TABLE]

Proof.

As the equilibrium solution maximizes the total entropy, we have $\overline{S}\geq S_{0}$; whence (5.1) is a dissipative solution. By the same token, it is a maximal solution and any other solution “larger” in the sense of $\succ$ has the same entropy, meaning it concides with the solution (5.1). ∎

5.2. Dafermos’ admissibility criteria

Finally, we discuss some implications of Dafermos’ admissibility criteria introduced in Section 1.1. As we have seen in the proof of Theorem 4.1, maximimality of the entropy production by dissipative solutions has been enforced by maximizing the integral

[TABLE]

modulo the cut–off function $\alpha$. Actually the same result could have been achieved by maximizing

[TABLE]

Thus, a natural question arises, namely, what would be the resulting solution if we maximized successively the entropy for a family $\lambda_{n}\to\infty$. The following result implies: if there is a solution such that it has maximal Laplace transform of the total entropy evaluated at every $\lambda_{n}$ from a sequence $\lambda_{n}\to\infty$, then it satisfies the entropy production criterium by Dafermos in the form (1.8), that is, it is maximal with respect to the partial ordering $\succ_{\mathcal{D}}$.

With this motivation in mind, we employ the following notation for $F,G\in L^{\infty}(0,\infty)\cap\mathrm{BV}_{\rm loc}([0,\infty))$: We say that $F\succ_{\mathcal{D}}G$ provided there exists $\delta>0$ such that $F(t+)\geqslant G(t+)$ for all $t\in(0,\delta)$; and we say that $F\succ_{\mathcal{F}}G$ provided there exists a sequence $(\tau_{n})_{n}$, $\tau_{n}>0$, $\tau_{n}\to 0$ such that $F(\tau_{n}+)\geqslant G(\tau_{n}+)$.

Lemma 5.2**.**

Let $\mathcal{G}\subset L^{\infty}(0,\infty)\cap\mathrm{BV}_{\rm loc}([0,\infty))$ and assume that there exists $F\in\mathcal{G}$ such that for some sequence $\lambda_{n}\to\infty$ it holds

[TABLE]

*Then for every $G\in\mathcal{G}\setminus\{F\}$ one of the following cases holds: either $F\succ_{\mathcal{D}}G$ or $F\sim_{\mathcal{F}}G$, meaning, in the latter case it holds simultaneously $F\succ_{\mathcal{F}}G$ and $G\succ_{\mathcal{F}}F$. Accordingly, $F$ is maximal with respect to $\succ_{\mathcal{D}}$. *

Proof.

Fix an arbitrary $G\in\mathcal{G}\setminus\{F\}$ and assume that the first claim $F\succ_{\mathcal{D}}G$ is not valid. That is, for every $\delta>0$ there exists $t\in(0,\delta)$ such that $F(t+)<G(t+)$ This directly implies $G\succ_{\mathcal{F}}F$. So it remains to show the converse statement, i.e. $F\succ_{\mathcal{F}}G$.

Case 1: If there is a $\delta>0$ such that $F(t+)<G(t+)$ for all $t\in(0,\delta)$ then we write

[TABLE]

[TABLE]

[TABLE]

where the right hand side is strictly positive provided $n$ was chosen sufficiently large. This is a contradiction with the maximality of the corresponding Laplace transforms of $F$.

Case 2: If that is not the case, then the functions $F,G$ oscillate around each other in the sense that for every $\delta>0$ there exists $t\in(0,\delta)$ such that $F(t+)\geqslant G(t+)$. Hence $F\succ_{\mathcal{F}}G$ and the second claim, namely, $F\sim_{\mathcal{F}}G$, is valid.

Finally, as already observed in Section 1.1, maximality with respect to $\succ_{\mathcal{F}}$ implies maximality with respect to $\succ_{\mathcal{D}}$. Thus, $F$ is maximal with respect to $\succ_{\mathcal{D}}$ and the proof is complete. ∎

In other words, if for some initial data the Euler system (1.1)–(1.4) possesses a solution whose total entropy has maximal Laplace transforms evaluated at a sequence $\lambda_{n}\to\infty$, then the Dafermos’ criterium is satisfied. Note that since we have a semigroup it is enough to test the criterium at the time $t=0$.

We now argue that such a solution can be always chosen by our selection process considering a suitable order of minimizers in the procedure described in [2, Section 5.1]. The latter one considers bounded continuous functionals $\beta:\mathcal{S}\rightarrow R$ on the phase space. The final selection $\mathbb{U}$ maximizes (or minimizes) $\beta(\mathbb{U})$ pointwise in time by selecting the maximizer (or minimizer) of the functionals

[TABLE]

from the solution set. In [2, Section 5.1], the sequence $(\lambda_{n})_{n}$ is chosen to be dense in $(0,\infty)$. However, the density of $(\lambda_{n})_{n}$ is only needed to apply Lerch’s theorem implying the uniqueness of the Laplace transform. In fact, it is sufficient that the family of functionals $I_{n}$ in (5.3) separates points, which can be achieved by choosing a suitable increasing sequence $\lambda_{n}\to\infty$. This follows from the following generalization of classical Lerch’s theorem. Consequently, it is enough that the total entropy of our selection maximizes (5.2) for $\lambda=\lambda_{n}$, $n\in\mathbb{N}$.

Lemma 5.3**.**

Let $\lambda>0$ and $(\zeta_{n})_{n}\subset(0,\infty)$ be a sequence with the following property: for all $n,m\in\mathbb{N}$ there exists $k\in\mathbb{N}$ such that $\zeta_{n}+\zeta_{m}=\zeta_{k}$. Then the set of functionals

[TABLE]

separates points.

Proof.

Due to the condition on the sequence $(\zeta_{n})_{n}$, the set of finite linear combinations of functions $(e^{-\zeta_{n}t})_{n}$ forms a subalgebra of $C_{0}([0,\infty))$ which separates points and vanishes nowhere. Hence by the locally compact version of the Stone–Weierstrass theorem it is dense in $C_{0}([0,\infty))$ and consequently dense in $L^{1}(0,\infty)$. The claim now follows since $L^{1}(0,\infty)$ is the predual of $L^{\infty}(0,\infty)$: let $F$ be such that all the functionals above are zero. Then for every $g\in L^{1}(0,\infty)$

[TABLE]

and consequently $e^{-\lambda t}F(t)=0$ for a.e. $t\in(0,\infty)$ hence $F=0$. ∎

From the above Lemma we see that a minimal sequence $(\lambda_{n})_{n}$ needed for Lerch’s theorem is an arithmetic sequence of the form $\lambda_{n}=\lambda+n\zeta$ where $\lambda,\zeta>0$. For completeness we note that the converse implication regarding the Dafermos’ criterium is not valid. More precisely, if a function satisfies the Dafermos’ criterium (1.8) then it will not necessarily be chosen by our selection. This observation is based on the following result.

Lemma 5.4**.**

Let $\mathcal{G}\subset L^{\infty}(0,\infty)\cap\mathrm{BV}_{\rm loc}([0,\infty))$ and assume that $F\in\mathcal{G}$ satisfies $F\succ_{\mathcal{D}}G$ for every $G\in\mathcal{G}$. Then there exists $\lambda_{0}=\lambda_{0}(F,G)>0$ such that for all $\lambda\geqslant\lambda_{0}$ it holds

[TABLE]

Proof.

We write

[TABLE]

[TABLE]

Hence there exists $\lambda_{0}>0$ (depending on $F,G$) such that the right hand side is strictly positive provided $\lambda\geqslant\lambda_{0}$. ∎

In other words, one would like to start the selection procedure given by an increasing sequence $(\lambda_{n})_{n}$ with $\lambda_{1}$ as large as possible, in order to guarantee that the solution satisfying the Dafermos’ criterium is selected. However, a uniform choice of $\lambda_{1}$ is not a priori known at the moment.