arxiv: v1 [math.ap] 29 May 2018

Transcription

1 Non-uniqueness of admissible weak solution to the Riemann problem for the full Euler system in D arxiv: v [math.ap] 9 May 08 Hind Al Baba Christian Klingenberg Ondřej Kreml Václav Mácha Simon Markfelder May 30, 08 Institute of Mathematics of the Czech Academy of Sciences Žitná 5, CZ-5 67 Praha, Czech Republic Department of Mathematics, Würzburg University Emil-Fischer-Str. 40, Würzburg, Germany Abstract: The question of well- and ill-posedness of entropy admissible solutions to the multidimensional systems of conservation laws has been studied recently in the case of isentropic Euler equations. In this context special initial data were considered, namely the D Riemann problem which is extended trivially to a second space dimension. It was shown that there exist infinitely many bounded entropy admissible weak solutions to such a D Riemann problem for isentropic Euler equations, if the initial data give rise to a D self-similar solution containing a shock. In this work we study such a D Riemann problem for the full Euler system in two space dimensions and prove the existence of infinitely many bounded entropy admissible weak solutions in the case that the Riemann initial data give rise to the D self-similar solution consisting of two shocks and possibly a contact discontinuity. Contents Introduction The D Riemann solution 4 3 Main results 5 4 Existence of infinitely many admissible weak solutions 7 4. Sufficient condition for non-uniqueness Algebraic equations and inequalities Proof of Theorem Reduction to special case Proof of Proposition Shock speeds Constants C i and γ i Subsolution and admissibility inequalities albaba@math.cas.cz klingen@mathematik.uni-wuerzburg.de kreml@math.cas.cz macha@math.cas.cz simon.markfelder@mathematik.uni-wuerzburg.de

