Elementary Wave Interactions for the Simplified Combustion Model in Magnetogasdynamics ()
1. Introduction
Magnetogasdynamics is very important in studying engineering physics [1] - [9] . It is difficult to investigate the governing equations of Magnetogasdynamics flows; the corresponding results are less than the conventional gas dynamics. When the velocity field and the magnetic field are everywhere orthogonal, the magnetogasdynamics flow is still important.
In [3] , Helliwell discussed the non-conducting inviscid gas at rest when ionization of the gas takes place across a shock wave and obtained that the magnetogasdynamic combustion wave has similar properties with the conventional gas dynamics. In [6] , Mareev investigated the above problem further and they could use the obtained results to study the hypersonic gas flow around a thin wedge in an axial magnetic field.
In [4] , Hu and Sheng constructed the unique Riemann solution of the one-dimensional inviscid flow
(1)
under the assumption
, where
,
, u,
and
are respectively the specific volume, pressure, velocity, transverse magnetic field and magnetic permeability.
is the specific total energy. and e is the specific internal energy. The Riemann problem of the conventional gas dynamics combustion models was investigated by many people ( [10] [11] [12] [13] , etc.). Zhang and Zheng [13] studied the Riemann problem of the conventional gas dynamics flow of combustible ideal gases
(2)
with an infinite rate of reaction which is described by
(3)
Under the proposed global entropy conditions, they constructed uniquely the Riemann solutions by the characteristic analysis. In [11] , we modified the above global entropy conditions and constructed the unique solution of the generalized Riemann problem for (2) and (3).
In [5] , we obtained uniquely the Riemann solutions for (1) and (3) under the modified entropy conditions in [11] with the following initial data
(4)
where
are arbitrary constants,
and
is a constant.
is the specific total energy, where q is the chemical binding energy. The temperature T satisfies Boyle and Gay-Lussac’s law:
.
is the ignition temperature. For polytropic gases, we know
and
, where
is the adiabatic exponent. For simplicity, we usually assume that R and
remain unchanged during the reaction. We also assume that the combustion process is exothermic, i.e. the energy used up in recombing the atoms to form the new molecules is smaller for the burnt gas than the binding energy of the unburnt gas [10] .
In the present paper, we are concerned with the wave interactions of the elementary waves of the Chapman-Jouguet (CJ) model (1) and (3). we can capture some interesting combustion phenomena by investigating the elementary wave interactions. For example, in most cases expect for Case 5. In Section 3, the combustion wave may be extinguished by the contact discontinuity or the shock wave. Especially, a detonation wave may be transformed into a deflagration wave by the contact discontinuity (see Case 1 and Case 2 in Section 3) and a deflagration wave may be transformed into a detonation wave by the contact discontinuity (see Case 4 in Section 3) or by the shock wave (see Case 6 in Section 3) which shows the transition between the detonation wave and the deflagration wave. In the case that when combustion waves interact with the contact discontinuity, our results are very different from the conventional gas dynamics combustion model (2) and (3) where there is no transition from the detonation (deflagration) wave to the deflagration (detonation) wave and the combustion wave can not be extinguished by the contact discontinuity.
This paper is organized as follows. In Section 2, we present the results of the Riemann problem for the CJ model (1), (3) with the initial values (4). In Section 3, the elementary wave interactions are considered case by case under the modified entropy conditions in [11] .
2. Preliminaries
As a preparation, we study the Riemann problem for the CJ model (1), (3) with the initial data (4) and we refer the detailed discussions to [4] [5] .
There are three eigenvalues of (1) which are
,
and
. If
and
, (1) is strictly hyperbolic. The characteristic fields
are genuinely nonlinear and
is linearly degenerate.
Considering the self-similar solution
, for any smooth solution we have
(5)
The forward or backward rarefaction waves
passing through the point
are
(6)
The Rankine-Hugoniot jump conditions at
are as follows
(7)
where
, etc.
The contact discontinuity J is given by
(8)
and J is a curve in the
space and the projection on the
plane is a straight line parallel to the p-axis. Denote J by
when
, and
when
.
If
in (7), we get the forward or backward shock waves
passing through the point
(9)
where
and
.
If
in (7), we obtain the combustion wave curve in the
plane
(10)
Draw two straight lines from
and they are tangent to the above curve. We call the tangent points A with
and B with
Chapman-Jouguet detonation (CJDT) and Chapman-Jouguet deflagration (CJDF), respectively. From the RH condition (7), we should disregard the curve between the points C and D. We call the curve between C and A weak detonation (WDT) and the curve above A strong detonation (SDT), the curve between D and B weak deflagration (WDF) and the curve below B strong deflagration (SDF), respectively (see Figure 1).
Figure 1. The combustion wave curve in the (τ, p).
From the known Jouguet’s rule in [13] , there are at most three different kinds of wave series that can be linked to the state
:
1)
or
(containing no combustion waves),
2)
or
(containing no DT waves),
3)
or
(containing no DF waves),
where
is the state at
with the ignition temperature
, and the symbol “+” means “followed by”. Notice that we let the temperature behind the pre-compressive shock wave which connects the state (l) and the ignition point (i) be the ignition point
, we just need to construct the deflagration wave curve which is the successor to the pre-compressive shock wave from the point (i).
In the
-plane we have
Denote
, where
denotes
or
or both of them, and
denotes
, here
and
.
Now we study the combustion wave curves in the
plane and construct the backward combustion wave curve
from the state
.
From (9), for the backward wave
we get
and from (10), we know
and
thus we obtain the backward combustion wave curve in the
plane (see Figure 2).
(11)
Now denote the backward DF and DT wave curve by
and
, respectively, where
Figure 2. The combustion wave curve in the plane (u, p).
Denote the backward wave curve
which can be linked to the state
, then
Similarly, we can construct the forward wave curve
that can be linked to the state
.
Since the image of J in
is a straight line which parallels to the τ-axis and the projection on the plane
is a point, J is a plane curve in
and the projection on the plane
is a straight line which parallels to the p-axis. Thus the Riemann prblem for (1) is much more complicated than that of the conventional gas dynamics.
When
, the gas on both sides are burnt, no combustion wave will occur.
When
and
are not both zero, there may exist more than one intersection points of
and
. Each intersection point corresponds to a unique Riemann solution. When the intersection point is unique, the solution is also unique, otherwise, in order to obtain the unique solution we select it under the following modified global entropy conditions (MGEC) ( [11] ):
We select the unique solution from the nine intersection points (at most) of the forward wave curves connecting (r) and the backward wave curves connecting (l) in the following order:
1) the solution with the propagating speed of combustion wave as low as possible;
2) the solution with the parameter
as small as possible, where
is defined as oscillation frequency of
between the set
and the set
;
3) the solution containing as many combustion wave as possible.
Case 1.
. In this case, the gas is unburnt on the left side, the gas is burnt on the right side, i.e.
,
. When there exists only one intersection point of
and
, we obtain the unique solution is a detonation wave solution
or
if
, or
or
if
.
When there are three intersection points of
and
(see Figure 3), from the modified global entropy condition A, we discard the intersection point of
and
. Denote the intersection point of
and
by
and the intersection point of
and
by
. Denote the temperature at the point
,
on
by
,
, respectively. The temperature at
on
is greater than
since the combustion process is exothermic.
Subcase 1.1.
.
Figure 3. ql > 0, qr = 0 and there are there interactions.
When
, then
,
, from the condition B, we select
and obtain a non-combustion wave solution
or
or
(Figure 4(a)).
When
, then
,
, from the condition C, we select
and obtain a combustion wave solution
or
(Figure 4(b)).
Subcase 1.2.
.
From the condition A, we discard the possible detonation DT wave solution and find that the possible Riemann solution is
or
or
or
or
. According to the modified global entropy conditions we obtain the unique Riemann solution as follows (see Figure 5).
1) When
,
, then
,
, from the condition C, we select
and obtain a combustion wave solution containing a DF (Figure 5(b)).
2) When
,
(
), then
,
, from the condition B, we select
and obtain a non-combustion wave solution (Figure 5(a)).
3) When
,
, then
,
, from the condition B, we select
and obtain a non-combustion wave solution (Figure 5(a)).
4) When
,
(
), then
,
, from the condition C, we select
and obtain a combustion wave solution containing a DF (Figure 5(b)).
Case 2.
and there are two intersection points of
and
(see Figure 6).
Subcase 2.1.
.
In this case, we select the point
or
and obtain the possible solutions
or
or
or
. Now we select the unique Riemann solution as follows.
When
, then
,
, from the condition C, we select
and obtain a combustion wave solution
or
(Figure 7(a)).
When
, then
,
, from the condition B, we select
and obtain a non-combustion wave solution
or
(Figure 7(b)).
Subcase 2.2.
.
There are two possibilities: one is that there exists only one intersection point of
and
and we obtain the unique Riemann solution is
Figure 4. Solutions in Subcase 1.1. (a) non-combustion wave; (b) combustion wave DF.
Figure 5. Solutions in Subcase 1.2. (a) non-combustion wave; (b) combustion wave DF.
Figure 6. ql > 0, qr = 0 and there are two interactions.
Figure 7. Solutions in Subcase 2.1. (a) combustion wave DT; (b) non-combustion wave.
or
, the other one is that there are three possible solutions which are the noncombution wave solution
or
or
, or the DF combustion wave solution
or
, or the DT combustion wave solution
or
. Similarly, according to the modified global entropy conditions we obtain the unique solution as follows (see Figure 8). From the global entropy condition A, we discard the DT combustion wave solution.
1) When
,
, then
,
, from the condition C, we select
and obtain a combustion wave solution containing a DF (Figure 8(b)).
2) When
,
(
), then
,
, from the condition B, we select
and obtain a non-combustion wave solution (Figure 8(a)).
Figure 8. Solutions in Subcase 2.2. (a) non-combustion wave; (b) combustion wave DF.
3) When
,
, then
,
, from the condition B, we select
and obtain a non-combustion wave solution (Figure 8(a)).
4) When
,
(
), then
,
, from the condition C, we select
and obtain a combustion wave solution containing a DF (Figure 8(b)).
Case 3.
and the gas on the both sides are unburnt. In this case ,we know that
,
.
If the intersection point of
and
is unique, the solution is
if
, or
if
. Otherwise, there are two possible subcases: one is that there is an intersection point of
and
, the other is that there is no intersection point of
and
.
Case 3.1. In the former subcase (see Figure 9), we discuss it in the following two subcases.
Subcase 3.1.1.
.
From the condition A, we just need to consider the intersection points 1, 2, 3, 4. We should select the unique solution from the four possible solutions (see Figure 10).
It is obvious that
for (a), and it holds that
for (b), (c) and (d). From the condition B, we select the intersection point 1 and obtain the unique non-combustion wave solution
or
or
.
Subcase 3.1.2.
. In a similar way as the above discussions in Subcase 3.1.1., we obtain that the unique Riemann solution is still the non-combustion wave solution
or
or
. The only difference is that here the contact discontinuity appears.
Case 3.2. In the latter subcase, there are only two possibilities:
intersects
only or
intersects
only. We just need to consider the former. If the intersection point is unique, the solution is
if
, or
if
, otherwise, there are at most three intersection points (see Figure 11).
Subcase 3.2.1.
.
From the condition A, the intersection point of
and
should be discarded. We denote the intersection point of
and
by
and denote the intersection point of
and
by
,
Figure 9. There is an intersection point of
and
.
Figure 10. The possible solutions in Subcase 3.1.1. (a) solution corresponds to point 1; (b) solution corresponds to point 2; (c) solution corresponds to point 3; (d) solution corresponds to point 4.
Figure 11.
intersects
only.
respectively. We denote the temperature at the point
,
on
by
,
, respectively (see Figure 12).
Since
we have
, then
,
. From the condition C, we select
and obtain a combustion wave solution
.
Subcase 3.2.2.
. In a similar way as the above discussions in Subcase 3.2.1., we obtain that the unique Riemann solution is still the combustion wave solution
. The only difference is that there exists the contact discontinuity in Subcase 3.2.1.
Based on the above analysis, we have the following result.
Theorem 2.1. There exists a unique piecewise smooth solution to the Riemann problem (1) and (3) with the initial data (4) under the modified global entropy conditions (MGEC).
3. Wave Interactions for the Combustion Model (1) and (3)
In this section, we discuss the wave interactions of the elementary waves for our combustion model (1) and (3). Consider the equations (1) and (3) with the following initial data
(12)
for arbitrary
.
In order to capture the interesting combustion phenomena, in the present paper we investigate the following kinds of interactions:
The interaction of a combustion wave
or
and a contact discontinuity:
,
,
,
;
The interaction of a combustion wave
or
and a shock wave:
,
.
In what follows, we construct the solutions of the wave interactions case by case.
Case 1.
In the
plane,
and
will encounter each other at a finite time and a new Riemann problem is formed with (l) and (r) as its left-hand side state and right-hand side state, respectively (Figure 13). We solve this new Riemann problem in the
plane. It is obvious that
,
and
.
From the analysis of the wave curves in the
plane, we know that
, and
.
There are two possible cases: there is only one intersection point of
and
, or there are three intersection points of
and
. From the arguments in the former section, we know the result is as follows:
When
, from the modified global entropy conditions we obtain that
(13)
where “
” means the result of the wave interaction.
Figure 13. The interaction of
and
.
When
, the results are given by
(14)
Case 2.
In this case (Figure 14), we know
, and for the contact discontinuity it holds that
and
. Similar with the discussions as the above case, we know that when
,
(15)
when
, the results are the same as the above except that J appears here.
Theorem 3.1. For Case 1 and Case 2, when a strong detonation combustion wave interacts with a contact discontinuity which is of a jump decrease (or increase) in density in the propagating direction of the strong detonation wave, the strong detonation combustion wave may be extinguished or cross the contact discontinuity at once, or be transformed into a deflagration wave. And the contact discontinuity may disappear or not.
Remark 3.2. The result of the interaction of the strong detonation wave and the contact discontinuity is very different from that of the corresponding case in the conventional gas dynamics combustion model (4) where there is no transition from the detonation wave to the deflagration wave and the combustion wave can not be extinguished by the contact discontinuity.
Case 3.
For this case, we know that
and
will encounter each other at a finite time in the
plane, and a new Riemann problem is formed with (l) and (r) as its left-hand side state and right-hand side state, respectively. Similarly, we solve the new Riemann problem in the
plane.
Figure 14. The interaction of
and
.
It is easily shown that
and
,
. We proceed as follows (Figure 15). Based on the analysis in the
plane, we notice specially that the temperature is just the ignition point on the wave front of the
, i.e. the state
where the symbol “
” means the ignition point corresponding to the state (m). Since there are three intersection points of
and
, from the global entropy conditions we know that when
,
(16)
when
,
(17)
Theorem 3.3. For this case, when a weak deflagration combustion wave interacts with a contact discontinuity which is of a jump decrease in density in the propagating direction of the weak deflagration wave, the weak deflagration combustion wave may be extinguished, or cross the contact discontinuity at once. And the contact discontinuity may disappear or not.
Case 4.
In this case (Figure 16), we know
, and for the contact discontinuity we know
and
. There are two possibilities: one is that there is only one intersection point of
and
, the other one is that there are three intersection points of
and
. Therefore we know the results are as follows:
When
, the results are as follows
(18)
when
, the results are the same as the above except that J appears here.
Figure 15. The interaction of
and
.
Figure 16. The interaction of
and
.
Theorem 3.4. For this case, when a weak deflagration combustion wave interacts with a contact discontinuity which is of a jump increase in density in the propagating direction of the weak deflagration wave, the weak deflagration combustion wave may be transformed into a detonation wave, or be extinguished, or cross the contact discontinuity at once. And the contact discontinuity may disappear or not.
Remark 3.5. The result of the interaction of the weak deflagration wave and the contact discontinuity is very different from that of the corresponding case in the conventional gas dynamics combustion model (2), where there is no transition from the deflagration wave to the detonation wave and the combustion wave can not be extinguished by the contact discontinuity.
Case 5.
In this case (Figure 17), we know for
it holds that
, and for
it holds that
, where
and
are respectively the velocity of
and
. Therefore,
will overtake
at a finite time and a new Riemann problem is formed. We still solve this problem in the
plane.
,
and the curves
are located below the curve
according to the condition A. From Lemma 2.6., we know
are located below
and the half straight line
. Considering that
and
from Lemma 2.6., we know that
is located above the curve
and the half
Figure 17. The interaction of
and
.
straight line
. Therefore,
does not intersect with
and the result of the interaction is expressed by
(19)
Theorem 3.6. For this case, when a shock wave overtakes a strong detonation combustion wave, it speeds up the strong detonation combustion wave and the contact discontinuity may appear or not.
Case 6.
In this case (Figure 18), similar with the discussions as the just above case, we know that
will overtake
at a finite time and a new Riemann problem is formed with (l) and (r) as its left-hand side state and right-hand side state, respectively. Notice that
,
and (r) is located above the curve
, we know the temperature at this point is higher than the ignition temperature, therefore when
, the results is described by
(20)
when
, the result is the same expect that the contact discontinuity appears here.
Theorem 3.7. For this case, when a shock wave overtakes a weak deflagration combustion wave, the weak deflagration combustion wave may be extinguished or be transformed into a detonation wave. And the contact discontinuity may appear or not.
After the discussions of the elementary wave interactions, we summarize the results as follows.
By investigating the kinds of wave interactions of the elementary waves, we can capture some interesting combustion phenomena. For example, the combustion wave may be extinguished by the contact discontinuity or by the shock wave. Especially, a detonation wave may be transformed into a weak deflagration wave coalescing with pre-compression shock wave by the contact
Figure 18. The interaction of
and
.
discontinuity, and a deflagration wave may be transformed into a detonation wave by the contact discontinuity or by the shock wave, which describe the transition between the detonation wave and the deflagration wave. Notice specially that the contact discontinuity play a very important role in magnetogasdynamic combustion model (1) and (3), which is very different from the conventional gas dynamic combustion model (2) and (3).
Since the reaction rate in our model is infinite which is an idealized hypothesis, while our model is still very important in application, we will investigate the initial value problem for the self-similar Zeldovich-von Neumann-Döring (ZND) model in magnetogasdynamic combustion with finite reaction rate in our coming works.
Supported
Partially by the Foundation for Young Scholars of Shandong University of Technology (No. 115024).