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 of Mathematics, Chinese Academy of Sciences, Beijing, 100080 Yuxi Zheng Department of Mathematics, The Pennsylvania State University, PA 1680 Abstract We present a characteristic decomposition of the potential flow equation in the self-similar plane. The decomposition allows for a proof that any wave adjacent to a constant state is a simple wave for the adiabatic Euler system. This result is a generalization of the well-known result on -d steady potential flow and a recent similar result on the pressure gradient system. Keywords: Simple waves, characteristic decomposition, Riemann invariants, pressure gradient equation, potential flow, shock waves, -D Riemann problem, gas dynamics. AMS subject classification: Primary: 35L65, 35J70, 35R35; Secondary: 35J65. 1 Introduction The one-dimensional wave equation with constant speed c has an interesting decomposition u tt c u xx = 0 (1) ( t + c x )( t c x )u = 0, () 1 Research partially supported by NSF of China with No. 103010, NSF from Beijing Municipality, Fok Ying Tong Educational Foundation, and the Key Program from Beijing Educational Commission with no. KZ0051008018. Research partially supported by NSF-DMS-0305497, 0305114. 1
or ( t c x )( t + c x )u = 0 (3) known from elementary text books. One can rewrite them as + u = 0, or + u = 0, (4) where ± = t ± c x. Sometimes, the same fact is written in Riemann invariants t R + c x R = 0, t S c x S = 0 (5) for the Riemann invariants R := ( t c x )u, S := ( t + c x )u. (6) For a pair of system of hyperbolic conservation laws u v ] t + f(u, v) g(u, v) ] x = 0 0 ], (7) it is known that a pair of Riemann invariants exist so that the system can be rewritten as { t R + λ 1 (u, v) x R = 0, (8) t S + λ (u, v) x S = 0, where (R, S) are the Riemann invariants and the λ s are the two eigenvalues of the system. These decompositions and Riemann invariants are useful in the construction of solutions, for example, the construction of D Alembert formula, and proof of development of singularities (5]). An example of the system is the system of isentropic irrotational steady two-dimensional Euler equations for compressible ideal gases { (c u )u x uv(u y + v x ) + (c v )v y = 0, u y v x = 0 (9) supplemented by Bernoulli s law c γ 1 + u + v = k, (10) where γ > 1 is the gas constant while k > 0 is an integration constant. This system has two unknowns (u, v), and by the existence of Riemann invariants, any solution adjacent to a constant state is a simple wave. A simple wave means a solution (u, v) that depends on a single parameter rather than the pair parameters (x, y). Since there is the lack of
explicit expressions, the concept of Riemann invariants plays a limited role in a much broader sense, e. g., to treat the full Euler equations. In recent years, the pressure gradient system u t + p x = 0 v t + p y = 0 E t + (up) x + (vp) y = 0, (11) where E = p + (u + v )/, has been known to have simple waves adjacent to a constant state (u, v, p) in the self-similar variable plane (ξ, η) = (x/t, y/t). This system has three equations and no Riemann invariants have been found. But the equation for the unknown variable p in the (ξ, η) plane (p ξ )p ξξ ξηp ξη + (p η )p ηη + (ξp ξ + ηp η ) allows for a decomposition p (ξp ξ + ηp η ) = 0 (1) + p = m + p, m + := r4 λ + p p r, (13) where (r, θ) denotes the polar coordinates of the (ξ, η) plane and + = θ + λ 1 + r ; = θ + λ 1 p r ; λ ± = ± r (r p). (14) For convenience of verification we state that the p equation in polar coordinates takes the form (p r )p rr + p r p θθ + p r p r + 1 p (rp r) rp r = 0. (15) The characteristics are defined by In addition, we know that dθ dr = λ ±. (16) for some nice factors n ±. These facts allow for expressions ± λ ± = n ± ± p (17) ( ± λ ± ) = ( ± λ ± )f ± (18) for some nice factors f ±. This decomposition leads directly to the fact that Proposition 1. A state adjacent to a constant state for the pressure gradient system must be a simple wave in which p is constant along characteristics of a plus (or minus) family. 3
These lead to the desire to consider the pseudo-steady isentropic irrotational Euler system which has three equations with source terms, (ρu) ξ + (ρv ) η = ρ, (ρu + p(ρ)) ξ + (ρuv ) η = 3ρU, (ρuv ) ξ + (ρv + p(ρ)) η = 3ρV, where (ξ, η) = (x/t, y/t), (U, V ) = (u ξ, v η) is the pseudo-velocity, and the pressure p = p(ρ) is the function of the density ρ. It turns out that we are unable to find explicit forms of the Riemann invariants, but decompositions similar to + λ = m λ hold for some m, presented in Section 4. We use the characteristic decomposition of Section 4 to establish in Section 5 that adjacent to a constant state a wave must be a simple wave for the pseudo-steady irrotational isentropic Euler system. A simple wave for this case is such that one family of wave characteristics are straight lines and the physical quantities velocity, speed of sound, pressure, and density are constant along the wave characteristics. Further, using the fact that entropy and vorticity are constant along the pseudo-flow characteristics (the pseudo-flow lines), our irrotational result extends to the adiabatic full Euler system, see Section 5. (19) Two-by-two system Consider a hyperbolic system in the Riemann invariants { t R + λ 1 x R = 0, t S + λ x S = 0. (0) So we find that where λ,r := R λ. We go on to find λ := ( t + λ x )λ = λ,r R, (1) and so 1 R = 1λ λ 1 λ λ 1 R, () 1 ( 1 λ,r λ ) = 1λ λ 1 λ λ 1 λ λ,r, (3) which implies 1 λ = ( 1 λ λ 1 + ) 1λ,R λ. (4) λ λ 1 λ,r 4
The elegant form is undermined by the dependence on the Riemann invariant R via the term λ,r. It is not very useful when the explicit form of the Riemann invariants are not known. But we think it is worth mentioning. For example, (4) can be applied directly to show that all characteristics in a wave adjacent to a constant state are straight and thus such a wave is a simple wave. 3 Steady Euler system Let us build explicitly the characteristic decomposition for the steady Euler system for isentropic irrotational flow (9)(10) in the absence of the explicit form of the Riemann invariants. The same technique can be extended to the case of pseudo-steady flows in Section 4. We write the system in the form u v The matrix has eigenvalues ] x + uv c u c v c u 1 0 ] u v λ ± = uv ± c (u + v c ) u c which are solutions to the characteristic equation We have the left eigenvectors where we have used the relation The characteristic form of the system is therefore ] ] which is equivalent to We then have ] y = 0. (5) (= dy ), (6) dx λ + uv c u λ + c v = 0. (7) c u l ± u v l ± = 1, λ ], (8) λ ± λ = c v c u. (9) x + λ ± l ± u v y = 0, (30) ± u + λ ± v = 0. (31) λ = x λ + λ y λ = u λ u + v λ v = ( u λ v λ /λ + ) u. (3) 5
We shall ignore the similar calculation for + λ + for simplicity of notation. Now that the term λ differs from u by a lower-order factor, we shall focus our attention on u. First we see that we can derive a second-order equation for u, i. e., or equivalently u xx u yy ( uv c v u y c u ) c v u x = 0, (33) x uv c u u xy + c v c u u yy = c v ( ( ) uv u c u c v y x We now compute the ordered derivative + u to find We find that ( c u ) c v ) u x x. (34) + u = u xx + (λ + + λ )u xy + λ + λ u yy + + λ u y. (35) + λ = + u( u λ v λ /λ ). (36) Thus we obtain (( ) c + u = v uv u c u c v y ( c u x c v )x u ) x +(u x + λ + u y )u y ( u λ v λ /λ ). (37) We notice that the above right-hand side is a quadratic form in (u x, u y ), once we substitute v x by u y. So we compute further. We use the Bernoulli s law to find (c ) x = (γ 1)(uu x + vu y ). (38) So we find ( uv c v )x ( c u c v )x = (c v ) vu x (c u v + γu ) + uu y (c + γv )], = 1 (c v ) uu x (c v u + γu γv ) vu y (c v + γv u γu )]. (39) We now compute the factor u λ λ 1 v λ. We use Bernoulli s law to find We use the characteristic equation (c ) u = (γ 1)u, (c ) v = (γ 1)v. (40) (c u )λ + uvλ + c v = 0 (41) 6
to obtain u λ = λ (γ+1)u vλ +(γ 1)u λ (c u )+uv, v λ = λ (γ 1)v uλ +(γ+1)v λ (c u )+uv. (4) We then simply compute to find u λ λ 1 v λ = γ + 1 λ (c u ) + uv Coming back to our equation for + u, we have (c u )(c v ) + u (uλ v) 3 c λ. (43) = u xu(c u v + γu γv ) +u x u y ( vu v 3 + 3γvu γv 3 + Q) (44) +u yu(c + γv ) + λ + Q], where we have introduced the notation where Q := (c u )(c v ) λ (c u ) + uv We then factorize the quadratic form to find finally where γ + 1 (uλ c v) 3 = (γ + 1)(c v )(uh vc) 3, (45) λ (c u )H(cH uv) H := u + v c. (46) + u = u(c u v + γu γv ) ( (c u )(c v x u + α y u) u, (47) ) α = u(c + γv ) + λ + Q λ u(c u v + γu γv ). (48) Proposition There holds the identity for α(u, v) given in (48) and m given by + u = m( x u + α y u) u (49) m = u(c u v + γu γv ). (50) (c u )(c v ) 7
We use the relation to go back to + λ. We find u = λ /( u λ v λ λ 1 + ) (51) u λ v λ λ 1 + = 4c + (γ 3)(u + v )]vλ (c + u ) + (c v )u] (c v )(c u )λ (c u ) + uv] So we have =: G. (5) or + ( 1 G λ ) = m α u λ G, (53) + λ = (m α u + + (ln G )) λ. (54) Proposition 3 There holds the identity + λ = m λ (55) for some m = m(u, v)( x u + β(u, v) y u). A similar identity holds for + λ +. We remark that in the application on simple waves, the equation for u is sufficient and the equations for λ ± are not needed. 4 Pseudo-steady Euler We consider the two-dimensional isentropic irrotational ideal flow in the self-similar plane (ξ, η) = (x/t, y/t). There holds the Bernoulli s law c γ 1 + U + V = ϕ, (56) where c is the speed of sound, (U, V ) = (u ξ, v η) are the pseudo-velocity, while (u, v) is the physical velocity, and ϕ is the pseudo-potential such that The equations of motion can be written as ϕ ξ = U, ϕ η = V. (57) { (c U )U ξ UV (U η + V ξ ) + (c V )V η = c + U + V, V ξ U η = 0. (58) 8
We can rewrite the equations of motion in a new form ] ] ] u UV c V + c U c U u v 1 0 v ξ η = 0 (59) to draw as much parallelism to the steady case as possible. We emphasize the mixed use of the variables (U, V ) and (u, v), i. e., (U, V ) is used in the coefficients while (u, v) is used in differentiation. This way we obtain zero on the right-hand side for the system. The eigenvalues are similar as before: dη dξ = Λ ± = UV ± c (U + V c ). (60) U c The left eigenvectors are And we have similarly l ± = 1, Λ ]. (61) ± u + Λ ± v = 0. (6) Our Λ ± now depend on more than (U, V ). But, let us regard Λ ± as a simple function of three variables Λ ± = Λ ± (U, V, c ) as given in (60). Thus we need to build differentiation laws for c. We can directly obtain ( ) ( ) c c + Uu ξ + V v ξ = 0, + Uu η + V v η = 0. (63) γ 1 γ 1 We have So we move on to compute ξ η ± c = (γ 1)(U ± u + V ± v). (64) ± Λ ± = U Λ ± ± U + V Λ ± ± V + c Λ ± ± c = U Λ ± ( ± u 1) + V Λ ± ( ± v Λ ± ) + c Λ ± ± c = U Λ ± ± u + V Λ ± ± v + c Λ ± ± c U Λ ± V Λ ± Λ ±. (65) We need to handle the term U Λ ± + V Λ ± Λ ±. We show it is zero. Recalling that (c U )Λ + UV Λ + c V = 0, (66) and regarding that Λ depends on the three quantities (U, V, c ) independently, we can easily find Λ U = Λ(UΛ V ) Λ(c U ) + UV, Λ UΛ V V = Λ(c U ) + UV. (67) 9
Thus Λ U + ΛΛ V = 0. (68) Therefore we end up with ± Λ ± = U Λ ± Λ 1 V Λ ± (γ 1) c Λ ± (U Λ 1 V ) ] ± u. (69) Thus, if one of the quantities (u, v, c, Λ ) is a constant along Λ, so are all the rest. So far the properties are very similar to the steady case. We derive an equation for u. We have a similar second-order equation for u ( UV u ηη = c V u η c U ) c V u ξ. (70) We have similarly ξ We compute + u = u ξξ + (Λ + + Λ )u ξη + Λ Λ + u ηη + + Λ u η = ( ) c V UV c U c V η ( ] c U ξ c V )ξ u ξ + + Λ u η. + Λ = U Λ + U + V Λ + V + c Λ + c = U Λ 1 Λ V Λ (γ 1) c Λ (U 1 Λ V )] + u ( U Λ + Λ + V Λ ). We continue to find { + u = c V u η V c U (c V ) ξ U(c V ) + V U ξ (c V ) V U((c ) ξ V V ξ )] u ξ ((c ) (c V ) ξ UU ξ )(c V ) (c U )((c ) ξ V V ξ )] } + + u u η U Λ 1 Λ V Λ (γ 1) c Λ (U 1 Λ V ) ] u η ( U Λ + Λ + V Λ ). We apply the rule U ξ = u ξ 1, V ξ = v ξ = u η to find { + u = c V u η u c U (c V ) η U(c V ) + V u ξ (c V ) +V U((γ 1)Uu ξ + (γ + 1)V u η )] u ξ ((γ + 1)Uu (c V ) ξ + (γ 1)V u η )(c V ) +(c U )((γ 1)Uu ξ + (γ + 1)V u η )]} + + u u η U Λ 1 Λ V Λ (γ 1) c Λ (U 1 Λ V ) ] V v ξ Uu ξ c U u η ( U Λ + Λ + V Λ ). (71) (7) (73) (74) 10
We note that there appear terms which are linear in the derivatives of (u, v) in addition to the pure quadratic form as in the steady case. The pure quadratic form is identical to the steady case, so we do not need to handle it further. The linear form can be handled as follows. First we use the derivatives (Λ U, Λ V ) to compute Then we have U Λ + Λ + V Λ = U Λ + 1 Λ c V c U V Λ = (UΛ V ) c U. (75) V v ξ Uu ξ c U u η ( U Λ + Λ + V Λ ) = U c U u. (76) Thus the linear form is also in the direction of Λ. Combining the steps we end up with Theorem 4. There holds where and + u = U(c U V + γu γv ) ( (c U )(c V ξ u + A η u) u U ) c U u, (77) A := U(c + γv ) + Λ + Q Λ U(c U V + γu γv ), (78) Q := (c U )(c V ) Λ (c U ) + UV γ + 1 c Λ (UΛ V ) 3. (79) We then have Theorem 5. There holds + ( Λ ) = m Λ. (80) Similarly there holds ( + Λ + ) = n + Λ +. 5 Application: Simple waves For a system of hyperbolic conservation laws in one-space dimension, a centered rarefaction wave is a simple wave, in which one family of characteristics are straight lines and the dependent variables are constant along a characteristic. See any text book on systems of conservation laws, e. g., Courant and Friedrichs ], p.59, and others 9]3]. Simple waves for the two-dimensional steady Euler system are similar, i. e., one family of characteristics are straight and the velocity are constant along the characteristics. 11
For the two-dimensional self-similar pressure gradient system, see 4], simple waves can be defined similarly, i. e., one nonlinear family of characteristics are straight and the pressure is constant along them. We note that we do not require the velocity to be constant. This way, by the characteristic decomposition, we find that a wave adjacent to a constant state is a simple wave. In the construction of solutions to the two-dimensional Riemann problem for the Euler system, see any of the sources 7]9]8], it is important to know how to construct solutions adjacent to a constant state in addition to the constructions of the interaction of rarefaction waves (6]), subsonic solutions, and transonic shock waves. Simple waves play a fundamental role in describing and building up solutions of flow problems. ( See pp.59 60, ]). In 1], Canić and Keyfitz explored simple waves further and generalized Courant and Friedrichs theorem by allowing the coefficients of a system (say (7)) to depend on the independent variables (t, x) linearly as well as the dependent variables (u, v). Here the characteristic decomposition 1 λ = m λ allows us to conclude that Theorem 6. Adjacent to a constant state in the self-similar plane of the potential flow system is a simple wave in which the physical variables (u, v, c) are constant along a family of characteristics which are straight lines. 5.1 Simple waves for full Euler Consider the adiabatic Euler system for an ideal fluid ρ t + (ρu) = 0, (ρu) t + (ρu u + pi) = 0, (ρe) t + (ρeu + pu) = 0, (81) where E := 1 u + e, where e is the internal energy. For a polytropic gas, there holds e = 1 p γ 1 ρ, where γ > 1. In the self-similar plane and for smooth solutions, the system takes the form 1 ρ sρ + u ξ + v η = 0, s u + 1p ρ ξ = 0, s v + 1p (8) ρ η = 0, 1 γp sp + u ξ + v η = 0, 1
where s := (u ξ) ξ + (v η) η, which we call pseudo-flow directions, as opposed to the other two characteristic directions, called (pseudo-)wave characteristics. We easily derive and s (pρ γ ) = 0, (83) ω t + (uω) x + (vω) y + ( p y ρ ) x ( p x ρ ) y = 0 (84) for the vorticity ω := v x u y. So entropy pρ γ is constant along the pseudo-flow lines. For a region Ω whose pseudo-flow lines come from a constant state, we see that the entropy is constant in the region. For isentropic region, vorticity has zero source of production since ( py ) ρ x ( px ) ρ y = 0. Thus vorticity tω = v ξ u η (setting t = 1 then) satisfies s (ω/ρ) = 0. (85) Hence, for a region whose pseudo-flow lines come from a constant state, the vorticity must be zero everywhere. So the region is irrotational and isentropic. Thus our formulas for the potential flow apply. We have Theorem 7. Adjacent to a constant state in the self-similar plane of the adiabatic Euler system is a simple wave in which the physical variables (u, v, c, p, ρ) are constant along a family of wave characteristics which are straight lines, provided that the region is such that its pseudo-flow characteristics extend into the state of constant. Acknowledgment Yuxi Zheng would like to thank the mathematics department at Capital Normal University for its hospitality during his visit when this work was done. Jiequan Li thanks Matania Ben-Artzi for his interest. References 1] Canić, S., Keyfitz, B. L.: Quasi-one-dimensional Riemann problems and their role in self-similar two-dimensional problems, Arch. Rational. Mech. Anal., 144 (1998), 33 58. ] Courant, R., and Friedrichs, K. O.: Supersonic flow and shock waves, Interscience, New York, 1948. 3] Dafermos, C.: Hyperbolic conservation laws in continuum physics (Grundlehren der mathematischen Wissenschaften), pp.443, Springer, 000. 4] Dai, Zihuan; and Zhang, Tong: Existence of a global smooth solution for a degenerate Goursat problem of gas dynamics, Arch. Rational Mech. Anal., 155 (000), 77 98. 13
5] Lax, P.: Development of singularities of solutions of nonlinear hyperbolic partial differential equations, J. Mathematical Phys., 5 (1964), 611 613. 6] Li, Jiequan: On the two-dimensional gas expansion for compressible Euler equations, SIAM J. Appl. Math., 6 (001), 831 85. 7] Li, Jiequan; Zhang, Tong; and Yang, Shuli: The two-dimensional Riemann problem in gas dynamics, Pitman monographs and surveys in pure and applied mathematics 98. Addison Wesley Longman limited, 1998. 8] Zhang, Tong; and Zheng, Yuxi: Conjecture on the structure of solution of the Riemann problem for two-dimensional gas dynamics systems, SIAM J. Math. Anal., 1 (1990), 593 630. 9] Zheng, Yuxi: Systems of conservation laws: Two-dimensional Riemann problems, 38 PNLDE, Birkhäuser, Boston, 001. 14