2 6 Proof of Theorem The case R < Solution to algebraic equations Positivity of p, ε and ε Admissibility inequalities The case R > The case R = Solution to algebraic equations Positivity of p, ε and ε Admissibility inequalities Introduction In this paper we consider the full compressible Euler system in the whole two-dimensional space, i.e. t ϱ div x ϱv = 0, t ϱv div x ϱv v x p = 0, [ ] t ϱ v ϱ eϱ, p div x ϱ v ϱ eϱ, p p v = 0,. where the density ϱ = ϱt, x R, the pressure p = pt, x R and the velocity v = vt, x R are unknown functions of the time t [0, and the position x = x, x R. Furthermore we consider an ideal gas, i.e. eϱ, p = c v p ϱ, where c v > 0 is a constant called the specific heat at constant volume. We recall that from a physical point of view c v = f, where f is the number of degrees of freedom. The system of equations. is complemented with the initial condition ϱ, v, p0, x = ϱ 0, v 0, p 0 x in R.. Further, we supplement the equations. with the entropy condition t ϱ sϱ, p div x ϱ sϱ, pv 0,.3 where sϱ, p is the entropy. For the ideal gas we have p c v sϱ, p = log ϱ cv. By admissible weak solution - we will also use the term weak entropy solution - we mean a triple of L -functions ϱ, v, p that fulfill.-.3 in the sense of distributions, see Section 3 for precise definitions. In this paper we consider Riemann initial data { ϱ 0, v 0, p 0 ϱ, v x :=, p if x < 0 ϱ, v, p if x > 0,.4 where ϱ ± R, v ± R and p ± R are constants. We write v = v, v T for the components of the velocity v. Throughout this paper we additionally assume v, = v, = 0, where v ±, are the first components of the velocities v ±.

3 Following the groundbreaking results of De Lellis and Székelyhidi [8, 9] concerning ill-posedness of the incompressible Euler equations, the question of well- and ill-posedness of admissible weak solutions to the isentropic compressible Euler equations was raised. In fact, De Lellis and Székelyhidi themselves showed in [9] the existence of initial data ϱ 0, v 0 L R for which there exist infinitely many admissible weak solutions such initial data is called wild initial data. In the context of hyperbolic conservation laws this result in particular shows that the notion of admissible weak solution, i.e. solution satisfying the system of conservation laws and an inequality arising from the companion law, is probably not suitable for the problem in the case of more than one space dimension. This result was further improved by Chiodaroli [3] and Feireisl [0] who showed that the illposedness is related to the irregularity of the velocity field, proving that there exist wild initial data ϱ 0, v 0 with ϱ 0 C R. Finally, De Lellis, Chiodaroli and Kreml [4] showed the existence of Lipschitz wild initial data. The method of the proof was motivated by the observation of Székelyhidi [7] that for the incompressible Euler system, vortex sheet initial data are wild. Hence, the authors in [4] studied the Riemann problem for the isentropic Euler system and found an example of Riemann initial data allowing for existence of infinitely many admissible weak solutions while being generated by a backwards rarefaction wave i.e. a compression wave. Since the Riemann problem is a basic building block in the theory of one-dimensional systems of conservation laws, it was only natural to study its properties in multi-dimensional case further. In fact, Chen and Chen [] already proved that if the solution to the Riemann problem consists only of rarefaction waves, such solution is unique in the class of multi-dimensional admissible weak solutions. The same was later proved also in []. On the non-uniqueness side, Chiodaroli and Kreml [5] proved that whenever the Riemann initial data give rise to the self-similar solution consisting of two shocks, there exist also infinitely many other bounded admissible weak solutions. Moreover they also showed that the self-similar solution may not be entropy rate admissible, meaning that some of the other solutions constructed may produce more total entropy. The study of the Riemann problem for the isentropic Euler system continued with results of Chiodaroli and Kreml [6] and Klingenberg and Markfelder [4] who independently showed ill-posedness in the case of Riemann initial data giving rise to the self-similar solution consisting of one rarefaction wave and one shock. More precisely in [6] the authors need a certain smallness condition for the initial data, whereas the result in [4] covers all cases of such Riemann initial data and also solves the case of Riemann initial data giving rise to self-similar solution consisting of a single shock. Finally, Březina, Chiodaroli and Kreml [] showed that the above mentioned results of illposedness of admissible weak solutions hold also in the presence of a contact discontinuity in the self-similar solution arising from the multidimensionality of the problem. However, the question, whether a single contact discontinuity is a unique admissible weak solution to the isentropic Euler system or not, is still an open problem; it is worth mentioning that this is actually a direct analog of the vortex sheet initial data considered by Székelyhidi in [7]. In view of results for the isentropic system a natural question arises, whether these results can be transferred to the case of full Euler equations. It was shown by Chen and Chen in [] that uniqueness of rarefaction waves actually holds also for the full system, see also [3]. On the nonuniqueness side, we note the result of Markfelder and Klingenberg [5], who showed that for the isentropic system there exist infinitely many solutions conserving the energy. This however can be viewed as a proof of existence of infinitely many solutions to the full system with a proper choice of c v ; all these solutions also conserve the entropy. In addition to that we want to mention the result by Feireisl et al. [], who proved existence of wild initial data for the full Euler equations. Their solutions conserve the entropy and have piecewise constant density and pressure. In this paper we prove non-uniqueness of admissible weak solutions for a larger class of Riemann initial data, namely for initial data that give rise to self-similar solution consisting of two shocks and possibly a contact discontinuity. The paper is structured as follows. In Section we recall some results for the one-dimensional Riemann problem for the full Euler system and present a condition on the initial states ϱ ±, v ±, p ± such that the self-similar solution consists of two shocks and possibly a contact discontinuity. In 3

4 Section 3 we present our two main results, Theorem 3.5 and Theorem 3.6. In Section 4 we provide the technical tools used in the proof, in particular we present a key lemma relating the existence of infinitely many admissible weak solutions to the existence of a single so-called admissible fan subsolution. Finally, in Sections 6 and 5 we present proofs of our two main results. The D Riemann solution The initial data. of the problem we study are one-dimensional, we have ϱ 0, v 0, p 0 x only with no dependence on the x variable. Therefore it is reasonable to search for solutions which share this property, i.e. for solutions satisfying ϱ, v, pt, x = ϱ, v, pt, x. Moreover since we assume v, = v, = 0, we search for solutions with the property v t, x 0. Hence we obtain the D Riemann problem t ϱ ϱv = 0, t ϱv ϱv p = 0, [ ] t ϱv ϱ eϱ, p ϱv ϱ eϱ, p p v = 0,. for the unknown functions ϱ, v, pt, x with initial data { ϱ, v, p0, x = ϱ 0, v, 0 p 0 ϱ, v x :=,, p if x < 0 ϱ, v,, p if x > 0.. The problem.-. is a model example of a one-dimensional Riemann problem for a system of conservation laws and has been studied extensively by several authors, e.g. by Smoller [6, Chapter 8 B]. It is well known that this problem admits a self-similar solution i.e. solution consisting of functions of a single variable x t belonging to the class of BV functions. This solution to.-. is unique in the class of BV self-similar functions and consists of three waves connected by constant states, where the first wave is either an admissible shock or a rarefaction wave, the second wave is a contact discontinuity and the third wave is again either an admissible shock or a rarefaction wave. In the rest of this paper we will call this solution the D Riemann solution. Table shows all the possibilities of the structure of the D Riemann solution. -wave -wave 3-wave -wave -wave 3-wave contact shock - contact shock - - rarefaction - contact rarefaction shock - - shock contact - shock - shock shock contact shock shock - rarefaction shock contact rarefaction rarefaction - - rarefaction contact - rarefaction - shock rarefaction contact shock rarefaction - rarefaction rarefaction contact rarefaction Table : All the 8 possibilities of the structure of the D Riemann solution It is easy to observe that the D Riemann solution is an admissible weak solution to the original D problem.-.. In this paper we will focus on the cases, where the D Riemann solution contains two admissible shocks and possibly a contact discontinuity. Therefore we recall the following proposition stating conditions for the initial data under which such D Riemann solution emerges. 4

5 Proposition.. Assume that either p p and v, v, < p p p p and v, v, < p p c v ϱ or p c vp c v ϱ. p c vp Then the D Riemann solution to problem.-. consists of a -shock, possibly a -contact discontinuity and a 3-shock. The intermediate states ϱ M, v M, p M and ϱ M, v M, p M are given by p M is the unique solution of p M p c v ϱ p c v p M p M p = v, v, ; ϱ p c v p M p M p v M, = v, c v ; ϱ p c vp M ϱ M = ϱ p c vp M p M c vp ; ϱ M = ϱ p c vp M p M c vp ; and the following properties hold. If ϱ M = ϱ M then the -contact discontinuity does not appear. p M > max{p, p }, ϱ M > ϱ and ϱ M > ϱ. The shock speeds denoted by σ, σ satisfy σ ± ϱ ± ϱ M± = ϱ ± v ±, ϱ M± v M,. The speed of the contact discontinuity is given by v M,. For the proof of Proposition. we refer to [6]. More details about self-similar solutions to D systems of conservation laws can be also found in [7]. 3 Main results In this section we state the definition of solutions to the problem.-. and our main results. Definition 3.. Let ϱ 0, v 0, p 0 L R. The triplet ρ, v, p L 0, R ; R R R is called a weak solution of.-. if the following system of identities is satisfied ϱ t ψ ϱv x ψ dx dt ϱ 0 xψ0, x dx = R R for all test functions ψ Cc [0, R, ϱv t ϕ ϱv v : x ϕ p div x ϕ dx dt ϱ 0 xv 0 x ϕ0, x dx = R R for all test functions ϕ Cc [0, R ; R, [ 0 R ϱ v ϱ eϱ, p t φ for all test functions φ C c [0, R. R ϱ v ϱ eϱ, p p ϱ0 x v 0 x ϱ 0 xeϱ 0 x, p 0 x ] v x φ dx dt φ0, x dx =

6 Definition 3.. Let ϱ 0, v 0, p 0 L R. A weak solution of.-. is called admissible or weak entropy solution, if the following inequality is satisfied for all non-negative test functions ϕ Cc [0, R ; R 0 : 0 R ϱ sϱ, p t ϕ ϱ sϱ, pv x ϕ dx dt R ϱ 0 xsϱ 0 x, p 0 xϕ0, x dx The solutions which we are going to construct will have a special structure, partially motivated by the structure of the D Riemann solution. Therefore we introduce first the notion of N-fan partition N =, and then N-fan solutions. Definition 3.3. Let N {, } and µ 0 < µ µ 0 < µ < µ real numbers. An N-fan partition of 0, R is a set of N open sets Ω, Ω, Ω Ω, Ω, Ω, Ω of the form Ω = {t, x : t > 0 and x < µ 0 t}; Ω i = {t, x : t > 0 and µ i t < x < µ i t}; Ω = {t, x : t > 0 and x > µ N t}. Definition 3.4. An admissible weak solution ϱ, v, p of the full Euler system.-. with Riemann initial data.4 is called an N-fan solution, if there exists N {, } and an N-fan partition of 0, R such that ϱ, p are constant on each set Ω,Ω i, Ω of the N-fan partition; v is constant on each set Ω,Ω i, Ω of the N-fan partition; ϱ, v, p = ϱ ±, v ±, p ± on the sets Ω ±. See also Fig.. t t ϱ, p = ϱ, p = const. ϱ, p = ϱ, p = const. ϱ, p = ϱ, p = const. ϱ, v, p = ϱ, v, p ϱ, v, p = ϱ, v, p ϱ, v, p = ϱ, v, p ϱ, v, p = ϱ, v, p x x a -fan solution b -fan solution Figure : Structure of an N-fan solution. On the leftmost and rightmost sets i. e. on Ω ±, the solution coincides with the initial states. On each set in the middle Ω i, i {, N}, the density ϱ and pressure p are constant, whereas the velocity v does not need to be constant. Now we are ready to state our main theorems. Theorem 3.5. Let c v > solution consists either and let the Riemann initial data.4 be such that the D Riemann of a -shock and a 3-shock or of a -shock, a -contact discontinuity and a 3-shock. Then there exist infinitely many admissible weak solutions to.-.4. These solutions are all -fan solutions. 6

7 Theorem 3.6. Let c v > 0. Let ϱ ±, p ± > 0 and V R be given. i If ϱ < ϱ or ϱ = ϱ and p > p let v, = V. Then there exists U = Uϱ ±, p ±, V such that for all v, < U there exist infinitely many admissible weak solutions to.-.4. These solutions are all -fan solutions. ii If ϱ > ϱ or ϱ = ϱ and p < p let v, = V. Then there exists U = Uϱ ±, p ±, V such that for all v, > U there exist infinitely many admissible weak solutions to.-.4. These solutions are all -fan solutions. Remark 3.7. Theorem 3.5 is the most general result presented here, covering all the cases of Riemann initial data giving rise to the self-similar solution containing two shocks. Solutions constructed in the proof of this Theorem are all -fan solutions. It might be of interest that at least if the difference between the velocities v, and v, is large enough, one can construct also another type of solutions, namely -fan solutions, this is stated in Theorem 3.6. We emphasize that the set of initial data in Theorem 3.6 is a strict subset of the set of initial data in Theorem 3.5. Finally, we refer the reader to Remark 6.7 for the reason why Theorem 3.6 does not cover the case ϱ = ϱ and p = p. Remark 3.8. Note that the general Theorem 3.5 requires the specific heat at constant volume coefficient c v to be strictly larger than, whereas the results of Theorem 3.6 hold for all positive c v. 4 Existence of infinitely many admissible weak solutions The content of this section shows the method to prove non-uniqueness of admissible weak solutions. As in the isentropic case see [4] and following works the idea is to work with a so called admissible fan subsolution - here a quintuple of piecewise constant functions satisfying certain system of partial differential equations. Then, Theorem 4. tells us that the existence of such an admissible fan subsolution implies existence of infinitely many weak entropy solutions to.-.4. In what follows we denote by S0 the space of symmetric traceless matrices and I denotes the identity matrix. Definition 4.. An admissible N-fan subsolution to the Euler system.-.3 with Riemann initial data.4 is a quintuple ϱ, v, U, C, p : 0, R R R S0 R R of piecewise constant functions, which satisfies the following properties:. There exists an N-fan partition Ω, Ω, Ω or Ω, Ω, Ω, Ω of 0, R and for i {, N} there exist constants ϱ i R, v i R, U i S0, C i R and p i R, such that ϱ, v, U, C, p = i {,} ϱ i, v i, v i v i v i N I, v i, p i Ωi ϱ i, v i, U i, C i, p i Ωi, where ϱ ±, v ±, p ± are the constants given by the initial condition.4.. For i {, N} the following inequality holds in the sense of definiteness: v i v i U i < C i I. i= 7

8 3. For all test functions ψ, ϕ, φ Cc [0, R, R R R the following identities hold: [ ϱ t ψ ϱ v x ψ ] dx dt ϱ 0 ψ0, dx = 0, 0 R R [ ϱ v t ϕ ϱ U : x ϕ p ] ϱ C div x ϕ dx dt ϱ 0 v 0 ϕ0, dx = 0, 0 R R [ ] 0 R ϱ C c v p t φ ϱ C c v p v x φ dx dt ϱ 0 v0 c v p 0 φ0, dx = For every non-negative test function ϕ Cc [0, R, R 0 the inequality [ ] ϱ sϱ, p t ϕ ϱ sϱ, p v x ϕ dx dt ϱ 0 sϱ 0, p 0 ϕ0, dx 0 0 R R is fulfilled. R 4. Sufficient condition for non-uniqueness The key tool in the proof of Theorems 3.5 and 3.6 is the following theorem relating existence of infinitely many admissible N-fan solutions to the existence of a single admissible N-fan subsolution Theorem 4.. Let ϱ ±, v ±, p ± be such that there exists an admissible N-fan subsolution ϱ, v, U, C, p to the Cauchy problem.-.4. Then there are infinitely many admissible weak N-fan solutions ϱ, v, p to.-.4 with the following properties: ϱ, p = ϱ, p, v = v ± on Ω ±, v = C i a.e. on Ω i, i =, N. The proof of Theorem 4. is based on a convex integration technique for the pressureless incompressible Euler equations which was introduced in seminal works of De Lellis and Székelyhidi [8, 9] and is summarized in the following Proposition. Proposition 4.3. Let ṽ, Ũ R S0 and C > 0 such that ṽ ṽ Ũ < C I. Furthermore let Ω R R open. Then there exist infinitely many maps v, U L R R, R S0 with the following properties.. v and U vanish outside Ω.. For all test functions ψ, ϕ Cc R R, R R it holds that v x ψ dx dt = 0, Ω v t ϕ U : x ϕ dx dt = 0. Ω 3. ṽ v ṽ v Ũ U = C I is fulfilled almost everywhere in Ω. For the proof of Proposition 4.3 we refer to [4, Lemma 3.7]. The proof of Theorem 4. using Proposition 4.3 is now quite straightforward. 8

9 Proof. Let ϱ, v, U, C, p an admissible N-fan subsolution. Furthermore for i =, N let v i, U i L R R, R S0 pairs of functions as in Proposition 4.3 with ṽ = v i, Ũ = U i, C = C i and Ω = Ω i. It suffices to prove that ϱ, v, p = ϱ, v N v i, p is an admissible weak solution to.-.4. The fact that it is indeed a N-fan solution is clear. Let ψ, ϕ, φ Cc [0, R, R R R be test functions. Then we get [ ϱ t ψ ϱ v x ψ ] dx dt ϱ 0 ψ0, dx 0 R R [ N ] = ϱ t ψ ϱ v v i x ψ dx dt ϱ 0 ψ0, dx 0 R i= R N = ϱ i v i x ψ dx dt = 0, Ω i 0 R = 0 R R i= [ ] ϱ v t ϕ ϱ v v : x ϕ p div x ϕ [ ϱ v t ϕ ϱ U : x ϕ N i= Ω i i= dx dt p ϱ C div x ϕ [ ] ϱ i v i t ϕ ϱ i U i : x ϕ dx dt = 0, [ 0 R ϱ v c v p ϱ 0 v0 c v p 0 φ0, dx R [ = 0 R ϱ C c v p t φ ϱ 0 v0 N c v p 0 φ0, dx ] t φ ϱ v c v p v x φ dx dt ϱ C c v p i= Ω i ϱ 0 v 0 ϕ0, dx R ] dx dt ϱ 0 v 0 ϕ0, dx R ] v x φ dx dt C i ϱ i c v p i v i x φ dx dt = 0. For every non-negative test function ϕ Cc [0, R, R 0 it holds [ ] ϱ sϱ, p t ϕ ϱ sϱ, p v x ϕ dx dt ϱ 0 sϱ 0, p 0 ϕ0, dx 0 R R [ ] = ϱ sϱ, p t ϕ ϱ sϱ, p v x ϕ dx dt ϱ 0 sϱ 0, p 0 ϕ0, dx 0 R R N ϱ i sϱ i, p i v i x ϕ dx dt 0. Ω i i= 4. Algebraic equations and inequalities Since the admissible N-fan subsolution consists of a quintuple of piecewise constant functions which satisfy a certain set of partial differential equations and inequalities, it is easy to observe that these differential constraints are equivalent to a system of algebraic equations and inequalities which arise from the Rankine-Hugoniot conditions on the interfaces between the sets of the N-fan partition. More precisely we have the following. 9

10 Proposition 4.4. Let ϱ ±, p ± R, v ± R be given. The constants µ 0, µ R resp. µ 0, µ, µ R and ϱ i, p i R, αi v i = R γi δ, U β i = i S i δ i γ 0, i and C i R for i =, N define an admissible N-fan subsolution to the Cauchy problem.-.4 if and only if they fulfill the following algebraic equations and inequalities: Order of the speeds: µ 0 < µ resp. µ 0 < µ < µ 4. Rankine Hugoniot conditions on the left interface: µ 0 ϱ ϱ = ϱ v, ϱ β 4. µ 0 ϱ v, ϱ α = ϱ v, v, ϱ δ 4.3 µ 0 ϱ v, ϱ β = ϱ v, C ϱ γ p p 4.4 µ 0 ϱ v C c v p ϱ c v p = 4.5 ϱ v C c v p v, ϱ c v p β Rankine Hugoniot conditions on the right interface: µ N ϱ N ϱ = ϱ N β N ϱ v, 4.6 µ N ϱ N α N ϱ v, = ϱ N δ N ϱ v, v, 4.7 CN µ N ϱ N β N ϱ v, = ϱ N γ N ϱ v, p N p 4.8 C N µ N ϱ N c v p N ϱ v c v p = 4.9 C N ϱ N c v p N β N ϱ v c v p v, If N = we additionally have Rankine Hugoniot conditions on the middle interface: µ ϱ ϱ = ϱ β ϱ β 4.0 µ ϱ α ϱ α = ϱ δ ϱ δ 4. C µ ϱ β ϱ β = ϱ γ C ϱ γ p p 4. C µ ϱ c C v p ϱ c v p = 4.3 C ϱ c C v p β ϱ c v p β Subsolution conditions for i =, N: C i αi βi > Ci Ci α i γ i β i γ i δ i α i β i >

11 Admissibility condition on the left interface: µ 0 ϱ sϱ, p ϱ sϱ, p ϱ sϱ, p β ϱ sϱ, p v, 4.6 Admissibility condition on the right interface: µ N ϱ sϱ, p ϱ N sϱ N, p N ϱ sϱ, p v, ϱ N sϱ N, p N β N 4.7 If N = we additionally have an admissibility condition on the middle interface: µ ϱ sϱ, p ϱ sϱ, p ϱ sϱ, p β ϱ sϱ, p β 4.8 Proof. There is really nothing to prove in this proposition, differential constraints for piecewise constant functions translate to Rankine-Hugoniot conditions in a standard way and the subsolution conditions state that the matrix Ci I v i v i U i is positive definite, which holds if and only if its trace and its determinant are both positive. As we already mentioned, in order to prove our main theorems, due to Theorem 4. it is always enough to find a single N-fan subsolution. Moreover, we are going to make a special ansatz for the subsolutions. Recall that we assume throughout this paper that v, = v, = 0. We will look for N-fan subsolutions satisfying α i = 0 & δ i = α i β i = 0, 4.9 i =, N. One can easily check that this choice implies that equations 4.3, 4.7 and 4. are trivially satisfied. We also note that in the case N = this is in fact a consequence of the set of Rankine-Hugoniot conditions, for more details see [5, Lemma 4.]. For simplicity of notation, we will use in the rest of the paper notation v ± in place of v ±, and similarly v M instead of v M,. 5 Proof of Theorem Reduction to special case We recall that the proof of Theorem 3.5 is finished as soon as we find a single admissible -fan subsolution. Since the structure of the -fan subsolution is similar to the D Riemann solution to the problem, we will look for the -fan subsolution as a suitable small perturbation of the D Riemann solution. Moreover, we claim that we can make the following choice without loss of generality. Proposition 5.. Let the assumption of Theorem 3.5 hold and assume furthermore that v M = 0, where v M is the second component of the velocity of the D Riemann solution see Proposition.. Then there exist infinitely many admissible weak -fan solutions to the Riemann problem.-.4. Proposition 5. implies Theorem 3.5 using the Galilean invariance of the Euler equations. Indeed, if we consider new initial data as follows: ϱ 0,new, v 0,new, p 0,new := ϱ 0, v 0 0, v M T, p the D Riemann solution to the problem.-.3 with 5. will fulfill v new M = 0. We apply Proposition 5. and find infinitely many admissible weak -fan solutions ϱ, v, p to the Riemann problem with initial data 5.. Then, again by Galilean invariance, ϱ, v 0, v M T, p are admissible weak -fan solutions to the Riemann problem with original initial data.4, and the proof of Theorem 3.5 is finished. What remains now is to prove Proposition 5..

12 5. Proof of Proposition 5. Let the assumptions of Proposition 5. be true. As we mentioned above, we will look for an admissible -fan subsolution as a suitable perturbation of the D Riemann solution. This perturbation will be quantified by a small parameter ε > 0. For convenience we define functions A, B, D : R R by Aε := ϱ ϱ M εϱ M ε ϱ ϱ ϱ M εϱ M ε ϱ ; Bε := ϱ ϱ ϱ M εϱ M ε v v p p Aε; Dε := v ϱ ϱ M εϱ M ε ϱ v ϱ ϱ M εϱ M ε ϱ. First of all we want to show some properties of the functions A and B. It is easy to deduce from Proposition. and the assumption v M = 0, that v = p M p c v ϱ p c v p M, 5. v = p M p c v ϱ p c v p M. 5.3 Since p M > max{p, p }, we have v < 0 < v. Furthermore we obtain from Proposition. and 5., 5.3 that ϱ M ϱ c v p M p = ϱ ϱ M ϱ p c v p M = v, p M p ϱ M ϱ c v p M p = ϱ ϱ M ϱ p c v p M = v. p M p This leads to v p p ϱ M ϱ ϱ ϱ M = v p M p p M p, and finally v p p ϱ M ϱ ϱ ϱ M = v p M p p M p, B0 = = ϱ ϱ ϱ M ϱ M v v p p ϱ ϱ M ϱm ϱ ϱ ϱ M ϱm ϱ [ v = ϱ ϱ ϱ M ϱ M p p ϱ M ϱ v p p ϱ M ϱ v v ϱ ϱ M ϱ ϱ M [ v p M p = ϱ ϱ ϱ M ϱ M ] p M p v v v > 0. p M p p M p Now, by continuity of the function B, there exists an ε max, > 0 such that Bε > 0 for all ε 0, ε max, ]. Because ϱ M > ϱ, there exists ε max, > 0 such that ϱ M ε ϱ > 0 for all ε 0, ε max, ]. Next we want to show that there is an ε max,3 > 0 such that Aε 0 for all ε 0, ε max,3 ]. To this end, let us first assume, that A0 0. Then, by continuity of the function A, there exists such an ε max,3 > 0. Now consider the case where A0 = 0. Then we obtain Aε = ϱ ϱ M εϱ M ε ϱ ϱ ϱ M εϱ M ε ϱ = ε ϱ ϱ ε ϱ ϱ ϱ ϱ ϱ M ϱ M A0 }{{} =0 = ε ϱ ϱ ε ϱ ϱ ϱ ϱ ϱ M ϱ M, ]

13 which has at most two zeros: If ϱ = ϱ then and if ϱ ϱ then Aε = 0 ε = 0; Aε = 0 ε = 0 or ε = ϱ ϱ ϱ ϱ ϱ M ϱ M ϱ ϱ. Hence there exists ε max,3 > 0 such that Aε 0 for all ε 0, ε max,3 ]. We set ε max := min{ε max,, ε max,, ε max,3 } and then we have Aε 0, ϱ M ε ϱ > 0, Bε > 0, ϱ M ε ϱ > 0, 5.4 for all ε 0, ε max ]. 5.. Shock speeds Next we define the functions µ 0, µ, µ : 0, ε max ] R by µ 0 ε := [ ϱm Dε ϱ ϱ ϱ M εv v ε ] ϱ M ε ϱ Bε ; Aε ϱ M ε ϱ µ ε := [ Dε ] ϱ M ε ϱ ϱ M ε ϱ Bε ; Aε µ ε := [ Dε ϱ ϱ ϱ M εv v ϱ M ε Aε ϱ ] M ε ϱ Bε. ϱ M ε ϱ Note first that the functions µ 0, µ, µ are well-defined because of 5.4. We claim that these functions define perturbations of the shock speeds σ, v M and σ of the D Riemann solution. More precisely we have Proposition 5.. It holds that lim µ 0ε = σ, ε 0 lim µ ε = v M, ε 0 lim µ ε = σ. ε 0 Proof. We start with the Rankine-Hugoniot conditions for the D Riemann solution σ ϱ ϱ M = ϱ v ϱ M v M ; 5.5 σ ϱ v ϱ M v M = ϱ v ϱ M v M p p M ; 5.6 σ ϱ M ϱ = ϱ M v M ϱ v ; 5.7 σ ϱ M v M ϱ v = ϱ M v M ϱ v p M p 5.8 and we obtain by eliminating σ, σ and p M that ϱ M ϱ ϱ v ϱ M v M ϱ ϱ M ϱ M v M ϱ v = ϱ v ϱ v ϱm v M ϱ M p p ϱ M ϱ ϱ ϱ M. In order to solve this equation for v M, we write it as follows where the constant A0 v M D0 vm E = 0, 5.9 E := p p ϱ M ϱ ϱ M ϱ v ϱ ϱ M ϱ M ϱ v ϱ ϱ M ϱ M ϱ 3

14 only depends on the initial states. Now we have to consider two cases, namely A0 = 0 and A0 0. Let us start with A0 = 0. Then we easily deduce that D0 = ϱ ϱ M ϱ M ϱ v v, which does not vanish because v < 0 < v and ϱ M > ϱ. Hence we get from 5.9, that v M = E D Next we want to compute lim ε 0 µ ε and compare it with 5.0. Keeping in mind that we are considering the case A0 = 0, we get D0 ϱ M ϱ ϱ M ϱ B0 = v v ϱ ϱ M ϱ M ϱ ϱ M ϱ ϱ M ϱ ϱ ϱ ϱ M ϱ M v v = v v ϱ ϱ M ϱ M ϱ ϱ M ϱ ϱ ϱ M v v = 0. Hence we can apply L Hospital s rule. We obtain lim ε 0 A ε = lim ϱ ϱ M ε ϱ ϱ ϱ M ε ϱ ϱ M ε ϱ ϱ M ε ϱ ε 0 = ϱ ϱ ϱ ϱ ϱ M ϱ M. A short calculation shows that this is non-zero: From A0 = 0 we can deduce that ϱ ϱ ϱ ϱ = ϱ M ϱ M ϱ M ϱ M. This means that ϱ ϱ and ϱ M ϱ M have the same sign, which implies ϱ ϱ ϱ M ϱ M 0. Since ϱ ϱ > 0, we have A 0 < 0, in particular A 0 0 Furthermore a long but straightforward computation yields [ lim Dε ] A 0 E ϱ M ε ϱ ϱ M ε ϱ Bε = ε 0 D0. Hence by L Hospital s rule we obtain lim µ ε = E ε 0 D0 and recalling 5.0, we deduce lim µ ε = ε 0 v M. Let us now consider the case A0 0. Then we obtain from 5.9 that v M = [ D0 ± ] D0 A0 A0 E. The correct sign in the equation above is because v M = 0 and D0 > 0, which easily follows from v < 0 < v. Furthermore it is simple to check that D0 A0 E = ϱ M ϱ ϱ M ϱ B0. Then it is easy to conclude µ 0 = v M. To finish the proof of Proposition 5. we have to show that lim µ ε = v M ε 0 lim µ 0ε = σ and lim µ ε = σ. It is straightforward to deduce that ε 0 ε 0 µ 0 ε = v ϱ M ε ϱ M ε ϱ µ ε v ; µ ε = v ϱ M ε ϱ M ε ϱ µ ε v. Alternatively, this can be verified by considering the admissibility criterion. implies that 4

15 On the other hand we get from 5.5 and 5.7, that σ = v ϱ M ϱ M ϱ vm v ; σ = v ϱ M ϱ M ϱ vm v. Hence we easily deduce lim ε 0 µ 0 ε = σ and lim ε 0 µ ε = σ. Because of σ < v M < σ and the continuity of the functions µ 0, µ, µ, we may assume that µ 0 ε < µ ε < µ ε for all ε 0, ε max ]. If this is not the case we redefine ε max to be a bit smaller that the smallest positive value of ε for which µ 0 ε < µ ε < µ ε is violated. 5.. Constants C i and γ i In order to proceed further we need to introduce a second positive parameter δ > 0. We define the functions C, C, γ, γ : 0, ε max ] 0, p M by [ C ε, δ := ϱ M ε µ 0 ε µ ε µ 0 ε c v p M δ p ϱ v ] µ εc v p M δ ϱ v c v p v ; [ C ε, δ := ϱ M ε µ ε µ ε µ ε c v p M δ p ϱ v ] µ εc v p M δ ϱ v c v p v ; [ γ ε, δ := ϱ M ε C ε, δ ϱ v p M δ p µ 0 ε ] ϱ M εµ ε ϱ v ; ϱ M ε [ γ ε, δ := ϱ M ε C ε, δ ϱ v p M δ p µ ε ] ϱ M εµ ε ϱ v. ϱ M ε Note that these functions are well-defined because of the arguments above. More precisely, it holds that µ 0 ε µ ε 0 and µ ε µ ε 0 for all ε 0, ε max ]. Proposition 5.3. It holds that lim C ε, δ = ; v M ε,δ 0 lim C ε, δ = ; v M ε,δ 0 lim γ vm ε, δ = ; ε,δ 0 lim γ vm ε, δ =. ε,δ 0 Proof. To prove this, we need the Rankine Hugoniot conditions of the D Riemann solution in the energy equation σ ϱ v c v p ϱ M v M c v p M = ϱ v c v p v ϱ M v M c v p M v M ; σ ϱ M v M c v p M ϱ v c v p = ϱ M v M c v p M v M ϱ v c v p v

16 We obtain that lim C ε, δ = v M and lim C ε, δ = v M by using Proposition 5. ε,δ 0 ε,δ 0 and 5., 5.. The fact that lim γ ε, δ = v M ε,δ 0 and lim γ ε, δ = v M ε,δ 0 can be shown analogously by using the Rankine Hugoniot conditions 5.6, 5.8. We continue the proof of the Proposition 5. by observing that the perturbations defined above indeed help to define an admissible -fan subsolution. Proposition 5.4. If there exists ε, δ 0, ε max ] 0, p M such that the following inequalities are fulfilled, then there exists an admissible -fan subsolution to the Cauchy problem.-.4. Order of the speeds: µ 0 ε < µ ε < µ ε 5.3 Subsolution conditions: C ε, δ µ ε > C ε, δ µ ε > C ε, δ C ε, δ γ ε, δ µ ε γ ε, δ > C ε, δ C ε, δ γ ε, δ µ ε γ ε, δ > Admissibility condition on the left interface: µ 0 ε ϱ M ε sϱ M ε, p M δ ϱ sϱ, p 5.8 ϱ M ε sϱ M ε, p M δ µ ε ϱ sϱ, p v Admissibility condition on the right interface: µ ε ϱ sϱ, p ϱ M ε sϱ M ε, p M δ ϱ sϱ, p v ϱ M ε sϱ M ε, p M δ µ ε 5.9 Proof. Let there be ε, δ 0, ε max ] 0, p M such that hold. In order to show that there exists an admissible -fan subsolution, we use Proposition 4.4, i.e. we define the constants appearing in Proposition 4.4 as follows: µ 0 := µ 0 ε; µ := µ ε; µ := µ ε; ϱ := ϱ M ε; ϱ := ϱ M ε; 0 v := v := ; p µ ε := p := p M δ; γ ε, δ 0 γ ε, δ 0 U := ; U 0 γ ε, δ := 0 γ ε, δ C := C ε, δ; C := C ε, δ. It is a matter of straightforward calculation to check that with this choice hold. ; 6

