Advances in Mathematical Physics
Volume 2016 (2016), Article ID 1783689, 12 pages
http://dx.doi.org/10.1155/2016/1783689

Research Article

Interactions of Delta Shock Waves for Zero-Pressure Gas Dynamics with Energy Conservation Law

Wei Cai¹ and Yanyan Zhang²

¹The Basic Department, The First Aeronautic Institute of the Air Force, Xinyang 464000, China
²College of Mathematics and Information Science, Xinyang Normal University, Xinyang 464000, China

Received 1 April 2016; Accepted 15 May 2016

Academic Editor: Ming Mei

Copyright © 2016 Wei Cai and Yanyan Zhang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We study the interactions of delta shock waves and vacuum states for the system of conservation laws of mass, momentum, and energy in zero-pressure gas dynamics. The Riemann problems with initial data of three piecewise constant states are solved case by case, and four different configurations of Riemann solutions are constructed. Furthermore, the numerical simulations completely coinciding with theoretical analysis are shown.

1. Introduction

As is well known, the system of zero-pressure gas dynamics consisting of conservation laws of mass and momentum, which is also called the transport equations, or Euler equations for pressureless fluids, has been extensively investigated since the 90s of 20th century. It is derived from Boltzmann equations [1] and the flux-splitting scheme of the full compressible Euler equations [2, 3] and can be used to describe the motion process of free particles sticking together under collision [4] and the formation of large-scale structures in the universe [5, 6].

However, we have to mention that, as having no pressure, the energy transport must be taken into account for the considered media. Therefore, it is very necessary to consider the conservation law of energy in zero-pressure gas dynamics. To this end, we study the one-dimensional zero-pressure gas dynamics governed by the conservation laws of mass, momentum, and energy:

where and represent the density and velocity, respectively, is the internal energy and assumed to be nonnegative, and is the internal energy per unit mass. The regions in the physical space where and are identified with the vacuum regions of the flow. Here, is considered as an independent variable just for convenience.

System (1) was early studied by Kraiko [7]. In contrast to the traditional zero-pressure gas dynamics system that contains only the conservation laws of mass and momentum, in order to construct the solution of (1) for arbitrary initial data, a new type of discontinuities which are different from the classical ones and carry mass, impulse, and energy are needed. In [8, 9], system (1) was further discussed. Some special integral identities were introduced to define the delta-shock solutions and construct the Rankine-Hugoniot relation for delta shock waves. Moreover, using these integral identities, the balance laws describing mass, momentum, and energy transport from the area outside the delta shock wave front onto its front were derived. What is more, the delta shock wave type solutions for multidimensional zero-pressure gas dynamics with the energy conservation law were defined in [10].

A delta shock wave is a generalization of an ordinary shock wave. Roughly speaking, it is a kind of discontinuity, on which at least one of the state variables may develop an extreme concentration in the form of a weighted Dirac delta function with the discontinuity as its support. It is more compressive than an ordinary shock wave and more characteristics enter the discontinuity line. Physically, the delta shock waves describe the process of formation of the galaxies in the universe and the process of concentration of particles. As for delta shock waves, there are numerous excellent papers, see [11–21] and so forth. Nevertheless, compared to these results, a distinctive feature for (1) is that the Dirac delta functions develop in both state variables and simultaneously, which is quite different from those aforementioned, in which only one state variable contains the Dirac delta function. In fact, the theory of delta shock waves with Dirac delta functions developing in both state variables has been established by Yang and Zhang [22, 23] for a class of nonstrictly hyperbolic systems of conservation laws.

In the past over two decades, the investigation of interactions of delta shock waves has been increasingly active. This is important not only because of their significance in practical applications but also because of their basic role as building blocks for the general mathematical theory of quasi-linear hyperbolic equations. And the results on interactions are also touchstones for the numerical schemes. Specifically, Sheng and Zhang [18] discussed the overtaking of delta shock waves and vacuum states in one-dimensional zero-pressure gas dynamics. By solving the two-dimensional Riemann problems for zero-pressure gas dynamics with three constant states, Cheng et al. [24] studied the interactions among delta-shock waves, vacuums, and contact discontinuities. In addition, with the help of a generalized plane wave solution, Yang [25] studied a type of generalized plane delta-shock wave for the -dimensional zero-pressure gas dynamics and investigated the overtaking of two plane delta shocks. For more works on the interactions of delta shock waves, we refer to [26–29] and so forth.

Motivated by the discussions above, in the present paper, we are concerned with the interactions among delta shock waves, vacuum states, and contact discontinuities in solutions. Therefore, we study the Riemann problem of (1) with initial data of three piecewise constant states as follows:

where , , are arbitrary constants and , are any two fixed points on -axis.

We will deal with the Riemann problem (1), (2) case by case along with constructing the solutions. For this purpose, it is necessary to consider whether two adjacent waves intersect and interact with each other when constructing the global solution. However, it is often not so easy to see whether two delta shock waves meet and how they interact with each other. Therefore, some technical treatments are needed.

This paper is arranged as follows. In Section 2, the delta shock solution of (1) is reviewed and a general case when the delta shock wave is emitted at the beginning with a nonzero initial data is considered. Section 3 discusses the interactions of the delta shock waves and vacuum states. The Riemann solutions of (1), (2) are constructed globally. Finally, four kinds of numerical simulations coinciding with the theoretical analysis are presented in Section 4.

2. Delta-Shock Solution

This section briefly reviews the delta shock solution of (1) with the initial data

and the detailed study can be found in [30].

System (1) has a triple eigenvalue and two right eigenvectors and . Since for , is linearly degenerate, which means that the elementary waves involve only contact discontinuities. The self-similar solution is constructed by two cases.

For the case , the solution containing two contact discontinuities and a vacuum state besides two constants is expressed as

where is a smooth function satisfying and .

For the case , the singularity of solutions must develop because of the overlap of characteristic lines. Therefore, the solution involving a delta shock wave is introduced.

Let be the delta shock solution of the form

and then the following generalized Rankine-Hugoniot relation holds

where . In order to ensure the uniqueness, the delta shock wave should satisfy the entropy condition

which means that the characteristics on both sides of the discontinuity are in-coming.

Under the entropy condition (7), by solving the ordinary differential equations (6) with the initial data , one has

For convenience, we now consider a special case when a delta shock wave is emitted at the beginning with the initial data

satisfying It yields from (6) and (9) that

One can check that the delta shock solution (10) satisfies the following:

 (1) is a monotone function of . (2)If , then

 While if , then  (3)

3. Interactions of Delta Shock Waves

In this section, we analyze the interactions of delta shock waves. To ensure that all the cases are covered completely, according to the relation among , our discussion is divided into four cases:

(1);(2);(3);(4).

Case 1 (). In this case, two delta shock waves and will be emitted from and , respectively, as shown in Figure 1.

Figure 1: .

According to what has been discussed in Section 2, these two delta shock waves are uniquely determined by

We have by entropy condition (6), which means that will overtake at a finite time. The intersection point is calculated by

which yields that

At the intersection , the new initial data are formed as follows:

satisfying .

In view of , a new delta shock wave will generate after interaction and we denote it with . The trajectory, velocity, and weights of can be uniquely obtained by solving the ordinary differential equations (6) with the initial date (16). The detail is omitted.

Thus, the result of interaction of two delta shock waves is still a single delta shock wave. This fact can be formulated as

Case 2 ( (when , the structure of solution is similar)). In this situation, a delta shock wave determined by (12) is emitted from and two contact discontinuities and with a vacuum in between are emitted from , as shown in Figure 2.
Since the propagating speed of satisfies , so must meet the contact discontinuity at , and a new delta shock wave forms, which is subjected to the generalized Rankine-Hugoniot relation

with the initial data

Figure 2: .

Therefore, by solving (18) and (19), we have

It is clear that will cross the vacuum region with a varying propagation speed. Noting that , so will penetrate over the whole vacuum region and then meet at a finite time. The intersection is determined by

At , a new initial value problem is formed and can be solved similar to Case  1. We denote the delta shock wave connecting two constant states and with after the interaction of and .

The conclusion of this case is that the delta shock wave will penetrate over the whole vacuum region between two contact discontinuities. This fact is expressed as

Case 3 ((when , the structure of solution is similar)). Similar to Case  2, there are a delta shock wave, two contact discontinuities, and a vacuum near on the -plane, as shown in Figure 3.
The delta shock wave collides with at first and a new delta shock wave generates. However, since , cannot penetrate over the vacuum region and finally has as its asymptote. This fact is symbolized as

Figure 3: .

Case 4 (). In this situation, both the contact discontinuities with a vacuum state in between are emitted from and , respectively. Noting that the contact discontinuities and own the same propagating speed, thus there is no collision of waves and the solution is expressed as

which is called a collisionless solution, as shown in Figure 4.

Figure 4: .

4. Numerical Simulations

In order to verify the validity of the interactions of delta shock waves and vacuum states mentioned in Section 3, we present some representative numerical simulations in this section. Many more numerical tests have been performed to make sure that what are presented are not numerical artifacts.

To discretize the system, we employ the second-order nonoscillatory central schemes [31] with cells and . In what follows, by taking and , we simulate the interaction of waves by four cases. For convenience, each situation will be simulated at two different times.

Case 1 (). We take the initial data as follows:

