arxiv: v3 [math.ap] 26 May 2015

Transcription

1 Delta shock wave for a 3 3 hyperbolic system of conservation laws arxiv: v3 [math.ap] 26 May 215 Richard De la cruz Juan Galvis Juan Carlos Juajibioy Leonardo Rendón August 27, 218 XV International Conference on Hyperbolic Problems: Theory, Numerics, Applications - HYP 214 Abstract We study the one-dimensional Riemann problem for a hyperbolic system of three conservation laws of Temple class. The system is a simplification of a recently propose system of five conservations laws by Bouchut and Boyaval that models viscoelastic fluids. An important issue is that the considered 3 3 system is such that every characteristic field is linearly degenerate. We study the Riemann problem for this system and under suitable generalized Rankine-Hugoniot relation and entropy condition, both existence and uniqueness of particular delta-shock type solutions are established. 1 Introduction In this work, we show existence and uniqueness of delta shock wave for the Suliciu relaxation system [1, 19] ρ t +ρu) x =, ρu) t +ρu 2 +s 2 v) x =, 1) ρv) t +ρuv +u) x =, with s = const. >. The concept of delta shock wave is a generalization of an ordinry shock wave. This generalization was introdduced by Korchinski in the year of 1977 in his PhD thesis [14]. He considered the Riemann problem for system u t u2) x =, v t uv) x =. 2) Universidad Pedagógica y Tecnológica de Colombia, Tunja, Colombia. richard.delacruz@uptc.edu.co Universidad Nacional de Colombia, Bogotá, Colombia. jcgalvisa@unal.edu.co Universidad Nacional de Colombia, Bogotá, Colombia. jcjuajibioyo@unal.edu.co Universidad Nacional de Colombia, Bogotá, Colombia. lrendona@unal.edu.co 1

2 Motivated by some numerical results, he constructed the unique Riemann solution using generalized delta functions to obtain singular shocks satisfying 2) in the sense of distributions. After that, in 1994, Tan, Zhang and Zheng established in [2] the existence, uniqueness and stability of delta shock waves for a viscous perturbation of system 2) where the first equation is replaced by u t + u 2) x = εu xx, with ε >. Other works dealing with delta-shock waves are due to Ercole [11] who in 2 obtained a delta shock solution as a limit of smooth solutions by the vanishing viscosity method. Sheng and Zhang [18] discussed the Riemann problem for pressureless gas system. In 25, Brenier [2] considered the Riemann problem for the Chaplygin gas system. Other works for the Chaplygin gas system can be found in [6, 12, 13]. Some more recent works on delta shocks for general hyperbolic conservation laws are due to Danilov and Mitrovic [7, 8], where they described delta shock wave generation from continuous initial data by using smooth approximations in the weak sense. The works mentioned above on the delta shock solutions are for particular cases of 2 2 systems. For particular 3 3 non-strictly hyperbolic systems with repeated eigenvalues, we mention the works of Panov and Shelkovich [16], Shelkovich [17] and Cheng [5]. For a strictly hyperbolic system, see [9] for partial results of delta shock solution. For new developments of delta shock waves, its applications and historical notes is [21]. Coming back to system 1), we can deduce that the associated eigenvalues are given by, λ 1 = u s ρ, λ 2 = uandλ 3 = u+ s ρ wherethecorrespondingriemanninvariants arer 1 = s 2 v su, R 2 = v + 1 ρ and R 3 = s 2 v +su. From the expressions for the eigenvalues and the Riemann invariants we obtain λ 1 = R 3 s sr 2, λ 2 = 1 R 3 R 1 ) and λ 3 = sr 2 R 1 s. We can see that system 1) is linearly degenerate. Recently, Lu et al. [15] showed existence of solutions for the Cauchy problem associated to the Suliciu relaxation system 1) with bounded initial data ρ,x),u,x),v,x)) = ρ x),u x),v x)), ρ x), x R, 3) subject to the following conditions: H1) The functions ρ, u and v satisfy c 1 u x) sv x) c 2, c 3 u x)+sv x) c 4, v x)+ 1 ε+ρ x) > c 5, where ε > is a small and the constants c i, i = 1,...,5, satisfy c 5 c 4 c 1 >. H2) The total variations of u x) sv x) and u x)+sv x) are bounded. The existence result for the Cauchy problem 1) 3) includes solutions in vacuum regions [15, Theorem 1]. When ρ x) ρ >, the condition H1 becomes: H1) The functions ρ, u and v satisfy c 1 u x) sv x) c 2, c 3 u x) + sv x) c 4, and v x)+ 1 ρ x) > c 5, where c i, i = 1,...,5, are constants satisfying c 5 c 4 c 1 >. In [9], we showed uniqueness of solutions for the Cauchy problem associated to the Suliciu relaxation system 1) 3) with ρ x) ρ >. Moreover, Theorem 2 in [9] ensures that when we consider v x) = 1 ρ x), problem 1) is reduced to following Chaplygin gas system ρt +ρu) x =, ρu) t +ρu 2 s2 ρ ) x =.