17 5..3 Subsolution and admissibility inequalities In order to finish the proof of Proposition 5., we have to find ε, δ 0, ε max ] 0, p M such that the conditions are satisfied. We start with noting that we already have 5.3 fulfilled for all ε 0, ε max ]. Let us now investigate the subsolution conditions We start with the terms in the first parenthesis in We obtain by using that δ 0, p M C ε, δ C ε, δ γ ε, δ = c v δ ϱ M ε µ ε δ ϱ M ε µ 0 ε µ ε C ε, 0 γ ε, 0 c v δ ϱ M ε µ ε p M ϱ M ε µ 0 ε µ ε C ε, 0 γ ε, 0 }{{} =:R ε γ ε, δ = c v δ ϱ M ε µ ε δ ϱ M ε µ ε µ ε C ε, 0 γ ε, 0 c v δ ϱ M ε µ ε p M ϱ M ε µ ε µ ε C ε, 0 γ ε, 0 }{{} =:R ε where Propositions 5. and 5.3 together with the fact that v M = 0 imply that lim R ε = 0; ε 0 lim R ε = 0. ε 0 Therefore R ε and R ε become arbitrary small if we choose ε small. Because of c v >, there exists ε δ 0, ε max ] for each δ 0, p M, such that C ε, δ γ ε, δ > 0 & C ε, δ γ ε, δ > 0 5. hold for all ε 0, ε δ. Similarly we handle terms in the second parenthesis in inequalities We obtain C ε, δ C ε, δ µ ε δ γ ε, δ = ϱ M ε C ε, 0 µ ε γ ε, 0 } {{} =:R 3ε 5.3 µ ε δ γ ε, δ = ϱ M ε C ε, 0 µ ε γ ε, } {{} =:R 4ε With the same arguments as above, we obtain that for each δ 0, p M there exists ε δ 0, ε max ] such that C i ε, δ µ ε γ i ε, δ > 0, i =, 5.5 hold for all ε 0, ε δ. Combining 5. and 5.5 we obtain 5.6 and 5.7 while summing together 5. and 5.5 we obtain 5.4 and 5.5. To finish the proof we have to show that we can achieve that in addition the admissibility conditions 5.8 and 5.9 hold. Note that in the limit ε, δ 0, 0 the admissibility conditions 5.8 and 5.9 turn into the admissibility conditions of the D Riemann solution according to Proposition 5.. Since the latter are fulfilled strictly, we can choose δ > 0 sufficiently small and also ε 0, min{ ε δ, ε δ} sufficiently small such that 5.8 and 5.9 hold. This finishes the proof of Proposition 5.. 7

