On the density of wild initial data for the compressible Euler system

Transcription

1 arxiv: v1 [math.ap] 31 Dec 218 On the density of wild initial data for the compressible Euler system Eduard Feireisl Christian Klingenberg Simon Markfelder January 1, 219 Institute of Mathematics of the Academy of Sciences of the Czech Republic Žitná 25, CZ Praha 1, Czech Republic and Institute of Mathematics, TU Berlin Strasse des 17.Juni, Berlin Department of Mathematics, Würzburg University Emil-Fischer-Str. 4, 9774 Würzburg, Germany Abstract We consider a class of wild initial data to the compressibleeuler 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. Keywords: compressible Euler system, convex integration, wild data, weak strong uniqueness Contents 1 Introduction Thermodynamics Smooth and weak solutions Main result The researchof E.F. leading to these results has received funding from the Czech Sciences Foundation (GAČR), Grant Agreement S. The Institute of Mathematics of the Academy of Sciences of the Czech Republic is supported by RVO:

2 2 Dissipative measure valued (DMV) solutions Definitions Generating a DMV solution Application of the weak strong uniqueness principle 1 4 Reachability 11 5 Reachability for the Riemann problem 14 A Wild solutions 18 1 Introduction The motion of a fluid in gas dynamics is described by means of the standard variables: the fluid density = (t,x), the macroscopic velocity u = u(t,x), and the absolute temperature ϑ = ϑ(t,x) satisfying the Euler system: t+div x (u) =, t (u)+div x (u u)+ x p(,ϑ) =, ( ) [( ) ] 1 1 t 2 u 2 +e(,ϑ) +div x 2 u 2 +e(,ϑ) u +div x (p(,ϑ)u) =. (1.1) We suppose the fluid is confined to a smooth bounded domain R N, N = 2,3, with the impermeable boundary: u n =. (1.2) The state of the system at the reference time t = is given by the initial conditions: (, ) =, u(, ) = u, ϑ(, ) = ϑ. (1.3) 2