3 The classical Riemann problem associated to the Suliciu relaxation system 1) 4) has been extensively studied, for instance in [1, 3, 4, 9]. In [1], it is shown uniqueness for the generalized Riemann problem for the Suliciu relaxation system. In the present paper, we show the existence and uniqueness of solutions for the Riemann problem associated with the Suliciu relaxation system with initial data ρ l,u l,v l ), if x <, ρ x),u x),v x)) = 4) ρ r,u r,v r ), if x >, in which the left and right constant states ρ l,u l,v l ) and ρ r,u r,v r ), with ρ l,ρ r >, satisfy the conditions H2 and λ 1 ρ l,u l,v l ) λ 3 ρ r,u r,v r ), i.e., they do not satisfy condition H1 globally. Under this assumptions we have the following situations: 1. ρ l,u l,v l ) and ρ r,u r,v r ) satisfy locally the condition H1, i.e., ρ l,u l,v l ) satisfy α 1 u l sv l α 2, α 3 u l sv l α 4 and v l + 1 ρ l > α 5 where α i, i = 1,...,5, are suitable constants satisfying α 5 α 4 α 1 >, ρ r,u r,v r ) satisfy β 1 u r sv r β 2, β 3 u r sv r β 4 and v r + 1 ρ r > β 5 where β i, i = 1,...,5, are suitable constants satisfying β 5 β 4 β 1 >. Let c 1 = minα 1,β 1 } and c 4 = maxα 1,β 1 }. Then there are c i, i = 1,...,4 such that } ul sv c 1 l c u r sv 2 r } ul +sv and c 3 l c u r +sv 4, r but is not possible to find a constant c 5 such that c 5 c 4 c 1 >. 2. Only ρ r,u r,v r ) satisfy locally the condition H1, i.e., ρ l,u l,v l ) satisfy β 1 u l sv l β 2, β 3 u l sv l β 4 and v l + 1 ρ l > β 5 where β i, i = 1,...,5, are suitable constants satisfying β 5 β 4 β Only ρ l,u l,v l ) satisfy locally the condition H1, i.e., ρ r,u r,v r ) satisfy β 1 u r sv r β 2, β 3 u r sv r β 4 and v r + 1 ρ r > β 5 where β i, i = 1,...,5, are suitable constants satisfying β 5 β 4 β Neither ρ l,u l,v l ) nor ρ r,u r,v r ) satisfy the local condition H1.

4 De la cruz et al. [9] studied only the first case. The idea used by the authors in [9] can be extended to the other cases. We considered all the cases for which existence and uniqueness of solutions can be guarantied fortheriemannproblemassociateothesuliciusystem. In[9]itisshownthatiftheinitialdata 3), with ρ x) ρ >, satisfies H1 and H2, then, the Riemann problem has a unique solution. It is also shown that H1 and H2 imply the Lax shock condition λ 1 ρ l,u l,v l ) < λ 3 ρ l,u l,v l ). Therefore, the Riemann problem for the Suliciu relaxation system has a classical solution in the region Γ 1 = ρ,u,v) : ρ > and λ 1 ρ l,u l,v l ) < λ 3 ρ,u,v)}. In this paper we show that the Suliciu relaxation system presents nonclassical solutions in the sense of delta shocks in the region Γ 2 = ρ,u,v) : ρ >,λ 1 ρ l,u l,v l ) λ 3 ρ,u,v) and u l u) 2 s 2 v l v)1/ρ 1/ρ l )}. Itremains tostudyif wecan obtain somekindof solution in theregion, ) R 2 Γ 1 Γ 2 ). This object of our current research. 2 Delta shock solutions Denote by BMR) the space of bounded Borel measures on R. the definition of a measure solution of Suliciu relaxation system in BMR) can be given as follows. Definition 2.1. A tripleρ, u, v) constitutes a measure solution to the Suliciu relaxation system, if it holds that a) ρ L, ),BMR)) C, ),H s R)), b) u L, ),L R)) C, ),H s R)), c) v L loc, ),L loc R)) C, ),H s R)), s >, d) u and v are measurable with respect to ρ at almost for all t, ), and I 1 = R φ t +uφ x )dρ =, I 2 = R uφ t +uφ x )dρ+ R s2 vφ x dx =, I 3 = R vφ t +uφ x )dρ+ R uφ xdx =, for all test function φ C R+ R). 5) Recall that atwo-dimensional weighted delta functionws)δ L supportedonasmoothcurvel parameterizedast = ts), x = xs)c s d)isdefinedby ws)δ L,φt,x) = d c ws)φts),xs))ds for all φ C R2 ). We next present the definition of a delta shock wave. Definition 2.2. A triple distribution ρ,u,v) is called a delta shock wave if it is represented in the form ρ l,u l,v l )t,x), x < xt), ρ,u,v)t,x) = wt)δx xt)),u δ t),gt)), x = xt), 6) ρ r,u r,v r )t,x), x > xt), and satisfies Definition 2.1, where ρ l,u l,v l )t,x) and ρ r,u r,v r )t,x) are piecewise smooth bounded solutions of the Suliciu relaxation system 1).

5 We set dx = u δt) since the concentration in ρ need to travel at the speed of discontinuity. Hence, we say that a delta shock wave 6) is a measure solution to the Suliciu relaxation system 1) if and only if the following relation holds, dxt) = u δ t), dwt) = [ρ]u δ t)+[ρu], dwt)u δ t) = [ρu]u δ t)+[ρu 2 +s 2 7) v], dwt)gt) = [ρv]u δ t)+[ρuv +u]. In fact, for any test function φ C R+ R), from 5), we obtain I 1 = φ t +uφ x )dρ = R I 2 = u δ t)[ρu]+[ρu 2 +s 2 v] dwt)u } δt), I 3 = u δ t)[ρv]+[ρuv +u] dwt)gt) } φ. u δ t)[ρ]+[ρu] dwt) and } φ, Relations 7) are called the generalized Rankine-Hugoniot condition. It reflects the exact relationship among the limit states on two sides of the discontinuity, the weight, propagation speed and the location of the discontinuity. In addition, to guarantee uniqueness, the delta shock wave should satisfy the admissibility entropy) condition λ 3 ρ r,u r,v r ) u δ t) λ 1 ρ l,u l,v l ). Now, the generalized Rankine-Hugoniot condition 7) is applied to the Riemann problem1) 4) with left and right constant states U =,u,v ) and U + = ρ +,u +,v + ), respectively, satisfying the condition H2, the fact λ 3 ρ +,u +,v + ) λ 1,u,v ) and u u + ) 2 s2 ρ + v v + ) s2 v v + ). 8) Therefore, the Riemann problem 1) 4) is reduced to solving 7) with initial data under entropy condition From 7) and 9), it follows that t =, x) =,w) =,g) =, 9) u + + s ρ + u δ t) u s. 1) wt) = [ρ]xt)+[ρu]t, wt)u δ t) = [ρu]xt)+[ρu 2 +s 2 v]t, and wt)gt) = [ρv]xt)+[ρuv +u]t. 11) Multiplying the first equation in 11) by u δ t) and then subtracting it from the second one, we obtain that [ρ]xt)u δ t) [ρu]u δ t)t [ρu]xt)+[ρu 2 +s 2 v]t =, that is, d [ρ] 2 x2 t) [ρu]xt)t+ [ρu2 +s 2 ) v] t 2 =, 2