18 6 Proof of Theorem 3.6 Unlike in the case of -fan solutions as in Theorem 3.5, we cannot use here the D Riemann solution and try to perturb it in a suitable way in order to find a -fan subsolution. On the other hand, since we don t have the middle interface, the number of equations and inequalities in Proposition 4.4 is lower than in the case N =. Since i could be only equal to in Proposition 4.4, for simplicity of notation we write α, β, γ, δ and C instead of α, β, γ, δ and C. We also recall that we set α = 0 and δ = αβ = 0, see 4.9. The subsolution conditions then simplify to C β > 0; 6. C C γ β γ > It is not difficult to observe see also [5, Lemma 4.3] that the necessary condition for to be satisfied is C γ > β which motivates us to introduce instead of C and γ new unknowns ε := C γ β ; 6.3 ε := C β ε. 6.4 The set of algebraic identities and inequalities from Proposition 4.4 then simplifies into Order of the speeds µ 0 < µ 6.5 Rankine Hugoniot conditions on the left interface: µ 0 ϱ ϱ = ϱ v ϱ β 6.6 µ 0 ϱ v ϱ β = ϱ v ϱ β ε p p 6.7 µ 0 cv p ϱ v c v p ϱ β ε ε = c v p ϱ v v c v p ϱ β ε ε β ; 6.8 Rankine-Hugoniot conditions on the right interface: µ ϱ ϱ = ϱ β ϱ v 6.9 µ ϱ β ϱ v = ϱ β ε ϱ v p p 6.0 µ cv p ϱ β ε ε c v p ϱ v = c v p ϱ β ε ε β c v p ϱ v v ; 6. Subsolution conditions: Admissibility conditions ε > 0 6. ε > 0 ; 6.3 µ 0 L L L v L β 6.4 µ L L L β L v

19 where for simplicity of notation we introduced L i as p c v L i = ϱ i sϱ i, p i = ϱ i log i ϱ cv, i =,,. i We observe that we have 6 equations and 5 inequalities for 7 unknowns µ 0, µ, ϱ, p, β, ε, ε. Moreover to these 5 inequalities we have to keep in mind other inequalities which have to be satisfied, namely ϱ > 0 and p > 0. We define the following quantities, functions of the Riemann initial data R = ϱ ϱ 6.6 A = ϱ v ϱ v 6.7 H = ϱ v ϱ v p p 6.8 u = v v 6.9 B = A RH = ϱ ϱ u ϱ ϱ p p. 6.0 Since we have six equations for seven unknowns it is reasonable to choose one unknown as a parameter and express other unknowns as functions of Riemann initial data and this parameter. We choose the density in the middle region ϱ as a parameter. Note that we may assume B > 0, this can be achieved by taking u large enough as in the assumptions of Theorem The case R < Solution to algebraic equations Recalling that R < 0 means ϱ < ϱ we express µ 0 and µ from equations 6.6, 6.7, 6.9, 6.0 as µ 0 = A R B ϱ ϱ 6. R ϱ ϱ µ = A R B ϱ ϱ 6. R ϱ ϱ and we observe that µ 0 < µ if ϱ > ϱ > ϱ. In what follows we therefore always assume that the set of possible values of ϱ is ϱ ϱ,. There exists also another solution µ 0, µ, which we do not use, see Remark 6.6 at the end of section Next we express β in two ways which will be useful later. Using 6.9 and 6.6 respectively we obtain β = ϱ ϱ v µ 0 ϱ ϱ ϱ ; 6.3 β = ϱ ϱ v µ ϱ ϱ ϱ. 6.4 It is technically more difficult to express p, ε and ε. We start with rewriting 6.0 as p = p ϱ v ϱ β ε µ ϱ ϱ. 6.5 Then, a lengthy yet straightforward computation yields the following equations as a consequence of 6.8, 6. β v ϱ ε ϱ ϱ ϱ ϱ β v ε ε = β v [ϱ ϱ p p c v ϱ p c v ϱ p ] ; 6.6 v βϱ ε ϱ ϱ ϱ ϱ v βε ε = v β [ϱ ϱ p p c v ϱ p c v ϱ p ]

20 At this point we assume that v β and v β are nonzero and refer to Lemma 6. for proof. We continue by expressing ε ε from 6.7 and plugging this to 6.6. We have ε ε = [ϱ ϱ p p c v ϱ p c v ϱ p ] v β ϱ ϱ v β ε ϱ ϱ 6.8 ϱ and consequently 6.6 becomes β v ϱ ϱ ϱ ϱ v β β v v β ϱ ϱ ϱ ϱ ε =ϱ [ϱ ϱ p p c v ϱ p ϱ p ] ϱ [ϱ ϱ p p c v ϱ p ϱ p ]. 6.9 Expressing the right hand side of 6.9 we end up with β v ϱ ϱ ϱ ϱ v β β v v β ϱ ϱ ϱ ϱ ε =ϱ p R ϱ ϱ p p c v ϱ ϱ p ϱ p Once again it is useful to introduce further notation to simplify resulting expressions. We define X = ϱ ϱ p p c v ϱ ϱ p ϱ p 6.3 Y = β v β v ϱ ϱ ϱ v β v β ϱ ϱ ϱ 6.3 Z = p ϱ v ϱ β µ ϱ ϱ This way, assuming Y 0 which will be justified in Lemma 6., we rewrite 6.30 as ε = ϱ p R X ϱ Y 6.34 and plugging this into 6.5 and assuming Y ϱ R 0 which will be justified in Lemma 6., too, we get p = Y Z X Y ϱ R, 6.35 this in turn yields ε = ϱ RZ X 6.36 ϱ Y ϱ R and finally ε is expressed using 6.8 ε = ϱ ϱ [ϱ ϱ p p c v ϱ p ϱ p ] v β ϱ ϱ ε v β ϱ 6.. Positivity of p, ε and ε Before we continue further, we remind the following useful expressions which can be derived from 6.6 and 6.9 β µ 0 = ϱ ϱ v µ v β = ϱ ϱ ϱ v µ µ β = ϱ ϱ µ v 6.40 β v = ϱ ϱ ϱ µ v. 6.4 In particular we observe that the signs of v β and β v are the same as v µ 0 and µ v respectively. 0