3 1.1 Thermodynamics The system (1.1) contains two thermodynamics functions: the presure p = p(, ϑ) and the internal energy e = e(,ϑ). In accordance with the Second law, we postulate the existence of the entropy s = s(,ϑ) related to p and e via Gibbs equation ( ) 1 ϑds = De+D p. If the functions, u, and ϑ are continuously differentiable, the relations (1.1) give rise to the entropy balance: t (s(,ϑ))+div x (s(,ϑ)u) =. (1.4) In the context of weak solutions considered in this paper, it is customary to relax(1.4) to inequality t (s(,ϑ))+div x (s(,ϑ)u). (1.5) The Euler system (1.1) can be written in the conservative variables: the density, the momentum m = u, and the energy E = 1 2 u 2 +e(,ϑ): t+div x m =, ( ) m m t m+div x + x p(,m,e) =, ( t E +div x E m ( )+div x p(,m,e) m ) =. (1.6) 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: p(, ϑ) >, e(, ϑ) ϑ > ; (1.7) (,m,e) S(,m,E) s(,m,e) is a strictly concave function of (,m,e). (1.8) 1.2 Smooth and weak solutions Smooth solutions are at least continuously differentiable on the set [,T), >, ϑ >, and satisfy ( ) pointwise. Smooth solutions are known to exist locally on a maximal interval [,T max ) as long as the initial data are smooth enough and satisfy the corresponding compatibility conditions on, see last part of Section 4 for details. Moreover, it is well known that T max < for a generic class of data. The weak solutions satisfy (1.1), (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 3

4 called admissible if the entropy inequality (1.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 ( ) admits infinitely many admissible weak solutions defined on a given time interval (, 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 (;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 (1.6) in the literature use such a splitting method. A vector field m can be decomposed by means of Helmholtz decomposition, i. e. it can be written as the sum of a solenoidal vector field H[m] and the gradient of a scalar field. More precisely, the Helmholtz projection operator H is defined as m = H[m]+H [m], where H = x Φ, x Φ = div x m, ( x Φ m) n =, Φ dx =. Accordingly, for m = v+ x Φ, v = H[m], the system (1.6) can be written in the form t+ x Φ =, [ ( m m t v+h div x 1 )] m 2 N I =, [ ( m m t ( x Φ)+H div x 1 )] ( ) m 2 1 N I m 2 + x + x p(,m,e) =, N ( t E +div x E m ( )+div x p(,m,e) m ) =. (1.9) To the best of our knowledge, there are only two convex integration ansatzes for obtaining wild solutions to the Euler system (1.6) in the literature. Let us begin with the first one, which is used in [9]. This ansatz is based on replacing equations (1.9) 2, (1.9) 3 by the system ( m m t v+div x 1 ) m 2 N I =, ( ) 1 m 2 t ( x Φ)+ x + x p(,m,e) =, N ( ) div x v =, ( ) t+ x Φ =. (1.1) Equation (1.1) 1 represents the volume preserving part of the motion, while (1.1) 2 may be seen as a wave equation governing the propagation of acoustic waves. Note that the fields v and x Φ 4

5 as well as the fluxes in (1.1) are orthogonal with respect to the L 2 scalar product, specifically, [ m m v x Φ dx =, 1 ] [( ) ] m 2 1 N I m 2 : +p(,m,e) I =, N where the latter is true since it is a product of a traceless matrix and a multiple of the identity matrix. Although (1.1) 1 is apparently overdetermined, it admits (infinitely many) weak solutions for any fixed and Φ, cf. Appendix A. Definition 1.1. A weak solution [, m, E] of the Euler system (1.6) is called wild solution if [,m,e] satisfy (1.1), with v = H[m], x Φ = H [m]. The corresponding initial data [,m,e ] are called wild initial data. Note that a technique similar to (1.1) 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], [1] among others. We will first focus on the ansatz (1.1) and afterwards extend our results to the wild solutions obtained via the Riemann problem in Section 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 { } L 1 +,s (;R N+2 ) = [,m,e] L 1 (;R N+2 ), E, s(,m,e) s >. Note that this is not very restrictive as any admissible weak solution satisfies the minimum principle s(,m,e)(t,x) ess inf y s((y),m (y),e (y)) a.a. in (,T), cf. [2, Section 2.1.1]. We say that a sequence of data {,n,u,n,ϑ,n } n=1, or, equivalently, {,n,m,n,e,n } n=1, (W)-converges to [,m,e ], if [,n,m,n,e,n ] (W) [,m,e ],,n >, s(,n,m,n,e,n ) s > ; 5

6 [,n,m,n,e,n ] [,m,e ] in L 1 (;R N+2 ); (1.11) theinitialdata[,n,m,n,e,n ]giverisetoasequenceofadmissibleweaksolutions[n,m n,e n ] satisfying T ( mn m n 1 ) m n 2 I : 2 n N xϕ dx dt as n for any ϕ Cc ((,T) ). n (1.12) Note that, in view of (1.1), H [ div x ( mn m n n 1 m n 2 N n )] I = whenever [n,m n,e n ] are wild solutions in the sense of Definition 1.1. In particular, (1.12) is trivially satisfied for any sequence of wild initial data in the sense of Definition 1.1. Finally, notethat thesequence ofsolutions {n,m n,e n } n=1 need not a prioriconverge strongly, and, consequently, (1.12) does not impose any similar restriction on a possible limit solution emanating from [,m,e ]. Definition 1.2. We say that a trio [,m,e ] is reachable if there exists a sequence of initial data {,n,m,n,e,n } n=1 such that [,n,m,n,e,n ] (W) [,m,e ]. 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 L 1 +,s (;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 p(,ϑ) = ϑ, and e(,ϑ) = c v ϑ, c v >. Equivalently, in the conservative variables, p = 1 ( E 1 ) m 2, ϑ = 1 ( E 1 ) m 2 >. (1.13) c v 2 c v 2 Here the kinetic energy is defined as The entropy reads m 2 if m = = if = m 2 otherwise. S(,m,E) = c v log 6 ( ) E 1 m 2 2. c v1+ 1 cv

7 Here is our main result. Theorem 1.3. Let s R be given. Let 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 L 1 +,s (;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 L 1 +,s (;R N+2 ). 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 [,m,e ] giving rise to a smooth solution [,m,e] defined on a maximal time interval [,T max ). Then we proceed in several steps: Assuming the data [,m,e ] are reachable we show that the associated sequence of solutions {n,m n,e n } n=1 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 Section 2. Next we realize that, thanks to the strong convergence required in (1.11), the limit DMV solution starts from the initial data represented by the parameterized family of Dirac masses, { δ(x),m (x),e (x)} x. Thanks to the general weak strong uniqueness principle, we conclude that the DMV solution coincides with the strong solution on the time interval [,T max ), more specifically, the DMV solution is represented by the parameterized family of Dirac masses, { δ(t,x),m(t,x),e(t,x) } (t,x) (,T max), see Section 3. In particular, the solutions [n,m n,e n ] converge strongly and we may assume n, m n m, E n E a.a. in (,T max ). (1.14) Thanks to the strong convergence (1.14), condition (1.12), and smoothness of the limit solution, we deduce that ( ) m m div x div x 1 ( ) m 2 N x = in (,T). see Section 3 for details. In particular, as the limit solution is continuous up to t =, ( ) m m div x div x 1 ( ) N m 2 x = in. (1.15) 7

8 In Section 4, we show that (1.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 Section 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 γ = 1+ 1 c v. We say that [,m,e] is an admissible weak solution to the Euler system ( ) if C weak ([,T];L γ ()),, m C weak ([,T];L 2γ γ+1 (;R N )),E C weak ([,T];L 1 ()), E, and the following holds: T for any ϕ C 1 c ([,T) ); T [ m t ϕ+ ( ) m m for any ϕ C 1 c([,t) ;R N ), ϕ n = ; T [ t ϕ+m x ϕ] dx dt = ϕ(, ) dx ] : x ϕ+p(,m,e)div x ϕ [ ( ) ] m E t ϕ+(e +p(,m,e)) x ϕ for any ϕ C 1 c([,t) ); dx dt = m ϕ(, ) dx dx dt = E ϕ(, ) dx (2.1) Z(s(,m,E ))ϕ(, ) dx T [Z(s(,m,E)) t ϕ+z(s(,m,e))m x ϕ] dx dt (2.2) for any ϕ C 1 c ([,T) ), ϕ, and Z(s) Z, Z (s). 8

9 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 { } P = [,m,e], m R N, E be the phase space associated to the Euler system. Definition 2.3 (DMV solution). A parametrized family of probability measures {V t,x } (t,x) (,T) on the set P, V L weak ( )((,T) ;prob[p]) is called a dissipative measure valued (DMV) solution to the compressible Euler system with the initial data {V,x } x, if the following holds: T for any ϕ C 1 c ([,T) ); [ V t,x ; t ϕ+ V t,x ;m x ϕ] dx dt = V,x ; ϕ(, ) dx T [ V t,x ;m t ϕ+ V t,x ; m m T = x ϕ : dµ C V,x ;m ϕ(, ) dx ] : x ϕ+ V t,x ;p(,m,e) div x ϕ dx dt for any ϕ C 1 c([,t) ;R N ), ϕ n =, where µ C is a (vectorial) signed measure on [,T] ; V τ,x ;E dx+d(τ) V,x ;E dx (2.3) for a.a. τ (,T), where D L (,T), D is called dissipation defect; T V,x ;Z(s(,m,E)) ϕ(, ) dx [ V t,x ;Z(s(,m,E)) ] t ϕ dx dt+ T for any ϕ C 1 c ([,T) ), ϕ, and Z(s) Z, Z (s) ; 9 [ V t,x ;Z(s(,m,E))m x ϕ] dx dt

10 τ µ C M([,τ] ) c D(t) dt for a.a. τ (,T). 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 (2.1) has been replaced by its integrated version (2.3). The reader may consult [2] for a thorough discussion of this new concept of solution. 2.2 Generating a DMV solution Suppose now that,n in L 1 (), m,n m in L 1 (), and E,n E in L 1 (),,n > a.a. in, s(,n,m,n,e,n ) s >, (2.4) where [,n,m,n,e,n ] are initial data of admissible weak solutions [n,m n,e n ] defined on (,T). Passing to a suitable subsequence, we may suppose that {n,m n,e n } n=1 generates a Young measure {V t,x } (t,x) (,T) ), V t,x prob[p]. As shown in [2, Section 2.1], the Young measure {V t,x } (t,x) (,T) ) is a DMV solution of the Euler system in the sense of Definition 2.3. Moreover, in view of (2.4), the initial measure {V,x } x reads V,x = δ(x),m (x),e (x) for a.a. x. 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 Section 2.2, we suppose now that [,n,m,n,e,n ] (W) [,m,e ], where [,m,e ] are now regular initial data generating a C 1 solution [,m,e] of the Euler system ( ) on a maximal time interval [,T max ). Without loss of generality, we may suppose that < T < T max. Applying the weak strong uniqueness principle we obtain In particular, we may assume V t,x = δ(t,x),m(t,x),e(t,x) for a.a. (t,x) (,T). (3.1) n, m n m, E n E a.a. in (,T) 1

11 for the associated weak solutions [n,m n,e n ]. Consequently, as [n,m n,e n ] satisfy (1.12), we get T [( ) m m : 2 x ϕ 1 ] m 2 N xϕ dx dt = for any ϕ Cc ((,T) ). Finally, as the limit solution is continuous up to the time t =, we may infer that in particular, [( ) m m : 2 x ϕ 1 ] m 2 x ϕ dx = for any ϕ Cc (). (3.2) N For C 2 initial data, relation (3.2) can be rewritten as a non linear differential equation ( ) ( ) m m 1 m 2 div x div x x = in. N 4 Reachability In this section we finish the proof of Theorem 1.3. Introducing w = m we may write (3.2) in a concise form [ w w : 2 x ϕ 1 ] N w 2 x ϕ dx = for any ϕ Cc (). (4.1) We consider solutions of (4.1) in the space L 2 (;R N ), { } S = w L 2 (;R N ) w solves (4.1) It is easy to check that: the set S is closed in L 2 (;R N ); if w S, then λw S for any λ R. Lemma 4.1. The set S is nowhere dense in L 2 (;R N ), meaning the (closure of the) set S does not contain any ball in L 2 (;R N ). Proof. Arguing by contradiction, we suppose that there is w L 2 (;R N ) and r > such that the ball centered at w with radius r is contained in S. Hence [ (w+w ) (w+w ) : 2 xϕ 1 ] N w+w 2 x ϕ dx = for any ϕ Cc (), (4.2) and any w, w L 2 (;R N ) r.. 11

12 Next, we show that this implies that the ball centered at with radius r is contained in S, too. To this end, we write [ (w+w ) (w+w ) : 2 xϕ 1 ] N w+w 2 x ϕ dx [ = w w : 2 xϕ 1 ] [ N w 2 x ϕ dx+ (w w +w w) : 2 xϕ 2 ] N w w x ϕ dx [ + w w : 2 x ϕ 1 ] N w 2 x ϕ dx. Thus using the fact that w solves (4.1) and (4.2) we conclude [ w w : 2 xϕ 1 ] [ N w 2 x ϕ dx+ (w w +w w) : 2 xϕ 2 ] N w w x ϕ for any ϕ C c () and any w, w L 2 (;R N ) r. dx = (4.3) If relation (4.3) holds for any w, w L 2 (;R N ) r, it must hold for λw, λ 1. Consequently, we get from (4.3), [ = (λ 2 λ) w w : 2 x ϕ 1 ] N w 2 x ϕ dx [ +λ w w : 2 x ϕ 1 ] [ N w 2 x ϕ dx+λ (w w +w w) : 2 x ϕ 2 ] N w w x ϕ dx [ =(λ 2 λ) w w : 2 xϕ 1 ] N w 2 x ϕ dx for any λ 1. We deduce [ w w : 2 xϕ 1 ] N w 2 x ϕ dx = for any ϕ C c (), and any w, w L 2 (;R N ) r. Thus the set S contains the ball of radius r > centered at but as it is invariant with respect to multiplication by any real constant, we conclude which is obviously false as soon as N 2. S = L 2 (;R N ), Lemma 4.2. The set of reachable data in the sense of Definition 1.2 is closed in L 1 ((,T) ). 12

13 Proof. Let [,n,m,n,e,n ] be a sequence of reachable data such that [,n,m,n,e,n ] [,m,e ] in L 1 ((,T) ;R N+2 ). To see that the limit is reachable, we have to find a generating sequence satisfying (1.11), (1.12). To this end, consider the generating sequences of data [,n,m,m,n,m,e,n,m ] satisfying [,n,m,m,n,m,e,n,m ] (W) [,n,m,n,e,n ] as m, meaning, in particular, T ( mn,m m n,m n,m 1 N ) m n,m 2 I : 2 xϕ dx dt as m for any ϕ Cc ((,T) ) n,m for the corresponding solutions n,m, m n,m and any fixed n. Now, let { ϕ k } k=1, ϕk C c ((,T) ) be a countable dense set in the Sobolev space W l,2 ((,T) ) C c ((,T) ) l,2 By the diagonal method, we can find a subsequence C 2 ([,T] ) as soon as l > N +5. (4.4) 2 [,(n,m)(j),m,(n,m)(j),e,(n,m)(j) ] [,m,e ] as j in L 1 ((,T) ;R N+2 ) such that T ( m(n,m)(j) m (n,m)(j) (n,m)(j) 1 N m (n,m)(j) 2 ) I : 2 x ϕk dx dt as j for any k. (4.5) (n,m)(j) Indeed, for given j 1, we find n(j) such that [,n(j),m,n(j),e,n(j) ] [,m,e ] < 1 j, and m = m(n(j)) such that T [,n(j),m(n(j)),m,n(j),m(n(j)),e,n(j),m(n(j)) ] [,n(j),m,n(j),e,n(j) ] < 1 j, ( mn(j),m(n(j)) m n(j),m(n(j)) 1 N n(j),m(n(j)) 13 m n(j),m(n(j)) 2 ) I n(j),m(n(j)) for all ϕ k C c : 2 x ϕk dx dt < 1 j ((,T) ), k j.

14 Finally, by virtue of (4.4), any 2 x ϕ can be uniformly approximated by 2 x ϕk, and, as the energy E n(j),m(j) are bounded, we conclude that (4.5) holds for any ϕ C c ((,T) ). Thus we have shown that whence [,m,e ] is reachable. [,n(j),m(j),m,n(j),m(j),e,n(j),m(j) ] (W) [,m,e ] as j ; 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 [,m,e ] giving rise to local-in-time regular solutions is dense in the space L 1 +,s (;R N+2 ). To this end, we record the following result by Schochet [12, Theorem 1] that asserts the existence of smooth solutions whenever: R N is a bounded domain with a sufficiently smooth boundary, say of class C ; the initial data [,E,ϑ,E,u,E ] belong to the class,e,ϑ,e W m,2 (), u,e W m,2 (;R N ),,E, ϑ,e > in, m > N; the compatibility conditions hold for k =,1,...,m. k tu,e n = In particular, any trio of initial data, [,m,e ] C (;R N+2 ), m =, = >, ϑ = ϑ > on a neighborhood of, (4.6) gives rise to a smooth solution defined on a maximal time interval [,T max ). Seeing that the set (4.6) is dense in L 1 +,s (;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 (1.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 (1.1). 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 (1.9). The first method [9] was explained above, see (1.1). The second method considers 14

15 initial data constant in each of the two half spaces (Riemann data), see [1] and also [6], [5], [1] for the isentropic case. Here instead of (1.1) equations (1.9) 2 and (1.9) 3 are replaced by ( m m t v+div x 1 ) m 2 ( ) N I U =, div x v =, ( ) 1 m 2 ( ) t ( x Φ)+div x (U)+ x + x p(,m,e) =, t+ x Φ =, N with a suitable matrix field U. The density, the acoustic potential x Φ as well as the field U are piecewise constant in (,T). In particular, there exists a partition of the set (,T) into a finite number of sectors Q i,i = 1,...m, such that the wild solutions satisfy Q i [ m m 1 N m 2 I ] : 2 xϕ dx dt = for any ϕ Cc (Q i ), i = 1,...,m. Consequently, we may accommodate the data reachable by this method replacing (1.12) by T ( mn m n 1 ) m n 2 I : 2 x N ϕ dx dt as n for any ϕ C c (Q), n n where Q (,T) is an open set such that Q = [,T]. The above observation motivates the following extension of the concept of reachability. We say that Q is a partition of the domain (,T) if Q (,T) is an open set; Q = [,T]. We say that a family of partitions 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 Q contains a subsequence such that (H) Q c as k for some Q Q, Q c n(k) where the symbol (H) denotes convergence in the Hausdorff metric and Q c n denotes the complement Q c n ([,T] )\Q n. An example of a closed partition set related to the convex integration method is { } Q = Q (,T) Q = ((,T) )\( M i=1 H i), H i a hyperplane in R N+1, where M is a given positive integer. Note that this is indeed a closed partition set since (,T) is a bounded subset of R N+1. The following property is well known, cf. [11]: Q c n (H) Q c and K Q is a compact set there exist n(k) such that K Q n for all n n(k). (5.1) 15

16 Next, we introduce the following generalization of (W) convergence. Let Q be a closed partition set. We say that a sequence of data {,n,m,n,e,n } n=1 (W[Q])-converges to [,m,e ], if [,n,m,n,e,n ] (W[Q]) [,m,e ],,n >, s(,n,m,n,e,n ) s > ; [,n,m,n,e,n ] [,m,e ] in L 1 (;R N+2 ); theinitialdata[,n,m,n,e,n ]giverisetoasequenceofadmissibleweaksolutions[n,m n,e n ] satisfying T ( mn m n 1 ) m n 2 I : 2 x N ϕ dx dt as n for any ϕ C c (Q) (5.2) for some Q Q. n n Definition 5.1. Let Q be a closed partition set. We say that a trio [,m,e ] is Q reachable if there exists a sequence of initial data {,n,m,n,e,n } n=1 such that [,n,m,n,e,n ] (W[Q]) [,m,e ]. Obviously any reachable data in the sense of Definition 1.2 are Q-reachable so the set of Q reachable data is always larger for any closed partition set. In order to adapt the arguments used in Sections 3, 4 we need the following result. Lemma 5.2. Let Q be a closed partition set. Then the set of Q reachable data is closed in L 1 ((,T) ;R N+2 ). Proof. Let [,n,m,n,e,n ] be a sequence of Q reachable data such that [,n,m,n,e,n ] [,m,e ] in L 1 ((,T) ;R N+2 ). Let [,n,m,m,n,m,e,n,m ] be the corresponding generating sequences, [,n,m,m,n,m,e,n,m ] (W[Q]) [,n,m,n,e,n ] as m. In particular, the corresponding solutions m,n, m m,n satisfy T ( mn,m m n,m 1 ) m n,m 2 I : 2 N xϕ dx dt as m for any ϕ Cc (Q n ) n,m n,m 16

17 for some Q n Q and any fixed n. As the partition set Q is closed, there exists a partition Q Q such that K Q a compact set there exist n(k) such that K Q n for all n n(k). (5.3) at least for a suitable subsequence (not relabeled). Let { ϕ k } k=1, ϕ C c (Q) be a countable dense subset of the Sobolev space W l,2 (Q) C c (Q) l,2 Using property (5.3), we can find a subsequence C 2 (Q) as soon as l > N +5. (5.4) 2 [,n(j),m(j),m,n(j),m(j),e,n(j),m(j) ] [,m,e ] as j in L 1 ((,T) ;R N+2 ) satisfying T ( mn(j),m(j) m n(j),m(j) n(j),m(j) 1 N ) m n(j),m(j) 2 I : 2 xϕ k dx dt as j for any k. (5.5) n(j),m(j) By virtue of (5.4), any 2 xϕ can be uniformly approximated by 2 xϕ k, and, as the energies E n(j),m(j) are bounded, we conclude that (5.5) holds for any ϕ Cc (Q). We may infer that [,n(j),m(j),m,n(j),m(j),e,n(j),m(j) ] (W[Q]) [,m,e ] as j ; whence [,m,e ] is Q reachable. Finally, we observe that any regular initial data [,m,e ] satisfy (3.2). To see this, we first observe, similarly to Section 3, that T [( ) m m : 2 x ϕ 1 ] m 2 N xϕ dx dt = for any ϕ Cc (Q). As, m are smooth, we get div x div x ( m m ) 1 N m 2 x = in Q. However, as Q = [,T] and the second derivatives of ρ, m are continuous, this implies ( ) m m div x div x 1 N m 2 x = in (,T). 17

18 In particular, we obtain (3.2) for the initial data. We have shown the following extension of Theorem 1.3. Theorem 5.3. Let s R be given. Let R N, N = 2,3 be a bounded smooth domain. Let Q be a closed partition set in (,T). Then the complement of the set of Q reachable data is an open dense set in L 1 +,s (;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]. A 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 and ϑ be piecewise constant strictly positive functions defined on. And let T > be given. Then there exists a constant Λ > such that for any Λ Λ there is u L (;R N ) such that the problem (1.1), (1.5), (1.8) admits infinitely many weak solutions in (,T), with the initial data, m = u, E = 1 2 u 2 +e(,ϑ ). In addition, these solutions enjoy the following properties: 1 m 2 N = Λ+p for a.a. (t,x) including t =, (A.1) div x m =, (t,x) = (x) for all t, S(t,x) = S (x) for all t, and the renormalized entropy inequality (2.2) holds as equality. (A.2) (A.3) (A.4) Note that here (1.1) is satisfied with Φ =. 18

19 References [1] Al Baba, H., Klingenberg, C., Kreml, O., Mácha, V., Markfelder, S.: Non-uniqueness of admissible weak solutions to the Riemann problem for the full Euler system in 2D. Preprint, submitted, 218, arxiv: [2] Březina, J., Feireisl, E.: Measure-valued solutions to the complete Euler system. Preprint, to appear in J. Math. Soc. Japan, 217, arxiv: [3] Chen, G.-Q., Frid, H., Li, Y.: Uniqueness and stability of Riemann solutions with large oscillation in gas dynamics. Comm. Math. Phys. 228(2), , 22 [4] Chiodaroli, E.: A counterexample to well-posedness of entropy solutions to the compressible Euler system. J. Hyperbolic Differ. Equ. 11(3), , 214 [5] Chiodaroli, E., De Lellis, C., Kreml, O.: Global ill-posedness of the isentropic system of gas dynamics. Comm. Pure Appl. Math. 68(7), , 215 [6] Chiodaroli, E., Kreml, O.: On the energy dissipation rate of solutions to the compressible isentropic Euler system. Arch. Ration. Mech. Anal. 214(3), , 214 [7] De Lellis, C., Székelyhidi Jr., L.: On admissibility criteria for weak solutions of the Euler equations. Arch. Ration. Mech. Anal. 195(1), , 21 [8] Feireisl, E.: Weak solutions to problems involving inviscid fluids. In: Mathematical Fluid Dynamics, Present and Future, pp , Springer Proceedings in Mathematics and Statistics 183, Springer-Verlag, Tokyo, 216 [9] Feireisl, E., Klingenberg, C., Kreml, O., Markfelder, S.: On oscillatory solutions to the complete Euler system. Preprint, submitted, 217, arxiv: [1] 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), , 218 [11] Pironneau, O.: Optimal shape design for elliptic systems. Springer series in Computational Physics, Springer-Verlag, New York, 1984 [12] Schochet, S.: The compressible Euler equations in a bounded domain: existence of solutions and the incompressible limit. Comm. Math. Phys. 14(1), 49 75, 1986 [13] Székelyhidi Jr., L., Wiedemann, E.: Young measures generated by ideal incompressible fluid flows. Arch. Ration. Mech. Anal. 26(1), ,