6 which is equivalent to [ρ]x 2 t) 2[ρu]xt)t+[ρu 2 +s 2 v]t 2 =. 12) From 12), one can find u δ t) := u δ is a constant and xt) = u δ t. Then, 12) can be rewritten [ρ]u 2 δ 2[ρu]u δ +[ρu 2 +s 2 v] =. 13) When [ρ] = ρ + =, the situation is very simple and one can easily calculate the solution u δ = u +u + 2 +s 2 [v] 2 [u], xt) = u δ t, 14) wt) = u u + )t, gt) = [ρuv+u] u δ [ρu], which obviously satisfies the entropy condition 1). From condition 8), s 2[v] 1 2 λ 1U ) λ 3 U + )) 2 < 1 2 [u]λ 1U ) λ 3 U + )) and u δ u s ) = 1 u + + s ) u s ) +s 2 [v] ). 2 [u] Similarly, we can deduce that u δ u + + s ) = 1 2 because u s ) u + + s ) +s 2 [v] ), [u] s 2[v] 1 2 λ 1U ) λ 3 U + )) 2 < 1 2 [u]λ 1U ) λ 3 U + )). When [ρ] = ρ +, the discriminant of the quadratic equation 13) is = 4[ρu] 2 4[ρ][ρu 2 +s 2 v] = ρ + [u] 2 s 2 [ρ][v] > and then we can find u δ = [ρu] /4 [ρ], S ± xt) = [ρu] /4 [ρ] t, ) = wt) = ± /4t, gt) = [ρu][ρv]±[ρv] /4+[ρ][ρuv+u] t. [ρ] /4 15) Next, with the help of the entropy condition 1), we will choose the admissible solution from 15). Observe that by the entropy condition and since the system is strictly hyperbolic, we have that u + s < u + < u + + s u s < u < u + s. ρ + ρ +

7 Observe that, [ρ]λ 1,u,v )+[ρu] = ρ + u s [ρ]λ 3 ρ +,u +,v + )+[ρu] = u + s ) u + s ) ρ + u + + s ρ + )) >, )) >, [ρ]λ 1,u,v )) 2 2[ρu]λ 1,u,v )+[ρu 2 ]+s 2 [v] = ρ + u u + s ) 2 + s2 +s 2 [v], [ρ]λ 3 ρ +,u +,v + )) 2 2[ρu]λ 3 ρ +,u +,v + )+[ρu 2 ]+s 2 [v] = u u + s ) 2 s2 +s 2 [v], ρ + ρ + then, for the solution given in 15) by S + ), we have u δ λ 1,u,v ) and u δ λ 3 ρ +,u +,v + ), which imply that the entropy condition 1) is valid. When λ 1,u,v ) = λ 3 ρ +,u +,v + ), we have trivially that λ 1,u,v ) = u δ = λ 3 ρ +,u +,v + ). Now, for the solution S ), when < ρ + we have u δ λ 3 ρ +,u +,v + ) = λ 3 U ) λ 3 U + ))+ [ρu] 2 [ρ][ρu 2 +s 2 v] [ρ] <, and when > ρ +, that u δ λ 1,u,v ) = ρ +λ 1 U ) λ 1 U + ))+ [ρu] 2 [ρ][ρu 2 +s 2 v] [ρ] >. showing that the solution S ) in 15) does not satisfy the entropy condition 1). Thus we have proved the following result. Theorem 2.1. Given left and right constant states ρ l,u l,v l ) and ρ r,u r,v r ), respectively, such that satisfy the condition H2, λ 1 ρ l,u l,v l ) λ 3 ρ r,u r,v r ) and 8), that is, u l u r ) 2 s 2 v l v r )/ρ r s 2 v l v r )/ρ l. Then, the Riemann problem 1) 4) admits a unique entropy solution in the sense of measures. This solution is of the form ρ l,u l,v l ), if x < u δ t, ρ,u,v)t,x) = wt)δx u δ t),u δ,gt)), if x = u δ t, ρ r,u r,v r ), if x > u δ t, where u δ, wt) and gt) are show in 14) for [ρ] = or S + ) in 15) for [ρ].

8 3 Numerical illustration In this section, we show numerical evidence of delta shock solution for the Suliciu relaxation system using the Lax-Friedrichs method. In the numerical test, with s = 1, we consider the initial data given by 9,5,14/5), if x <, ρ,u,v )x) = 1,3,2), if x >, the spatial discretization parameter for N = 178 points and a constant CFL = The numerical results at final time t =.1 is displayed in Figure 1. The exact solution at time t is 9,5,14/5), if x < u δ t, ρ,u,v)t,x) = Atδx u δ t),u δ,bt), if x = u δ t, 1,3,2), if x > u δ t, with u δ = , A = and B = Conclusions In this paper we obtain delta-shock waves for the Riemann problem associated to the Soliciu relaxation system. This work complements and extend the partial results concerning deltashock waves for the Riemann problem for the Suliciu relaxation system that were obtained in [9]. We mention that in [1, 3] the classical Riemann problem for the Suliciu relaxation system was solved. In our analysis we consider all possible cases for the existence and uniqueness of delta-shock solution for the Suliciu relaxation system and by the Lax-Friedrichs method we show numerical evidence of such solutions. We identified particular regions where we can present classical solutions and delta shock waves. However, we observed that there are regions in which the Riemann problem for this system does not have classical solution neither delta shock waves, which encourages us to seek solutions in some different sense to be able to analyze this missing regions. References [1] F. Bouchut. Nonlinear Stability of Finite Volume Methods for Hyperbolic Conservation Laws and Well-Balanced Schemes for Sources. Front. Math.,. Birkhäuser Verlag, Basel, 24. [2] Y. Brenier. Solutions with concentration to the Riemann problem for one-dimensional Chaplygin gas dynamics. J. Math. Fluid Mech., 7:S326 S331, 25. [3] G. Carbou, B. Hanouzet, and R. Natalini. Semilinear behavior for totally linearly degenerate hyperbolic systems with relaxation. J. Differential Equations., 2461): , 29. [4] C. Chalons and F. Coquel. Navier-Stokes equations with several independent pressure laws and explicit predictor corrector schemes. Numer. Math., 113): , 25.