21 Lemma 6.. For u = v v sufficiently large it holds Proof. We have v µ 0 = ϱ u R = ϱ u R v µ 0 > 0 & µ v > ϱ ϱ u R ϱ ϱ p p ϱ ϱ ϱ ϱ ϱ ϱ u ϱ ϱ p p ϱ ϱ ϱ ϱ R Now we distinguish two cases. First, if p p, then we easily observe recall ϱ > ϱ that However if p < p we get v µ 0 > R ϱ u ϱ ϱ u > v µ 0 > R ϱ u ϱ ϱ u ϱ ϱ p p 6.45 and the expression on the right hand side can be made positive assuming u is large enough. Similarly we proceed with quantity µ v. We recall that Lemma 6. justifies the assumptions v β and v β nonzero we made in the process of deriving 6.35 and We continue by analyzing whether one can find ϱ such that p, ε and ε are all positive. This part starts with rewriting the expression 6.3 for Y. We have Y = ϱ ϱ ϱ v β v β ϱ ϱ ϱ β v β v v β β v = ϱ ϱ ϱ ϱ v µ 0 µ v = ϱ ϱ v µ 0 v µ 0 ϱ ϱ µ v µ v, 6.46 where the last equality holds although obviously β / {µ 0, µ }. In what follows our strategy is to express certain quantities as functions of ϱ and u while treating ϱ ±, p ± and v as data. This way we obtain Y = ϱ R B ϱ u ϱ ϱ ϱ ϱ v ϱ ϱ ϱ u Bϱ ϱ ϱ ϱ ϱ u Bϱ ϱ Lemma 6.. For u = v v large enough it holds Y < 0 and Y ϱ R < 0 for all ϱ > ϱ. Proof. To prove Lemma 6. we examine the limit of Y as u. First we observe B ϱ ϱ ϱ ϱ ϱ ϱ u ϱ ϱ ϱ ϱ ϱ > 6.48 and therefore the second term in the square brackets of 6.47 is positive. We continue by proving

22 that the last term converges to zero. We write v Bϱ ϱ ϱ ϱ ϱ u Bϱ ϱ ϱ ϱ ϱ u = v u Bϱ ϱ u ϱ ϱ ϱ Bϱ ϱ ϱ u ϱ ϱ Bϱ ϱ ϱ u ϱ Bϱ ϱ ϱ u ϱ ϱ ϱ = v u Bϱ ϱ u ϱ ϱ ϱ Bϱ ϱ ϱ u ϱ ϱ 6.49 and observe that both terms in the numerator of the fraction converge to distinct quantities and the denominator has finite limit as u. Therefore the whole fraction has a finite nonzero limit and the whole term therefore can be made arbitrarily small by choosing u sufficiently large. Altogether we have proved that the expression in the square brackets in 6.47 is positive at least for u sufficiently large and since R < 0 we conclude that Y < 0 for large u. Also we observe that Y has a finite limit as u. The same two properties obviously hold also for Y ϱ R. Remark 6.3. At this point we emphasize the need to take v as fixed and by taking u large enough recover v. This procedure would not work the other way round, i.e. taking v fixed. If we would consider v fixed, the expression 6.47 would become Y = ϱ R B ϱ ϱ ϱ u ϱ ϱ v It is not difficult to figure out that lim u Bϱ ϱ ϱ ϱ ϱ u Bϱ ϱ ϱ ϱ ϱ u < 0, 6.5 B ϱ ϱ ϱ u ϱ ϱ which yields Y positive and this would cause troubles later when studying the sign of ε. X Now let us turn our attention to expressions X and Z. We recall here the definition 6.3 of X = ϱ ϱ p p c v ϱ ϱ p ϱ p. 6.5 We see that X does not depend on u and may be positive or negative depending on values of ϱ ±, p ± and ϱ. Concerning Z we have from 6.33 Z = p ϱ v ϱ β µ ϱ ϱ = p ϱ ϱ ϱ ϱ v µ 6.53 and we see that Z is always positive. Moreover we plug in the expression 6. to obtain Z = p ϱ ϱ ϱ ϱ v µ = p ϱ ϱ ϱ ϱ R = p ϱ ϱ ϱ u ϱ R B ϱ ϱ ϱ u ϱ ϱ ϱ ϱ Rp p u ϱ ϱ ϱ ϱ ϱ. 6.54

23 Observing that the quantity in parenthesis on the last line of 6.54 has finite nonzero limit as u we conclude that Z grows as u. We are now ready to study the signs of p and ε. Lemma 6.4. For u = v v large enough it holds p > 0 and ε > 0 for all ϱ > ϱ. Proof. Since we have from 6.35 p = Y Z X Y ϱ R, 6.55 Y Z X Y ϱ R has we easily see that for large u the leading term is Y ϱ which grows like R u, whereas a finite limit as u. Since both Y and Y ϱ R are negative at least for large u and Z > 0, we conclude that p > 0 for u large enough. Similarly we have from 6.36 ε = ϱ RZ X ϱ Y ϱ R and using similar arguments as above for p we see that the term positive whereas RZ Y ϱ R 6.56 grows like u and is X ϱ Y ϱ R has finite limit as u and therefore ε > 0 for u large enough. Next we study the sign of ε. We have Lemma 6.5. For ρ > c v ϱ ϱ and u = v v large enough it holds ε > 0. ϱ 4c v ϱ ϱ 6.57 Proof. We start from 6.37 ε = ϱ ϱ [ϱ ϱ p p c v ϱ p ϱ p ] ε = p c v ϱ ϱ p ϱ c v ϱ ϱ ϱ ϱ ϱ ε ϱ ϱ R ϱ ϱ ϱ Bϱ ϱ u ϱ ϱ v β ϱ ϱ v β ϱ ϱ Rv ϱ B ϱ ϱ ϱ ϱ ϱ u Keeping in mind that both p and ε grow as u we identify the leading terms of ε as u. We have ε = Y Z ϱ c v ϱ RZ ϱ Y ϱ R ϱ ϱ Y ϱ R ϱ ϱ ϱ ϱ R Bϱ ϱ u ϱ ϱ Z = Y ϱ c v ϱ ϱ R R ϱ ϱ Y ϱ R ϱ ϱ ϱ =: Bϱ ϱ u ϱ ϱ l.o.t. l.o.t. Z E l.o.t., 6.59 ϱ Y ϱ R where l.o.t. stands for lower order terms with respect to u Z. We know that ϱ Y ϱ R is negative so we continue by studying the remaining part E of the leading term of ε and to make ε > 0 we want to show that E < 0. 3

24 Recall that from we already know that lim Y = ϱ R 6.60 u ϱϱ ϱ ϱ ϱ ϱ and therefore a straightforward calculation yields lim E = R ϱ ϱ c v ϱ c v ϱ u ϱ 6.6 ϱϱ ϱ ϱ ϱ ϱ Since we have R < 0, we want the term in square parenthesis of 6.6 to be positive. It is not difficult to show that this can be done by choosing ϱ large enough with respect to ϱ ±, more specifically if then lim u E < 0. ϱ > c v ϱ ϱ ϱ 4c v ϱ ϱ Admissibility inequalities Finally we have to check whether the admissibility inequalities are satisfied. start by plugging the expression 6.3 into 6.4 and we obtain We v µ 0 ϱ L ϱ L We already know that v µ 0 > 0 at least for u large enough, see Lemma 6.. Therefore we need to ensure that ϱ L ϱ L Using the definition of L i the inequality 6.64 is equivalent to p cv ϱ cv pcv ϱ cv We proceed similarly with 6.5 where we plug in the expression 6.4 for β to get µ v ϱ L ϱ L Again by Lemma 6. we know that µ v > 0 at least for u large enough and thus we need which is equivalent to p cv ϱ cv Combining 6.65 and 6.68 we end up with p cv ϱcv max ϱ L ϱ L, 6.67 pcv ϱ cv { p c v ϱ cv. 6.68, } p cv ϱ cv However we already know that p grows as u, so regardless of the choice of ϱ we can always ensure that 6.69 is satisfied by choosing u large enough. We have showed that for any ϱ satisfying 6.57 there exists some uϱ such that if u > uϱ then we can construct a -fan subsolution. Then the point i of Theorem 3.6 in the case R < 0 is proved by defining U = Uϱ ±, p ±, v := v inf uρ, 6.70 ρ where the inf is taken among ρ such that 6.57 holds. 4

25 Remark 6.6. There exists also another solution to equations 6.6, 6.7, 6.9, 6.0, namely µ 0 = A R B ϱ ϱ 6.7 R ϱ ϱ µ = A R B ϱ ϱ 6.7 R ϱ ϱ and here it holds µ 0 < µ for ϱ < ϱ < ϱ. However this solution is not convenient to work with, since it violates Lemma 6.. More precisely, one of the inequalities 6.4 holds with the opposite sign. This causes troubles in the analysis of the admissibility inequalities , because these then yield one upper bound and one lower bound for p instead of two lower bounds. Then naturally the argument with taking u sufficiently large fails. 6. The case R > 0 The proof in the case R > 0 follows the same steps. The expressions for µ 0, µ stay the same and again we see that µ 0 < µ if ϱ > ϱ > ϱ. Since we treat v as fixed in this case, instead of 6.47 we get Y = ϱ R B ϱ ϱ ϱ u ϱ ϱ v Bϱ ϱ ϱ ϱ ϱ u Bϱ ϱ ϱ ϱ ϱ u Similarly as in the case R < 0 we show that the term in the square brackets is positive at least for u sufficiently large, so in this case we conclude that Y > 0 and obviously also Y ϱ R > 0. Since nothing changes in terms X and Z, this yields that p and ε are both positive at least for u sufficiently large. The analysis of ε changes again with respect to the fact that v is now fixed instead of v. Instead of 6.58 we obtain ε = p c v ϱ ϱ ϱ ϱ p ϱ c v ϱ and hence ε = Z ϱ Y ϱ R Y ϱ c v ϱ ϱ ϱ ϱ ε ϱ ϱ ϱ Rv ϱ B ϱ ϱ ϱ ϱ ϱ u 6.74 Z ϱ R l.o.t. =: l.o.t ϱ Y ϱ RẼ Since the limit of Y as u changes to lim Y = ϱ R, 6.76 u ϱ ϱ ϱ ϱ ϱ ϱ a direct calculation yields lim u Ẽ = R It is not difficult to obtain that lim u Ẽ > 0 if c v ϱ ϱ c v ϱ ϱ ϱ ϱ ϱ ϱ ϱ ϱ > ϱ c v ϱ ϱ 4c v ϱ ϱ and thus choosing such ϱ we always find u large enough such that ε > 0. Since nothing changes regarding the admissibility inequalities, we conclude the proof of point ii of Theorem 3.6 in the case R > 0 by the same arguments as in the case of point i in the case R < 0. 5

26 6.3 The case R = Solution to algebraic equations This case has to be treated separately since even the expression for the speeds of interfaces µ 0, µ are considerably different. We have µ 0 = ϱ u p p v v 6.79 ϱ ϱ ϱ u µ = ϱ u p p v v 6.80 ϱ ϱ ϱ u and we immediately see that since we examine the case where u > 0 we are forced to assume ϱ > ϱ = ϱ in order to ensure µ 0 < µ. The relations stay the same in the case R = 0, however the expressions for quantities X and Y simplify a little bit to X = ϱ p p c v ϱ ϱ p p = ϱ p p c v ϱ ϱ 6.8 v β Y = ϱ ϱ ϱ v β β v 6.8 β v and consequently we have and ε = X ϱ Y 6.83 p = Z X Y, 6.84 with Z defined as in Finally, ε is given by Positivity of p, ε and ε We claim that Lemma 6. holds also in the case R = 0. Indeed, it is even easier here to see that v µ 0 = ϱ u p p 6.85 ϱ ϱ ϱ u µ v = ϱ u p p 6.86 ϱ ϱ ϱ u and both expressions on the right hand sides are positive at least for u large enough. Next, we observe that the sign of X depends on the sign of p p, namely sign X = sign p p 6.87 and in particular this shows that this presented construction simply does not work in the case ϱ = ϱ and p = p, because in this case we end up with ε = 0. Also, knowing the sign of X and having in mind the expression 6.83 for ε we see that in order to ensure ε > 0 we need to have sign Y = sign X = sign p p Therefore we study now the sign of Y in a similar manner as in the proof of Lemma 6.. We start this part by claiming that the expression 6.46 still holds and can be written as v µ 0 Y = ϱ ϱ µ v v µ 0 µ v 6

27 Using 6.79, 6.80, keeping v fixed and expressing v = v u we end up after a simple calculation with Y = 4ϱ ϱ v ϱ u ϱ ϱ ϱ ϱp p ϱ ϱ u3 ϱ ϱ ϱ p p ϱ u We observe that in the limit u we get Y 4ϱ ϱ ϱ < 0. Therefore we conclude that ε > 0 for u large enough if we keep v fixed and if p > p. Similarly if we keep v fixed, express v = v u we end up instead of 6.90 with Y = 4ϱ ϱ ϱ u ϱ ϱ v ϱ ϱp p ϱ ϱ u3 ϱ ϱ ϱ p p ϱ u 6.9 and in the limit u we get Y 4ϱ ϱ ϱ > 0. Therefore we conclude that ε > 0 for u large enough if we keep v fixed and if p < p. Concerning the sign of p we have similarly as in 6.54 Z = ϱ ϱ u 4ϱ ϱ p p ϱ ϱ p p 4ϱ ϱ u, 6.9 thus we see that Z > 0 and Z u as u. Therefore using 6.84 we conclude p > 0. Remark 6.7. From the considerations above it is also clear why we do not state anything about the case ϱ = ϱ and p = p in Theorem 3.6. In this case X = 0 and the only way how one could choose ε > 0 and satisfy 6.30 would be to have also Y = 0. However, even if p = p we observe that Y has nonzero finite limit as u and hence this construction fails in this case. We finish this section with the study of the sign of ε. Lemma 6.8. For and u = v v large enough it holds ε > 0. ρ > c v ϱ 6.93 Proof. Recalling 6.37 we have ε = v β ϱ ϱ [ϱ ϱ p p c v ϱ p ϱ p ] ε ϱ ϱ v β ϱ ϱ c v ϱ c v ϱ ϱ v β ϱ ϱ = p p ε ϱ ϱ ϱ ϱ v β ϱ We already know, that p grows as u as u. The second term on the right hand side of 6.94 does not depend on u at all. Concerning the third term, we know that ε has a finite limit as u and now we are going to show that the same property holds also for the term vβ v β ϱ ϱ ϱ Plugging in 6.80 we obtain. Indeed, by a simple calculation using 6.4 we have v β ϱ ϱ = ϱ v µ v β ϱ ϱ v µ ϱ ϱ v µ v µ = ϱ ϱ = ϱ ϱ = ϱ ϱ ϱ u ϱ ϱ ϱ u ϱ ϱ ϱ u ϱ ϱ 3v v p p ϱ u 4v u p p ϱ u 4v 3u p p ϱ u