^{1}

^{*}

^{1}

We investigate the elementary wave interactions for the simplified combustion model in magnetogasdynamics. Under the modified entropy conditions, we construct the unique solution and observe some interesting phenomena ; such as , the combustion wave may be extinguished by the contact discontinuity or the shock wave. Especially, the transition between the detonation wave and the deflagration wave is also captured.

Magnetogasdynamics is very important in studying engineering physics [

In [

In [

{ τ t − u x = 0 , u t + ( p + B 2 2 μ ) x = 0 , ( E + B 2 τ 2 μ ) t + ( p u + B 2 u 2 μ ) x = 0 , (1)

under the assumption B = k ρ , where τ > 0 , p ≥ 0 , u, B ≥ 0 and μ are respectively the specific volume, pressure, velocity, transverse magnetic field and magnetic permeability. E = e + u 2 2 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 ( [

{ u t + p x = 0 , τ t − u x = 0 , E t + ( u p ) x = 0 , (2)

with an infinite rate of reaction which is described by

q ( x , t ) = { 0 , if sup 0 ≤ y ≤ t T ( x , y ) > T i ; q ( x , 0 ) , otherwise . (3)

Under the proposed global entropy conditions, they constructed uniquely the Riemann solutions by the characteristic analysis. In [

In [

( τ , p , u , q ) ( x , 0 ) = ( τ ± , p ± , u ± , q ± ) , ± x > 0 , (4)

where τ ± > 0 , p ± , u ± are arbitrary constants,

q ± = { 0 , if T ± > T i , 0 or q 0 , if T ± ≤ T i ,

and q 0 > 0 is a constant. E = e + u 2 2 + q is the specific total energy, where q is the chemical binding energy. The temperature T satisfies Boyle and Gay-Lussac’s law: p τ = R T . T i is the ignition temperature. For polytropic gases, we know e = e ( T ) and E = u 2 2 + p τ γ − 1 + q , where γ > 1 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 [

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 [

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 [

There are three eigenvalues of (1) which are λ 1 = − ( p − e p B B τ μ + e τ e p ) 1 2 , λ 2 = 0 and λ 3 = ( p − e p B B τ μ + e τ e p ) 1 2 . If e p > 0 and e τ + p > 0 , (1) is strictly hyperbolic. The characteristic fields λ 1,3 are genuinely nonlinear and λ 2 is linearly degenerate.

Considering the self-similar solution ( τ , p , u ) ( ζ ) ( ζ = x t ) , for any smooth solution we have

{ ζ d τ = − d u , ζ d u = d ( p + B 2 u 2 μ ) , ζ d ( E + B 2 u 2 μ τ ) = d ( u p + B 2 u 2 μ u ) . (5)

The forward or backward rarefaction waves R → ← passing through the point ( τ 0 , p 0 , u 0 ) are

{ p τ γ = p 0 τ 0 γ , u = u 0 ± ∫ p 0 p γ p τ + B 2 τ μ γ p d p . (6)

The Rankine-Hugoniot jump conditions at ζ = σ are as follows

{ σ [ τ ] = − [ u ] , σ [ u ] = [ p + B 2 u 2 μ ] , σ [ E + B 2 u 2 μ τ ] = [ u p + B 2 u 2 μ u ] , (7)

where [ τ ] = τ r − τ l , etc.

The contact discontinuity J is given by

[ u ] = [ p + B 2 2 μ ] = 0 , (8)

and J is a curve in the ( τ , p , u ) space and the projection on the ( p , u ) plane is a straight line parallel to the p-axis. Denote J by J < when p l < p r , τ l < τ r , and J > when p l > p r , τ l > τ r .

If [ q ] = 0 in (7), we get the forward or backward shock waves S → ← passing through the point ( τ 0 , p 0 , u 0 )

{ ( p + θ 2 p 0 + θ 2 ( 3 B 2 2 μ + B 0 2 2 μ ) ) τ = ( p 0 + θ 2 p + θ 2 ( 3 B 0 2 2 μ + B 2 2 μ ) ) τ 0 , u = u 0 ± ( p + B 2 2 μ − p 0 − B 0 2 2 μ ) − τ − τ 0 p + B 2 2 μ − p 0 − B 0 2 2 μ , (9)

where θ 2 = γ − 1 γ + 1 and B 0 = k τ 0 .

If [ q ] ≠ 0 in (7), we obtain the combustion wave curve in the ( τ , p ) plane

D u ( 0 ) : ( τ − θ 2 τ 0 ) ( p + θ 2 ( p 0 + B 0 2 2 μ + 3 B 2 2 μ ) ) = ( 1 − θ 4 ) τ 0 p 0 + θ 2 τ 0 μ [ B 0 2 ( 3 − θ 2 ) + B 2 ( 1 − 3 θ 2 ) ] + 2 θ 2 q 0 . (10)

Draw two straight lines from ( τ 0 , p 0 ) and they are tangent to the above curve. We call the tangent points A with τ < τ 0 and B with τ > τ 0 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

From the known Jouguet’s rule in [

1) S u ( l ) or R u ( l ) (containing no combustion waves),

2) ( i ) + W D F ( i ) or ( i ) + C J D F ( i ) + R ( C J D F ( i ) ) (containing no DT waves),

3) S D T ( l ) or C J D T ( l ) + R ( C J D T ( l ) ) (containing no DF waves),

where i ≡ i ( l ) ≡ ( u i , p i , τ i , q ) is the state at S ( l ) with the ignition temperature T i , 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 T i , 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 ( τ , p ) -plane we have

R u ( l ) : p τ γ = p − τ − γ , ( 0 < p < p l ) ,

S u ( l ) : ( τ − θ 2 τ l ) ( p + θ 2 ( p l + B l 2 2 μ + 3 B 2 2 μ ) ) = ( 1 − θ 4 ) τ l p l + θ 2 τ l μ [ B l 2 ( 3 − θ 2 ) + B 2 ( ( 1 − 3 θ 2 ) ] , ( p > p l ) ,

S D T ( l ) : ( τ − θ 2 τ l ) ( p + θ 2 ( p l + B l 2 2 μ + 3 B 2 2 μ ) ) = ( 1 − θ 4 ) τ l p l + θ 2 τ l μ [ B l 2 ( 3 − θ 2 ) + B 2 ( 1 − 3 θ 2 ) ] + 2 θ 2 q 0 , ( p > p A ) ,

W D F ( i ) : ( τ − θ 2 τ i ) ( p + θ 2 ( p i + B i 2 2 μ + 3 B 2 2 μ ) ) = ( 1 − θ 4 ) τ i p i + θ 2 τ i μ [ B i 2 ( 3 − θ 2 ) + B 2 ( 1 − 3 θ 2 ) ] + 2 θ 2 q 0 , ( ( p D ) i < p < p i ) ,

R ( C J D T ( l ) ) : p τ γ = p A τ C γ , ( p < p A ) ,

R ( C J D F ( l ) ) : p τ γ = ( p B ) i ( τ B ) i γ , ( p < ( p B ) i ) .

Denote W ( l ) : = W S ( l ) ∪ W D ( l ) , where W S ( l ) denotes ( W S ( l ) , q l = 0 ) or ( W S ( l ) , q l > 0 ) or both of them, and W D ( l ) denotes W D T ( l ) ∪ W D F ( l ) , here W D T ( l ) : = S D T ( l ) ∪ C J D T ( l ) ∪ R ( C J D T ( l ) ) and W D F ( l ) : = W D F ( i ) ∪ C J D F ( i ) ∪ R ( C J D T ( i ) ) .

Now we study the combustion wave curves in the ( u , p ) plane and construct the backward combustion wave curve W ← ( l ) from the state ( l ) = ( τ l , p l , u l , q l ) .

From (9), for the backward wave S ← τ ( Q 0 τ ) we get

u = u 0 − ( p + B 2 2 μ − p 0 − B 0 2 2 μ ) ( τ 0 − τ ) , u < u 0 ,

and from (10), we know

τ = ( 1 − θ 4 ) τ 0 p 0 + θ 2 τ 0 μ [ B 0 2 ( 3 − θ 2 ) + B 2 ( 1 − 3 θ 2 ) ] + 2 θ 2 q 0 p + θ 2 ( p 0 + B 0 2 2 μ + 3 B 2 2 μ ) + θ 2 τ 0 ,

and

τ − τ 0 = ( 1 − θ 2 ) τ 0 ( p 0 − p ) + θ 2 τ 0 μ ( B 0 2 − B 2 ) + 2 θ 2 q 0 p + θ 2 ( p 0 + B 0 2 2 μ + 3 B 2 2 μ ) ,

thus we obtain the backward combustion wave curve in the ( u , p ) plane (see

D ← τ ( 0 ) : u = u 0 − ( p + B 2 2 μ − p 0 − B 0 2 2 μ ) ( 1 − θ 2 ) τ 0 ( p − p 0 ) + θ 2 τ 0 μ ( B 2 − B 0 2 ) − 2 θ 2 q 0 p + θ 2 ( p 0 + B 0 2 2 μ + 3 B 2 2 μ ) . (11)

Now denote the backward DF and DT wave curve by W ← D F ( l ) and W ← D T ( l ) , respectively, where

W ← D F ( l ) : = W D F → ( i s ) ∪ C J D F → ( i s ) ∪ R ← ( C J D F → ( i s ) ) ,

W ← D T ( l ) : = S D T → ( l ) ∪ C J D T → ( l ) ∪ R ← ( C J D T → ( l ) ) .

Denote the backward wave curve W ← ( l ) which can be linked to the state ( l ) = ( τ l , p l , u l , q l ) , then

W ← ( l ) : = W ← S ( l ) ∪ W ← D ( l ) ,

W ← S ( l ) : = ( W ← S ( l ) , q l = 0 ) or ( W ← S ( l ) , q l > 0 ) ,

W ← D ( l ) : = W ← D F ( l ) ∪ W ← D T ( l ) .

Similarly, we can construct the forward wave curve W → ( r ) that can be linked to the state ( r ) = ( τ r , p r , u r , q r ) .

Since the image of J in ( τ , p , u ) is a straight line which parallels to the τ-axis and the projection on the plane ( u , p ) is a point, J is a plane curve in ( τ , p , u ) and the projection on the plane ( u , p ) 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 q l = q r = 0 , the gas on both sides are burnt, no combustion wave will occur.

When q l and q r are not both zero, there may exist more than one intersection points of W ← ( l ) and W → ( r ) . 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) ( [

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 T ( ξ ) between the set { ξ ∈ R 1 : T ( ξ ) ≤ T i } and the set { ξ ∈ R 1 : T ( ξ ) > T i } ;

3) the solution containing as many combustion wave as possible.

Case 1. q l > 0 , q r = 0 . In this case, the gas is unburnt on the left side, the gas is burnt on the right side, i.e. W ← ( l ) = W ← S ( l ) ∪ W ← D F ( l ) ∪ W ← D T ( l ) , W → ( r ) = W → S ( r ) . When there exists only one intersection point of W ← ( l ) and W → ( r ) , we obtain the unique solution is a detonation wave solution D T → + R → or S → if p l τ l γ = p r τ r γ , or D T → + J + R → or S → if p l τ l γ ≠ p r τ r γ .

When there are three intersection points of W ← ( l ) and W → S ( r ) (see

Subcase 1.1. p l τ l γ = p r τ r γ .

When T r ≤ T i , then β ( ∗ S ) = 0 , β ( ∗ D F ) = 2 , from the condition B, we select ∗ S and obtain a non-combustion wave solution R ← or S ← + R → or S → (

When T r > T i , then β ( ∗ S ) = 1 , β ( ∗ D F ) = 1 , from the condition C, we select ∗ D F and obtain a combustion wave solution D T → + R → or S → (

Subcase 1.2. p l τ l γ ≠ p r τ r γ .

From the condition A, we discard the possible detonation DT wave solution and find that the possible Riemann solution is R ← or S ← + J + R → or S → or D T → + J + R → or S → . According to the modified global entropy conditions we obtain the unique Riemann solution as follows (see

1) When T r > T i , T D F > T i , then β ( ∗ S ) = 1 , β ( ∗ D F ) = 1 , from the condition C, we select ∗ D F and obtain a combustion wave solution containing a DF (

2) When T r > T i , T D F ≤ T i ( ⇒ T S ≤ T i ), then β ( ∗ S ) = 1 , β ( ∗ D F ) = 3 , from the condition B, we select ∗ S and obtain a non-combustion wave solution (

3) When T r ≤ T i , T S ≤ T i , then β ( ∗ S ) = 0 , β ( ∗ D F ) = 2 , from the condition B, we select ∗ S and obtain a non-combustion wave solution (

4) When T r ≤ T i , T S > T i ( ⇒ T D F > T i ), then β ( ∗ S ) = 2 , β ( ∗ D F ) = 2 , from the condition C, we select ∗ D F and obtain a combustion wave solution containing a DF (

Case 2. q l > 0 , q r = 0 and there are two intersection points of W ← ( l ) and W → S ( r ) (see

Subcase 2.1. p l τ l γ = p r τ r γ .

In this case, we select the point ∗ S or ∗ D T and obtain the possible solutions S ← + R → or S → or D T → + R → or S → . Now we select the unique Riemann solution as follows.

When T r > T i , then β ( ∗ S ) = 1 , β ( ∗ D T ) = 1 , from the condition C, we select ∗ D T and obtain a combustion wave solution D T → + R → or S → (

When T r ≤ T i , then β ( ∗ S ) = 0 , β ( ∗ D T ) = 2 , from the condition B, we select ∗ S and obtain a non-combustion wave solution S ← + R → or S → (

Subcase 2.2. p l τ l γ ≠ p r τ r γ .

There are two possibilities: one is that there exists only one intersection point of W ← ( l ) and W → ( r ) and we obtain the unique Riemann solution is D T → + J > + R →

or S → , the other one is that there are three possible solutions which are the noncombution wave solution S ← or R ← + J + R → or S → , or the DF combustion wave solution D T → + J + R → or S → , or the DT combustion wave solution D T → + J < + R → or S → . Similarly, according to the modified global entropy conditions we obtain the unique solution as follows (see

1) When T r > T i , T D F > T i , then β ( ∗ S ) = 1 , β ( ∗ D F ) = 1 , from the condition C, we select ∗ D F and obtain a combustion wave solution containing a DF (

2) When T r > T i , T D F ≤ T i ( ⇒ T S ≤ T i ), then β ( ∗ S ) = 1 , β ( ∗ D F ) = 3 , from the condition B, we select ∗ S and obtain a non-combustion wave solution (

3) When T r ≤ T i , T S ≤ T i , then β ( ∗ S ) = 0 , β ( ∗ D F ) = 2 , from the condition B, we select ∗ S and obtain a non-combustion wave solution (

4) When T r ≤ T i , T S > T i ( ⇒ T D F > T i ), then β ( ∗ S ) = 2 , β ( ∗ D F ) = 2 , from the condition C, we select ∗ D F and obtain a combustion wave solution containing a DF (

Case 3. q l > 0 , q r > 0 and the gas on the both sides are unburnt. In this case ,we know that W ← ( l ) = W ← S ( l ) ∪ W ← D F ( l ) ∪ W ← D T ( l ) , W → ( r ) = W → S ( r ) ∪ W → D F ( r ) ∪ W → D T ( r ) .

If the intersection point of W ← ( l ) and W → ( r ) is unique, the solution is D T → + D T ← if p l τ l γ = p r τ r γ , or D T → + J + D T ← if p l τ l γ ≠ p r τ r γ . Otherwise, there are two possible subcases: one is that there is an intersection point of W ← S ( l ) and W → S ( r ) , the other is that there is no intersection point of W ← S ( l ) and W → S ( r ) .

Case 3.1. In the former subcase (see

Subcase 3.1.1. p l τ l γ = p r τ r γ .

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

It is obvious that β = 0 for (a), and it holds that β = 2 for (b), (c) and (d). From the condition B, we select the intersection point 1 and obtain the unique non-combustion wave solution R ← or S ← + R → or S → .

Subcase 3.1.2. p l τ l γ ≠ p r τ r γ . 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 R ← or S ← + J + R → or S → . The only difference is that here the contact discontinuity appears.

Case 3.2. In the latter subcase, there are only two possibilities: W ← ( l ) intersects W → D T ( r ) only or W → ( + ) intersects W ← D T ( l ) only. We just need to consider the former. If the intersection point is unique, the solution is D T → + D T ← if p − τ − γ = p + τ + γ , or D T → + J + D T ← if p − τ − γ ≠ p + τ + γ , otherwise, there are at most three intersection points (see

Subcase 3.2.1. p l τ l γ ≠ p r τ r γ .

From the condition A, the intersection point of W ← D T ( l ) and W → D T ( r ) should be discarded. We denote the intersection point of W → D T ( r ) and W ← S ( l ) by S ∗ and denote the intersection point of W → D T ( + ) and W ← D F ( l ) by D ∗ F ,

respectively. We denote the temperature at the point S ∗ , D ∗ F on W → D T ( r ) by T S , T D F , respectively (see

Since T S > T i we have T D F > T i , then β ( S ∗ ) = 2 , β ( D ∗ F ) = 2 . From the condition C, we select D ∗ F and obtain a combustion wave solution D F → + D T ← .

Subcase 3.2.2. p l τ l γ = p r τ r γ . 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 D F → + D T ← . 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).

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

( B , ρ , u ) ( x , 0 ) = { ( B l , ρ l , u l ) , − ∞ < x ≤ x 1 , ( B m , ρ m , u m ) , x 1 < x ≤ x 2 , ( B r , ρ r , u r ) , x 2 < x < ∞ , (12)

for arbitrary x 1 , x 2 ∈ R .

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 S D T → or W D F → and a contact discontinuity: J > S D T → , J < S D T → , J > W D F → , J < W D F → ;

The interaction of a combustion wave S D T → or W D F → and a shock wave: S D T → S → , W D F → S → .

In what follows, we construct the solutions of the wave interactions case by case.

Case 1. J > S D T →

In the ( x , t ) plane, J > and S D T → 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 (

From the analysis of the wave curves in the ( u , p ) plane, we know that W ← ( l ) = W ← S ( l ) ∪ W ← D F ( l ) ∪ W ← D T ( l ) , and W → ( r ) = W → S ( r ) .

There are two possible cases: there is only one intersection point of W ← D T ( l ) and W → S ( r ) , or there are three intersection points of W ← ( l ) and W → S ( r ) . From the arguments in the former section, we know the result is as follows:

When p l τ l γ = p r τ r γ , from the modified global entropy conditions we obtain that

J > S D T → → { S D T → + R → or S → , S ← or R ← + R → or S → , D F → + R → or S → , (13)

where “ → ” means the result of the wave interaction.

When p l τ l γ ≠ p r τ r γ , the results are given by

J > S D T → → { S D T → + J + R → or S → , S ← or R ← + J + R → or S → , D F → + J + R → or S → . (14)

Case 2. J < S D T →

In this case (

J < S D T → → { S ← or R ← + R → or S → , S D T → + R → or S → , D F → + R → or S → , (15)

when p l τ l γ ≠ p r τ r γ , 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. J > W D F →

For this case, we know that J > and W D F → will encounter each other at a finite time in the ( x , t ) 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 ( u , p ) plane.

It is easily shown that u l = u m and p l > p m , τ l > τ m . We proceed as follows (

J > W D F → → { S ← or R ← + R → or S → , D F → + R → or S → , (16)

when p l τ l γ ≠ p r τ r γ ,

J > W D F → → { S ← or R ← + J + R → or S → , D F → + J + R → or S → . (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. J < W D F →

In this case (

When p l τ l γ = p r τ r γ , the results are as follows

J < D F → → { D T → + R → or S → , S ← or R ← + R → or S → , D F → + R → or S → , (18)

when p l τ l γ ≠ p r τ r γ , the results are the same as the above except that J appears here.

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. S D T → S →

In this case (

straight line p = p m ( u > u m ) . Therefore, W → S ( r ) does not intersect with W ← S ( l ) ∪ W ← D F ( l ) and the result of the interaction is expressed by

S D T → S → → { S D T → + R → or S → , p l τ l γ = p r τ r γ , S D T → + J + R → or S → , p l τ l γ ≠ p r τ r γ . (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. W D F → S →

In this case (

W D F → S → → { W D F → + R → or S → , R → or S → , D T → + R → or S → , (20)

when p l τ l γ ≠ p r τ r γ , 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

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.

Partially by the Foundation for Young Scholars of Shandong University of Technology (No. 115024).

The authors declare that there is no conflict of interest regarding the publication of this paper.

Liu, Y.J. and Sun, W.H. (2019) Elementary Wave Interactions for the Simplified Combustion Model in Magnetogasdynamics. Journal of Applied Mathematics and Physics, 7, 983-1000. https://doi.org/10.4236/jamp.2019.74066