arxiv: v1 [math.ap] 24 Dec 2018

Transcription

1 arxiv:8.0997v [math.ap] 4 Dec 08 Non uniqueness of admissible weak solutions to the compressible Euler equations with smooth initial data Elisabetta Chiodaroli Ondřej Kreml Václav Mácha Sebastian Schwarzacher Institute of Mathematics, Czech Academy of Sciences Žitná 5, 5 67 Praha, Czech Republic Dipartimento di Matematica, Università di Pisa Via F. Buonarroti /c, 567 Pisa, Italy Abstract: We consider the isentropic Euler equations of gas dynamics in the whole twodimensional space and we prove the existence of a C initial data which admit infinitely many bounded admissible weak solutions. Taking advantage of the relation between smooth solutions to the Euler system and to Burgers equation we construct a smooth compression wave which collapses into a perturbed Riemann state at some time instant T > 0. In order to continue the solution after the formation of the discontinuity, we apply the theory developed by De Lellis and Székelyhidi in [] [] in order to construct infinitely many solutions. We introduce the notion of an admissible generalized fan subsolution to be able to handle data which are not piecewise constant and we reduce the argument to the finding of a single generalized subsolution. Introduction In this paper, we study the Cauchy problem for the isentropic Euler equations of gas dynamics in two space dimensions, i.e. we look for a density ρt,x and a velocity vt,x satisfying t ρ+div x ρv = 0. t ρv+div x ρv v+ x pρ = 0. in 0,T R together with the initial conditions ρ,v0, = ρ 0,v 0 in R..3 The system.-. is a paradigmatic example of a hyperbolic system of conservation laws. The mathematical theory of such systems of partial differential equations is still far from being completely understood in more than one space dimension. It is however well known and observed already in the case of the most simple example of a scalar law, the Burgers equation, that singularities appear even in case of a smooth initial data and that weak solutions with no additional properties are not unique. The notion of admissibility in the form of entropy inequalities was used for one-dimensional systems and multi-dimensional scalar equations to restore uniqueness of weak solutions at least in proper classes of solutions; for more details about this nowadays classical theory we refer for example to the monograph [9]. elisabetta.chiodaroli@unipi.it kreml@math.cas.cz macha@math.cas.cz schwarz@math.cas.cz

