The 2-d isentropic compressible Euler equations may have infinitely many solutions which conserve energy

Transcription

1 The -d isentropic compressible Euler equations may have infinitely many solutions which conserve energy Simon Markfelder Christian Klingenberg September 15, 017 Dept. of Mathematics, Würzburg University, Germany Abstract We consider the D isentropic compressible Euler equations. It was shown in [] that there exist Riemann initial data as well as Lipschitz initial data for which there exist infinitely many weak solutions that fulfill an energy inequality. In this note we will prove that there is Riemann initial data for which there exist infinitely many weak solutions that conserve energy, i.e. they fulfill an energy equality. As in the aforementioned paper we will also show that there even exists Lipschitz initial data with the same property. Contents 1 Introduction and main result Proof of theorem Definitions Sufficient condition for non-uniqueness Algebraic equations and inequalities Proof Proof of corollary Conclusion 10 1

2 1 Introduction and main result For hyperbolic conservation laws the quest for showing that the initial value problem is well-posed is more than a century old. The focus has been on the compressible Euler equations in particular. In one space dimensions the seminal work of Glimm [7] and DiPerna [4] gave rise to the hope that this goal might be achievable. In the two dimensional case progress has been made to this end for particular self similar solutions, see Chen and Feldman [1]. In all these cases an admissibility condition, which mimics the nd law of thermodynamics, needs to be invoked. Thus results for the equations of two-dimensional isentropic compressible gas dynamics, that prove that with these admissibility conditions the initial value problem may lead to infinitely many admissible weak solutions see [6] have been met with quite a surprise. Even in situations, where one can solve the initial value problem explicitly via an admissible solution, one could prove that in addition to this standard solution infinitely many other admissible weak solutions exist [], [3] and [8]. The hunt is open for finding new admissibility conditions proving this line of research a fluke. The admissibility criterion used so far is entropy dissipation. For the isentropic Euler equations energy takes the role of entropy, so this means that one was seeking weak solutions that dissipate energy. So could one rid oneself of these infinitely many extra solution by demanding the more physical criterion that the solutions conserve energy? In this paper we show that this is not the case. We can show that there exist Lipschitz initial data that may lead to infinitely many solutions of the two dimensional Euler equations all of which conserve energy. To this end we consider the Cauchy problem for the -d isentropic compressible Euler system, i.e. t ϱ + div x ϱ v = 0, t ϱ v + div x ϱ v v + x pϱ = 0, ϱ0, x = ϱ 0 x, v0, x = v 0 x where the unknowns, more precisely the density ϱ = ϱt, x R + and the velocity v = vt, x R, are functions of time t [0, and position x = x 1, x R. We will focus both on Riemann and on Lipschitz initial data ϱ 0, v 0, see below. Furthermore we consider the polytropic pressure law, more specifically pϱ = ϱ. However the results can be extended to more general pressure laws 1 pϱ = K ϱ γ, where K > 0 and γ 1 are constants. In this paper the aim is to find energy conserving weak solutions, which are defined as follows. Definition 1. energy conserving weak solution A weak solution to 1 is energy conserving if the energy equality [ ] 1 1 t ϱ v + P ϱ + div x ϱ v + P ϱ + pϱ v = 0, 3

3 holds in the weak sense, where the pressure potential is given by P ϱ = ϱ. In the case of more general pressures the pressure potential reads P ϱ = K γ 1 ϱγ in the case γ > 1 and P ϱ = K ϱ logϱ for γ = 1. In the first part of the paper we consider Riemann initial data, i.e. { ϱ, v ϱ, v0, x = ϱ 0, v 0 x := if x < 0 ϱ +, v + if x > 0, 4 where ϱ ± R + and v ± R are constants. In the past the system 1 with initial data 4 was often discussed concerning uniqueness of so-called admissible weak solutions defined as follows. Definition. admissible weak solution A weak solution to 1 is admissible if the energy inequality [ ] 1 1 t ϱ v + P ϱ + div x ϱ v + P ϱ + pϱ v 0. 5 holds in the weak sense. It was shown by exploiting the convex integration method [5], [6] that for some initial states ϱ ±, v ± R + R there exist infinitely many admissible weak solutions, see [3], [] and [8]. These infinitely many solutions that are produced by the convex integration method are usually called wild solutions. In this note we want to show that there are initial states ϱ ±, v ± R + R for which there are infinitely many weak solutions that fulfill the energy equation 3 rather than the energy inequality 5. In other words we will prove that there are not only infinitely many admissible weak solutions but also infinitely many energy conserving weak solutions. The following theorem is our main result: Theorem 3. Let pϱ = ϱ. There exist initial states ϱ ±, v ± R + R such that there are infinitely many energy conserving weak solutions to 1, 4. This theorem will be proved by methods first presented in [3] and []. As a consequence of theorem 3 we can show the following: Corollary 4. Let pϱ = ϱ. There exist Lipschitz continuous initial data ϱ 0, v 0 such that there are infinitely many energy conserving weak solutions to 1. Proof of theorem 3.1 Definitions We proceed as in [3] and therefore recall the definition of a fan partition. 3

4 Definition 5. fan partition, see [3, Definition 4] Let µ 0 < µ 1 real numbers. A fan partition of 0, R consists of three open sets Ω, Ω 1, Ω + of the form Ω = {t, x : t > 0 and x < µ 0 t}, Ω 1 = {t, x : t > 0 and µ 0 t < x < µ 1 t}, Ω + = {t, x : t > 0 and x > µ 1 t}. Furthermore we define S := {M R M symmetric} and S 0 := {M S trm = 0}. We also recall the definition of a fan subsolution, where we slightly adjust [3, Defnitions 5 and 6] for our needs. More precisely we have an equality in 6 see below in contrast to the inequality in [3, Definition 6]. Definition 6. energy conserving fan subsolution, cf. [3, Definitions 5 and 6] An energy conserving fan subsolution to the Euler system 1 with initial condition 4 is a triple ϱ, v, u : 0, R R + R S0 of piecewise constant functions, which satisfies the following properties: 1. There exists a fan partition of 0, R and constants ϱ 1 R +, v 1 R and u 1 S0, such that ϱ, v, u = ϱ i, v i, v i v i v i Id 1 Ωi + ϱ 1, v 1, u 1 1 Ω1, i {,+} where ϱ ±, v ± are the given initial states.. There is a constant C 1 R + such that in the sense of definiteness. v 1 v 1 u 1 < C 1 Id 3. For all test functions ψ, φ Cc [0, R, R R the following identities hold: ϱ t ψ + ϱ v x ψ dx dt + ϱ 0 x ψ0, x dx = 0, 0 R R [ v ϱ v t φ + ϱ v 1Ω Ω + + u 1 1 Ω1 : D x φ 0 R + pϱ + 1 ] ϱ 1 C 1 1 Ω1 div x φ dx dt + ϱ 0 x v 0 x φ0, x dx = 0. R 4

5 4. For every non-negative test function ϕ Cc [0, R, R + 0 the equation [ P ϱ R ϱ v 1 Ω Ω + + C 1 1 Ω1 t ϕ + P ϱ + pϱ + 1 ] ϱ v 1 Ω Ω + + C 1 1 Ω1 v x ϕ dx dt is fulfilled. + R P ϱ 0 x + ϱ 0 x v 0x ϕ0, x dx = 0 6. Sufficient condition for non-uniqueness We obtain a slightly adjusted version of [3, Proposition 3.1]. Theorem 7. cf. [3, Proposition 3.1] Let ϱ ±, v ± be such that there exists an energy conserving fan subsolution ϱ, v, u to 1, 4. Then there are infinitely many energy conserving weak solutions ϱ, v to 1, 4 with the following properties: ϱ = ϱ, vt, x = vt, x for almost all t, x Ω Ω +, vt, x = C 1 for almost all t, x Ω 1. Proof. For the proof we refer to the proofs of [3, Proposition 3.1] and [, Proposition 3.6]. Note that the energy conservation 6 of the fan subsolution ensures that the energy equation 3 holds, i.e. that the weak solutions are energy conserving..3 Algebraic equations and inequalities As in [3, Proposition 4.1], definition 6 can be translated into a system of algebraic equations and inequalities. Note again that we get here energy equations instead of energy inequalities. Proposition 8. cf. [3, Proposition 4.1] Let ϱ ±, v ± be given. The constants µ 0, µ 1 R, ϱ 1 R +, v1 v 1 = 1 R u1, u v 1 = 11 u 1 1 S0 1 u 1 1 u 1 11 and C 1 R + define an energy conserving fan subsolution to 1, 4 if and only if they fulfill the following algebraic equations and inequalities: 5

6 Order of the speeds: µ 0 < µ 1 7 Rankine Hugoniot conditions on the left interface: µ 0 ϱ ϱ 1 = ϱ v ϱ 1 v 1 8 µ 0 ϱ v 1 ϱ 1 v 1 1 = ϱ v 1 v ϱ 1 u µ 0 ϱ v ϱ 1 v 1 = ϱ v + ϱ 1 u pϱ pϱ 1 ϱ 1 C 1 10 Rankine Hugoniot conditions on the right interface: µ 1 ϱ 1 ϱ + = ϱ 1 v 1 ϱ + v + 11 µ 1 ϱ 1 v 1 1 ϱ + v + 1 = ϱ 1 u 1 1 ϱ + v + 1 v + 1 µ 1 ϱ 1 v 1 ϱ + v + = ϱ 1 u 1 11 ϱ + v + + pϱ 1 pϱ + + ϱ 1 C 1 13 Subsolution condition: v1 1 + v1 < C 1 14 C1 C1 v u 1 11 v 1 u 1 11 u 1 1 v 1 1 v 1 > 0 15 Energy equation on the left interface: v C 1 µ 0 P ϱ + ϱ P ϱ 1 ϱ 1 = P ϱ + pϱ v P ϱ 1 + pϱ 1 v 1 + ϱ v v ϱ 1 v 1 C 1 16 Energy equation on the right interface: C 1 µ 1 P ϱ 1 + ϱ 1 P ϱ v + + ϱ + = P ϱ 1 + pϱ 1 v 1 P ϱ + + pϱ + C 1 v + + ϱ 1 v 1 ϱ v + + v + 17 Proof. The above proposition can be proved analogously to [, Proposition 5.1]. 6

7 .4 Proof Finally we prove theorem 3. Proof. Let the initial data be given by We set ϱ = 1 ϱ + = 4 0 v = 0 v + = 0 µ 0 = 7 ϱ 1 = u 1 = µ 1 = 0 0 v 1 = C 1 = Simple computations show that the equations and inequalities of proposition 8 hold. Hence there exists an energy conserving fan subsolution and according to theorem 7 infinitely many energy conserving weak solutions to 1 with initial data 4, 18. The reader might wonder where the values in 19 come from. Hence we want to expose the way we reached to them. Let us begin by looking for weak solutions to 1 with initial data 4, 18 that are admissible instead of energy conserving. This means that the energy inequality 5 holds rather than the energy equation 3. To this end, as shown in [], [3] and [8], one has to look for admissible fan subsolutions. Admissible fan subsolutions are similar to energy conserving fan subsolutions as defined in definition 6, with the only difference that the energy equation 6 turns into an energy inequality. One ends up with algebraic equations and inequalities, see e.g. [, Proposition 5.1], like those presented in proposition 8 with the difference that again the energy equations 16, 17 are replaced by energy inequalities. Since there are six equations for eight unknowns, the idea in [3] and also in [8] was to choose two unknowns as parameters and express the other six unknowns as functions of the two parameters. For convenience one replaces the unknowns u 1 11 and C 1 by δ 1 = C 1 v 1 u 1 11, δ = C 1 v u 1 11 and chooses ϱ 1, δ as parameters. For the special initial data 18 one obtains the following proposition, shown by the authors in [8]. 7

8 Proposition 9. see [8, Theorem 5.] There exists an admissible fan subsolution to 1 with initial data 4, 18 if there exist constants ϱ 1, δ R + that fulfill ϱ < ϱ 1 < ϱ +, 0 δ 1 ϱ 1 > 0, 1 v 1 ϱ 1 v pϱ + pϱ 1 ϱ 1 P ϱ ϱ P ϱ 1 ϱ ϱ 1 δ 1 ϱ 1 ϱ 1 v 1 ϱ 1 + v δ 1 ϱ 1 + δ ϱ ϱ 1 v 1 ϱ 1 v, ϱ ϱ 1 v + v 1 ϱ 1 pϱ 1 + pϱ + ϱ + P ϱ 1 ϱ 1 P ϱ + ϱ 1 ϱ + δ 1 ϱ 1 ϱ 1 v + + v 1 ϱ 1 + δ 1 ϱ 1 + δ ϱ 1 ϱ + v + v 1 ϱ 1, ϱ 1 ϱ + where we define the functions 1 v 1 ϱ 1 := ϱ v ϱ + ϱ 1 ϱ + v + ϱ 1 ϱ ϱ 1 ϱ ϱ + [ + ϱ ϱ + pϱ pϱ + ] ϱ + ϱ v v + ϱ 1 ϱ ϱ + ϱ and δ 1 ϱ 1 := pϱ 1 pϱ + ϱ ϱ 1 ϱ ϱ ϱ 1 ϱ 1 ϱ ϱ + + v v + [ + ϱ ϱ + pϱ pϱ + ] ϱ+ ϱ ϱ + ϱ v v + 1. ϱ 1 ϱ 5 Note that these functions are well-defined for ϱ < ϱ 1 < ϱ + and for the initial states 18. Here the conditions δ > 0 and 1 ensure the subsolution conditions, 0 guarantees the correct order of the speeds and the inequalities and 3 correspond to the energy inequalities for the left and the right interface, respectively see proof of [8, Theorem 5.]. Next we are going to plot the regions in the ϱ 1 -δ -plane where the inequalities 0-3 are fulfilled, see figure 1. Proposition 9 claims that each point ϱ, δ lying in the region shown in figure 1 d corresponds to an admissible fan subsolution. Let us now return to this paper s aim, namely to find energy conserving solutions. The procedure presented above will lead to a similar claim as stated in proposition 9 with the difference 8

9 δ 1.5 δ ρ 1 a Condition ρ 1 b Condition δ 1.5 δ ρ 1 c Condition ρ 1 d Condition 0-3 Figure 1: Regions in the ϱ 1 -δ -plane where the conditions of proposition 9 hold that the inequalities, 3 are replaced by equations. Therefore, in order to find an energy conserving fan subsolution, we have to find a point ϱ 1, δ in the region shown in figure 1 c which in addition lies on the boundary of the regions shown in figures 1 a and b. Hence it suffices to find the apex in figure 1 d. The values 19 correspond exactly to this apex. 3 Proof of corollary 4 Proof. In order to prove corollary 4 we proceed as in the proof of [, Corollary 1.]. We solve the system 1 with initial data 4, 18 backwards in time. It is an easy observation that if 9

10 ϱt, x, vt, x is a solution to 1 then ϱt, x, ṽt, x = ϱ t, x, v t, x solves 1, too. Hence in order to solve 1 with initial data 4, 18 backwards in time, it suffices to solve 1 forward in time, where now the initial states are switched { ϱ+, v ϱ, v0, x = ϱ 0, ṽ 0 x := ϱ 0, v 0 x = + if x < 0 ϱ, v if x > 0, 6 and ϱ ±, v ± given in 18. By well-known methods one can show that there is a self-similar solution ϱ, ṽ to 1, 6 which consists of one rarefaction, see [8, Proposition 1.3] and the references therein. Note that rarefaction solutions are Lipschitz in the spatial variable x and they conserve the energy in the sense of 3. Hence we can simply set the initial data to ϱ 0, v 0 x := ϱ, ṽt = 1, x, which is Lipschitz and leads to Riemann data 4, 18 for t = 1. Consequently according to theorem 3 there exist infinitely many energy conserving weak solutions to 1 with Lipschitz initial data as above. All these solutions coincide for t [0, 1]. 4 Conclusion In this note we showed that even if we insist that solutions conserve energy for the two-dimensional isentropic Euler equations one can find Lipschitz initial data which gives rise to infinitely many weak solutions. The proof makes use of the techniques shown in [8]. Thus if one wants to find a criterion to rule out these solutions one has to look elsewhere. Acknowledgement Both authors were partially funded by the European Research Council under the European Unions Seventh Framework Programme FP7/ / ERC Grant Agreement We thank Eduard Feireisl for hosting us. 10

11 References [1] Chen, Gui-Qiang; Chen, Jun; Feldman, Mikhail: Transonic flows with shocks past curved wedges for the full Euler equations. Discrete Contin. Dyn. Syst [] Chiodaroli, E., De Lellis, C., Kreml, O.: Global ill-posedness of the isentropic system of gas dynamics. Comm. Pure Appl. Math. 687, [3] Chiodaroli, E., Kreml, O.: On the energy dissipation rate of solutions to the compressible isentropic euler system. Arch. Ration. Mech. Anal. 143, [4] DiPerna, Ronald J. Convergence of approximate solutions to conservation laws. Archive for Rational Mechanics and Analysis :. [5] De Lellis, C., Székelyhidi Jr., L.: The euler equations as a differential inclusion. Ann. of Math. 1703, [6] De Lellis, C., Székelyhidi Jr., L.: On admissibility criteria for weak solutions of the euler equations. Arch. Ration. Mech. Anal. 1951, [7] Glimm, J.: Solutions in the large for nonlinear hyperbolic systems of equations. Communications on pure and applied mathematics [8] Markfelder, S., Klingenberg, C.: The Riemann problem for the multidimensional isentropic system of gas dynamics is ill-posed if it contains a shock. accepted in Arch. Ration. Mech. Anal., see also arxiv: