Singularity formation for compressible Euler equations

Size: px
Start display at page:

Download "Singularity formation for compressible Euler equations"

Transcription

1 Singularity formation for compressible Euler equations Geng Chen Ronghua Pan Shengguo Zhu Abstract In this paper, for the p-system and full compressible Euler equations in one space dimension, we provide an equivalent and a sharp condition on initial data, respectively, under which the classical solution must break down in finite time. Moreover, we provide time-dependent lower bounds on density for arbitrary classical solutions for these two equations. Our results have no restriction on the size of solutions. 21 Mathematical Subject Classification: 76N15, 35L65, 35L67. Key Words: Singularity formation, compressible Euler equations, p- system, conservation laws, large data. 1 Introduction In this paper, we consider the initial value problem for the compressible Euler equations in Lagrangian coordinates in one space dimension, ( 1 2 u2 + e) τ t u x =, (1.1) u t + p x =, (1.2) t + (u p) x =, (1.3) School of Mathematics, Georgia Institute of Technology, Atlanta, GA 3332 USA (gchen73@math.gatech.edu). School of Mathematics, Georgia Institute of Technology, Atlanta, GA 3332 USA (panrh@math.gatech.edu). Department of Mathematics, Shanghai Jiao Tong University, Shanghai 224, P.R.China (zhushengguo@sjtu.edu.cn). Supported in part by National Natural Science Foundation of China under grant , Natural Science Foundation of Shanghai under grant 14ZR14231 and China Scholarship Council. 1

2 where ρ is the density, τ = ρ 1 is the specific volume, p is the pressure, u is the velocity, e is the specific internal energy, t R + is the time and x R is the spatial coordinate. This model is used to describe the gas dynamics. We assume that the gas is ideal polytropic, so and p = K e S cv τ γ with adiabatic gas constant γ > 1, (1.4) e = pτ γ 1, where S is the entropy, K and c v are positive constants, see [7] or [24]. For most gases, 1 < γ < 3. For C 1 solutions, it follows that (1.3) is equivalent to the entropy equation : S t =. (1.5) When the entropy is constant, the flow is isentropic, then (1.1) and (1.2) become a closed system, known as the p-system: with τ t u x =, (1.6) u t + p x =, (1.7) p = K τ γ, (1.8) where, without loss of generality, we still use K to denote the constant in pressure. Compressible Euler equations and p-system are two of the most important models for hyperbolic conservations laws u t + f(u) x =, (1.9) where u = u(x, t) R n is the unknown vector and f : R n R n is the nonlinear flux. System (1.9) typically admits discontinuity, i.e. shock wave, even when initial data are C. The lack of regularity is the major difficulty in analyzing these systems. Now the well-posedness of small total variation solutions for (1.9) including Compressible Euler equations and p-system are fairly well understood [1, 8]. However, large data results, which means results without restriction on the size of solutions, are still very limited. The main purpose of this paper is to study the breakdown of large data classical solutions for both p-system and full Euler equations, which is related to the formation of shock wave. The previous large data results by 2

3 Lax in [1] for p-system and the first author, R. Young and Q. Zhang in [6] for full Euler equations do not include the most practical case 1 < γ < 3 in gas dynamics. After resolving the case when 1 < γ < 3 in this paper, we will give a complete picture on the mechanism of breakdown of classical solutions for both p-system and full Euler equations. Singularity formation for hyperbolic conservation laws has been studied by a big amount of articles. A survey on the history of these articles can be found in [8]. In one space dimension, the singularity formation for small data solution, i.e. solution around a constant state, has been well understood, where we refer the reader to [9, 1, 12, 13, 18]. For large data problem, scalar conservation law has been fairly completely understood [8]. For uniformly strictly hyperbolic system with two unknowns, Lax provided singularity formation result [1] in His result can be directly applied to p-system (1.6) (1.8) when γ 3. However Lax s proof does not cover the p-system when 1 < γ < 3, because the system might lose its uniformly strict hyperbolicity as density goes to zero in infinite time. Indeed, a Riemann problems connecting two extreme sides of two interacting strong rarefaction waves generates vacuum instantaneously when t >, [24]. Smoothing out this data implies the existence of a C 1 -solution such that inf (x,t) ρ(x, t) as t +. We can also see the difficulty of proving singularity formation from the Riccati equation established by Lax. The p-system satisfies some Riccati equation y = ay 2, (1.1) where y(t) denotes some gradient variable, a is a positive function on density and the derivative is along a characteristic direction. To prove the singularity formation when y() <, i.e. initial data have compression, one needs to show that a dt =, (1.11) where the integral is along a characteristic. For uniformly strictly hyperbolic system or small data problem or large data problem of p-system with γ 3, the leading coefficient in the Riccati equation is uniformly away from zero hence (1.11) is clearly correct. However, for p-system with γ (1, 3), coefficient a is vanishing as density approaches zero. In this paper, when 1 < γ < 3, we establish a time-dependent lower bound on density, using which we prove (1.11) and then the singularity formation when initial data have compression together with Lax s decomposition. Combing with existing results when γ 3, we prove our main 3

4 theorem for p-system: Theorem 2.4. This theorem can be understood as: 1. For p-system (1.6) (1.8) with smooth initial data and γ > 1, classical solution breaks down if and only if the initial data are forward or backward compressive somewhere. Remark 1.1. This theorem gives a complete picture of the mechanism of singularity formation for isentropic gas. Here a wave is compression when gradient variable s x or r x is negative somewhere, where s and r are some Riemann invariants in forward and backward directions respectively which will be specified later. Furthermore, p-system with general pressure law satisfies a similar result, which is given in Theorem 2.9. Then we consider the non-isentropic Euler equations (1.1) (1.4). Variation of entropy makes the extension of large data result from (1.6) (1.8) to (1.1) (1.4) highly nontrivial. For example, in p-system, Riemann invariants are constant along characteristics, while this is not true for full Euler equations. And Riccati equations for full Euler equations are in a more complicated form comparing to (1.1). In [6], the first author, R. Young and Q. Zhang first resolved the singularity formation when γ 3, when initial entropy has finite total variation. In this case, the leading coefficient in the Riccati equation does not vanish when density approaches zero, on the contrary, the leading coefficient vanishes as density approaches infinity. By studying the propagation of Riemann invariants, the authors in [6] established a uniform upper bound on density for smooth solution when initial entropy has finite total variation and then prove the singularity formation result by analyzing the Riccati equations found in [2, 17]. However, in the case when 1 < γ < 3, one meets a more essential difficulty comparing to the case when γ 3, similar as p-system, related to the loss of uniformly strict hyperbolicity near vacuum. In fact, the leading coefficient in the Riccati equation is also vanishing when density approaches zero. One needs some lower bound estimates on density. In this paper, we establish a time-dependent lower bound on density for any smooth solutions in Corollary 3.7, which we believe is the first one for non-isentropic gas dynamics to the limit of our knowledge. Then we prove singularity formation result for 1 < γ < 3. Our result can be understood as (see precise statement at Theorems 3.5 and 3.1): 4

5 2. For compressible Euler equations (1.1) (1.4) with smooth initial data and γ > 1, classical solution breaks down if the initial data satisfy some condition, describing compression, somewhere. Remark 1.2. In Subsection 3.5, we find smooth stationary solutions without compression showing that the conditions on initial data for singularity formation provided in Theorems 3.5 is a sharp condition, under which classical solution must break down in finite time. Shock-free examples in [5, 17, 25] will also give us some clues why this condition is sharp. Furthermore, the result for full Euler equations is consistent with the one for p-system: when the initial entropy oscillation is weak enough, classical solution must break down in finite time even when there are only weak initial compressions. Our results are also correct for Euler equations in Eulerian coordinates, whose classical solution is equivalent to the one in Lagrangian coordinates, c.f. [8, 24]. Note that the singularities we study are corresponding to shock formation, see Remark 2.6 for an explanation. At the time of the blowup, the L norm of u, τ, ρ, e and S are all finite. Our conclusions are all large data results, i.e. there are no restriction on the size of the solutions. Other large data singularity formation results for Euler equations in one space dimension can be found at [4, 5]. Other research works considering the rate on which density approaches zero for specially solution can be found at [14] for p-system and [5] for full system. And there are several papers discussing the difficulty for the well-posedness of Euler equations caused by the vacuum [3, 19]. For compressible Euler equations in multiple space dimensions, there are also some singularity formation results [6, 2, 22, 23]. However, different from the Euler equations in one space dimension, the mechanism of singularity formation in multiple space dimension is still far from fully understood. This paper is divided into three sections. In Section 2, we review Lax s original result in [1], then prove our singularity formation result for p- system when 1 < γ < 3. In Section 3, we prove our results for the compressible Euler equations. 2 Singularity formation for p-system In this section, we consider singularity formation for p-system (1.6) (1.8). The proof of our main theorem (Theorem 2.4) is based on the study of Lax s characteristic decomposition established first for general hyperbolic 5

6 system with two unknowns in [1]. To make this paper self-contained, we first review Lax s decomposition in Subsection 2.1. Then in Subsection 2.2, we apply Lax s decomposition to p-system with γ-law pressure and prove a singularity formation result when γ 3 using similar argument as Lax used for general strictly hyperbolic system with two unknowns. Unfortunately, we could not directly use Lax s argument when 1 < γ < 3, because density might approach zero as time goes to infinity. To overcome this difficulty, we need to find a time-dependent lower bound on density, which will be done in Subsection 2.3. This finally directs to the proof of Theorem 2.4. Finally, in Subsection 2.4, we extend the result for p-system with γ-law pressure to p-system with general pressure. 2.1 Lax s result for system with two unknowns This part is basically taken from Lax s paper [1] in Consider a system of two first-order partial differential equations u t + f x =, v t + g x =, (2.1) where f and g are functions of u and v. Carrying out the differentiation in (2.1), we obtain u t + A u x =, (2.2) where u = ( u v ) ( ) fu f and A = v. Suppose that this system is strictly hyperbolic, i.e. the matrix A has real and distinct eigenvalues λ and µ for relevant values of u and v. Use l λ and l µ to denote the left eigenvectors corresponding to eigenvalues λ and µ, respectively. Multiplying (2.2) by l λ and l µ respectively, we have where we denote g u l λ u =, l µ u =, = t + λ x, = t + µ x. Suppose there exists an integrating factor φ such that w = φ λ l λ u =. (2.3) g v 6