9 [5] H. Cheng. Riemann problem for one-dimensional system of conservation laws of mass, momentum and energy in zero-pressure gas dynamics. Differ. Equ. Appl., 44): , 212. [6] H. Cheng and H. Yang. Riemann problem for the isentropic relativistic Chaplygin Euler equations. Z. Angew. Math. Phys., 63:429 44, 212. [7] V. G. Danilov and D. Mitrovic. Delta shock wave formation in the case of triangular hyperbolic system of conservation laws. J. Differential Equations., 24512): , 28. [8] V. G. Danilov and V. M. Shelkovich. Delta-shock wave type solution of hyperbolic systems of conservation laws. Q. Appl. Math., 633):41 427, 25). [9] R. De la cruz, J. Galvis, J. Juajibioy, and L. Rendón. Delta shock wave for the Suliciu Relaxation System. Adv. Math. Phys., Art. ID :11 pp, 214. [1] R. De la cruz and J. Juajibioy. On the Riemann problem for a hyperbolic system of Temple class. In G. Tost and O. Vasilieva, editors, Analysis, modelling, optimization, and numerical techniques, Springer Proceedings in Mathematics & Statistics, Vol. 121, 215. [11] G. Ercole. Delta-shock waves as self-similar viscosity limits. Q. Appl. Math., LVIII1): , 2. [12] L. Guo, W. Sheng, and T. Zhang. The two-dimensional Riemann problem for isentropic Chaplygin gas dynamic system. Comm. Pure Appl. Math., 9: , 21. [13] D.-X. Kong and C. Wei. Formation and propagation of singularities in one-dimensional Chaplygin gas. J. Geometry. Phys., 8:58 7, 214. [14] D. Korchinski. Solution of a Riemann problem for a 2 2 system of conservation laws possessing no classical weak solution. PhD thesis, Adelphi University, [15] Y.-G. Lu, C. Klingenberg, L. Rendón, and D.-Y. Zheng. Global solutions for a simplified shallow elastic fluids model. Abstr. Appl. Anal., ID 92248:5 pages, 214. doi: /214/ [16] E. Panov and V. M. Shelkovich. δ -Shock waves as a new type of solutions to systems of conservation laws. J. Differential Equations., 2281):49 86, 26. [17] V. M. Shelkovich. The Riemann problem admitting δ-,δ -shocks, and vacuum states the vanishing viscosity approach). J. Differential Equations., 2312):459 5, 26. [18] W. Sheng and T. Zhang. The Riemann problem for transportation equation in gas dynamics. Technical Report 654, American Mathematical Society, [19] I. Suliciu. On modelling phase transition by means of rate-type constitutive equations. Shock wave structure. Int. J. Engng. Sci., Vol. 28No. 8):pp , 199. [2] D. Tan, T. Zhang, and Y. Zheng. Delta shock waves as limits of vanishing viscosity for hyperbolic systems of conservation laws. J. Differential Equations., 1121):1 32, [21] H. Yang and Y. Zhang. New developments of delta shock waves and its applications in systems of conservation laws. J. Differential Equations., 25211): , 212.

10 a) b) a) b) a) b) Figure 1: Numerical solution. Top: a) Theinitial data ρ x) = ρ,x). b) Numerical solution of ρ at time t =.1. Middle: a) The initial data u x) = u,x). b) Numerical solution of u at time t =.1. Bottom: a) The initial data v x) = v,x). b) Numerical solution of v at time t =.1