2 This theory in particular motivates us to look for weak solutions i.e. solutions in the sense of distributions which satisfy in addition the entropy inequality. In the case of the system.-. this is rather a form of energy inequality, since the only existing mathematical entropy is the physical total energy of the system. Therefore, we call a weak solution admissible, if it satisfies the energy inequality t ρeρ+ρ +div v x ρeρ+pρ+ρ v v 0,.4 in the sense of distributions on [0,T R. In.4, eρ denotes the internal energy that is related to the pressure through the relation pr = r e r..5 In their groundbreaking work [], [] on the existence of infinitely many solutions to the incompressible Euler equations in multiple space dimensions, De Lellis and Székelyhidi also noted that this theory can be applied to the compressible isentropic system as well and proved in [] the existence of initial data ρ 0,v 0 L R for which there exists infinitely many admissible weak solutions to.-.3. This result was further improved by the first author in [5] and by Feireisl in [5], proving the existence of initial data ρ 0 C R,v 0 L R for which there exists infinitely many admissible weak solutions, locally in time and globally in time respectively. We call such data allowing for infinitely many admissible weak solutions wild data. Quite surprisingly in [6], the first two authors together with De Lellis proved that even in the class of Lipschitz initial data there exists an example of wild data. The proof in [6] relied substantially on the convex integration method developed by De Lellis and Székelyhidi [], [] for the incompressible Euler equations combined with a clever analysis of the Riemann problem for.-.. A building block in the construction carried out in [6] are piecewise constant subsolutions to.-. which were also inspired by the work [9] of Székelyhidi on irregular solutions of the incompressible Euler equations with vortex-sheet data. The wild initial data in [6] are chosen so to generate a compression wave collapsing to a shock giving then rise to a particular Riemann datum at a certain time instant t = T starting from which infinitely many admissible solutions forward in time can be constructed. Such a solution ρ,v is indeed Lipschitz on the time interval 0, T but not smoother, since the first derivatives are not continuous on the boundaries of the wave. A natural question therefore arises, whether a similar property can be achieved also for even more regular initial data. In this work, we provide the positive answer. Even so the answer itself is hardly surprising the technical machinery that was necessary to be built in order to produce these counterexamples is rather involving; indeed, we provide technical tools in the proof which may well be of further use. Our main theorem reads as follows. Theorem.. Let pρ = ρ and let T > 0. There exist initial data ρ 0,v 0 C R and δ 0 > 0, such that there exists infinitely many bounded admissible weak solutions to the Euler system.-.3 on the time interval 0,T +δ 0. These solutions all coincide with a unique C solution to the Euler system.-.3 on the time interval 0,T. Even if C data are shown to allow for infinitely many bounded admissible solutions, such solutions still display a high level of irregularity, i.e. they are only bounded and highly oscillatory. An open question for the compressible Euler equations concerns the level of regularity of solutions displaying such a non standard behavior as in Theorem.. We note that the related question concerning regularity of solutions needed for the energy conservation was recently studied in [6]. This issue has been now extensively tackled for the incompressible Euler equations leading to the proof of the famous Onsager conjecture. For details see [4], [8] for one side of the conjecture and [3], [], [], [3], [4], [0], [7] and [8] for the other side. The paper is structured as follows. In Section we introduce the notion of smooth compression wave, which is one of the building blocks in our proof. In Section 3 we define subsolutions in such

3 a way that the proof of Theorem. is reduced to finding a single subsolution so to apply the theory of De Lellis and Székelyhidi. In our case we generalize the definition of subsolution given in [6] since the assumption that the subsolutions consist of piecewise constant functions is too restrictive for our construction. The main ingredient to the proof of the Proposition 3.6 which states the relation between the existence of infinitely many admissible solutions and the existence of a single properly defined subsolution is Lemma 3. which is an easy generalization of [6, Lemma 3.7]. The main part of the proof of Theorem. is contained in Section 4. Remark.. In this paper, it is useful to work both with the velocity v and the momentum m = ρv. We emphasize here that we never run into troubles arising from possible presence of vacuum regions using either of these notions, since we always work with strictly positive densities. Smooth compression wave In this section we recall some standard facts about the one-dimensional rarefaction waves solving the Riemann problem for the Euler system and introduce the smooth compression wave. Since we wanttokeepthissectionrelatedtothetwo-dimensionalcase, wedenotev = v,v, m = m,m and we work here with the system for unknowns ρ,v in the space-time domain t,x R + R, i.e. t ρ+ x ρv = 0. t ρv + x ρv +pρ = 0. ρ,v 0, = ρ 0,v 0..3 It is well known that the characteristic speeds of this system are and the functions λ = v p ρ, λ = v + p ρ.4 ρ p r w = v + dr, 0 r ρ p r w = v dr.5 0 r are, respectively, - and -Riemann invariants. The classical theory, see for example [9, Theorem 7.6.6], yields that every i-riemann invariant is constant along any i-rarefaction wave. It is then not difficult to observe that for example in the case of a -rarefaction wave, under the condition w const, the system.-.3 reduces to a simple Burgers equation for λ t λ +λ x λ = 0.6 λ 0, = λ 0..7 As it is well known, the solution to the Burgers equation can be obtained by the method of characteristics as long as the characteristic curves do not intersect each other. The characteristic curves are defined as x = λ 0 rt+r,.8 for r R, and the solution λ to the Burgers equation.6 is constant along these curves, taking values λ 0 r. The classical theory also yields, that for the Riemann problem λ 0 x = { λ for x < 0 λ + for x > 0,.9 with λ < λ +, the problem.6-.7 admits a Lipschitz solution in the form of a rarefaction wave λ for x < λ t x λ t,x = for λ t < x < λ + t.0 t λ + for x > λ + t. 3

4 Reversing time and space we obtain a solution to the Burgers equation called the compression wave. Namely for initial data λ for x < λ T λ 0 x = x for λ T < x < λ + T. T for x > λ + T λ + for some T > 0 and λ > λ +, the solution to.6-.7 on the time interval 0,T is Lipschitz and has the form λ for x < λ T t λ t,x = x for λ T t < x < λ + T t. T t for x > λ + T t, λ + in time t = T the solution has the form λ T,x = { λ for x < 0 λ + for x > 0,.3 and for t > T the solution consists of a shock. We introduce the smooth compression wave as follows. Consider initial data in the form λ for x < λ T ζ f 0 x for λ T ζ < x < λ T +ζ λ 0 x = x for λ T +ζ < x < λ + T ζ.4 T f+ 0 x for λ + T ζ < x < λ + T +ζ for x > λ + T +ζ, λ + where ζ > 0 is small and f± 0 are smooth strictly monotone functions with strictly monotone first derivatives such that λ 0 C R. The solution to the Burgers equation.6-.7 then has the following form on the time interval 0,T λ for x < λ T t ζ f t,x for λ T t ζ < x < λ T t+ζ T t T λ t,x = x for λ T t+ζ T t < x < λ + T t ζ T t.5 T t T T f + t,x for λ + T t ζ T t < x < λ + T t+ζ T for x > λ + T t+ζ λ + and in particular at time t = T a discontinuity appears λ for x < ζ f T,x for ζ < x < 0 λ T,x = f + T,x for 0 < x < ζ λ + for x > ζ..6 Here again the functions f ± t, are smooth and strictly monotone with strictly monotone first derivatives for all t 0,T] and λ C 0,T R. It is also easy to observe that the left derivative of f T, and the right derivative of f + T, at point x = 0 are equal to. 4

5 Choosing the initial data.3 in such a way that w defined in.5 is constant in R and λ defined in.4 takes the form.4 we obtain a smooth solution to the Euler system.-.3 on the time interval 0,T in the form w const on [0,T] R and λ given by.5,.6. Such solution satisfies the energy equality.4 with an equality sign. The main part of the proof of Theorem. then consists of proving that the solution on the time interval [0,T] can be prolonged in a non-unique way. In particular we will take values of ρ,v at time t = T as the initial data and prove the existence of infinitely many admissible weak solutions starting from this initial data. 3 Subsolutions Our aim in this section is to properly define admissible subsolutions. Unlike in previous works on the construction of non unique admissible solutions, here it is not enough to work with piecewise constant subsolutions. First, we recall the crucial lemma based on the convex integration technique developed for the incompressible Euler equations which allows to pass from a single subsolution to infinitely many solutions, see [6, Lemma 3.7]. In what follows, the symbol S0 stands for the space of symmetric matrices with zero trace and I denotes the identity matrix. Lemma 3.. Let ṽ,ũ R S0 and let C > 0 be such that ṽ ṽ ũ < C I. For any open set Ω R R there are infinitely many maps v,u L R R,R S0 such that i v and u vanish identically outside Ω ii div x v = 0 and t v +div x u = 0 iii ṽ +v ṽ +v ũ+u = C I a.e. in Ω. We need to generalize this lemma in order to be able to work with piecewise continuous subsolutions for the compressible Euler equations instead of piecewise constant ones. We can prove the following Lemma. Lemma 3.. Let Ω R R be an open set. Let ρ, m,ũ, C CΩ,R + R S 0 R + such that m m ρ ũ < C I 3. pointwise in Ω in the sense of positive definiteness. Let us assume furthermore that C is a bounded function on Ω. Then, there are infinitely many maps m,u L R R,R S 0 such that i m and u vanish identically outside Ω. ii div x m = 0 3. t m+div x u = iii It holds a.e. in Ω m+m m+m ρ ũ+u = C I 3.4 In order to prove Lemma 3., we follow the strategy outlined in [6] for the proof of Lemma 3.7 therein, but we need to borrow some ingredients also from [5]. 5

6 3. Proof of Lemma 3. We define X 0 to be the space of m,u Cc Ω,R S0 which satisfy ii and the pointwise inequality m+m m+m ρ ũ + u < C I. We then take the closure X of X 0 in the L weak topology. Since C is a bounded continuous function on Ω, X is a bounded weakly closed subset of L. This implies that the L weak topology is metrizable on X, thus producing a complete metric space X,d. Observe that any element in X satisfies i and ii. Our aim is to show that on a residual set in the sense of Baire category iii holds. Following [6], we define for any N N\{0} the map I N as follows: to m,u we associate the corresponding restrictions of these maps to N,N B N 0. We then consider I N as a map from X,d to Y, where Y is the space L N,N B N 0,R S 0 endowed with the strong L topology. Arguing as in [, Lemma 4.5] it is easily seen that I N is a Baire- map and hence, from a classical theorem in Baire category, its points of continuity form a residual set in X. We claim that CLAIM if m,u is a point of continuity of I N, then iii holds a.e. on N,N B N 0. CLAIM implies then iii for those maps at which all I N are continuous which is also a residual set. The proof of CLAIM is achieved as in [, Lemma 4.6] showing that: Property If m,u X 0, then there is a sequence m k,u k X 0 converging weakly to m,u for which lim inf m+m k L k Γ m+m L Γ +β C ρ dxdt m+m L Γ, Γ where Γ = N,N B N 0 and β depends only on Γ. Indeed assuming that Property holds, fix then a point m,u X where I N is continuous and assume by contradiction that iii does not hold on Γ. By definition of X there is a sequence m k,u k X 0 converging weakly to m,u. Since the latter is a point of continuity for I N, we then have that m k m strongly in L Γ. We apply Property to each m k,u k and find a sequence {m k,j,u k,j } such that lim inf m+m k,j L j Γ m+m k L Γ +β C ρ dxdt m+m k L Γ Γ and m k,j,u k,j m k,u k. A standard diagonal argumentthen allowsto conclude the existence of a sequence m k,jk,u k,jk which converges weakly to m,u and such that lim inf m+m k,jk L k Γ m+m L Γ +β C ρ dxdt m+m L Γ > m+m L Γ. Γ However this contradicts the assumption that m,u is a point of continuity for I N. Thus, everything is reduced to proving Property, i.e. to constructing the sequence m k,u k. For the proof of Property we refer the reader to [5, Lemma 4.5] whose statement is analogous; the only two differences are that on the one hand the role of C is played in [5] by the function χ which depends on time only and on the other hand here ρ depends not only on the space variable as in [5], but also on time. However these differences are not substantial for the proof of [5, Lemma 4.5] and hence of Property, which can be carried out similarly under the current setting. 3. From subsolutions to solutions As we stated in the previous section, we assume here that the initial data for the isentropic Euler equations take the form ρ 0,v 0 = { ρ 0,v 0 for x < 0 ρ 0 +,v 0 + for x > 0, 3.5 6

7 where ρ 0 ±,v0 ± are C functions on,0 and 0, respectively such that there exists ζ > 0 and constants R ±,V ± R + R such that ρ 0,v0 = R,V for x < ζ ρ 0 +,v0 + = R +,V + for x > ζ. Moreover we use the momentum formulation of the Euler equations instead of using the velocity, i.e. we introduce m = ρv and rewrite the Euler system.-.3 to 3.6 t ρ+div x m = m m t m+div x + x pρ = ρ The admissibility or energy inequality.4 is rewritten as t ρeρ+ m +div x ρeρ+pρ+ m ρ ρ ρ,m,0 = ρ 0,m m ρ In order to define subsolutions we will work with in this paper we start with a generalized notion of fan partition. Definition 3.3. A generalized fan partition of 0, R consists of three open sets P,P,P + such that P = {t,x : t > 0 and x < ν t} 3. P = {t,x : t > 0 and ν t < x < ν + t} 3. P + = {t,x : t > 0 and ν + t < x }, 3.3 where ν ± t is a couple of continuous functions satisfying ν t < ν + t for all t > 0 and ν 0 = ν + 0 = 0. Definition 3.4. A generalized fan subsolution to the compressible Euler equations is a quadruple of piecewise continuous functions ρ,m,u,c : 0, R R +,R,S0,R + satisfying the following properties. i There exists a generalized fan partition P,P,P + of 0, R such that ρ,m,u,c = ρ i,m i,u i,c i χ Pi, 3.4 i=,,+ where ρ i,m i,u i,c i are continuous functions on P i, u ± = m± m± ρ ± and m ± ρ ± I, C ± = m± ρ ± lim ρ ±,m ± t,x = ρ 0 t 0 + ±,m0 ± x for ±x > ii It holds m m u < C ρ I 3.6 pointwise in P in the sense of positive definiteness. iii The quadruple ρ, m, u, C solves the following system of partial differential equations in the sense of distributions t ρ+div x m = t m+div x u+ x pρ+ C =

8 Definition 3.5. A generalized fan subsolution ρ, m, u, C to the compressible Euler equations is called admissible if it satisfies that i there exists Kt : 0, R + such that for all x,t P it holds eρ + pρ + C x,t = Kt, 3.9 ρ ρ ii the following inequality is satisfied in the sense of distributions t ρeρ+ C +div x ρeρ+pρ+ C m ρ The importance of assumptions in the definitions of the admissible generalized fan subsolution is revealed in the following Proposition. Proposition 3.6. Let p be any C function and ρ 0,m 0 be such that there exists at least one admissible generalized fan subsolution ρ, m, u, C to the Euler equations with initial data 3.9. Then there exists infinitely many bounded admissible weak solutions ρ, m to such that ρ = ρ. Proof. WeuseLemma3.intheregionP with ρ, m,ũ, C = ρ,m,u,c andachieveinfinitely many maps m,u. We set ρ = ρ, m = m + m, u = u + u and C = C. We need to show that ρ, m is an admissible weak solution to the Euler equations First of all we observe that due to 3.4 we know that in the region P it holds u = m m ρ m ρ I 3. C = m ρ, 3. whereas we already know that hold in regions P and P + due to the definition of the subsolution. Combining 3. with 3.7 we easily observe that ρ, m satisfy 3.7. Similarly combining 3.3 with 3.8 we get t m+div x u+ x pρ+ C = and plugging in 3. and 3. we get 3.8. Finally, we observe that using 3.9, 3. and the fact that m is supported in P we have div x ρeρ+pρ+ C m = div x ρ eρ +pρ + C m = Ktdiv x m = 0. ρ ρ 3.4 Therefore it is enough to sum 3.0 with 3.4 to obtain 3.0. The proof of Proposition 3.6 is concluded by observing that as τ 0+ the Lebesgue measure of P {t = τ} converges to zero and therefore the attainment of the initial conditions follows from The existence of a subsolution In this section we study further the admissible generalized fan subsolution which we will look for. 8

9 4. Basic outline We start with examining some properties of solution to the Burgers equation starting from the initial data of the form we are interested in. More precisely, following the discussion in Section and in particular the form of the solution to the Burgers equation at a time instant when the smooth compression wave collapses.6, we consider now λ 0 such that λ for x < ζ f λ 0 x x for ζ < x < 0 = 4. f + x for 0 < x < ζ for x > ζ where f ± are differentiable functions with the following properties λ + lim f x = λ ξ x 0 f + x = λ + +ξ lim x 0+ λ ξ > λ + +ξ lim f x = lim f + x = x 0 x where ξ = ζ/t. Lemma 4.. Let λ 0 has the form 4. such that 4. holds. Then there exist functions s ± t such that s 0 = s + 0 = 0, s t < s + t for all t > 0 and the admissible solution λ t,x to the Burgers equation.6-.7 is C on the sets x < s t and x > s + t. Proof. Since λ > λ +, the admissible solution to the Burgers equation contains a shock. The shock curve satisfies the Rankine-Hugoniot conditions sλ R λ L = λ R λ L 4.3 yielding s = λ R +λ L, 4.4 where s is the shock speed and λ R,λ L are values on the right and on the left side of the shock respectively. Even though the precise position of the shock curve in spacetime is obviously related to the functions f ±, it follows by the continuity of the related quantities that the shock speed will always belong to the interval s λ + +λ ξ, λ + +λ +ξ 4.5 and therefore it is enough to set s ± t = λ + +λ ±ξt. 4.6 Our goal is to find an admissible generalized fan subsolution to the compressible Euler equations with the initial data generated by a smooth compression wave. We will search for a subsolution with the property that ρ,m,u and C do not depend on x and are therefore functions only of t and x. In particular we have to make the following steps i Find functions ν t, ν + t describing the generalized fan partition. 9

10 ii Construct solutions to the Euler equations in the regions P, P + with the initial data generated by a smooth compression wave. These C solutions will be constructed by the method of characteristics for the Burgers equation.6 for λ in regions P and P + while keeping w constant. This will be possible provided ν < s and ν + > s + for all t > 0. iii In order to the equations and the inequality 3.0 to be satisfied on the discontinuities given by x = ν ± t, the appropriate Rankine-Hugoniot conditions need to be satisfied on these two interfaces. iv We need to ensure that in the region P the equations are satisfied, the inequalities 3.6, 3.0 are satisfied and the condition 3.9 holds. In the following we will investigate the points above. The actual construction is to some extent reverse; by first defining ρ 0,m 0 close to a certain piecewise constant Riemann state and then deducing from these functions the fan partition. This construction can be found at the end of the paper. 4. A general set of conditions We set and we introduce the following notation m = α,β 4.7 γ γ u = 4.8 γ γ ν ± t = d dt ν ±t 4.9 ρ ν ±,m ν ±t := ρ ±,m ± t, ν ± t 4.0 f L t := ft, ν t 4. f R t := ft, ν + t 4. for any function f defined in P. Let us now write down the set of Rankine-Hugoniot conditions described in the point iii above. For each time t > 0 we have Rankine-Hugoniot conditions on the left interface: ν ρ ν ρ L = m ν β L 4.3 ν m ν α L = mν m ν ρ ν ν m ν β L = mν ρ ν Rankine-Hugoniot conditions on the right interface: γ L 4.4 +γ L +pρ ν pρ L C L ; 4.5 ν + ρ R ρ ν + = β R m ν ν + α R m ν + = γ R mν + mν + ρ ν + ν + β R m ν + = γ R mν + ρ ν pρ R pρ ν ++ C R ; 4.8 0

11 admissibility condition on the left interface: ν ρ ν eρ ν ρ L eρ L +ν m ν ρ ν C L [ ρ ν eρ ν +pρ ν mν ρ ν ρ L eρ L +pρ L β L ρ L ] + m ν m ν ρ ν β L admissibility condition on the right interface: ν + ρ R eρ R ρ ν +eρ ν C R ++ν + m ν + ρ ν + [ ρ R eρ R +pρ R β ] R ρ ν ρ +eρ ν ++pρ ν + mν + C R R ρ ν + β R m ν + + ρ R C L ρ L m ν + ρ ν + ; Finally we write down the conditions mentioned in the point iv above. We recall that we search for subsolution independent of x, therefore we have from the subsolution condition 3.6 transforms to t ρ + x β = 0 4. t β + x pρ + C γ = 0, 4. α +β < ρ C 4.3 C α C +γ ρ β γ γ αβ > 0, ρ ρ 4.4 we rewrite 3.9 as ρ eρ +pρ + C t,x = ρ t,x Kt 4.5 and we can formulate 3.0 as t ρ eρ + C + x ρ eρ +pρ + C β ρ 4.3 Simplifications and ansatz We continue with several observations which were already introduced in previous works, for more details see [7, Section 4]. First, our initial data ρ 0,m 0 will be chosen in such a way that m 0 = v0 = m0 + = v This directly implies m = m + 0 and we will therefore look for subsolution with similar property, namely α 0 and γ 0. This choice implies that 4.4 and 4.7 are trivially satisfied and simplify to β < ρ C 4.7 C +γ C β γ > ρ

12 The necessary condition for to be satisfied is C γ > β ρ see [7, Lemma 4.3], which motivates us to introduce and further simplify to ε := C γ β ρ 4.9 ε := C β ρ ε 4.30 ε > ε > Now we rewrite the remaining equations and inequalities using ε, ε instead of C and γ and plugging in also our choice of the pressure law pρ = ρ, which yields eρ = ρ. We get Rankine-Hugoniot conditions on the left interface: ν ρ ν ρ L = m ν β L 4.33 ν m ν β L = mν ρ ν Rankine-Hugoniot conditions on the right interface: +ρ ν β L ρ L ρ L ε L; 4.34 ν + ρ R ρ ν + = β R m ν ν + β R m ν + = β R +ρ R ρ +ε R mν + R ρ ν ρ ν + ; admissibility condition on the left interface: ν ρ ν ρ L+ν m ν ρ ν ρ ν m ν ρ L β L + m ν 3 β L ε L +ε L ρ L ρ ν β3 L ρ L β Lε L +ε L ρ L admissibility condition on the right interface: β ν + ρ R ρν + +ν R + + ε R +ε R mν +R ρ R ρ ν + ρ R β R ρ ν +m ν β R ρ + β Rε R +ε R mν + 3 R ρ R ρ ν + differential equations in P : ; 4.37 ; 4.38 t ρ + x β = t β + x ρ + β +ε = 0; 4.40 ρ admissibility in P : ρ + β +ε +ε t,x = ρ t,x Kt 4.4 ρ t ρ + β +ε +ε +Kt x β ρ

13 Using 4.4 we moreover rewrite 4.40 and 4.4 to t β + x ρ K 3ρ ε = 0; 4.43 t ρ K ρ +Kt x β At this point we introduce the following ansatz for the part of the subsolution supported in the region P. We will look for the admissible generalized fan subsolution having the following properties. ρ is constant in the whole P ; K is constant in the whole P ; β depends only on t and is independent of x. With this ansatz, 4.39 is satisfied trivially, 4.44 is satisfied trivially as an equation and 4.43 simplifies into t β x ε = Since we assume that β is independent of x, the same has to hold also for x ε, which implies that ε has to be linear in the x -variable and 4.45 therefore becomes where l = lt = ν + t ν t = t β = ε R ε L, 4.46 l t 0 ν + t ν tdt > Moreover we rewrite once again, this time using 4.4 to replace ε by K. Note that as a consequence of our ansatz, ε remains in this set of equaitons and inequalities the only function defined in P with different values of ε L and ε R. In all other functions we can skip the subscripts L and R, since these functions are constant for fixed t. We obtain Rankine-Hugoniot conditions on the left interface: ν ρ ν ρ = m ν β 4.48 ν m ν β = mν ρ ν +ρ ν ρ K +3ρ +ε L; 4.49 Rankine-Hugoniot conditions on the right interface: ν + ρ ρ ν + = β m ν ν + β m ν + = ρ K 3ρ ε R mν + ρ ν + ρ ν + ; 4.5 Admissibility condition on the left interface: ν ρ ν +ρ + mν ρ ν ρ K ρ ν mν + mν 3 ρ ν βk; 4.5 Admissibility condition on the right interface: ν + ρ K ρ ρ ν + mν +R ρ ν βk ρ ν+m ν+ mν ρ ν

14 4.4 Solution for Riemann initial data As it was already shown in [6], there exist infinitely many admissible weak solutions starting from Riemann initial data generated by a compression wave. We will use the same result as a starting point, however we present a proof of existence of little bit different solutions than those proved to exist in [6]. Here we assume that the initial data have the form ρ 0,v 0 = { R,V for x < 0 R +,V + for x > Lemma 4.. Let pρ = ρ. Let ρ 0,v 0 have the form 4.54 with R =, R + = 4, V = 0, 8, V + = 0,0. Then there exist infinitely many bounded admissible weak solutions to the Euler equations.-.3. Proof. First, it is easy to observe, that such initial data are generated by a compression wave, see [6, Lemma 6.]. We will show that there exists a piecewise constant admissible fan subsolution which yields the existence of infinitely many bounded admissible weak solutions using Proposition 3.6. Such subsolution is given by the following set of numbers ν = ν + = ρ = β = ε = ε = K = It is a matter of computation to check that equations are satisfied and inequalities are satisfied as strict inequalities. The equation 4.46 is of course satisfied trivially since we work now with piecewise constant functions. 4.5 Generalization to nonconstant ε Let us recall again that the initial data we consider at this point in the proof have the form and are generated by a smooth compression wave which approximates the compression wave from Lemma 4.. This means that the Riemann invariant w defined in.5 is constant the value of this constant is actually w = 4 and the initial characteristic speed λ 0 x defined in.4 has the form In particular taking ζ > 0 small enough, we know that the values of ρ 0 ±,m0 ± are close to R ±,M ±. Therefore the same property holds also for ρ ν ±,m ν ±, since the values of ρ ±,m ± are propagated along characteristics of the Burgers equation We proceed with the construction of the subsolution as follows. We fix the values ρ = K = For fixed t > 0 we solve the equations to obtain ν ±, β and ε L as a functions of the parameter ε := ε L ε R. We denote M ± = R ± V ±. 4

15 We denote R = ρ ν ρ ν R = R R + = A = m ν m ν A = M M + = H = mν ρ ν mν + ρ ν + B = A RH = ρ ν ρ ν + +ρ ν ρ ν m ν ρ ν mν + ρ ν + ρ ν ρ ν +ρ ν ρ ν M B = R R + M + R R + R R R R + = and we emphasize that for the initial conditions under consideration R, A, B takes values in the neighborhood of R,A,B, i.e. we know that R R δ,r+ δ A A δ,a+ δ B B δ,b + δ 4.64 for some small δ = δζ > 0. We observe that B Rε < 0 provided ε < ε and therefore we express the solutions to as functions of the values ρ ν ±,mν ± and parameter ε : ν = A R + B +Rε ρν + R ρ ν 4.65 ν + = A R R B +Rε ρν ρ ν + β = mν ρ ν + m ν + ρ ν R ε L = ν ρν ρν mν ρ ν + R 4.66 B +Rε ρ ν ρν We also note that due to continuous dependence of various formulas on data the fact that the admissibility inequalities and subsolution inequalities are satisfied as strict inequalities for values R ±,M ± and ε = 0 implies that these inequalities will be satisfied also for ρ ν ±,m ν ± and ε > 0 provided δ and ε are sufficiently small. 4.6 Solution to 4.46 We have to construct the subsolution in such a way that the equation 4.46, which can be written as d dt βt = ε t lt, 4.69 is satisfied. In order to solve this equation we choose to emphasize the dependence of all functions on ε, hence obtaining βt = βt,ε t lt = lt,ε t To simplify notation, in the rest of this section we skip writing bars over β and l and treat them as functions of two variables, t and ε. In particular we have d dt βt,ε t = t βt,ε t+ ε βt,ε t d dt ε t 4.7 5

16 and the equation 4.69 becomes ε βt,ε t d dt ε t = ε t lt,ε t tβt,ε t. 4.7 We complement 4.7 with the natural initial condition We denote ε 0 = ft,ε t := ε βt,ε t 4.74 gt,ε t := t βt,ε t Note that the existence of solution to 4.7 cannot be solved by the Picard-Lindelöf theorem because in particular the function gt,ε develops a singularity at t = 0. Below we provide an existence proof based on the Banach fixed point argument with a special emphasis put upon the possible singularities of gt,ε. For T > 0 we consider a Banach space X T = {h C[0,T],h0 = 0} with a norm h XT = sup t [0,T] ft. For r > 0 we denote B XT r = {h X T, h XT r}. We consider an operator F on X T defined as Fδ = ε where ε is a solution to ft,δ t d dt ε t = ε t lt,δ t gt,δ t, ε 0 = The solution to 4.76 can be written as ε t = t 0 { gτ,δ τf τ,δ τexp as far as the integral on the right hand side exists. t τ } l s,δ sf s,δ sds dτ Lemma 4.3. Let there be a constant c > 0 such that gt,δ t < c t for some α <. Then α there exists T > 0 such that F : B XT B XT. B R B R Proof. It suffices to show that the integral on the right hand side of 4.77 is finite for, say, t = and that it converges to 0 as t 0. It is easy to observe that for δ B B XT R there exist constants c > 0 and c > 0 such that lt,δ t c t,c t. Similarly expressing ft,δ t from 4.67 we have ρ ν ft,δ t = ρν + B +Rδ ρ ν ρν and thus ft,δ t c,c. Thus, we have from 4.77 t ε t c 0 { t τ α exp c τ } c sds dτ = c t c 0 τ c τ α t c c dτ = c c + α t α It immediately follows that ε 0 as t 0. Furthermore, we may choose T > 0 in order to have c c c c + α t α c c + α T α B R. 6

17 Lemma 4.4. Let the assumptions of Lemma 4.3 be satisfied and let there be a constant C > 0 such that gt,δ t gt,δ t C δ δ L 0,T t α Then there exists T 0,T ] such that F : B R for all t < T and δ,δ B X T B B XT R B B XT R is a contraction. Proof. Take δ,δ B X T B R h i = ht,δ i for h = f,g,l and i =,. We have d dt ε ε = ε ε l f ε and let ε = Fδ and ε = Fδ. We also use a notation l f l f g f g f. 4.8 We use a notation L = l f, L = l f, G = g f and G = g f. From 4.8 we deduce t ε ε t = ε L L G G t τexp L sds dτ We have ε t ε t t 0 0 ε τ L L τ+ G G τ τ c We will estimate the differences L L and G G in terms of δ δ. First, it is easy to observe from 4.78 Similarly we have and therefore lt,δ t = τ t c ft,δ t ft,δ t C δ δ L 0,T. 4.8 t Plugging together 4.8 and 4.84 we get and together with B +Rδ ρν + R ρ ν + B +Rδ ρν ρ ν + dτ 4.83 lt,δ t lt,δ t C δ δ L 0,T t L τ L τ f f l f f + l l f l l C δ δ L 0,T τ 4.85 ε τ L τ L τ C δ δ L 0,T τα Using the assumption 4.80 we get also G τ G τ g f f f f + g g f C δ δ L 0,T τ α 4.87 Alltogether we end up with ε t ε t C δ δ L 0,T t 0 t c Now it is sufficient to choose T 0,T ] such that C T α <. τ c α dτ C δ δ L 0,T t α. As a consequence of the previous two lemmas and the Banach fixed point theorem we get Corollary 4.5. Under the assumptions of Lemmas 4.3 and 4.4 there exists T > 0 and a solution ε C0,T to such that for some positive constant C. ε t Ct α

18 4.7 Construction of the initial data We start the final part of the proof of Theorem. with the following lemma. Lemma 4.6. Let λ t,x be a solution to the Burgers equation.6-.7 with initial data λ 0 x in the form Let ν < s < s + < ν + and ν < ds dt < ds+ dt < ν +, where s ± is given by Lemma 4., and let λ t, ν t λ +ξ λ t, ν + t λ + ξ lim = c t 0+ t lim = c t 0+ t for some positive constants c ±. Then lim x 0 for some positive constants c ±. f x λ +ξ x = c lim x 0+ f + x λ + ξ x = c Proof. Let t > 0. We show the proof for ν t and f x. The solution to the Burgers equation.6-.7 is given by the method of characteristics. The characteristic lines are given by ft,x = x+tλ 0 x 4.9 and the solution is given by We introduce the inverse function g as λ t,ft,x = λ 0 x. 4.9 ft,gt,y = y 4.93 and note that this function is well defined in points y = ν ± t due to the assumptions of Lemma 4.6. From 4.93 we obtain y gt,y = x ft,gt,y = +tλ 0 gt,y 4.94 and in particular Similarly we obtain y gt, ν t = +tλ 0 gt, ν t 4.95 λ 0 t gt, ν t = gt, ν t +tλ 0 gt, ν t Denoting λ t := λ t, ν t we have λ t = λ 0 gt, ν t 4.97 and therefore it is enough to prove that the mapping gt, ν t behaves like a linear function in the neighborhood of zero. We thus denote and using 4.95 and 4.96 we get Gt := gt, ν t 4.98 d dt Gt = ν t λ 0 Gt +tλ 0 Gt

19 Moreover we also have d t = dt λ d dt λ0 Gt = λ0 Gt t gt, ν t+ν t y gt, ν t λ 0 Gt ν t λ 0 Gt = +tλ 0, Gt 4.00 which allows us to express λ 0 Gt = Plugging in 4.0 into 4.99 we end up with d dt λ t ν t λ 0 Gt t d dt λ t. 4.0 d dt Gt = ν t λ 0 Gt t d dt λ t. 4.0 We note that close to t = 0, λ 0 Gt = λ t takes values close to λ ξ, whereas due to the assumptions ν t < d dt s t = λ +λ + ξ and therefore ν t λ 0 Gt is bounded away from zero. Finally, due to 4.89 we have near t = 0 and thus Gt behaves like a linear function near t = 0. t d dt λ t C t 4.03 The final argument of our proof consists of using Lemma 4.6 as follows. We first prescribe functions ρ ν ±,mν ± t in a suitable way, such that 4.64 is satisfied. With these functions given, we solve the Rankine-Hugoniot conditions together with the differential equation 4.46 as it is described in Sections 4.5 and 4.6, in particular we obtain the curves ν ± t defining the generalized fan partition by 4.65 and Then we use the method of characteristics to map the given functions ρ ν ±,m ν ±t to the initial data ρ 0 ±,m 0 ±x and Lemma 4.6 provides a proof that the behaviour of ρ 0 ±,m0 ± around zero is the same as the behaviour of the prescribed functions ρ ν ±,m ν ± around zero. We use the relation between the solution λ to the Burgers equation.6 and the solution ρ, m to the Euler equations Recalling.4 and.5 we get ρ = w λ 8 m = w λ λ w We know that ρ,m is a solution to the Euler system if w is constant and λ is a continuous solution to the Burgers equation. Using 4.04 we also easily observe that if λ t, ν ± t λ 0±+c ± t then ρ ν ± ρ ±0±+C ± t and m ν ± ρ ±0±+C ± t and vice versa, if ρ ν ±,mν ± behave like t then so does λ t, ν ± t. This motivates us to prescribe at this point ρ ν ± t as follows. +ξ Ct / for 0 < t < δ ρ ν t = smooth for δ < t < δ 4.05 for t > δ and 4 ξ +Ct / for 0 < t < δ ρ ν + t = smooth for δ < t < δ for t > δ 9

20 for some positive constant C and small ξ, δ and δ > δ, such that ρ ν ± are monotone with monotone first derivatives. We set w 4, so m ν ± have to be defined as m ν ± t = 8ρ ν ± ρ t ν± t This choice of functions ρ ν ± and m ν ± ensures that there exists a constant c > 0 such that t ρ ν + tρ ν + + tm ν + tm ν + c It is then not very difficult to conclude from 4.67 and 4.75 that both gt,δ t C t α and 4.80 are satisfied with α = and thus, as stated in Corollary 4.5, there exists T and a solution ε C0,T to This means, that on 0,T we constructed functions ν ± t, βt, ε L t and ε t with ε R t = ε L t ε t satisfying 4.46 together with the Rankine-Hugoniot conditions with values of ρ and K given by Moreover, using the continuity argument, we know that the admissibility inequalities and subsolution inequalities are satisfied as strict inequalities for ε = 0 and values of ρ ν ±,m ν ± close to R ±,M ± and therefore there exists time δ 0 > 0 such that all these inequalities are satisfied as well with ε t constructed above. Therefore we constructed an admissible generalized fan subsolution on time interval 0,δ 0 with the initial data given by mapping the values of ρ ν ±,m ν ± from points x = ν ± t to t = 0 using the method of characteristics for the Burgers equation. Lemma 4.6 provides a proof that properties 4. are satisfied. Theorem. is proved. Acknowledgment O. Kreml, V. Mácha and S. Schwarzacher were supported by the GAČR Czech Science Foundation project GJ7-0694Y in the general framework of RVO: E. Chiodaroli was supported by the Italian National Grant FFABR 07. References [] Buckmaster, T.: Onsager s conjecture almost everywhere in time. Communications in Mathematical Physics, [] Buckmaster, T., De Lellis, C., Isett, P., Székelyhidi, L. Jr.: Anomalous dissipation for /5- holder Euler flows. Annals of Mathematics, [3] Buckmaster, T., De Lellis, C., Székelyhidi, L. Jr.: Dissipative Euler flows with Onsager-critical spatial regularity. Comm. Pure Appl. Math., [4] Buckmaster, T., De Lellis, C., Székelyhidi, L. Jr., Vicol, V.: Onsager s conjecture for admissible weak solutions. Comm. Pure Appl. Math., [5] Chiodaroli, E.: A counterexample to well posedness of entropy solutions to the compressible Euler system, J. Hyperbolic Differ. Equ [6] Chiodaroli, E., De Lellis, C., Kreml, O.: Global ill-posedness of the isentropic system of gas dynamics, Comm. Pure Appl. Math [7] Chiodaroli, E., Kreml, O.: On the energy dissipation rate of solutions to the compressible isentropic Euler system, Arch. Rational Mech. Anal [8] Constantin, P., E, W., Titi, E. S.: Onsager s conjecture on the energy conservation for solutions of Euler s equation, Comm. Math. Phys., t 0

21 [9] Dafermos, C. M.: Hyperbolic conservation laws in continuum physics, vol. 35 of Grundleheren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Third edition. Springer-Verlag, Berlin 00. [0] Daneri, S., Székelyhidi, L. Jr.: Non-uniqueness and h-principle for Hölder-continuous weak solutions of the Euler equations, Archive for Rational Mechanics and Analysis, [] De Lellis, C., Székelyhidi, L.: The Euler equations as a differential inclusion, Ann. of Math [] De Lellis, C., Székelyhidi, L.: On admissibility criteria for weak solutions of the Euler equations, Arch. Rational Mech. Anal [3] De Lellis, C., Székelyhidi, L. Jr.: Dissipative continuous Euler flows. Invent. Math., [4] Eyink, G. L.: Energy dissipation without viscosity in ideal hydrodynamics I. Fourier analysis and local energy transfer, Phys. D, [5] Feireisl, E.: Maximal dissipation and well-posedness for the compressible Euler system, J. Math. Fluid Mech [6] Feireisl, E., Gwaizda, P., Świerczewska-Gwiazda, A., Wiedemann, E.: Regularity and energy conservation for the compressible Euler equations., Arch. Ration. Mech. Anal [7] Isett, P.: A proof of Onsager s conjecture. arxiv: To appear in Ann. of Math., 06. [8] Isett, P.: Nonuniqueness and existence of continuous, globally dissipative Euler flows. ArXiv e-prints, October 07. [9] Székelyhidi, L.: Weak solutions to the incompressible Euler equations with vortex sheet initial data, C. R. Math. Acad. Sci. Paris