7 Such φ always exists at least locally. Similarly, we have z = for some functions w(u, v) and z(u, v), which are called Riemann invariants along characteristics with characteristic speeds λ and µ, respectively. For general hyperbolic systems with two unknowns, there always exist two Riemann invariants for different families, if we restrict our consideration to the small data solution, i.e. solution which is around a constant state. Furthermore, for p-system, there exist Riemann invariants for both characteristic families for arbitrary solutions. Then we differentiate w = in (2.3) on x, we have Also by (2.3), so w tx + λw xx + λ w w 2 x + λ z w x z x =. (2.4) = z = z (λ µ)z x, Substitute (2.5) into (2.4) and denote then we have α + λ w α 2 + z x = z λ µ. (2.5) α := w x, Denote by h a function of w and z which satisfies Using w = in (2.3), we have Substitute this into (2.6) gives h z = λ z λ µ z α =. (2.6) λ z λ µ. h = h w w + h z z = Multiplying (2.7) by e h and denoting λ z λ µ z. α + λ w α 2 + h α =. (2.7) α := e h α, 7

8 we have with α = a(t) α 2 (2.8) a(t) := e h λ w. This Riccati equation gives us a framework for studying the singularity formation and global existence of classical solutions for hyperbolic system with two unknowns. In fact, we could formally solve gradient variable α along a characteristic with speed λ: 1 α(t) = 1 α() + t a(t) dt where the integral is taken on a characteristic with speed λ. For simplicity, suppose that a(t) is always non-zero, which is also satisfied by the solution of p-system if initially a. Without loss of generality, we only consider the case that a(t) >. If α() <, i.e. initial solution is compressive somewhere in the λ direction, then α(t) breaks down under an extra condition that t a(t) dt =. (2.9) In Lax s original proof of singularity formation result, he only consider the hyperbolic system with uniformly strict hyperbolicity, i.e. characteristic speeds λ and µ are uniformly away from each other. In this case, a(t) has a positive lower bound hence (2.9) is automatically correct. If we restrict our consideration to small data problems, the function a(t) also has positive lower bound if it initially has one. However, for large data problem, a(t) does not in general have positive lower bound, even for one of the most important example of (2.1): p-system (1.6) (1.8) with 1 < γ < 3. Hence, Lax s result only covers special hyperbolic systems with two unknowns when we consider large data problems, such as p-system (1.6) (1.8) with γ 3. To resolve this issue for p-system with 1 < γ < 3, in this paper, we establish a time dependent lower bound on a(t) which helps proving (2.9) hence directs to a singularity formation result. Detail is given in Subsection

9 2.2 Application of Lax s result to p-system We apply Lax s decomposition and singularity formation result to the Cauchy problem of (1.6) (1.8) with smooth initial data u(x, ) and τ(x, ). We use the following coordinates, c.f. [2]. Denote η := τ c dτ = 2 Kγ γ 1 where the nonlinear Lagrangian sound speed c is It follows that τ γ 1 2 >, (2.1) c := p τ = K γ τ γ+1 2. (2.11) τ = K τ η 2 γ 1, p = K p η 2γ γ 1, (2.12) c = p τ = K c η γ+1 γ 1, where K τ, K p and K c are positive constants given by K τ := so that also ( 2 Kγ ) 2 γ 1, K p := K Kτ γ, and K c := Kγ K γ+1 2 τ, (2.13) γ 1 K p = γ 1 2γ K c and K τ K c = γ 1 2. (2.14) In this paper, we always use K with some subscripts to denote positive constants. We will not notify the reader again if there is no ambiguity. The forward and backward characteristics are described by dx dt = c and dx dt = c, and we denote the corresponding directional derivatives along these by + := t + c x and := t c x, respectively. Furthermore, the Riemann invariants are r := u η and s := u + η, (2.15) which satisfy + s = and r =, (2.16) 9

10 respectively. Then we denote gradient variables y := η γ+1 s x and q := η γ+1 r x, and show y and q satisfy following Riccati equations: Lemma 2.1. For C 1 solutions of (1.6)-(1.8), we have where + y = a 2 y 2, (2.17) q = a 2 q 2, (2.18) a 2 := K c γ+1 η 3 γ. (2.19) This lemma is an easy corollary of Lax s decomposition (2.8), where the detail calculation can be found in [2]. Proposition 2.2. (A corollary from [1]) Assume that initial data u(x, ) and τ(x, )) are C 1, u(x, ), τ(x, ), s x (x, ) and r x (x, ) are uniformly bounded above, and τ(x, ) is uniformly away from zero. When γ 3, classical solution of (1.6) (1.8) breaks down if s x (x, ) < or r x (x, ) <. (2.2) Proof. We will show that if s x (x, ) < or r x (x, ) < for some x, then singularity forms in finite time. Without loss of generality, we assume that s x (x, ) <, then y(x, ) < for some x. Then we denote the forward characteristic passing (x, ) as x + (t). By (2.17), where recall 1 y(x + (t), t) = 1 y(x, ) + a 2 = K c γ+1 η 3 γ. t a 2 dt, (2.21) By (2.15) and (2.16), η is uniformly bounded above, so a 2 has a positive constant lower bound when γ 3, hence right hand side of (2.21) approaches zero in finite time, which means singularity happens in finite time. However, when 1 < γ < 3, the function a 2 does not have constant positive lower bound, because the density has no constant positive lower bound. For example, in the interaction between two strong rarefaction simple waves, the density approaches zero as time goes to infinity. 1

11 2.3 Breakdown of classical solutions in p-system when 1 < γ < 3 In this section, we prove the singularity formation for the Cauchy problem in p-system when 1 < γ < 3, by providing a time dependent positive lower bound on density. Before showing the main theorem, we first give a lemma: Lemma 2.3. For C 1 solutions of (1.6)-(1.8), we have a priori bounds { } y(x, t) max, sup{y(x, )} x =: Y and { q(x, t) max, sup x } {q(x, )} =: Q. Proof. This Lemma is easily proved by Lemma 2.1. Then we give the first main theorem in this paper. Theorem 2.4. Assume that initial data (u(x, ), τ(x, )) are C 1, u(x, ), τ(x, ), s x (x, ) and r x (x, ) are uniformly bounded above, and τ(x, ) is uniformly away from zero. Global-in-time classical solution of (1.6) (1.8) with γ > 1 exists if and only if s x (x, ) and r x (x, ), for any x. (2.22) Remark 2.5. At a point (x, t), the solution is said to be forward rarefactive (resp. compressive) if s x (x, t) (resp. s x (x, t) < ); the solution is said to be backward rarefactive (resp. compressive) if r x (x, t) (resp. r x (x, t) < ). Hence the theorem can be understood as that classical global-in-time solution of p-system exists if and only if the initial data are nowhere compressive. If (2.22) is not satisfied, gradient blowup happens in finite time. Proof. Step (1). Sufficiency. This part can be proved by standard local existence and global a priori C 1 estimates argument, c.f. [11], under the help of the lower bound of density provided in [14]. To make this paper self-contained, we give the sketch of proof. One can first prove the local-in-time existence of classical solutions for (1.6) (1.8), c.f. [11], where the life-span depends on the C 1 -norm of u and τ and the positive lower bound on τ in the initial data. Then to extend the local-intime existence of classical solution to global-in-time existence, we only have 11

12 to get global a priori C 1 bounds on u and τ and the positive lower bound on τ for classical solutions. In fact, the upper bound on density and u can be easily found by studying (2.16). Furthermore, Lemma 2.1 and (2.22) tell us that y and q are always nonnegative, while Lemma 2.3 gives us the upper bounds on y and q. Finally, the classical solution also satisfies the time-dependent density lower bound in [14] when (2.22) is satisfied by the weak-strong uniqueness, c.f. [8]. Step (2). Necessity. We only have to consider the case 1 < γ < 3, in which a 2 does not have positive lower bound. In fact, it is enough to show that t lim a 2 dt =. and Hence t Via the definition of c in (2.12), (2.16) and Lemma 2.3, we have so by (2.15), s t = cs x = K c η γ+1 y K c η γ+1 Y, r t = cr x = K c η γ+1 q K c η γ+1 Q. (s r) t K c η γ+1 (Y + Q), η t K c 2 η γ+1 (Y + Q). Dividing the above inequality by η γ+1, then integrating both sides on t, we have 2(γ 1) γ 3 (η(x, t)) + 3 γ 2(γ 1) γ 3 (η(x, )) 3 γ Hence, when 1 < γ < 3, using (2.1) we have τ(x, t) { K [ τ 3 γ K c 4 (x, ) + 2 } 4 3 γ (Y + Q)t]. K c (Y + Q)t. 2 So { [ 3 γ a 2 K K τ K c 4 (x, ) + 2 } 1 (Y + Q)t]. 12

13 Hence, a 2 dt =, since τ 3 γ 4 (x, ) and Y + Q are bounded. So we prove singularity formation when 1 < γ < 3. The proof of the theorem is completed. Remark 2.6. Finally, we give a remark why the singularity in Theorem 2.4 is in fact a shock wave satisfying Lax entropy condition. First, at (x, t ) where the first singularity formation happens, there are some characteristics in the same family interacting with each other. Let s prove it by contradiction. Assume there are no characteristics in the same family interacting with each other at or before time t, then C 1 solution exists when t [, t ], which is contradict to the singularity formation at time t. t (x *, t * ) l1 l 2 S- S+ x Figure 1: Shock formation Hence, without loss of generality, we could find two characteristics l 1 and l 2 interacting with each other (see Figure 1) at (x, t ), and s x < in the region between l 1 and l 2 near (x, t ) on (x, t)-plane, because s x when (x, t) approaches (x, t ). Hence the solution is discontinuous when singularity forms because lim s(x, x x t ) < lim s(x, x x + t ). Finally we check the Lax entropy condition. For smooth solutions before 13

14 blowup, by (1.6) (1.7) and (2.1) (2.16), cs x = s t = u t + η t = p x + η t = η. So c + when (x, t) approaches (x, t ), hence the solution is discontinuous when singularity forms, and the Lax entropy condition is satisfied on the discontinuity. 2.4 p-system with general pressure law In this subsection, we consider the p-system (1.6) (1.7) with general C 3 pressure p(τ) satisfying p τ <, p ττ > (2.23) and lim p(τ) =, lim τ p(τ) = and pτ dτ < (2.24) τ 1 where condition (2.23) is dictated by physics when one uses (1.6) (1.7) to model gas dynamics, c.f. [21]. Furthermore, we assume that 1 pτ dτ = (2.25) which includes the γ-law pressure case. Applying Lax s decomposition to this case we easily have the following proposition, where the detail calculation can be found in [4] Proposition 2.7. [4] The smooth solutions of (1.6) (1.7) satisfy where y = a(τ)y 2, (2.26) q = a(τ)q 2, (2.27) a(τ) := p ττ 4( p τ ) 5 4 >, (2.28) and y := c s x, q := c r x, (2.29) 14

15 with Lagrangian wave speed c c(v) = p τ, and Riemann invariants s := u + 1 τ c(τ)dτ and r := u 1 τ c(τ)dτ. Furthermore, s t + cs x = and r t cr x =. (2.3) Secondly, we have the following lemma as same as Lemma (2.3). Lemma 2.8. For C 1 solutions of (1.6)-(1.7) and (2.23), we have a priori bounds { } y(x, t) max, sup{y(x, )} =: Y and x { } q(x, t) max, sup{q(x, )} x =: Q. Then we could state our theorem for the general pressure law case. Theorem 2.9. Assume that initial data (u(x, ), τ(x, )) are C 1, u(x, ), τ(x, ), s x (x, ) and r x (x, ) are uniformly bounded above, and τ(x, ) is uniformly away from zero. The pressure satisfies (2.23) (2.25). Furthermore, assume there exists some positive constant A, such that for any τ >, (5 + A)(p ττ ) 2 4p τ p τττ. (2.31) Then global-in-time classical solution of (1.6) (1.7) exists if and only if s x (x, ) and r x (x, ), for any x. (2.32) Remark 2.1. It is clear that for the singularity formation, conditions (2.24) (2.25) are not necessary. Condition (2.31) is a fairly mild assumption because the constant A can be arbitrarily large. For example, the γ-law pressure p = kτ γ with γ > satisfies conditions (2.31) and (2.23), and the pressure p = kτ γ with γ > 1 satisfies conditions (2.31) and (2.23) (2.25). Proof. We first remark on a fact for later references. By (2.3), (2.25), it is easy to see that τ has a uniform lower bound. We denote that τ min = min τ(x, t) (x,t) R R + 15

16 which is a positive constant only depending on the maximum values of s(x, ) and r(x, ). Then it is easy to check that c has a uniform upper bound by (2.23). If condition (2.32) holds, the global existence could be proved in an entirely similar way as we introduced in Theorem 2.4 together with the lower bound on density provided in [14]. If condition (2.32) fails, by a similar argument in Theorem 2.4, in order to prove singularity formation in finite time, we only have to show Hence it is sufficient to show that a(τ(x(t), t) dt =. (2.33) a(τ(x, t)) = 4( p τ ) K 1 + K 2 t (2.34) p ττ for some positive K 1 and K 2. Then we prove (2.34). By Proposition 2.7, we have then by Lemma 2.8 we have ( τ Hence τ(x,t) τ min τ min ( p τ (τ)) 1 4 ( p τ (τ)) 1 4 dτ 1 2 (y + q) = 1 c(sx + r x ) = c u x = c τ t, 2 ( dτ )t τ = ) c(τ) dτ = 1 τ min t 2 (y + q) 1 (Y + Q) 2 τ(x,) τ min ( p τ (τ)) 1 4 dτ (Y +Q)t K 3+K 4 t (2.35) for some positive constants K 3 and K 4. Compare (2.34) and (2.35), and use τ > τ min >, then it is easy to see that in order to get (2.34) we only have to show that ( 4( pτ ) 5 4 p ττ ) τ A( p τ (τ)) 1 4, (2.36) for some positive constant A. In fact, inequality (2.36) means that the left hand side of (2.34) grows slower than the left hand side of (2.35) as τ increases. It is easy to see that (2.36) is correct because of (2.31). Hence we finish the proof of this Theorem. 16

17 3 Compressible Euler equations In this section, we consider the Cauchy problem of compressible Euler equations (1.1) (1.4) with given smooth initial data (u(x, ), τ(x, ), S(x, )). We will establish a time-dependent lower bound on density for smooth solutions, then prove the main singularity formation results: Theorems 3.5 and 3.1. We will provide our result through several steps. We first introduce the coordinates and equations which are basically from [2], then review the uniform upper bound on density and velocity established in [6], and finally provide the lower bound on density and prove Theorem 3.5 when initial entropy has finite total variation and Theorem 3.1 when entropy has infinite total variation. Finally, in Subsection 3.5, we explain why the condition on initial data for singularity formation in our Theorems are the best we could expect. 3.1 Equations and coordinates We use the coordinates used in [2]. And all detail calculation in this subsection can be found in [2] or [6]. Define new variables (m, η) for (S, τ), by m := e S 2cv > (3.1) and η := τ c m dτ = 2 Kγ γ 1 where the nonlinear Lagrangian sound speed c is τ γ 1 2 >, (3.2) c := p τ = K γ τ γ+1 2 e S 2cv. (3.3) Without confusions, we still use η and c for full Euler equations. In fact, functions η and c in p-system equals to η and c for full Euler equations when m = 1, respectively. Similarly, in this section, if we use the same letter as in p-system to denote a function, this function is extended form the corresponding function in p-system, which should have no ambiguity. It follows that τ = K τ η 2 γ 1, p = K p m 2 η 2γ γ 1, (3.4) c = c(η, m) = K c m η γ+1 γ 1. 17

18 In these coordinates, for C 1 solutions, equations (1.1) (1.3) are equivalent to η t + c m u x =, (3.5) u t + m c η x + 2 p m m x =, (3.6) m t =, (3.7) where the last equation comes from (1.5), which is equivalent to (1.3), c.f. [8, 24]. Note that, while the solution remains C 1, m = m(x) is given by the initial data and can be regarded as a stationary quantity. We denote the Riemann invariants by r := u m η, s := u + m η. (3.8) Different from the isentropic case (m constant), for general non-isentropic flow, s and r vary along characteristics. The forward and backward characteristics are described by dx dt = c and dx dt = c, (3.9) and we denote the corresponding directional derivatives along these by + := t + c x and := t c x, respectively. Using (3.8), (3.4) and (2.14), equations (3.5) and (3.6) give Then we introduce gradient variables + s = 1 c m x (s r), 2γ m (3.1) r = 1 c m x (s r). 2γ m (3.11) y := m 3(3 γ) 2(3γ 1) η γ+1 ((u + m η) x 2 3γ 1 m x η) and q := m 3(3 γ) 2(3γ 1) η γ+1 ((u m η) x + 2 3γ 1 m x η), (3.12) and derive Riccati type equations for their evolution: Lemma 3.1. [2] For C 1 solutions of (1.1) (1.4), we have + y = a a 2 y 2, q = a a 2 q 2, (3.13) 18

19 where a := Kc γ Furthermore, [ γ 1 3γ 1 m m xx (3γ+1)(γ 1) m 2 ] (3γ 1) 2 x m 3(3 γ) 2(3γ 1) η 3(γ+1) +1, γ+1 a 2 := K c m 3(3 γ) 2(3γ 1) η 3 γ. (3.14) y or q iff u x or τ x. (3.15) The proof of this lemma can be found in [2]. The decomposition (3.13), which is generalized from the Lax s decomposition in [1] for hyperbolic system with two unknowns, was first provided by [17], and then found in [2] independently through another explanation. 3.2 Uniform upper bound on density In this part, we review a result on the uniform upper bounds for u and ρ provided by G. Chen, R. Young and Q. Zhang in [6], for later references. Assume that the initial entropy S(x) is C 1 and has finite total variation, so that V := 1 2c v while also, by (3.1), + S (x) dx = + m (x) m(x) dx <, (3.16) < M L < m( ) < M U, (3.17) for some constants M L and M U. Also, we assume that ρ > and u are bounded above initially. Hence, there exist positive constants M s and M r, such that, in the initial data, s ( ) < M s and r ( ) < M r. (3.18) In this section, we always assume (3.16) (3.18). We define two useful constants by N 1 := M s + V M r + V (V M s + V 2 M r ) e V 2, N 2 := M r + V M s + V (V M r + V 2 M s ) e V 2, where V := V 2γ, which clearly depend only on the initial data. By below proposition in [6], u and ρ are shown to be uniformly bounded above. 19

20 Proposition 3.2. [6] Assume system (1.1) (1.4), with provided initial data satisfying (3.16) (3.18), has a C 1 solution when t [, T ), then one has the uniform bounds u(x, t) N 1 + N 2 1 M U 2 2γ and ρ(x, t) N 1 + N 2 2 M L 1 2γ 1, (3.19) where T can be any positive number or infinity. And the bounds are independent of T. 3.3 Singularity formation when entropy has finite total variation For full Euler equations, the major difficulty we need to overcome in this paper is still how to find the time dependent lower bound on density when 1 < γ < 3. We first give a lemma to define a positive constant N for later use. Lemma 3.3. Assume the initial data u(x, ) and τ(x, ) are C 1 and uniformly bounded above, the initial entropy S(x, ) is C 2 and has bounded variation. Furthermore, suppose there is a positive constant M such that the initial entropy satisfies m (x) < M. Then, for C 1 solution of system (1.1) (1.4), there exists a positive constant N depending only on the initial data, such that a a 2 N if a. (3.2) Furthermore, if y or q is larger than N or less than N at (x, t), then + y = a a 2 y 2 < or q = a a 2 q 2 < at (x, t), (3.21) respectively. Proof. Clearly, the assumptions in Proposition 3.2 are also satisfied, by the assumptions in this lemma. Then by Proposition 3.2 and (3.2), we know η has a global upper bound depending only on the initial data, denoted via E U. By (3.14), it is easy to calculate that, if a, a = a 2 2 γ(γ+1)(3γ 1) ( m mxx 3γ+1 ) γ+1 3γ 1 m2 x η +1 m 3(3 γ) 2(3γ 1), (3.22) 2

21 which implies the uniform bound a /a 2 N, where 3γ 1 3γ 5 2 γ(γ+1)(3γ 1) N := M 3γ 1 EU ML, 1 < γ 5/3, 3γ 1 3γ 5 2 γ(γ+1)(3γ 1) M 3γ 1 EU MU, γ 5/3, (3.23) where M L and M U are defined in (3.17). By (3.2) and (3.13), (3.21) is clearly correct. Next we prove a key lemma. Lemma 3.4. For C 1 solutions of (1.1) (1.4) with initial data satisfying (3.16) (3.18) and m xx is uniformly bounded above, we have { } y(x, t) max N, sup{y(x, )} x =: Ȳ, { } and q(x, t) max N, sup{q(x, )} x =: Q. Proof. Without loss of generality, we only prove the inequality for y. Then the inequality for q will be proved in an entirely same way. We prove the inequality for y by contradiction. Assume that { } y(x, t ) > max N, sup{y(x, )} x. (3.24) We use Γ(t) with t [ t, t ] to denote the largest connected piece of forward characteristic x + (t) containing (x, t ) as its upper endpoint, such that y(x + (t), t) y(x, t ) for any points (x + (t), t) on Γ(t). Since { } y(x + (), ) max N, sup{y(x, )} x < y(x, t ), hence t >, then by the definition of Γ(t), we have y(x + ( t ), t ) = y(x, t ). 21

22 On the other hand, by the note before this Lemma, we know d dt y(x+ (t), t) <, when t = t, hence t < t and there exists some ˆt ( t, t ) such that So y(x + (ˆt ), ˆt ) > y(x, t ). y(x + ( t ), t ) = y(x, t ) < y(x + (ˆt ), ˆt ). However, this is impossible because y is decreasing on t on Γ with t [ t, t ] by (3.21) and (3.24). Hence we find a contradiction. So the lemma is proved. This lemma is consistent with Lemma 2.3 for p-system. In fact, in p- system, N =. Finally we prove our main theorem for full compressible Euler equations. Theorem 3.5. Assume that the initial data u(x, ) and τ(x, ) are C 1, u(x, ), τ(x, ), y(x, ) and q(x, ) are uniformly bounded above, τ(x, ) is uniformly away from zero, and S(x, ) is C 2 and has bounded variation. Furthermore, suppose there is a positive constant M such that the initial entropy satisfies m (x) < M. Then, for system (1.1) (1.4), there exists a positive constant N defined in (3.23) which depends only on the initial data, such that, if the initial data satisfy { } inf y(, ), q(, ) < N, (3.25) then u x and/or τ x blow up in finite time. Remark 3.6. When γ 3, this theorem has already been proved in [6]. To make the paper self-contained, we will still give the proof of this case. Proof. Suppose that (3.25) holds. Without loss of generality, we can assume that inf y < N. In fact the case when inf q < N can be proved in an entirely same way. Then there exist ε > and x such that y(x, ) < (1 + ε) N. (3.26) 22

23 Now considering the solution y(t) on the forward characteristic Γ x (t) starting at (x, ). Along this characteristic Γ x (t), by Lemma 3.3 and recall (3.26), we have + y(t) < and y(t) < (1 + ε) N for any t, which together with (3.2) implies that for all t, hence by (3.13) and a 2 >, Integrating it, we get a a 2 + y(t) = a a 2 y 2 (t) < 1 y(t) 1 y() + y 2 (t) (1 + ε) 2 <, t ( (1+ε) 2 ) a 2 y 2 (t) <. where the integral is along the forward characteristic. ( 1 1 (1+ε) 2 ) a 2 dt, (3.27) Case when γ 3. By (3.14) and Proposition 3.2, a 2 is positive and bounded above, so the right hand side of (3.27) approaches zero in finite time. This implies that y(t) approaches in finite time, so that τ x and/or u x blow up. This is a known result in [6] as discussed in Remark 3.6. Case when 1 < γ < 3. To prove the singularity formation in finite time, we only have to show a 2 dt =, by establishing a similar time dependent lower bound on density as the one for p-system. In fact, by (3.1) (3.11) and Lemma 3.4, s t = c s x + 1 c m x (s r) 2γ m = K c m 3γ+7 2(3γ 1) η γ+1 2 y K c 3γ 1 m xmη 2γ γ γ K 2 η γ+1 Ȳ K 3 η γ γ 23 c m x m (s r), c m x m (s r)

24 where to get positive constants K 2 and K 3, we need to use the uniform upper bound on density and bounds on initial entropy. Similarly r t = c r x + 1 c m x (s r) 2γ m = K c m 3γ+7 2(3γ 1) η γ+1 2 q K c 3γ 1 m xmη 2γ γ γ K 2 η γ+1 Q + K3 η γ γ c m x m (s r). Clearly, Ȳ and Q are both finite constants. Then we have So c m x m 2m η t = (s r) t [K 2 (Ȳ + Q) + 2K 3 ] η γ+1. η t K 4 [K 2 (Ȳ + Q) + 2K 3 ] η γ+1 (s r) because m defined in (3.1) has positive lower bound. Similar as in p-system, we have when 1 < γ < 3, τ(x, t) {K 5 [τ 3 γ 4 (x, ) + K4 ( K2 (Ȳ + Q) + 2K 3 ) t ]} 4 3 γ. Hence a 2 K 6 { τ 3 γ 4 (x, ) + K4 ( K2 (Ȳ + Q) + 2K 3 ) t } 1, (3.28) which gives t lim a 2 dt =. t Hence, singularity forms in finite time. Corollary 3.7. Suppose all assumptions in Theorem 3.5 hold. Then when 1 < γ < 3, τ(x, t) {K 5 [τ 3 γ 4 (x, ) + K4 ( K2 (Ȳ + Q) + 2K 3 ) t ]} 4 3 γ See the definitions of positive constants of K 2 K 5, Ȳ and Q in the proof of Theorem

25 3.4 Singularity formation when entropy has infinite total variation In this section, we consider the singularity formation when the entropy has infinite total variation on the whole real line but finite total variation on any bounded interval, i.e. entropy is local BV, which includes interesting cases such as periodic solutions. Since when we prove (3.19) for ρ(x, t) and u(x, t) in [6], we only have to consider the domain of dependence of the point (x, t). Hence, we can get similar upper bounds on density and velocity as those in (3.19) in below remark on a domain of dependence, when initial entropy is only local BV. Remark 3.8. Consider the domain of dependence of a point (X x,z, T x,z ), in Figure 2, which is denoted by Ω x,z, where [x, z] is the intersection interval of the domain of dependence and the initial line t =. t (Xx,z, Tx,z) Ωx,z x z x Figure 2: A domain of dependence Ω x,z We still assume (3.17) (3.18). But instead of (3.16), we suppose that S is C 1 and local BV, which shows V x,z := 1 2c τ z x S ( x) d x = z x m ( x) m( x) d x <. (3.29) By the same proof in [6], it is easy to get, for any point ( x, t) in Ω x,z with u( x, t) N 1x,z + N 2x,z 1 M U 2 2γ and ρ( x, t) N 1x,z + N 2x,z 2 N 1x,z := M s + V x,z M r + V x,z ( V x,z M s + V 2 x,z M r ) e V 2 x,z, N 2x,z := M r + V x,z M s + V x,z ( V x,z M r + V 2 x,z M s ) e V 2 x,z, M L 1 2γ 1 (3.3) 25

26 where V x,z := V x,z 2γ. Remark 3.9. Still under the assumptions and notations in Remark 3.8, we introduce some new notations for later reference. We denote the maximum density value in Ω x,z to be ρ x,z which satisfies ρ x,z N 1x,z + N 2x,z 1 M 2γ 1 L. 2 Denote the maximum η value in Ω x,z to be E U x,z which satisfies E U x,z 2 Kγ γ 1 ρ γ 1 2 x,z (3.31) by (3.2). Denote the time at the upper vertex of Ω x,z as T x,z which satisfies 1 2 T x,z γ+1 γ 1 z x K c E U x,z M U, by (3.9). And if further assume that m ( x) < M for any x [x, z], then 3γ 1 3γ 5 ( 2 a a 2 )x,z N γ(γ+1)(3γ 1) x,z := M 3γ 1 E U x,z ML, 1 < γ 5 3, 3γ 1 3γ 5 2 γ(γ+1)(3γ 1) M 3γ 1 E U x,z MU, γ 5 3, (3.32) ( ) where a a denotes the maximum value of a 2 x,z a 2 in Ω x,z when a. ( If a is always negative in Ω x,z, define a a 2 )x,z =. When 1 < γ < 3, assuming initial density has uniform lower bound, by a similar argument as for (3.28), we have in Ω x,z a 2 K 7x,z (1 + K 8x,z t) 1, (3.33) where K 7x,z > and K 8x,z > depend on the initial bounds in (3.17) (3.18), lower bound on initial density and initial bound V x,z defined in (3.29) for entropy. Theorem 3.1. Assume the initial data u(x, ) and τ(x, ) are C 1, u(x, ), τ(x, ), y(x, ) and q(x, ) are uniformly bounded above, τ(x, ) is uniformly away from zero, and S(x, ) is C 2 and local BV. Furthermore, suppose there is a positive constant M such that the initial entropy satisfies m (x) < M. 26

27 Then, for system (1.1) (1.4), if there exists some interval (x, z) such that the initial data satisfy y(x, ) N x,z (1 + B x,z ), (3.34) where N x,z is defined in (3.32) and B x,z satisfies B x,z (2 + B x,z ) (1 + B x,z ) > ( K c(γ+1) N 3(3 γ) x,zt x,z M U 3 γ 2(3γ 1) E U x,z ) 1, γ 3, ( K7x,z K 8x,z N x,z ln(1 + K 8x,z T x,z )) 1, 1 < γ < 3, (3.35) then u x and/or τ x blow up in finite time. See Remarks 3.8, 3.9 and Figure 2 for definitions of notations. Symmetric result holds in the backward direction for q. Remark The right hand side of (3.35) only depends on the initial data. For any given local BV initial entropy, condition (3.35) will be satisfied when B x,z is large enough, i.e. y(x, ) is negative enough. This means that singularity forms in finite time when the initial compression is strong enough somewhere. One sufficient condition on B x,z such that (3.35) is satisfied is that B x,z > ( K c(γ+1) N 3(3 γ) x,zt x,z M U 2(3γ 1) E U x,z ( K7x,z K 8x,z N x,z ln(1 + K 8x,z T x,z ) ) 3 γ 1, γ 3, (3.36) ) 1, 1 < γ < 3, This result is consistent with Theorem 3.5. In fact, when the initial entropy has finite total variation, T x, = while N x,, K 7x, and K 8x, are all finite, so B x, can be arbitrarily small. Hence, if y(x, ) < N x,, then blowup happens in finite time. Proof. We only consider the solution in the domain of dependence Ω x,z, and prove that singularity formation happens in Ω x,z. More precisely, we will show that y goes to negative infinity along the forward characteristic starting from the point (x, ) before T x,z. We still use y(t) to denote this characteristic. From now on, we restrict our consideration in Ω x,z. Then we can use all bounds in Remark 3.9. By (3.13), (3.32) and (3.34), we know + y(t) < and y(t) N x,z (1 + B x,z ), for any < t < T x,z. 27

28 Still by (3.32), when < t < T x,z, so a a 2 y 2 (t) (1 + B x,z ) 2 <, + y(t) = a a 2 y 2 (t) < B x,z(2 + B x,z ) (1 + B x,z ) 2 a 2 y 2 (t) <. Integrating it, we get 1 y(t) 1 y() + B x,z(2 + B x,z ) (1 + B x,z ) 2 t a 2 dt. (3.37) Hence the blowup happens when the right hand side of (3.37) equals to zero, i.e. when 1 y() = B t x,z(2 + B x,z ) (1 + B x,z ) 2 a 2 dt. (3.38) Now to complete the proof of the Theorem 3.1, the only thing left to show is that blowup happens before T x,z. By (3.34) we only have to show that 1 B x,z(2 + B x,z ) N x,z (1 + B x,z ) Finally, we prove (3.39) case by case. Case when γ 3. By (3.31), in Ω x,z, γ+1 a 2 K c M 3(3 γ) U T x,z a 2 dt. (3.39) 3 γ 2(3γ 1) E U x,z. Hence by (3.35), clearly (3.39) is satisfied. Case when 1 < γ < 3. By (3.33), we know to prove (3.39), we only have to show that 1 B x,z(2 + B x,z ) N x,z (1 + B x,z ) T x,z K 7x,z (1 + K 8x,z t) 1 dt, (3.4) which is clearly correct by (3.35). Hence we complete the proof of this theorem. 28

29 3.5 Further discussion Finally we give an example to show the sharpness of the Theorem 3.5. We consider a global smooth stationary solution which satisfies u =, m = m(x) and τ = τ(x) p x (x) =. (3.41) Then we provide the profile of m(x), satisfiing all conditions in Theorem 3.5, by which we can get τ(x) using (3.41) and (3.4). First, by (3.41) and (3.4), q(x) = y(x) = γ 1 γ(3γ 1) m xm 3(γ 3) 2(3γ 1) η γ (3.42) We note that N is the best estimate we could have now for the upper bound of a a 2 = 2 ( γ(γ+1)(3γ 1) m mxx 3γ+1 ) γ+1 3γ 1 m2 x η +1 m 3(3 γ) 2(3γ 1), (3.43) by (3.22), when m m xx 3γ+1 3γ 1 m2 x. Comparing (3.42) with (3.43), we see if γ 1 γ(3γ 1) m x = 2 γ(γ+1)(3γ 1) ( m mxx 3γ+1 ) 3γ 1 m2 x which is equivalent to or 6γ 2 +3γ+1 γ(3γ 1) m2 x = 2mm xx S xx = 5γ+1 4c vγ(3γ 1) S2 x, (3.44) where S is the entropy satisfying (3.1), then either y(x) = q(x) = a (x) a 2 (x). a (x) a 2 (x) It is clear that we could find a positive solution S(x) of (3.44) in the region x [1, 2], once we set S(1) to be a number large enough. Then choose the function S(x) to be almost constant in either (, 1) or (2, ). Finally we get τ(x) by (3.41) and (3.4). And it is easy to see that this solution is a global smooth solution. Now we already found a smooth stationary solution, satisfying either y(x) = a (x) a 2 (x) or q(x) = a (x) a 2 (x) for each x [1, 2]. On the other hand, 29 or

30 y(x), q(x) and a (x) a 2 (x) are all almost zero in the set outside x (1 ε, 2+ε), for a very small constant ε, by (3.42) (3.43) and m x is almost zero in that set. In conclusion, inf ( y, q ) min a (x,t) R R + a 2. Till now, N, provided by the initial data, is the best estimate for min (x,t) R R + a a 2. If one could not improve this estimate, then the condition (3.25) provides a sharp condition, under which classical solution must break down in finite time. Finally, when the maximum oscillation of initial entropy approaches zero, M approaches zero, then N and N x,z tend to zero. When m is constant, N = and N x,z =. Hence, Theorems 3.5 and 3.1 for non-isentropic Euler equations are consistent with Theorem 2.4 for p-system: when the initial entropy oscillation is weak enough, classical solution must break down in finite time even there are only weak initial compressions. Acknowledgement We sincerely appreciate Professor Alberto Bressan for his very helpful suggestions and discussions when we wrote this paper. References [1] A. Bressan, Hyperbolic systems of conservation laws: the 1-dimensional Cauchy problem, (Oxford Univ. Press, Oxford 2). [2] G. Chen, Formation of singularity and smooth wave propagation for the non-isentropic compressible Euler equations, J. Hyperbolic Differ. Equ., 8:4 (211), [3] G. Chen, H. K. Jenssen, No TVD fields for 1-D isentropic gas flow, Comm. Partial Differential Equations, 38:4 (213), [4] G. Chen and R. Young, Smooth solutions and singularity formation for the inhomogeneous nonlinear wave equation, J. Differential Equations, 252:3 (212),

31 [5] G. Chen and R. Young, Shock free solutions for the compressible Euler equations, to appear in Arch. Ration. Mech. Anal.. [6] G. Chen, R. Young and Q. Zhang, Shock formation in the compressible Euler equations and related systems, J. Hyperbolic Differ. Equ., 1:1 (213), [7] R. Courant and K. O. Friedrichs, Supersonic flow and shock waves, (Wiley-Interscience, New York, 1948). [8] C. Dafermos, Hyperbolic conservation laws in continuum physics, Third edition, (Springer-Verlag, Heidelberg 21). [9] F. John, Formation of singularities in one-dimensional nonlinear wave propagation, Comm. Pure Appl. Math., 27 (1974) [1] P. Lax, Development of singularities of solutions of nonlinear hyperbolic partial differential equations, J. Mathematical Phys., 5:5 (1964) [11] T.-T. Li, Global classical solutions for quasilinear hyperbolic systems, (Wiley, New york, 1994). [12] T.-T. Li, Y. Zhou and D.-X. Kong, Weak linear degeneracy and global classical solutions for general quasilinear hyperbolic systems,comm. Partial Differential Equations, 19:7-8 (1994) [13] T.-T. Li, Y. Zhou and D.-X. Kong, Global classical solutions for general quasilinear hyperbolic systems with decay initial data, Nonlinear Analysis, Theory, Methods & Applications, 28:8 (1997) [14] L. Lin, On the vacuum state for the equations of isentropic gas dynamics, J. Math. Anal. Appl., 121:2 (1987), [15] L. Lin, Vacuum states and equidistribution of the random sequence for Glimm s scheme, J. Math. Anal. Appl., 124:1, (1987), [16] L. Lin; Y. Zheng, Existence and nonexistence of global smooth solutions for quasilinear hyperbolic systems. Chinese Ann. Math. Ser. B 9:3 (1988), [17] L. Lin, H. Liu and T. Yang, Existence of globally bounded continuous solutions for nonisentropic gas dynamics equations, J. Math. Anal. Appl., 29 (1997)

32 [18] T.-P. Liu, The development of singularities in the nonlinear waves for quasi-linear hyperbolic partial differential equations, J. Differential Equations, 33 (1979) [19] T.-P. Liu and J. Smoller, On the vacuum state for the isentropic gas dynamics equations. Adv. in Appl. Math. 1:4 (198), [2] T. Makino, S. Ukai and S. Kawashima, Sur la solution à support compact de equations d Euler compressible, Japan J Appl Math, 33 (1986) [21] R. Menikoff and B. J. Plohr The Riemann problem for fluid flow of real materials, Reviews of Modern Physics, 61:1 (1989), [22] M. A. Rammaha, Formation of singularities in compressible fluids in two-space dimensions, Proc. Amer. Math. Soc., 17:3 (1989) [23] T. Sideris, Formation of singularities in three-dimensional compressible fluids, Comm. Math. Phys., 11 (1985) [24] J. Smoller, Shock waves and reaction-diffusion equations, (Springer- Verlag, New York 1982). [25] B. Temple and R. Young, A paradigm for time-periodic sound wave propagation in the compressible Euler equations, Methods Appl. Anal., 16:3 (29)

Shock formation in the compressible Euler equations and related systems

Shock formation in the compressible Euler equations and related systems Shock formation in the compressible Euler equations and related systems Geng Chen Robin Young Qingtian Zhang Abstract We prove shock formation results for the compressible Euler equations and related systems

More information

Author(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)

Author(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1) Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL

More information

Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations

Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations Jiequan Li 1 Department of Mathematics, Capital Normal University, Beijing, 100037 Tong Zhang Institute

More information

Convergence to the Barenblatt Solution for the Compressible Euler Equations with Damping and Vacuum

Convergence to the Barenblatt Solution for the Compressible Euler Equations with Damping and Vacuum Arch. Rational Mech. Anal. 176 (5 1 4 Digital Object Identifier (DOI 1.17/s5-4-349-y Convergence to the Barenblatt Solution for the Compressible Euler Equations with Damping Vacuum Feimin Huang, Pierangelo

More information

Applications of the compensated compactness method on hyperbolic conservation systems

Applications of the compensated compactness method on hyperbolic conservation systems Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,

More information

L 1 stability of conservation laws for a traffic flow model

L 1 stability of conservation laws for a traffic flow model Electronic Journal of Differential Equations, Vol. 2001(2001), No. 14, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu ftp ejde.math.unt.edu (login:

More information

Simple waves and characteristic decompositions of quasilinear hyperbolic systems in two independent variables

Simple waves and characteristic decompositions of quasilinear hyperbolic systems in two independent variables s and characteristic decompositions of quasilinear hyperbolic systems in two independent variables Wancheng Sheng Department of Mathematics, Shanghai University (Joint with Yanbo Hu) Joint Workshop on

More information

Piecewise Smooth Solutions to the Burgers-Hilbert Equation

Piecewise Smooth Solutions to the Burgers-Hilbert Equation Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang

More information

Generic Singularities of Solutions to some Nonlinear Wave Equations

Generic Singularities of Solutions to some Nonlinear Wave Equations Generic Singularities of Solutions to some Nonlinear Wave Equations Alberto Bressan Deartment of Mathematics, Penn State University (Oberwolfach, June 2016) Alberto Bressan (Penn State) generic singularities

More information

EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem

EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem Electronic Journal of Differential Equations, Vol. 207 (207), No. 84, pp. 2. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS

More information

Global existence of smooth solutions to two-dimensional compressible isentropic Euler equations for Chaplygin gases

Global existence of smooth solutions to two-dimensional compressible isentropic Euler equations for Chaplygin gases Global existence of smooth solutions to two-dimensional compressible isentropic Euler equations for Chaplygin gases De-Xing Kong a Yu-Zhu Wang b a Center of Mathematical Sciences, Zhejiang University Hangzhou

More information

On the Front-Tracking Algorithm

On the Front-Tracking Algorithm JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 7, 395404 998 ARTICLE NO. AY97575 On the Front-Tracking Algorithm Paolo Baiti S.I.S.S.A., Via Beirut 4, Trieste 3404, Italy and Helge Kristian Jenssen

More information

The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method

The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method Alexis Vasseur, and Yi Wang Department of Mathematics, University of Texas

More information

u-= (u, v), x>o, j u, (u,, v,), x<o, U(X 0) (1) (1), A A2 only when u 0, in which case A 0. THE RIEMANN PROBLEM NEAR A HYPERBOLIC SINGULARITY II*

u-= (u, v), x>o, j u, (u,, v,), x<o, U(X 0) (1) (1), A A2 only when u 0, in which case A 0. THE RIEMANN PROBLEM NEAR A HYPERBOLIC SINGULARITY II* SIAM J. APPL. MATH. Vol. 48, No. 6, December 1988 (C) 1988 Society for Industrial and Applied Mathematics 006 THE RIEMANN PROBLEM NEAR A HYPERBOLIC SINGULARITY II* E. ISAACSONS" AN[) B. TEMPLE:I: Abstract.

More information

Decay rate of the compressible Navier-Stokes equations with density-dependent viscosity coefficients

Decay rate of the compressible Navier-Stokes equations with density-dependent viscosity coefficients South Asian Journal of Mathematics 2012, Vol. 2 2): 148 153 www.sajm-online.com ISSN 2251-1512 RESEARCH ARTICLE Decay rate of the compressible Navier-Stokes equations with density-dependent viscosity coefficients

More information

Instability of Finite Difference Schemes for Hyperbolic Conservation Laws

Instability of Finite Difference Schemes for Hyperbolic Conservation Laws Instability of Finite Difference Schemes for Hyperbolic Conservation Laws Alberto Bressan ( ), Paolo Baiti ( ) and Helge Kristian Jenssen ( ) ( ) Department of Mathematics, Penn State University, University

More information

u-- (u, v), x < 0, v,+1/2{bu:+2uv},,=o u--(u,v), (1) (1) and note that A1 A2 only when u 0 in which case A 0.

u-- (u, v), x < 0, v,+1/2{bu:+2uv},,=o u--(u,v), (1) (1) and note that A1 A2 only when u 0 in which case A 0. SIAM J. APPL. MATH. Vol. 48, No. 6, December 1988 1988 Society for Industrial and Applied Mathematics OO7 THE RIEMANN PROBLEM NEAR A HYPERBOLIC SINGULARITY III* E. ISAACSON AND B. TEMPLE Abstract. This

More information

Structurally Stable Singularities for a Nonlinear Wave Equation

Structurally Stable Singularities for a Nonlinear Wave Equation Structurally Stable Singularities for a Nonlinear Wave Equation Alberto Bressan, Tao Huang, and Fang Yu Department of Mathematics, Penn State University University Park, Pa. 1682, U.S.A. e-mails: bressan@math.psu.edu,

More information

Hyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan

Hyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global

More information

Blowup phenomena of solutions to the Euler equations for compressible fluid flow

Blowup phenomena of solutions to the Euler equations for compressible fluid flow J. Differential Equations 1 006 91 101 www.elsevier.com/locate/jde Blowup phenomena of solutions to the Euler equations for compressible fluid flow Tianhong Li a,, Dehua Wang b a Department of Mathematics,

More information

Waves in a Shock Tube

Waves in a Shock Tube Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid

More information

On the Cauchy Problems for Polymer Flooding with Gravitation

On the Cauchy Problems for Polymer Flooding with Gravitation On the Cauchy Problems for Polymer Flooding with Gravitation Wen Shen Mathematics Department, Penn State University. Email: wxs27@psu.edu November 5, 2015 Abstract We study two systems of conservation

More information

WEAK ASYMPTOTIC SOLUTION FOR A NON-STRICTLY HYPERBOLIC SYSTEM OF CONSERVATION LAWS-II

WEAK ASYMPTOTIC SOLUTION FOR A NON-STRICTLY HYPERBOLIC SYSTEM OF CONSERVATION LAWS-II Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 94, pp. 1 14. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu WEAK ASYMPTOTIC

More information

Various lecture notes for

Various lecture notes for Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes

More information

A RIEMANN PROBLEM FOR THE ISENTROPIC GAS DYNAMICS EQUATIONS

A RIEMANN PROBLEM FOR THE ISENTROPIC GAS DYNAMICS EQUATIONS A RIEMANN PROBLEM FOR THE ISENTROPIC GAS DYNAMICS EQUATIONS KATARINA JEGDIĆ, BARBARA LEE KEYFITZ, AND SUN CICA ČANIĆ We study a Riemann problem for the two-dimensional isentropic gas dynamics equations

More information

UNIFORM DECAY OF SOLUTIONS FOR COUPLED VISCOELASTIC WAVE EQUATIONS

UNIFORM DECAY OF SOLUTIONS FOR COUPLED VISCOELASTIC WAVE EQUATIONS Electronic Journal of Differential Equations, Vol. 16 16, No. 7, pp. 1 11. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu UNIFORM DECAY OF SOLUTIONS

More information

Critical Thresholds in a Relaxation Model for Traffic Flows

Critical Thresholds in a Relaxation Model for Traffic Flows Critical Thresholds in a Relaxation Model for Traffic Flows Tong Li Department of Mathematics University of Iowa Iowa City, IA 52242 tli@math.uiowa.edu and Hailiang Liu Department of Mathematics Iowa State

More information

Presenter: Noriyoshi Fukaya

Presenter: Noriyoshi Fukaya Y. Martel, F. Merle, and T.-P. Tsai, Stability and Asymptotic Stability in the Energy Space of the Sum of N Solitons for Subcritical gkdv Equations, Comm. Math. Phys. 31 (00), 347-373. Presenter: Noriyoshi

More information

On the Dependence of Euler Equations on Physical Parameters

On the Dependence of Euler Equations on Physical Parameters On the Dependence of Euler Equations on Physical Parameters Cleopatra Christoforou Department of Mathematics, University of Houston Joint Work with: Gui-Qiang Chen, Northwestern University Yongqian Zhang,

More information

Conical Shock Waves for Isentropic Euler System

Conical Shock Waves for Isentropic Euler System Conical Shock Waves for Isentropic Euler System Shuxing Chen Institute of Mathematical Research, Fudan University, Shanghai, China E-mail: sxchen@public8.sta.net.cn Dening Li Department of Mathematics,

More information

NONCLASSICAL SHOCK WAVES OF CONSERVATION LAWS: FLUX FUNCTION HAVING TWO INFLECTION POINTS

NONCLASSICAL SHOCK WAVES OF CONSERVATION LAWS: FLUX FUNCTION HAVING TWO INFLECTION POINTS Electronic Journal of Differential Equations, Vol. 2006(2006), No. 149, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) NONCLASSICAL

More information

OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES

OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES RENJUN DUAN Department of Mathematics, City University of Hong Kong 83 Tat Chee Avenue, Kowloon, Hong Kong,

More information

M. HERTY, CH. JÖRRES, AND B. PICCOLI

M. HERTY, CH. JÖRRES, AND B. PICCOLI EXISTENCE OF SOLUTION TO SUPPLY CHAIN MODELS BASED ON PARTIAL DIFFERENTIAL EQUATION WITH DISCONTINUOUS FLUX FUNCTION M. HERTY, CH. JÖRRES, AND B. PICCOLI Abstract. We consider a recently [2] proposed model

More information

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

The 2-d isentropic compressible Euler equations may have infinitely many solutions which conserve energy 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,

More information

K. Ambika and R. Radha

K. Ambika and R. Radha Indian J. Pure Appl. Math., 473: 501-521, September 2016 c Indian National Science Academy DOI: 10.1007/s13226-016-0200-9 RIEMANN PROBLEM IN NON-IDEAL GAS DYNAMICS K. Ambika and R. Radha School of Mathematics

More information

Gas Dynamics Equations: Computation

Gas Dynamics Equations: Computation Title: Name: Affil./Addr.: Gas Dynamics Equations: Computation Gui-Qiang G. Chen Mathematical Institute, University of Oxford 24 29 St Giles, Oxford, OX1 3LB, United Kingdom Homepage: http://people.maths.ox.ac.uk/chengq/

More information

Radon measure solutions for scalar. conservation laws. A. Terracina. A.Terracina La Sapienza, Università di Roma 06/09/2017

Radon measure solutions for scalar. conservation laws. A. Terracina. A.Terracina La Sapienza, Università di Roma 06/09/2017 Radon measure A.Terracina La Sapienza, Università di Roma 06/09/2017 Collaboration Michiel Bertsch Flavia Smarrazzo Alberto Tesei Introduction Consider the following Cauchy problem { ut + ϕ(u) x = 0 in

More information

On Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations

On Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations On Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations G. Seregin, V. Šverák Dedicated to Vsevolod Alexeevich Solonnikov Abstract We prove two sufficient conditions for local regularity

More information

Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3

Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso

More information

Compressible Euler equations with second sound: asymptotics of discontinuous solutions. Beixiang Fang Reinhard Racke

Compressible Euler equations with second sound: asymptotics of discontinuous solutions. Beixiang Fang Reinhard Racke Universität Konstanz Compressible Euler equations with second sound: asymptotics of discontinuous solutions Beixiang Fang Reinhard Racke Konstanzer Schriften in Mathematik Nr. 306, August 0 ISSN 430-3558

More information

Answers to Problem Set Number 04 for MIT (Spring 2008)

Answers to Problem Set Number 04 for MIT (Spring 2008) Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,

More information

Hyperbolic Problems: Theory, Numerics, Applications

Hyperbolic Problems: Theory, Numerics, Applications AIMS on Applied Mathematics Vol.8 Hyperbolic Problems: Theory, Numerics, Applications Fabio Ancona Alberto Bressan Pierangelo Marcati Andrea Marson Editors American Institute of Mathematical Sciences AIMS

More information

Rarefaction wave interaction for the unsteady transonic small disturbance equations

Rarefaction wave interaction for the unsteady transonic small disturbance equations Rarefaction wave interaction for the unsteady transonic small disturbance equations Jun Chen University of Houston Department of Mathematics 4800 Calhoun Road Houston, TX 77204, USA chenjun@math.uh.edu

More information

STRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN PROBLEM FOR A NON-STRICTLY HYPERBOLIC SYSTEM WITH FLUX APPROXIMATION

STRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN PROBLEM FOR A NON-STRICTLY HYPERBOLIC SYSTEM WITH FLUX APPROXIMATION Electronic Journal of Differential Equations, Vol. 216 (216, No. 126, pp. 1 16. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu STRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN

More information

Shock Waves in Plane Symmetric Spacetimes

Shock Waves in Plane Symmetric Spacetimes Communications in Partial Differential Equations, 33: 2020 2039, 2008 Copyright Taylor & Francis Group, LLC ISSN 0360-5302 print/1532-4133 online DOI: 10.1080/03605300802421948 Shock Waves in Plane Symmetric

More information

On finite time BV blow-up for the p-system

On finite time BV blow-up for the p-system On finite time BV blow-up for the p-system Alberto Bressan ( ), Geng Chen ( ), and Qingtian Zhang ( ) (*) Department of Mathematics, Penn State University, (**) Department of Mathematics, University of

More information

Scalar conservation laws with moving density constraints arising in traffic flow modeling

Scalar conservation laws with moving density constraints arising in traffic flow modeling Scalar conservation laws with moving density constraints arising in traffic flow modeling Maria Laura Delle Monache Email: maria-laura.delle monache@inria.fr. Joint work with Paola Goatin 14th International

More information

Notes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow

Notes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.

More information

APPROXIMATE PERIODIC SOLUTIONS for the RAPIDLY ROTATING SHALLOW-WATER and RELATED EQUATIONS

APPROXIMATE PERIODIC SOLUTIONS for the RAPIDLY ROTATING SHALLOW-WATER and RELATED EQUATIONS 1 APPROXIMATE PERIODIC SOLUTIONS for the RAPIDLY ROTATING SHALLOW-WATER and RELATED EQUATIONS BIN CHENG Department of Mathematics University of Michigan Ann Arbor, MI 48109 E-mail:bincheng@umich.edu EITAN

More information

1 Math 241A-B Homework Problem List for F2015 and W2016

1 Math 241A-B Homework Problem List for F2015 and W2016 1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let

More information

Relaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni

Relaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)

More information

Non-linear Scalar Equations

Non-linear Scalar Equations Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here

More information

BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED

BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED TAOUFIK HMIDI AND SAHBI KERAANI Abstract. In this note we prove a refined version of compactness lemma adapted to the blowup analysis

More information

On the Three-Phase-Lag Heat Equation with Spatial Dependent Lags

On the Three-Phase-Lag Heat Equation with Spatial Dependent Lags Nonlinear Analysis and Differential Equations, Vol. 5, 07, no., 53-66 HIKARI Ltd, www.m-hikari.com https://doi.org/0.988/nade.07.694 On the Three-Phase-Lag Heat Equation with Spatial Dependent Lags Yang

More information

On the minimum of certain functional related to the Schrödinger equation

On the minimum of certain functional related to the Schrödinger equation Electronic Journal of Qualitative Theory of Differential Equations 2013, No. 8, 1-21; http://www.math.u-szeged.hu/ejqtde/ On the minimum of certain functional related to the Schrödinger equation Artūras

More information

Hyperbolic Systems of Conservation Laws

Hyperbolic Systems of Conservation Laws Hyperbolic Systems of Conservation Laws III - Uniqueness and continuous dependence and viscous approximations Alberto Bressan Mathematics Department, Penn State University http://www.math.psu.edu/bressan/

More information

REGULARITY THROUGH APPROXIMATION FOR SCALAR CONSERVATION LAWS

REGULARITY THROUGH APPROXIMATION FOR SCALAR CONSERVATION LAWS SIAM J. MATH. ANAL. c 1988 Society for Industrial and Applied Mathematics Vol. 19, No. 4, pp. 1 XX, July 1988 003 REGULARITY THROUGH APPROXIMATION FOR SCALAR CONSERVATION LAWS BRADLEY J. LUCIER Abstract.

More information

HAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM

HAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM Georgian Mathematical Journal Volume 9 (2002), Number 3, 591 600 NONEXPANSIVE MAPPINGS AND ITERATIVE METHODS IN UNIFORMLY CONVEX BANACH SPACES HAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM

More information

A PARADIGM FOR TIME-PERIODIC SOUND WAVE PROPAGATION IN THE COMPRESSIBLE EULER EQUATIONS

A PARADIGM FOR TIME-PERIODIC SOUND WAVE PROPAGATION IN THE COMPRESSIBLE EULER EQUATIONS METHODS AND APPLICATIONS OF ANALYSIS. c 2009 International Press Vol. 16, No. 3, pp. 341 364, September 2009 005 A PARADIGM FOR TIME-PERIODIC SOUND WAVE PROPAGATION IN THE COMPRESSIBLE EULER EQUATIONS

More information

Hyperbolic Gradient Flow: Evolution of Graphs in R n+1

Hyperbolic Gradient Flow: Evolution of Graphs in R n+1 Hyperbolic Gradient Flow: Evolution of Graphs in R n+1 De-Xing Kong and Kefeng Liu Dedicated to Professor Yi-Bing Shen on the occasion of his 70th birthday Abstract In this paper we introduce a new geometric

More information

Hyperbolic Conservation Laws Past and Future

Hyperbolic Conservation Laws Past and Future Hyperbolic Conservation Laws Past and Future Barbara Lee Keyfitz Fields Institute and University of Houston bkeyfitz@fields.utoronto.ca Research supported by the US Department of Energy, National Science

More information

Lecture Notes on Hyperbolic Conservation Laws

Lecture Notes on Hyperbolic Conservation Laws Lecture Notes on Hyperbolic Conservation Laws Alberto Bressan Department of Mathematics, Penn State University, University Park, Pa. 16802, USA. bressan@math.psu.edu May 21, 2009 Abstract These notes provide

More information

LOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S PARABOLIC-HYPERBOLIC COMPRESSIBLE NON-ISOTHERMAL MODEL FOR LIQUID CRYSTALS

LOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S PARABOLIC-HYPERBOLIC COMPRESSIBLE NON-ISOTHERMAL MODEL FOR LIQUID CRYSTALS Electronic Journal of Differential Equations, Vol. 017 (017), No. 3, pp. 1 8. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu LOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S

More information

UNIQUENESS OF SELF-SIMILAR VERY SINGULAR SOLUTION FOR NON-NEWTONIAN POLYTROPIC FILTRATION EQUATIONS WITH GRADIENT ABSORPTION

UNIQUENESS OF SELF-SIMILAR VERY SINGULAR SOLUTION FOR NON-NEWTONIAN POLYTROPIC FILTRATION EQUATIONS WITH GRADIENT ABSORPTION Electronic Journal of Differential Equations, Vol. 2015 2015), No. 83, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu UNIQUENESS OF SELF-SIMILAR

More information

Mem. Differential Equations Math. Phys. 36 (2005), T. Kiguradze

Mem. Differential Equations Math. Phys. 36 (2005), T. Kiguradze Mem. Differential Equations Math. Phys. 36 (2005), 42 46 T. Kiguradze and EXISTENCE AND UNIQUENESS THEOREMS ON PERIODIC SOLUTIONS TO MULTIDIMENSIONAL LINEAR HYPERBOLIC EQUATIONS (Reported on June 20, 2005)

More information

Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation

Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation Advances in Dynamical Systems and Applications ISSN 973-532, Volume 6, Number 2, pp. 24 254 (2 http://campus.mst.edu/adsa Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation

More information

On a Lyapunov Functional Relating Shortening Curves and Viscous Conservation Laws

On a Lyapunov Functional Relating Shortening Curves and Viscous Conservation Laws On a Lyapunov Functional Relating Shortening Curves and Viscous Conservation Laws Stefano Bianchini and Alberto Bressan S.I.S.S.A., Via Beirut 4, Trieste 34014 Italy. E-mail addresses: bianchin@mis.mpg.de,

More information

GLOBAL EXISTENCE OF SOLUTIONS TO A HYPERBOLIC-PARABOLIC SYSTEM

GLOBAL EXISTENCE OF SOLUTIONS TO A HYPERBOLIC-PARABOLIC SYSTEM PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 35, Number 4, April 7, Pages 7 7 S -99396)8773-9 Article electronically published on September 8, 6 GLOBAL EXISTENCE OF SOLUTIONS TO A HYPERBOLIC-PARABOLIC

More information

BLOW-UP OF SOLUTIONS FOR A NONLINEAR WAVE EQUATION WITH NONNEGATIVE INITIAL ENERGY

BLOW-UP OF SOLUTIONS FOR A NONLINEAR WAVE EQUATION WITH NONNEGATIVE INITIAL ENERGY Electronic Journal of Differential Equations, Vol. 213 (213, No. 115, pp. 1 8. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu BLOW-UP OF SOLUTIONS

More information

SEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT

SEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT Electronic Journal of Differential Equations, Vol. 2012 (2012), No. 139, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SEMILINEAR ELLIPTIC

More information

c 2009 Society for Industrial and Applied Mathematics

c 2009 Society for Industrial and Applied Mathematics SIAM J. MATH. ANAL. Vol. 41, No. 1, pp. 1 25 c 2009 Society for Industrial and Applied Mathematics EVOLUTION OF DISCONTINUITY AND FORMATION OF TRIPLE-SHOCK PATTERN IN SOLUTIONS TO A TWO-DIMENSIONAL HYPERBOLIC

More information

Blow up of solutions for a 1D transport equation with nonlocal velocity and supercritical dissipation

Blow up of solutions for a 1D transport equation with nonlocal velocity and supercritical dissipation Blow up of solutions for a 1D transport equation with nonlocal velocity and supercritical dissipation Dong Li a,1 a School of Mathematics, Institute for Advanced Study, Einstein Drive, Princeton, NJ 854,

More information

NONTRIVIAL SOLUTIONS FOR SUPERQUADRATIC NONAUTONOMOUS PERIODIC SYSTEMS. Shouchuan Hu Nikolas S. Papageorgiou. 1. Introduction

NONTRIVIAL SOLUTIONS FOR SUPERQUADRATIC NONAUTONOMOUS PERIODIC SYSTEMS. Shouchuan Hu Nikolas S. Papageorgiou. 1. Introduction Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 34, 29, 327 338 NONTRIVIAL SOLUTIONS FOR SUPERQUADRATIC NONAUTONOMOUS PERIODIC SYSTEMS Shouchuan Hu Nikolas S. Papageorgiou

More information

Hyperbolic conservation laws and applications Schedule and Abstracts

Hyperbolic conservation laws and applications Schedule and Abstracts Hyperbolic conservation laws and applications Schedule and Abstracts The Graduate Center, CUNY 365 Fifth Avenue New York, NY 10016 Science Center, Room 4102 Thursday, April 26th, 2012 9:30am till 5:30pm

More information

Hysteresis rarefaction in the Riemann problem

Hysteresis rarefaction in the Riemann problem Hysteresis rarefaction in the Riemann problem Pavel Krejčí 1 Institute of Mathematics, Czech Academy of Sciences, Žitná 25, 11567 Praha 1, Czech Republic E-mail: krejci@math.cas.cz Abstract. We consider

More information

Propagation of Singularities

Propagation of Singularities Title: Name: Affil./Addr.: Propagation of Singularities Ya-Guang Wang Department of Mathematics, Shanghai Jiao Tong University Shanghai, 200240, P. R. China; e-mail: ygwang@sjtu.edu.cn Propagation of Singularities

More information

S chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.

S chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1. Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable

More information

Blowup for Hyperbolic Equations. Helge Kristian Jenssen and Carlo Sinestrari

Blowup for Hyperbolic Equations. Helge Kristian Jenssen and Carlo Sinestrari Blowup for Hyperbolic Equations Helge Kristian Jenssen and Carlo Sinestrari Abstract. We consider dierent situations of blowup in sup-norm for hyperbolic equations. For scalar conservation laws with a

More information

c 2007 Society for Industrial and Applied Mathematics

c 2007 Society for Industrial and Applied Mathematics SIAM J MATH ANAL Vol 38, No 5, pp 474 488 c 007 Society for Industrial and Applied Mathematics OPTIMAL TRACING OF VISCOUS SHOCKS IN SOLUTIONS OF VISCOUS CONSERVATION LAWS WEN SHEN AND MEE REA PARK Abstract

More information

BLOW-UP OF SOLUTIONS FOR VISCOELASTIC EQUATIONS OF KIRCHHOFF TYPE WITH ARBITRARY POSITIVE INITIAL ENERGY

BLOW-UP OF SOLUTIONS FOR VISCOELASTIC EQUATIONS OF KIRCHHOFF TYPE WITH ARBITRARY POSITIVE INITIAL ENERGY Electronic Journal of Differential Equations, Vol. 6 6, No. 33, pp. 8. ISSN: 7-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu BLOW-UP OF SOLUTIONS FOR VISCOELASTIC EQUATIONS OF KIRCHHOFF

More information

Simultaneous vs. non simultaneous blow-up

Simultaneous vs. non simultaneous blow-up Simultaneous vs. non simultaneous blow-up Juan Pablo Pinasco and Julio D. Rossi Departamento de Matemática, F.C.E y N., UBA (428) Buenos Aires, Argentina. Abstract In this paper we study the possibility

More information

On Non-degeneracy of Solutions to SU(3) Toda System

On Non-degeneracy of Solutions to SU(3) Toda System On Non-degeneracy of Solutions to SU3 Toda System Juncheng Wei Chunyi Zhao Feng Zhou March 31 010 Abstract We prove that the solution to the SU3 Toda system u + e u e v = 0 in R v e u + e v = 0 in R e

More information

Euler Equations: local existence

Euler Equations: local existence Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u

More information

LARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION

LARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION Differential Integral Equations Volume 3, Numbers 1- (1), 61 77 LARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION Robert Glassey Department of Mathematics, Indiana University

More information

Low Froude Number Limit of the Rotating Shallow Water and Euler Equations

Low Froude Number Limit of the Rotating Shallow Water and Euler Equations Low Froude Number Limit of the Rotating Shallow Water and Euler Equations Kung-Chien Wu Department of Pure Mathematics and Mathematical Statistics University of Cambridge, Wilberforce Road Cambridge, CB3

More information

z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.

z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables. Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These

More information

Global unbounded solutions of the Fujita equation in the intermediate range

Global unbounded solutions of the Fujita equation in the intermediate range Global unbounded solutions of the Fujita equation in the intermediate range Peter Poláčik School of Mathematics, University of Minnesota, Minneapolis, MN 55455, USA Eiji Yanagida Department of Mathematics,

More information

GLOBAL EXISTENCE FOR THE ONE-DIMENSIONAL SEMILINEAR TRICOMI-TYPE EQUATIONS

GLOBAL EXISTENCE FOR THE ONE-DIMENSIONAL SEMILINEAR TRICOMI-TYPE EQUATIONS GLOBAL EXISTENCE FOR THE ONE-DIMENSIONAL SEMILINEAR TRICOMI-TYPE EQUATIONS ANAHIT GALSTYAN The Tricomi equation u tt tu xx = 0 is a linear partial differential operator of mixed type. (For t > 0, the Tricomi

More information

ANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS INVOLVING A NEW FRACTIONAL DERIVATIVE WITHOUT SINGULAR KERNEL

ANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS INVOLVING A NEW FRACTIONAL DERIVATIVE WITHOUT SINGULAR KERNEL Electronic Journal of Differential Equations, Vol. 217 (217), No. 289, pp. 1 6. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS

More information

Ray equations of a weak shock in a hyperbolic system of conservation laws in multi-dimensions

Ray equations of a weak shock in a hyperbolic system of conservation laws in multi-dimensions Proc. Indian Acad. Sci. (Math. Sci.) Vol. 126, No. 2, May 2016, pp. 199 206. c Indian Academy of Sciences Ray equations of a weak in a hyperbolic system of conservation laws in multi-dimensions PHOOLAN

More information

On Asymptotic Variational Wave Equations

On Asymptotic Variational Wave Equations On Asymptotic Variational Wave Equations Alberto Bressan 1, Ping Zhang 2, and Yuxi Zheng 1 1 Department of Mathematics, Penn State University, PA 1682. E-mail: bressan@math.psu.edu; yzheng@math.psu.edu

More information

On the simplest expression of the perturbed Moore Penrose metric generalized inverse

On the simplest expression of the perturbed Moore Penrose metric generalized inverse Annals of the University of Bucharest (mathematical series) 4 (LXII) (2013), 433 446 On the simplest expression of the perturbed Moore Penrose metric generalized inverse Jianbing Cao and Yifeng Xue Communicated

More information

ON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN

ON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN ON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN XIAN LIAO Abstract. In this work we will show the global existence of the strong solutions of the inhomogeneous

More information

1. Introduction. Consider a strictly hyperbolic system of n conservation laws. 1 if i = j, 0 if i j.

1. Introduction. Consider a strictly hyperbolic system of n conservation laws. 1 if i = j, 0 if i j. SIAM J MATH ANAL Vol 36, No, pp 659 677 c 4 Society for Industrial and Applied Mathematics A SHARP DECAY ESTIMATE FOR POSITIVE NONLINEAR WAVES ALBERTO BRESSAN AND TONG YANG Abstract We consider a strictly

More information

Stability of Mach Configuration

Stability of Mach Configuration Stability of Mach Configuration Suxing CHEN Fudan University sxchen@public8.sta.net.cn We prove the stability of Mach configuration, which occurs in moving shock reflection by obstacle or shock interaction

More information

Math 660-Lecture 23: Gudonov s method and some theories for FVM schemes

Math 660-Lecture 23: Gudonov s method and some theories for FVM schemes Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,

More information

Non-radial solutions to a bi-harmonic equation with negative exponent

Non-radial solutions to a bi-harmonic equation with negative exponent Non-radial solutions to a bi-harmonic equation with negative exponent Ali Hyder Department of Mathematics, University of British Columbia, Vancouver BC V6TZ2, Canada ali.hyder@math.ubc.ca Juncheng Wei

More information

A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.

A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1. A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion

More information