In this paper, we study the following equations
where , T are the density, velocity, pressure, temperature respectively. is the ignition temperature and q is the chemical binding energy. The total
energy , where e is the internal energy. The state equation is given by and . The process of combustion is exothermic . The early study for Chaplygin gas can be found in  .
In , the authors investigated the two-dimensional Riemann problem for isentropic Chaplygin gas dynamic system according to different combinations of four elementary waves and they delivered a complete classification to the above problem.
In , the authors studied two new types of self-similar solutions to the Chaplygin gas model in two space dimensions and they constructed the solutions to the interaction of two rarefaction simple waves.
The authors  investigated several two-dimensional Riemann problems for Chaplygin gas model and proved the existence of the unique solution.
In , the authors studied the appearance of delta shock wave and vacuum state in the vanishing pressure limit of Riemann solutions to the non-isentropic generalized Chaplygin gas equations.
The authors  constructed the solutions of Riemann problem for the isentropic relativistic Euler equations with the extended Chaplygin gas and studied the asymptotic limits of solutions to the Riemann problem for the relativistic Euler equations.
We use usually two mathematical models to investigate the combustion phenomena. One is the Chapman-Jouguet (CJ) model, and the other is the Zeldovich-von Neumann-Döring (ZND) model  . In , the authors obtained the partial combustion solutions for the CJ model. In , the authors investigated the following CJ model
here and is the specific volume, . It is followed that there is
unique solution for (2). In  we got uniquely the solution of the generalized Riemann problem for (2).
In  the authors studied the following selfsimilar ZND model
and constructed its unique solution under some hypothesis. In , the authors determined when the temperatures along the burning solutions of (3) are higher than the ignition temperature which is the necessary condition to study the limit behavior as the reaction velocity tends to infinity.
In , we obtained the Riemann solution of (1) with the following initial values
In this paper, we investigate the initial value problem for (1) with the following initial values
In , we considered the generalized Riemann problem for (1) and found some interesting phenomena.
It is shown that the structures of the Riemann solutions can retain their forms after perturbation for most of the cases, while for some other cases, the perturbation can make a big difference. It is observed that although there is no combustion wave in the corresponding Riemann solution, the combustion wave may occur after perturbation. It shows the instability of the unburnt gas.
The paper is arranged as follows. In Section 2, we list the preliminaries. We obtain the unique solution for the initial value problem (1) and (5) in Section 3. Section 4 gives our main result.
In this section we give some preliminaries    .
The characteristic roots for (1) are given by
it follows that (1.1) is strictly hyperbolic. The right characteristic vector of is
From , all characteristic fields are linearly degenerate.
The rarefaction wave curves (or ) are
and the shock wave curves (or ) are
The contact discontinuity J is
Further, the delta shock satisfies the entropy condition
The non-combustion wave curves can be shown in the phase plane (Figure 1).
From the R-H condition
Figure 1. The non-combustion wave curves. (i) The backward wave curves; (ii) The forward wave curves.
we find that
The forward combustion wave curve in (Figure 2) is
where or . When , it is the detonation wave curve and when , it is the deflagration wave curve .
When , the forward wave curve is
When , , it is studied in . Otherwise, we proceed as follows.
Case 2.1. When .
The backward wave curve is , and the forward wave curve is (Figure 3).
Figure 2. The combustion wave curve in (u, p).
Figure 3. The wave curves in Case 2.1.
Subcase 2.1.1 .
To guarantee the unique solution of the Riemann problem, we give the global entropy conditions:
We select the unique solution from the four intersection points (at most) in the following order:
1) the solution with the parameter as small as possible, where is defined as oscillation frequency of between the set and the set ;
2) the solution containing the combustion wave as many as possible.
For simplicity, is the intersection point of and and is the intersection point of and or The temperature is respectively , at the point , (Figure 4).
1) When , , then , , due to (ii), we pick out and get the combustion wave solution DF or DT (Figure 4(i)).
2) When , ( ), then , , due to (i), we pick out and get the noncombustion wave solution (Figure 4(ii)).
3) When , , then , , due to (i), we pick out and get the noncombustion wave solution (Figure 4(ii)).
4) When , ( ), then , , due to (ii), we pick out and get the combustion wave solution DF or DT (Figure 4(i)).
Subcase 2.1.2 .
It is shown that we can get the delta shock wave solution   .
Figure 4. The analysis in Subcase 2.1.1. (i) combustion wave solution; (ii) non-combustion wave solution.
Case 2.2. When . We go on as follows.
We know that , and (Figure 5).
Subcase 2.2.1 When (Figure 6).
Since we have for Case (1), for Case (2)-(4), we pick out the intersection point a and get the non-combustion wave solution.
Subcase 2.2.2 When . We get the delta shock wave solution for this case similarly.
Theorem 2.1. Under the given global entropy conditions, we obtain uniquely the Riemann solution for (1) and (4).
3. Elementary Wave Interactions for (1) with (5)
For simplicity, we just consider the case that the delta shock wave does not appear. We consider the four kinds of wave interactions:
Figure 5. The wave curves in Case 2.2.
Figure 6. The discussions in Subcase 2.2.1. Case (i) solution of the point a; Case (ii) solution of the point b; Case (iii) solution of the point c; Case (iv) solution of the point d.
A combustion wave and a shock wave , a combustion wave and a shock wave , a combustion wave and a rarefaction wave , a combustion wave and a rarefaction wave .
Case 3.1. and (Figure 7)
We know for , and for , where and are the velocity of and respectively. , . From the global entropy conditions, we know that the solution for this case is
Theorem 3.1. For the wave interaction between the shock wave and the detonation wave, we find that it may extinguish the combustion wave which shows that the unburnt gas is unstable. And after the elementary wave interaction process, we observe that the contact discontinuity may emerge.
Case 3.2. and
In this case (Figure 8), it follows that the shock wave will overtake the combustion wave and we need to investigate the new initial value problem. Considering the fact that , and due to , we get that the results from the global entropy conditions are described by
Figure 7. The interaction of and .
Figure 8. The interaction of and .
Theorem 3.2. In this situation, the result of the wave interaction between the shock wave and the deflagration wave is that the deflagration combustion wave may be extinguished. Furthermore, we find that the contact discontinuity may emerge after the elementary wave interaction process.
Case 3.3. and
For this situation (Figure 9), it follows that will interact with , and the new initial value problem is formed. Due to , and considering that , again from the global entropy conditions, we know that the result of the wave interaction process can be given by
Theorem 3.3. In this situation, the result of the wave interaction between the rarefaction wave and the detonation wave is that it may extinguish the combustion wave. It shows that the unburnt gas is unstable. Further, we find that the contact discontinuity may emerge after the elementary wave interaction process.
Case 3.4. and
In this case (Figure 10), we know that the rarefaction wave will interact with the combustion wave. Then a new initial value problem is formed. Considering the fact and , due to the global entropy conditions we get that the result of the wave interaction can be described by
Figure 9. The interaction of and .
Figure 10. The interaction of and .
Theorem 3.4. In this situation, the result of the wave interaction between the rarefaction wave and the deflagration wave is the deflagration wave may be extinguished. Furthermore, we find that the contact discontinuity may emerge after the elementary wave interaction process.
Based on the above analysis, we conclude our main results.
Theorem 4.1. There exists unique solution for the initial value problem (1) and (5). By studying the concrete elementary wave interactions, we observe interesting combustion phenomenon. The combustion process may be extinguished which reveals that the unburnt gas is unstable.
For simplicity, we suppose that the reaction rate in the combustion model (1) is infinite, Of course, it is still an important model to study the combustion phenomena in our many applications in every way. In our further study, we will study the self-similar ZND model which has a finite reaction rate.
Supported by the Foundation for Young Scholars of Shandong University of Technology (No. 115024).
 Brenier, Y. (2005) Solutions with Concentration to the Riemann Problem for One-Dimensional Chaplygin Gas Equations. Journal of Mathematical Fluid Mechanics, 7, S326-S331.
 Guo, L.H., Sheng, W.C. and Zhang, T. (2010) The Two-Dimensional Riemann Problem for Isentropic Chaplygin Gas Dynamics System. Communications on Pure & Applied Analysis, 9, 431-458.
 Lai, G., Sheng, W.C. and Zheng, Y.X. (2011) Simple Waves and Pressure Delta Waves for a Chaplygin Gas in Two-Dimensions. Discrete & Continuous Dynamical Systems, 31, 489-523.
 Ding, Q.X. and Guo, L.H. (2019) The Vanishing Pressure Limit of Riemann Solutions to the Non-Isentropic Euler Equations for Generalized Chaplygin Gas. Advances in Mathematical Physics, 2019, Article ID: 5253717.
 Zhang, Y.F. and Sun, M.N. (2020) Concentration Phenomenon of Riemann Solutions for the Relativistic Euler Equations with the Extended Chaplygin Gas. Acta Applicandae Mathematicae, 170, 539-568.
 Tan, D.C. and Zhang, T. (1992) Riemann Problem for the Selfsimilar ZND Model in Gas Dynamical Combustion. Journal of Differential Equations, 95, 331-369.
 Hsu, C.H. and Lin, S.S. (1997) Some Qualitative Properties of the Riemann Problem in Gas Dynamical Combustion. Journal of Differential Equations, 140, 10-43.