The numerical results are presented by Figures 5–7.

(a)

(b)

Figure 5: Numerical results of at (a) and (b) for Data 1.

(a)

(b)

Figure 6: Numerical results of at (a) and (b) for Data 1.

(a)

(b)

Figure 7: Numerical results of at (a) and (b) for Data 1.

We observe from Figures 5–7 that when , two delta shock waves appear at and , respectively. As increases, they will overtake each other and finally unify into a new delta shock wave at .

Case 2 (). We choose the following initial data

The numerical results are shown in Figures 8–10.

(a)

(b)

Figure 8: Numerical results of at (a) and (b) for Data 2.

(a)

(b)

Figure 9: Numerical results of at (a) and (b) for Data 2.

(a)

(b)

Figure 10: Numerical results of at (a) and (b) for Data 2.

From Figures 8–10, we can clearly see that, at , a delta shock wave and two contact discontinuities with a vacuum state in between are emitted from and , respectively. However, at , the delta shock wave penetrates over the whole vacuum region, and a new delta shock wave generates.

Case 3 (). The initial data are

The numerical results are shown in Figures 11–13.

(a)

(b)

Figure 11: Numerical results of at (a) and (b) for Data 3.

(a)

(b)

Figure 12: Numerical results of at (a) and (b) for Data 3.

(a)

(b)

Figure 13: Numerical results of at (a) and (b) for Data 3.

Figures 11–13 imply that a delta shock wave is emitted from , and two contact discontinuities with a vacuum in between are emitted from at . But the delta shock wave can not penetrate over the whole vacuum region even though time is on the increase. In this process, the region of vacuum state keeps expanding.

Case 4 (). We select the initial data to be

The numerical results are presented by Figures 14–16.

(a)

(b)

Figure 14: Numerical results of at (a) and (b) for Data 4.

(a)

(b)

Figure 15: Numerical results of at (a) and (b) for Data 4.

(a)

(b)

Figure 16: Numerical results of at (a) and (b) for Data 4.

From Figures 14–16, we observe that both the contact discontinuities with a vacuum state in between are emitted from and at , respectively. As time goes on, the vacuum state keeps continuously expanding and never disappears.

To sum up, all of the above numerical results clearly reveal the interactions of delta shock waves and vacuum states discussed in Section 3. We also indicate that because of the occurrence of singularity as the weighted Dirac delta functions, some oscillations appear in the numerical experiments as shown in Figures 8 and 11. It may be a challenge for numerical schemes when delta shock waves develop in solutions.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (No. 11501488), the Scientific Research Foundation of Xinyang Normal University (No. 0201318), and Nan Hu Young Scholar Supporting Program of XYNU.

References

F. Bouchut, “On zero pressure gas dynamics,” in Advances in Kinetic Theory and Computing, vol. 22, pp. 171–190, World Scientific, Singapore, 1994.

R. K. Agarwal and D. W. Halt, “A modified CUSP scheme in wave/particle split form for unstructured grid Euler flows,” in Frontiers of Computational Fluid Dynamics, pp. 155–163, World Scientific, Singapore, 1994.

Y. Li and Y. Cao, “Second order large particle difference method,” Science China Mathematics A, vol. 22, pp. 171–190, 1985.

Y. Brenier and E. Grenier, “Sticky particles and scalar conservation laws,” SIAM Journal on Numerical Analysis, vol. 35, no. 6, pp. 2317–2328, 1998.

S. F. Shandarin and Ya. B. Zeldovich, “The large-scale structure of the universe: turbulence, intermittency, structures in a self-gravitating medium,” Reviews of Modern Physics, vol. 61, no. 2, pp. 185–220, 1989.

E. Weinan, Y. G. Rykov, and Y. G. Sinai, “Generalized variational principles, global weak solutions and behavior with random initial data for systems of conservation laws arising in adhesion particle dynamics,” Communications in Mathematical Physics, vol. 177, no. 2, pp. 349–380, 1996.

A. Kraiko, “Discontinuity surfaces in medium without self-pressure,” Prikladnaia Matematika i Mekhanika, vol. 43, pp. 539–549, 1979.

B. Nilsson, O. S. Rozanova, and V. M. Shelkovich, “Mass, momentum and energy conservation laws in zero-pressure gas dynamics and δ-shocks: II,” Applicable Analysis, vol. 90, no. 5, pp. 831–842, 2011.

B. Nilsson and V. M. Shelkovich, “Mass, momentum and energy conservation laws in zero-pressure gas dynamics and delta-shocks,” Applicable Analysis, vol. 90, no. 11, pp. 1677–1689, 2011.

S. Albeverio, O. S. Rozanova, and V. M. Shelkovich, “Transport and concentration processes in the multidimensional zero-pressure gas dynamics model with the energy conservation law,” http://arxiv.org/abs/1101.5815

