Abstract
We consider a class of “wild” initial data to the compressible Euler system that give rise to infinitely many admissible weak solutions via the method of convex integration. We identify the closure of this class in the natural \(L^1\)-topology and show that its complement is rather large, specifically it is an open dense set.
Similar content being viewed by others
1 Introduction
The motion of a fluid in gas dynamics is described by means of the standard variables: the fluid density \(\varrho = \varrho (t,x)\), the macroscopic velocity \({{\mathbf {u}}}= {{\mathbf {u}}}(t,x)\), and the absolute temperature \(\vartheta = \vartheta (t,x)\) satisfying the Euler system:
We suppose the fluid is confined to a smooth bounded domain \(\varOmega \subset R^N\), \(N=2,3\), with the impermeable boundary:
The state of the system at the reference time \(t = 0\) is given by the initial conditions:
1.1 Thermodynamics
The system (1) contains two thermodynamics functions: the pressure \(p = p(\varrho , \vartheta )\) and the internal energy \(e = e(\varrho , \vartheta )\). In accordance with the Second law, we postulate the existence of the entropy \(s = s(\varrho , \vartheta )\) related to p and e via Gibbs’ equation
If the functions \(\varrho \), \({{\mathbf {u}}}\), and \(\vartheta \) are continuously differentiable, the relations (1) give rise to the entropy balance:
In the context of weak solutions considered in this paper, it is customary to relax (4) to inequality
The Euler system (1) can be written in the conservative variables: the density \(\varrho \), the momentum \({{\mathbf {m}}} = \varrho {{\mathbf {u}}}\), and the energy \(E = \frac{1}{2} \varrho |{{\mathbf {u}}}|^2 + \varrho e(\varrho , \vartheta )\):
Although we divide by \(\varrho \) when passing to the formulation in conservative variables, we can still treat vacuum by requiring the momentum to be zero in vacuum zones, see Sect. 2.
The forthcoming analysis leans essentially on the hypothesis of thermodynamics stability. The latter can be formulated either in the standard variables:
or in the conservative variables:
1.2 Smooth and weak solutions
Smooth solutions are at least continuously differentiable on the set \([0,T) \times \overline{\varOmega }\), \(\varrho > 0\), \(\vartheta > 0\), and satisfy (1–3) pointwise. Smooth solutions are known to exist locally on a maximal interval \([0,T_\mathrm{max})\) as long as the initial data are smooth enough and satisfy the corresponding compatibility conditions on \(\partial \varOmega \), see last part of Sect. 4 for details. Moreover, it is well known that \(T_\mathrm{max} < \infty \) for a “generic” class of data.
The weak solutions satisfy (1), (2) in the sense of distributions, the conservative variables are weakly continuous in time so that the initial conditions are well defined. The weak solutions are called admissible if the entropy inequality (5) holds in the weak sense. The existence of global in time weak (admissible) solutions for given initial data is an open problem. Recently, however, the method of convex integration developed in the context of the incompressible Euler system by De Lellis and Székelyhidi [7] has been adapted to identify a class of initial data for which the problem (1–3) admits infinitely many admissible weak solutions defined on a given time interval (0, T), [9]. Similar results have been obtained also for the associated Riemann problem in [1].
The initial data for which the problem admits local smooth solution will be termed regular, the data giving rise to infinitely many admissible weak solutions are called wild. Our goal is to identify the class of regular data that can be obtained as limits of wild data. Note that for the incompressible Euler system, the wild data (velocities) are dense in the Lebesgue space \(L^2(\varOmega ; R^N)\), see Székelyhidi and Wiedemann [13].
In the present paper, we restrict ourselves to the class of wild solutions resulting from a rather general splitting method which goes back to [8]. To the best of our knowledge, all convex integration results for the Euler system (6) in the literature use such a splitting method.
A vector field \({{\mathbf {m}}}\) can be decomposed by means of Helmholtz decomposition, i. e. it can be written as the sum of a solenoidal vector field \({{\mathbf {H}}}[{{\mathbf {m}}}]\) and the gradient of a scalar field.
More precisely, the Helmholtz projection operator \({{\mathbf {H}}}\) is defined as
Accordingly, for \({{\mathbf {m}}} = {{\mathbf {v}}} + \nabla _x\varPhi \), \({{\mathbf {v}}} = {{\mathbf {H}}} [{{\mathbf {m}}}]\), the system (6) can be written in the form
To the best of our knowledge, there are only two convex integration ansatzes for obtaining wild solutions to the Euler system (6) in the literature. Let us begin with the first one, which is used in [9]. This ansatz is based on replacing Eqs. (9)\(_2\) and (9)\(_3\) by the system
Equation (10)\(_1\) represents the volume preserving part of the motion, while (10)\(_2\) may be seen as a wave equation governing the propagation of acoustic waves. Note that the fields \({{\mathbf {v}}}\) and \(\nabla _x\varPhi \) as well as the fluxes in (10) are orthogonal with respect to the \(L^2\) scalar product, specifically,
where the latter is true since it is a product of a traceless matrix and a multiple of the identity matrix.
Although (10)\(_1\) is apparently overdetermined, it admits (infinitely many) weak solutions for any fixed \(\varrho \) and \(\varPhi \), cf. “Appendix A”.
Definition 1.1
A weak solution \([\varrho , {{\mathbf {m}}}, E]\) of the Euler system (6) is called wild solution if \([\varrho , {{\mathbf {m}}}, E]\) satisfy (10), with \({{\mathbf {v}}} = {{\mathbf {H}}}[{{\mathbf {m}}}]\), \(\nabla _x\varPhi = {{\mathbf {H}}}^\perp [{{\mathbf {m}}}]\). The corresponding initial data \([\varrho _0, {{\mathbf {m}}}_0, E_0]\) are called wild initial data.
As indicated above, Theorem A.1 yields solutions which are wild in the sense of Definition 1.1. This guarantees that the set of wild solutions, as well as the set of wild initial data are non-empty. However Definition 1.1 is more general, which means that there might be wild solutions (and hence wild initial data) in the sense of Definition 1.1 that cannot be obtained by Theorem A.1.
Note that a technique similar to (10) has been used also for the simplified isentropic Euler system by Chiodaroli [4].
The second convex integration ansatz, that is available in the literature, is based on the analysis of the corresponding Riemann problem, see [1, 6, 10] among others. We will first focus on the ansatz (10) and afterwards extend our results to the wild solutions obtained via the Riemann problem in Sect. 5.
1.3 Main result
We are ready to formulate our main result. To avoid the situation when the temperature approaches absolute zero, we restrict ourselves to the phase space
Note that this is not very restrictive as any admissible weak solution satisfies the minimum principle
cf. [2, Section 2.1.1].
We say that a sequence of data \(\left\{ \varrho _{0,n}, {{\mathbf {u}}}_{0,n}, \vartheta _{0,n} \right\} _{n=1}^\infty \), or, equivalently,
\(\left\{ \varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n} \right\} _{n=1}^\infty \), (W)-converges to \([\varrho _0, {{\mathbf {m}}}_0, E_0]\),
if
-
$$\begin{aligned} \varrho _{0,n}> 0,\quad \ s(\varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n}) \ge s_0 > -\infty ; \end{aligned}$$
-
$$\begin{aligned}{}[ \varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n} ] \rightarrow [\varrho _0, {{\mathbf {m}}}_0, E_0] \quad \ \hbox {in}\ L^1(\varOmega ; R^{N+2}); \end{aligned}$$(11)
-
the initial data \([\varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n}]\) give rise to a sequence of admissible weak solutions \([\varrho _n, {{\mathbf {m}}}_n, E_n]\) satisfying
$$\begin{aligned} \begin{aligned} \int _0^T \int _{\varOmega } \left( \frac{{{\mathbf {m}}}_n \otimes {{\mathbf {m}}}_n}{\varrho _n} - \frac{1}{N} \frac{ |{{\mathbf {m}}}_n|^2}{\varrho _n} {\mathbb {I}} \right) : \nabla _x^2 \varphi \ {\mathrm{d}} {x} \ {\mathrm{d}} t \rightarrow 0 \quad \ \hbox {as}\ n \rightarrow \infty \\ \quad \hbox {for any}\ \varphi \in C^\infty _c((0,T) \times \varOmega ). \end{aligned} \end{aligned}$$(12)
Note that, in view of (10),
whenever \([\varrho _n, {{\mathbf {m}}}_n, E_n]\) are wild solutions in the sense of Definition 1.1. In particular, (12) is trivially satisfied for any sequence of wild initial data in the sense of Definition 1.1.
Finally, note that the sequence of solutions \(\left\{ \varrho _{n}, {{\mathbf {m}}}_{n}, E_{n} \right\} _{n=1}^\infty \) need not a priori converge strongly, and, consequently, (12) does not impose any similar restriction on a possible limit solution emanating from \([\varrho _0, {{\mathbf {m}}}_0, E_0]\).
Definition 1.2
We say that a trio \([\varrho _0, {{\mathbf {m}}}_0, E_0]\) is reachable if there exists a sequence of initial data \(\left\{ \varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n} \right\} _{n=1}^\infty \) such that
Our goal is to show that the set of reachable data is “small” in the sense that its complement is an open dense set in \(\texttt {L}^1_{+,s_0} (\varOmega ; R^{N+2})\). For the sake of simplicity, we focus on the equations of gas dynamics, where the pressure is given by Boyle–Mariotte law
Equivalently, in the conservative variables,
Here the kinetic energy is defined as
The entropy reads
Here is our main result.
Theorem 1.3
Let \(s_0 \in R\) be given. Let \(\varOmega \subset R^N\), \(N=2,3\) be a bounded smooth domain.
Then the complement of the set of reachable data is an open dense set in \(\texttt {L}^1_{+,s_0} (\varOmega ; R^{N + 2})\).
Corollary 1.4
The complement of the set of wild initial data (in the sense of Definition 1.1) contains an open dense set in \(\texttt {L}^1_{+,s_0} (\varOmega ; R^{N + 2})\).
Proof
As mentioned above, wild initial data are reachable and hence the complement of the set of wild initial data contains the complement of the set of reachable data. Thus Theorem 1.3 yields the claim. \(\square \)
The main part of the paper is devoted to the proof of Theorem 1.3. Our strategy is to identify a large (dense) set of regular data for the Euler system that are not reachable. To see this, we consider initial data \([\varrho _0, {{\mathbf {m}}}_0, E_0]\) giving rise to a smooth solution \([\varrho , {{\mathbf {m}}}, E]\) defined on a maximal time interval \([0, T_{\mathrm{max}})\). Then we proceed in several steps:
-
Assuming the data \([\varrho _0, {{\mathbf {m}}}_0, E_0]\) are reachable we show that the associated sequence of solutions \(\{ \varrho _n, {{\mathbf {m}}}_n, E_n \}_{n=1}^\infty \) generates a Young measure that can be identified with a dissipative measure valued (DMV) solution to the Euler system in the sense of [2], see Sect. 2.
-
Next we realize that, thanks to the strong convergence required in (11), the limit DMV solution starts from the initial data represented by the parameterized family of Dirac masses,
$$\begin{aligned} \left\{ \delta _{\varrho _0(x),{{\mathbf {m}}}_0(x), E_0(x) } \right\} _{x \in \varOmega }. \end{aligned}$$Thanks to the general weak–strong uniqueness principle, we conclude that the DMV solution coincides with the strong solution on the time interval \([0, T_\mathrm{max})\), more specifically, the DMV solution is represented by the parameterized family of Dirac masses,
$$\begin{aligned} \left\{ \delta _{\varrho (t,x),{{\mathbf {m}}}(t,x), E(t,x) } \right\} _{(t,x) \in (0, T_{\mathrm{max}}) \times \varOmega }, \end{aligned}$$see Sect. 3. In particular, the solutions \([\varrho _n, {{\mathbf {m}}}_n, E_n]\) converge strongly and we may assume
$$\begin{aligned} \varrho _n \rightarrow \varrho , \ {{\mathbf {m}}}_n \rightarrow {{\mathbf {m}}},\ E_n \rightarrow E \ \quad \hbox {a.a. in}\ (0,T_{\mathrm{max}}) \times \varOmega . \end{aligned}$$(14) -
Thanks to the strong convergence (14), condition (12), and smoothness of the limit solution, we deduce that
$$\begin{aligned} \mathrm{div}_x\mathrm{div}_x\left( \frac{ {{\mathbf {m}}} \otimes {{\mathbf {m}}} }{\varrho } \right) - \frac{1}{N} \varDelta _x\left( \frac{|{{\mathbf {m}}}|^2}{\varrho } \right) = 0 \ \quad \hbox {in}\ (0,T) \times \varOmega . \end{aligned}$$see Sect. 3 for details. In particular, as the limit solution is continuous up to \(t= 0\),
$$\begin{aligned} \mathrm{div}_x\mathrm{div}_x\left( \frac{ {{\mathbf {m}}}_0 \otimes {{\mathbf {m}}}_0 }{\varrho _0} \right) - \frac{1}{N} \varDelta _x\left( \frac{|{{\mathbf {m}}}_0|^2}{\varrho _0} \right) = 0 \ \quad \hbox {in}\ \varOmega . \end{aligned}$$(15) -
In Sect. 4, we show that (15) can be satisfied for a very narrow class of the initial data. In particular, we complete the proof of Theorem 1.3.
A related result for the wild solutions obtained via the Riemann problem is shown in Sect. 5.
2 Dissipative measure valued (DMV) solutions
2.1 Definitions
First we recall the standard definition of the admissible weak solutions to the Euler system:
Definition 2.1
(Admissible weak solution) Let \(\gamma = 1 + \frac{1}{c_v}\). We say that \([\varrho , {{\mathbf {m}}}, E]\) is an admissible weak solution to the Euler system (1–3) if
and the following holds:
-
$$\begin{aligned} \int _0^T \int _{\varOmega } \left[ \varrho \partial _t \varphi + {{\mathbf {m}}}\cdot \nabla _x\varphi \right] \ {\mathrm{d}} {x} \ {\mathrm{d}} t = - \int _{\varOmega } \varrho _0 \varphi (0, \cdot ) \ {\mathrm{d}} {x} \end{aligned}$$
for any \(\varphi \in C^1_c([0,T) \times \overline{\varOmega })\);
-
$$\begin{aligned} \begin{aligned}&\int _0^T \int _{\varOmega } \left[ {{\mathbf {m}}}\cdot \partial _t\varvec{\varphi }+ \left( \frac{{{\mathbf {m}}}\otimes {{\mathbf {m}}}}{\varrho } \right) : \nabla _x\varvec{\varphi }+ p \left( \varrho , {{\mathbf {m}}}, E \right) \mathrm{div}_x\varvec{\varphi }\right] \ {\mathrm{d}} {x} \ {\mathrm{d}} t \\&\quad = - \int _{\varOmega } {{\mathbf {m}}}_0 \cdot \varvec{\varphi }(0, \cdot ) \ {\mathrm{d}} {x} \end{aligned} \end{aligned}$$
for any \(\varvec{\varphi }\in C^1_c([0,T) \times \overline{\varOmega }; R^N)\), \(\varvec{\varphi }\cdot {{\mathbf {n}}}|_{\partial \varOmega } = 0\);
-
$$\begin{aligned} \int _0^T \int _{\varOmega } \left[ E \partial _t \varphi + \left( E + p(\varrho , {{\mathbf {m}}}, E) \right) \left( \frac{{{\mathbf {m}}}}{\varrho } \right) \cdot \nabla _x\varphi \right] \ {\mathrm{d}} {x} \ {\mathrm{d}} t = - \int _{\varOmega } E_0 \varphi (0, \cdot ) \ {\mathrm{d}} {x} \end{aligned}$$(16)
for any \(\varphi \in C^1_c ([0,T) \times \overline{\varOmega })\);
-
$$\begin{aligned} \begin{aligned} -&\int _{\varOmega } \varrho _0 Z \left( s \left( \varrho _0, {{\mathbf {m}}}_0, E_0 \right) \right) \varphi (0, \cdot ) \ {\mathrm{d}} {x} \\&\ge \int _0^T \int _{\varOmega } \left[ \varrho Z \left( s \left( \varrho , {{\mathbf {m}}}, E \right) \right) \partial _t \varphi + Z \left( s \left( \varrho , {{\mathbf {m}}}, E \right) \right) {{\mathbf {m}}} \cdot \nabla _x\varphi \right] \ {\mathrm{d}} {x} \ {\mathrm{d}} t \end{aligned} \end{aligned}$$(17)
for any \(\varphi \in C^1_c([0,T) \times \overline{\varOmega })\), \(\varphi \ge 0\), and \(Z(s) \le \overline{Z}\), \(Z'(s) \ge 0\).
Remark 2.2
Here the entropy inequality is satisfied in the renormalized sense, see e.g. [3].
To carry out the programme delineated at the end of the preceding section, we introduce the concept of dissipative measure valued (DMV) solution, see [2]. Let
be the phase space associated to the Euler system.
Definition 2.3
(DMV solution) A parametrized family of probability measures \(\{ {\mathcal {V}}_{t,x} \}_{(t,x) \in (0,T) \times \varOmega }\) on the set \({\mathcal {P}}\),
is called a dissipative measure–valued (DMV) solution to the compressible Euler system with the initial data \(\{ {\mathcal {V}}_{0,x} \}_{x \in \varOmega }\), if the following holds:
-
$$\begin{aligned} \int _0^T \int _{\varOmega } \left[ \left\langle {\mathcal {V}}_{t,x}; \varrho \right\rangle \partial _t \varphi + \left\langle {\mathcal {V}}_{t,x}; {{\mathbf {m}}}\right\rangle \cdot \nabla _x\varphi \right] \ {\mathrm{d}} {x} \ {\mathrm{d}} t = - \int _{\varOmega } \left\langle {\mathcal {V}}_{0,x}; \varrho \right\rangle \varphi (0, \cdot ) \ {\mathrm{d}} {x} \end{aligned}$$
for any \(\varphi \in C^1_c([0,T) \times \overline{\varOmega })\);
-
$$\begin{aligned} \begin{aligned}&\int _0^T \int _{\varOmega } \bigg [ \left\langle {\mathcal {V}}_{t,x}; {{\mathbf {m}}}\right\rangle \cdot \partial _t\varvec{\varphi }+ \left\langle {\mathcal {V}}_{t,x} ; \frac{ {{\mathbf {m}}}\otimes {{\mathbf {m}}}}{\varrho } \right\rangle : \nabla _x\varvec{\varphi }\\&\qquad + \left\langle {\mathcal {V}}_{t,x}; p \left( \varrho , {{\mathbf {m}}}, E \right) \right\rangle \mathrm{div}_x\varvec{\varphi }\bigg ] \ {\mathrm{d}} {x} \ {\mathrm{d}} t \\&\quad = \int _0^T \int _{\varOmega } \nabla _x\varvec{\varphi }: \mathrm{d}\mu _C - \int _{\varOmega } \left\langle {\mathcal {V}}_{0,x}; {{\mathbf {m}}} \right\rangle \cdot \varvec{\varphi }(0, \cdot ) \ {\mathrm{d}} {x} \end{aligned} \end{aligned}$$
for any \(\varvec{\varphi }\in C^1_c([0,T) \times \overline{\varOmega }; R^N)\), \(\varvec{\varphi }\cdot {{\mathbf {n}}}|_{\partial \varOmega } = 0\), where \(\mu _C\) is a (vectorial) signed measure on \([0,T] \times \overline{\varOmega }\);
-
$$\begin{aligned} \int _{\varOmega } \left\langle {\mathcal {V}}_{\tau ,x}; E \right\rangle \ {\mathrm{d}} {x} + {\mathcal {D}}(\tau ) \le \int _{\varOmega } \left\langle {\mathcal {V}}_{0,x}; E \right\rangle \ {\mathrm{d}} {x} \end{aligned}$$(18)
for a.a. \(\tau \in (0,T)\), where \({\mathcal {D}} \in L^\infty (0,T)\), \({\mathcal {D}} \ge 0\) is called dissipation defect;
-
$$\begin{aligned} \begin{aligned}&- \int _{\varOmega } \left\langle {\mathcal {V}}_{0,x}; \varrho Z \left( s \left( \varrho , {{\mathbf {m}}}, E \right) \right) \right\rangle \varphi (0, \cdot ) \ {\mathrm{d}} {x} \\&\quad \ge \int _0^T \int _{\varOmega } \left[ \left\langle {\mathcal {V}}_{t,x} ; \varrho Z \left( s \left( \varrho , {{\mathbf {m}}}, E \right) \right) \right\rangle \right] \partial _t \varphi \ {\mathrm{d}} {x} \ {\mathrm{d}} t \\&\qquad + \int _0^T \int _{\varOmega } \left[ \left\langle {\mathcal {V}}_{t,x} ; Z \left( s \left( \varrho , {{\mathbf {m}}}, E \right) \right) {{\mathbf {m}}} \right\rangle \cdot \nabla _x\varphi \right] \ {\mathrm{d}} {x} \ {\mathrm{d}} t \end{aligned} \end{aligned}$$
for any \(\varphi \in C^1_c([0,T) \times \overline{\varOmega })\), \(\varphi \ge 0\), and \(Z(s) \le \overline{Z}\), \(Z'(s) \ge 0\);
-
$$\begin{aligned} \left\| \mu _C \right\| _{{\mathcal {M}}([0, \tau ] \times \varOmega ) } \le c \int _0^\tau {\mathcal {D}}(t) \ {\mathrm{d}} t \ \quad \hbox {for a.a.}\ \tau \in (0,T). \end{aligned}$$
The reader will have noticed that Definition 2.3 is not a mere measure–valued variant of Definition 2.1. In particular, the energy equation (16) has been replaced by its integrated version (18). The reader may consult [2] for a thorough discussion of this new concept of solution. Note in particular, that (18) implies
see [2, Remark 2.7]. Hence DMV solutions allow for treating vacuum.
2.2 Generating a DMV solution
Suppose now that
where \([\varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n}]\) are initial data of admissible weak solutions \([\varrho _n, {{\mathbf {m}}}_n, E_n]\) defined on \((0,T) \times \varOmega \).
Passing to a suitable subsequence, we may suppose that \(\{ \varrho _n, {{\mathbf {m}}}_n, E_n \}_{n=1}^\infty \) generates a Young measure
As shown in [2, Section 2.1], the Young measure \(\left\{ {\mathcal {V}}_{t,x} \right\} _{(t,x) \in (0,T) \times \varOmega )}\) is a DMV solution of the Euler system in the sense of Definition 2.3. Moreover, in view of (19), the initial measure \(\{ {\mathcal {V}}_{0,x} \}_{x \in \varOmega }\) reads
3 Application of the weak–strong uniqueness principle
The weak-strong uniqueness principle (Theorem 3.3 in [2]) asserts that a DMV solution coincides with the strong solution emanating from the same initial data at least on the life span of the latter. Evoking the situation from Sect. 2.2, we suppose now that
where \([\varrho _0, {{\mathbf {m}}}_0, E_0]\) are now regular initial data generating a \(C^1-\)solution \([\varrho , {{\mathbf {m}}}, E]\) of the Euler system (1–3) on a maximal time interval \([0, T_{\mathrm{max}})\). Without loss of generality, we may suppose that \(0< T < T_{\mathrm{max}}\).
Applying the weak strong uniqueness principle we obtain
In particular, we may assume
for the associated weak solutions \([\varrho _n, {{\mathbf {m}}}_n, E_n]\). Consequently, as \([\varrho _n, {{\mathbf {m}}}_n, E_n]\) satisfy (12), we get
Finally, as the limit solution is continuous up to the time \(t = 0\), we may infer that in particular,
For \(C^2\) initial data, relation (21) can be rewritten as a non–linear differential equation
4 Reachability
In this section we finish the proof of Theorem 1.3. Introducing \({{\mathbf {w}}} = \frac{ {{\mathbf {m}}}_0 }{\sqrt{\varrho _0}}\) we may write (21) in a concise form
We consider solutions of (22) in the space \(L^2(\varOmega ; R^N)\),
It is easy to check that:
-
the set \({\mathcal {S}}\) is closed in \(L^2(\varOmega ; R^N)\);
-
if \({{\mathbf {w}}} \in {\mathcal {S}}\), then \(\lambda {{\mathbf {w}}} \in {\mathcal {S}}\) for any \(\lambda \in R\).
Lemma 4.1
The set \({\mathcal {S}}\) is nowhere dense in \(L^2(\varOmega ; R^N)\), meaning the (closure of the) set \({\mathcal {S}}\) does not contain any ball in \(L^2(\varOmega ; R^N)\).
Proof
Arguing by contradiction, we suppose that there is \({{\mathbf {w}}}_0 \in L^2(\varOmega ; R^N)\) and \(r > 0\) such that the ball centered at \({{\mathbf {w}}}_0\) with radius r is contained in \({\mathcal {S}}\). Hence
Next, we show that this implies that the ball centered at 0 with radius r is contained in \({\mathcal {S}}\), too. To this end, we write
Thus using the fact that \({{\mathbf {w}}}_0\) solves (22) and (23) we conclude
If relation (24) holds for any \({{\mathbf {w}}}\), \(\Vert {{\mathbf {w}}} \Vert _{L^2(\varOmega ; R^N)} \le r\), it must hold for \(\lambda {{\mathbf {w}}}\), \(0 \le \lambda \le 1\). Consequently, we get from (24),
We deduce
Thus the set \({\mathcal {S}}\) contains the ball of radius \(r > 0\) centered at 0 but as it is invariant with respect to multiplication by any real constant, we conclude
which is obviously false as soon as \(N \ge 2\). \(\square \)
Lemma 4.2
The set of reachable data in the sense of Definition 1.2 is closed in \(L^1((0,T) \times \varOmega )\).
Proof
Let \([\varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n}]\) be a sequence of reachable data such that
To see that the limit is reachable, we have to find a generating sequence satisfying (11), (12). To this end, consider the generating sequences of data \([\varrho _{0,n,m}, {{\mathbf {m}}}_{0,n,m}, E_{0,n,m}]\) satisfying
meaning, in particular,
for the corresponding solutions \(\varrho _{n,m}\), \({{\mathbf {m}}}_{n,m}\) and any fixed n. Now, let
be a countable dense set in the Sobolev space
By the diagonal method, we can find a subsequence
such that
Indeed, for given \(j \ge 1\), we find n(j) such that
and \(m = m(n(j))\) such that
Finally, by virtue of (25), any \(\nabla ^2_x \varphi \) can be uniformly approximated by \(\nabla ^2_x \varphi ^k\), and, as the energy \(E_{n(j), m(j)}\) are bounded, we conclude that (26) holds for any \(\varphi \in C^\infty _c((0,T) \times \varOmega )\). Thus we have shown that
whence \([\varrho _0, {{\mathbf {m}}}_0, E_0]\) is reachable.
\(\square \)
We are ready to complete the proof of Theorem 1.3. The set of reachable data being closed (cf. Lemma 4.2), its complement is open. In view of Lemma 4.1, we have only to show that the set of the smooth initial data \([\varrho _0, {{\mathbf {m}}}_0, E_0]\) giving rise to local-in-time regular solutions is dense in the space \(\texttt {L}^1_{+, s_0}(\varOmega ; R^{N+2})\). To this end, we record the following result by Schochet [12, Theorem 1] that asserts the existence of smooth solutions whenever:
-
\(\varOmega \subset R^N\) is a bounded domain with a sufficiently smooth boundary, say \(\partial \varOmega \) of class \(C^\infty \);
-
the initial data \([\varrho _{0,E}, \vartheta _{0,E}, {{\mathbf {u}}}_{0,E}]\) belong to the class
$$\begin{aligned} \varrho _{0, E}, \vartheta _{0, E} \in W^{m,2}(\varOmega ), \ \quad {{\mathbf {u}}}_{0, E} \in W^{m,2}(\varOmega ; R^N),\quad \ \varrho _{0, E}, \ \vartheta _{0,E} > 0 \ \hbox {in}\ \overline{\varOmega }, \end{aligned}$$where \( m > N\);
-
the compatibility conditions
$$\begin{aligned} \partial ^k_t {{\mathbf {u}}}_{0,E} \cdot {{\mathbf {n}}}|_{\partial \varOmega } = 0 \end{aligned}$$hold for \(k=0,1,\dots ,m\).
In particular, any trio of initial data,
gives rise to a smooth solution defined on a maximal time interval \([0, T_{\mathrm{max}})\). Seeing that the set (27) is dense in \(\texttt {L}^1_{+, s_0}(\varOmega ; R^{N+2})\), we have completed the proof of Theorem 1.3.
Remark 4.3
Note that the result is basically independent of the choice of boundary conditions. Similarly, it can be easily extended to more general equation of state satisfying (8), including the isentropic case. Of course, the wild solution class is restricted only to those solutions that can be obtained by the splitting of the incompressible and acoustic part as specified in (10).
5 Reachability for the Riemann problem
To the best of our knowledge there are only two methods available in the literature that give rise to infinitely many solutions of the compressible Euler system, which are both based on a splitting as in (9). The first method [9] was explained above, see (10). The second method considers initial data constant in each of the two half spaces (Riemann data), see [1] and also [5, 6, 10] for the isentropic case. Here—instead of (10)—Eqs. (9)\(_2\) and (9)\(_3\) are replaced by
with a suitable matrix field \({\mathbb {U}}\). The density \(\varrho \), the acoustic potential \(\nabla _x\varPhi \) as well as the field \({\mathbb {U}}\) are piecewise constant in \((0,T) \times \varOmega \). In particular, there exists a partition of the set \((0,T) \times \varOmega \) into a finite number of sectors \(Q_i, i=1, \dots m,\) such that the wild solutions satisfy
Consequently, we may accommodate the data reachable by this method replacing (12) by
where \(Q \subset (0,T) \times \varOmega \) is an open set such that \(\overline{Q} = [0,T] \times \overline{\varOmega }\).
The above observation motivates the following extension of the concept of reachability. We say that Q is a partition of the domain \((0,T) \times \varOmega \) if
-
\(Q \subset (0,T) \times \varOmega \) is an open set;
-
\(\overline{Q} = [0,T] \times \overline{\varOmega }.\)
We say that a family of partitions \({\mathcal {Q}}\) is a closed partition set if it is closed with respect to the Hausdorff complementary topology. More specifically, any sequence \(\{ Q_n \}_{n=1}^\infty \subset {\mathcal {Q}}\) contains a subsequence such that
where the symbol \({\mathop {\rightarrow }\limits ^{(H)}}\) denotes convergence in the Hausdorff metric and \(Q^c_n\) denotes the complement \(Q^c_n \equiv ([0,T] \times \overline{\varOmega } ) \setminus Q_n\).
An example of a closed partition set related to the convex integration method is
where M is a given positive integer. Note that this is indeed a closed partition set since \((0,T)\times \varOmega \) is a bounded subset of \(R^{N+1}\).
The following property is well known, cf. [11]:
Next, we introduce the following generalization of \((W)-\)convergence. Let \({\mathcal {Q}}\) be a closed partition set. We say that a sequence of data \(\left\{ \varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n} \right\} _{n=1}^\infty \) \((W[{\mathcal {Q}}])\)-converges to \([\varrho _0, {{\mathbf {m}}}_0, E_0]\),
if
-
$$\begin{aligned} \varrho _{0,n}> 0,\quad \ s(\varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n}) \ge s_0 > -\infty ; \end{aligned}$$
-
$$\begin{aligned}{}[ \varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n} ] \rightarrow [\varrho _0, {{\mathbf {m}}}_0, E_0] \ \quad \hbox {in}\ L^1(\varOmega ; R^{N+2}); \end{aligned}$$
-
the initial data \([\varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n}]\) give rise to a sequence of admissible weak solutions \([\varrho _n, {{\mathbf {m}}}_n, E_n]\) satisfying
$$\begin{aligned} \int _0^T \int _{\varOmega } \left( \frac{{{\mathbf {m}}}_n \otimes {{\mathbf {m}}}_n}{\varrho _n} - \frac{1}{N} \frac{ |{{\mathbf {m}}}_n|^2}{\varrho _n} {\mathbb {I}} \right) : \nabla _x^2 \varphi \ {\mathrm{d}} {x} \ {\mathrm{d}} t \rightarrow 0 \ \quad \hbox {as}\ n \rightarrow \infty \ \hbox {for any}\ \varphi \in C^\infty _c(Q) \nonumber \\ \end{aligned}$$(29)for some \(Q \in {\mathcal {Q}}\).
Definition 5.1
Let \({\mathcal {Q}}\) be a closed partition set. We say that a trio \([\varrho _0, {{\mathbf {m}}}_0, E_0]\) is \({\mathcal {Q}}-\)reachable if there exists a sequence of initial data \(\left\{ \varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n} \right\} _{n=1}^\infty \) such that
Obviously any reachable data in the sense of Definition 1.2 are \({\mathcal {Q}}\)-reachable so the set of \({\mathcal {Q}}-\)reachable data is always larger for any closed partition set.
In order to adapt the arguments used in Sects. 3 and 4 we need the following result.
Lemma 5.2
Let \({\mathcal {Q}}\) be a closed partition set.
Then the set of \({\mathcal {Q}}-\)reachable data is closed in \(L^1((0,T) \times \varOmega ; R^{N+2})\).
Proof
Let \([\varrho _{0,n}, {{\mathbf {m}}}_{0,n}, E_{0,n}]\) be a sequence of \({\mathcal {Q}}-\)reachable data such that
Let \([\varrho _{0,n,m}, {{\mathbf {m}}}_{0,n,m}, E_{0,n,m}]\) be the corresponding generating sequences,
In particular, the corresponding solutions \(\varrho _{m,n}\), \({{\mathbf {m}}}_{m,n}\) satisfy
for some \(Q_n \in {\mathcal {Q}}\) and any fixed n.
As the partition set \({\mathcal {Q}}\) is closed, there exists a partition \(Q \in {\mathcal {Q}}\) such that
at least for a suitable subsequence (not relabeled).
Let
be a countable dense subset of the Sobolev space
Using property eq. (30), we can find a subsequence
satisfying
By virtue of (31), any \(\nabla ^2_x \varphi \) can be uniformly approximated by \(\nabla ^2_x \varphi ^k\), and, as the energies \(E_{n(j), m(j)}\) are bounded, we conclude that (32) holds for any \(\varphi \in C^\infty _c(Q)\). We may infer that
whence \([\varrho _0, {{\mathbf {m}}}_0, E_0]\) is \({\mathcal {Q}}-\)reachable. \(\square \)
Finally, we observe that any regular initial data \([\varrho _0, {{\mathbf {m}}}_0, E_0]\) satisfy (21). To see this, we first observe, similarly to Sect. 3, that
As \(\varrho \), \({{\mathbf {m}}}\) are smooth, we get
However, as \(\overline{Q} = [0,T] \times \overline{\varOmega }\) and the second derivatives of \(\rho \), \({{\mathbf {m}}}\) are continuous, this implies
In particular, we obtain (21) for the initial data.
We have shown the following extension of Theorem 1.3.
Theorem 5.3
Let \(s_0 \in R\) be given and \(\varOmega \subset R^N\), \(N=2,3\) be a bounded smooth domain. Let \({\mathcal {Q}}\) be a closed partition set in \((0,T) \times \varOmega \).
Then the complement of the set of \({\mathcal {Q}}-\)reachable data is an open dense set in \(\texttt {L}^1_{+,s_0} (\varOmega ; R^{N + 2})\).
As already pointed out, Theorem 5.3 accommodates the limits of the wild initial data constructed via the Riemann problem with M fans in the sense of [5].
References
Al Baba, H., Klingenberg, C., Kreml, O., Mácha, V., Markfelder, S.: Nonuniqueness of admissible weak solutions to the Riemann problem for the full Euler system in two dimensions. SIAM J. Math. Anal. 52(2), 1729–1760 (2020)
Březina, J., Feireisl, E.: Measure-valued solutions to the complete Euler system. J. Math. Soc. Jpn. 70(4), 1227–1245 (2018)
Chen, G.-Q., Frid, H., Li, Y.: Uniqueness and stability of Riemann solutions with large oscillation in gas dynamics. Commun. Math. Phys. 228(2), 201–217 (2002)
Chiodaroli, E.: A counterexample to well-posedness of entropy solutions to the compressible Euler system. J. Hyperbolic Differ. Equ. 11(3), 493–519 (2014)
Chiodaroli, E., De Lellis, C., Kreml, O.: Global ill-posedness of the isentropic system of gas dynamics. Commun. Pure Appl. Math. 68(7), 1157–1190 (2015)
Chiodaroli, E., Kreml, O.: On the energy dissipation rate of solutions to the compressible isentropic Euler system. Arch. Ration. Mech. Anal. 214(3), 1019–1049 (2014)
De Lellis, C., Székelyhidi Jr., L.: On admissibility criteria for weak solutions of the Euler equations. Arch. Ration. Mech. Anal. 195(1), 225–260 (2010)
Feireisl, E.: Weak solutions to problems involving inviscid fluids. In: Mathematical Fluid Dynamics, Present and Future, Springer Proceedings in Mathematics and Statistics, vol. 183, Springer, Tokyo, pp. 377–399 (2016)
Feireisl, E., Klingenberg, C., Kreml, O., Markfelder, S.: On oscillatory solutions to the complete Euler system. J. Differ. Equ. 269(2), 1521–1543 (2020)
Klingenberg, C., Markfelder, S.: The Riemann problem for the multidimensional isentropic system of gas dynamics is ill-posed if it contains a shock. Arch. Ration. Mech. Anal. 227(3), 967–994 (2018)
Pironneau, O.: Optimal shape design for elliptic systems. In: Springer Series in Computational Physics, Springer, New York (1984)
Schochet, S.: The compressible Euler equations in a bounded domain: existence of solutions and the incompressible limit. Commun. Math. Phys. 104(1), 49–75 (1986)
Székelyhidi Jr., L., Wiedemann, E.: Young measures generated by ideal incompressible fluid flows. Arch. Ration. Mech. Anal. 206(1), 333–366 (2012)
Acknowledgements
Open Access funding provided by Projekt DEAL. 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.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by C.De Lellis.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Wild solutions
Wild solutions
The existence of weak solutions for the Euler system in higher space dimensions is an outstanding open problem. Nevertheless, the recently developed adaptation of the convex integration technique to the incompressible Euler system provided certain results also for the compressible case. We report the following result proved in [9].
Theorem A.1
Let the initial data \(\varrho _0\) and \(\vartheta _0\) be piecewise constant strictly positive functions defined on \(\varOmega \). And let \(T > 0\) be given.
Then there exists a constant \(\varLambda _0 > 0\) such that for any \(\varLambda \ge \varLambda _0\) there is \({{\mathbf {u}}}_0 \in L^\infty (\varOmega ; R^N)\) such that the problem (1), (5), (8) admits infinitely many weak solutions in \((0,T) \times \varOmega \), with the initial data
In addition, these solutions enjoy the following properties:
and the renormalized entropy inequality (17) holds as equality.
Note that here (10) is satisfied with \(\varPhi = 0\).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Feireisl, E., Klingenberg, C. & Markfelder, S. On the density of “wild” initial data for the compressible Euler system. Calc. Var. 59, 152 (2020). https://doi.org/10.1007/s00526-020-01806-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s00526-020-01806-5