1. Introduction
In this paper, we study the following equations
(1)
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 [1]. The early study for Chaplygin gas can be found in [2] [3].
In [4], 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 [5], 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 [6] investigated several two-dimensional Riemann problems for Chaplygin gas model and proved the existence of the unique solution.
In [7], 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 [8] 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 [1] [9]. In [1], the authors obtained the partial combustion solutions for the CJ model. In [10], the authors investigated the following CJ model
(2)
here
and
is the specific volume,
. It is followed that there is
unique solution for (2). In [11] we got uniquely the solution of the generalized Riemann problem for (2).
In [12] the authors studied the following selfsimilar ZND model
(3)
and constructed its unique solution under some hypothesis. In [13], 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 [14], we obtained the Riemann solution of (1) with the following initial values
(4)
In this paper, we investigate the initial value problem for (1) with the following initial values
(5)
In [15], 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.
2. Preliminaries
In this section we give some preliminaries [14] [15] [16] [17].
The characteristic roots for (1) are given by
(6)
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
(7)
and the shock wave curves
(or
) are
(8)
The contact discontinuity J is
(9)
Suppose
(10)
(11)
(12)
When
,
when
,
Further, the delta shock
satisfies the entropy condition
(13)
The non-combustion wave curves can be shown in the phase plane
(Figure 1).
From the R-H condition
(i) (ii)
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
(14)
where
or
. When
, it is the detonation wave curve and when
, it is the deflagration wave curve [1].
When
, the forward wave curve
is
here
When
,
, it is studied in [16]. 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).
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
[14] [15] [16].
(i)(ii)
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 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.
4. Conclusions
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.
Funding
Supported by the Foundation for Young Scholars of Shandong University of Technology (No. 115024).