D. J. Korchinski, Solution of a Riemann problem for a 2×2 system of conservation laws possessing no classical weak solution [Ph.D. thesis], Adelphi University, Garden City, NY, USA, 1977.

B. L. Keyfitz and H. C. Kranzer, “A viscosity approximation to a system of conservation laws with no classical Riemann solution,” in Nonlinear Hyperbolic Problems, vol. 1402, pp. 185–197, Springer, Berlin, Germany, 1990.

H. C. Kranzer and B. L. Keyfitz, “A strictly hyperbolic system of conservation laws admitting singular shocks,” in Nonlinear Evolution Equations That Change Type, vol. 27 of The IMA Volumes in Mathematics and Its Applications, pp. 107–125, Springer, New York, NY, USA, 1990.

B. L. Keyfitz and H. C. Kranzer, “Spaces of weighted measures for conservation laws with singular shock solutions,” Journal of Differential Equations, vol. 118, no. 2, pp. 420–451, 1995.

P. Le Floch, “An existence and uniqueness result for two nonstrictly hyperbolic systems,” in Nonlinear Evolution Equations that Change Type, vol. 27, pp. 126–138, Springer, New York, NY, USA, 1990.

F. Huang and Z. Wang, “Well posedness for pressureless flow,” Communications in Mathematical Physics, vol. 222, no. 1, pp. 117–146, 2001.

V. Shelkovich, “The Riemann problem admitting -shocks, and vacuum states (the vanishing viscosity approach),” Journal of Differential Equations, vol. 231, no. 2, pp. 459–500, 2006.

W. Sheng and T. Zhang, “The Riemann problem for the transportation equations in gas dynamics,” Memoirs of the American Mathematical Society, vol. 137, no. 654, pp. 1–77, 1999.

D. C. Tan and T. Zhang, “Two-dimensional Riemann problem for a hyperbolic system of nonlinear conservation laws. I. Four- cases,” Journal of Differential Equations, vol. 111, no. 2, pp. 203–254, 1994.

D. C. Tan, T. Zhang, and Y. X. Zheng, “Delta-shock waves as limits of vanishing viscosity for hyperbolic systems of conservation laws,” Journal of Differential Equations, vol. 112, no. 1, pp. 1–32, 1994.

H. Yang, “Riemann problems for a class of coupled hyperbolic systems of conservation laws,” Journal of Differential Equations, vol. 159, no. 2, pp. 447–484, 1999.

H. Yang and Y. Zhang, “New developments of delta shock waves and its applications in systems of conservation laws,” Journal of Differential Equations, vol. 252, no. 11, pp. 5951–5993, 2012.

H. Yang and Y. Zhang, “Delta shock waves with Dirac delta function in both components for systems of conservation laws,” Journal of Differential Equations, vol. 257, no. 12, pp. 4369–4402, 2014.

H. Cheng, W. Liu, and H. Yang, “Two-dimensional Riemann problems for zero-pressure gas dynamics with three constant states,” Journal of Mathematical Analysis and Applications, vol. 343, no. 1, pp. 127–140, 2008.

H. Yang, “Generalized plane delta-shock waves for -dimensional zero-pressure gas dynamics,” Journal of Mathematical Analysis and Applications, vol. 260, no. 1, pp. 18–35, 2001.

V. G. Danilov and V. M. Shelkovich, “Dynamics of propagation and interaction of -shock waves in conservation law systems,” Journal of Differential Equations, vol. 211, no. 2, pp. 333–381, 2005.

M. Nedeljkov and M. Oberguggenberger, “Interactions of delta shock waves in a strictly hyperbolic system of conservation laws,” Journal of Mathematical Analysis and Applications, vol. 344, no. 2, pp. 1143–1157, 2008.

C. Shen and M. Sun, “Interactions of delta shock waves for the transport equations with split delta functions,” Journal of Mathematical Analysis and Applications, vol. 351, no. 2, pp. 747–755, 2009.

H. Cheng, “Delta shock waves for a linearly degenerate hyperbolic system of conservation laws of Keyfitz-Kranzer type,” Advances in Mathematical Physics, vol. 2013, Article ID 958120, 10 pages, 2013.

H. Cheng, “Riemann problem for one-dimensional system of conservation laws of mass, momentum and energy in zero-pressure gas dynamics,” Differential Equations & Applications, vol. 4, no. 4, pp. 653–664, 2012.

G.-S. Jiang and E. Tadmor, “Nonoscillatory central schemes for multidimensional hyperbolic conservation laws,” SIAM Journal on Scientific Computing, vol. 19, no. 6, pp. 1892–1917, 1998.