﻿ Elementary Wave Interactions for the Simplified Combustion Model in Magnetogasdynamics

Journal of Applied Mathematics and Physics
Vol.07 No.04(2019), Article ID:92168,18 pages
10.4236/jamp.2019.74066

Elementary Wave Interactions for the Simplified Combustion Model in Magnetogasdynamics

Yujin Liu*, Wenhua Sun

School of Mathematics and Statistics, Shandong University of Technology, Zibo, China    Received: March 22, 2019; Accepted: April 26, 2019; Published: April 29, 2019

ABSTRACT

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.

Keywords:

Riemann Problem, Wave Interaction, Hyperbolic Conservation Laws, Combustion, Magnetogasdynamics 1. Introduction

Magnetogasdynamics is very important in studying engineering physics  -  . 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  , 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  , 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  , Hu and Sheng constructed the unique Riemann solution of the one-dimensional inviscid flow

$\left\{\begin{array}{l}{\tau }_{t}-{u}_{x}=0,\\ {u}_{t}+{\left(p+\frac{{B}^{2}}{2\mu }\right)}_{x}=0,\\ {\left(E+\frac{{B}^{2}\tau }{2\mu }\right)}_{t}+{\left(pu+\frac{{B}^{2}u}{2\mu }\right)}_{x}=0,\end{array}$ (1)

under the assumption $B=k\rho$ , where $\tau >0$ , $p\ge 0$ , u, $B\ge 0$ and $\mu$ are respectively the specific volume, pressure, velocity, transverse magnetic field and magnetic permeability. $E=e+\frac{{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 (     , etc.). Zhang and Zheng  studied the Riemann problem of the conventional gas dynamics flow of combustible ideal gases

$\left\{\begin{array}{l}{u}_{t}+{p}_{x}=0,\\ {\tau }_{t}-{u}_{x}=0,\\ {E}_{t}+{\left(up\right)}_{x}=0,\end{array}$ (2)

with an infinite rate of reaction which is described by

$q\left(x,t\right)=\left\{\begin{array}{l}0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{ }\text{if}\text{\hspace{0.17em}}\underset{0\le y\le t}{\mathrm{sup}}T\left(x,y\right)>{T}_{i};\\ q\left(x,0\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{ }\text{otherwise}.\text{ }\end{array}$ (3)

Under the proposed global entropy conditions, they constructed uniquely the Riemann solutions by the characteristic analysis. In  , we modified the above global entropy conditions and constructed the unique solution of the generalized Riemann problem for (2) and (3).

In  , we obtained uniquely the Riemann solutions for (1) and (3) under the modified entropy conditions in  with the following initial data

$\left(\tau ,p,u,q\right)\left(x,0\right)=\left({\tau }^{±},{p}^{±},{u}^{±},{q}^{±}\right),±x>0,$ (4)

where ${\tau }^{±}>0,{p}^{±},{u}^{±}$ are arbitrary constants,

${q}^{±}=\left\{\begin{array}{l}0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{ }\text{if}\text{ }{T}^{±}>{T}_{i},\\ 0\text{ }\text{or}\text{ }{q}_{0},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{ }\text{if}\text{ }{T}^{±}\le {T}_{i},\end{array}$

and ${q}_{0}>0$ is a constant. $E=e+\frac{{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\tau =RT$ . ${T}_{i}$ is the ignition temperature. For polytropic gases, we know $e=e\left(T\right)$ and $E=\frac{{u}^{2}}{2}+\frac{p\tau }{\gamma -1}+q$ , where $\gamma >1$ is the adiabatic exponent. For simplicity, we usually assume that R and $\gamma$ 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  .

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   .

There are three eigenvalues of (1) which are ${\lambda }_{1}=-{\left(\frac{p-{e}_{p}\frac{B{B}_{\tau }}{\mu }+{e}_{\tau }}{{e}_{p}}\right)}^{\frac{1}{2}}$ , ${\lambda }_{2}=0$ and ${\lambda }_{3}={\left(\frac{p-{e}_{p}\frac{B{B}_{\tau }}{\mu }+{e}_{\tau }}{{e}_{p}}\right)}^{\frac{1}{2}}$ . If ${e}_{p}>0$ and ${e}_{\tau }+p>0$ , (1) is strictly hyperbolic. The characteristic fields ${\lambda }_{1,3}$ are genuinely nonlinear and ${\lambda }_{2}$ is linearly degenerate.

Considering the self-similar solution $\left(\tau ,p,u\right)\left(\zeta \right)\left(\zeta =\frac{x}{t}\right)$ , for any smooth solution we have

$\left\{\begin{array}{l}\zeta \text{d}\tau =-\text{d}u,\\ \zeta \text{d}u=\text{d}\left(p+\frac{{B}^{2}u}{2\mu }\right),\\ \zeta \text{d}\left(E+\frac{{B}^{2}u}{2\mu }\tau \right)=\text{d}\left(up+\frac{{B}^{2}u}{2\mu }u\right).\end{array}$ (5)

The forward or backward rarefaction waves $\stackrel{←}{\stackrel{\to }{R}}$ passing through the point $\left({\tau }_{0},{p}_{0},{u}_{0}\right)$ are

$\left\{\begin{array}{l}p{\tau }^{\gamma }={p}_{0}{\tau }_{0}^{\gamma },\\ u={u}_{0}±{\int }_{{p}_{0}}^{p}\frac{\sqrt{\gamma p\tau +\frac{{B}^{2}\tau }{\mu }}}{\gamma p}\text{d}p.\end{array}$ (6)

The Rankine-Hugoniot jump conditions at $\zeta =\sigma$ are as follows

$\left\{\begin{array}{l}\sigma \left[\tau \right]=-\left[u\right],\\ \sigma \left[u\right]=\left[p+\frac{{B}^{2}u}{2\mu }\right],\\ \sigma \left[E+\frac{{B}^{2}u}{2\mu }\tau \right]=\left[up+\frac{{B}^{2}u}{2\mu }u\right],\end{array}$ (7)

where $\left[\tau \right]={\tau }_{r}-{\tau }_{l}$ , etc.

The contact discontinuity J is given by

$\left[u\right]=\left[p+\frac{{B}^{2}}{2\mu }\right]=0,$ (8)

and J is a curve in the $\left(\tau ,p,u\right)$ space and the projection on the $\left(p,u\right)$ plane is a straight line parallel to the p-axis. Denote J by $\stackrel{<}{J}$ when ${p}_{l}<{p}_{r},{\tau }_{l}<{\tau }_{r}$ , and $\stackrel{>}{J}$ when ${p}_{l}>{p}_{r},{\tau }_{l}>{\tau }_{r}$ .

If $\left[q\right]=0$ in (7), we get the forward or backward shock waves $\stackrel{←}{\stackrel{\to }{S}}$ passing through the point $\left({\tau }_{0},{p}_{0},{u}_{0}\right)$

$\left\{\begin{array}{l}\left(p+{\theta }^{2}{p}_{0}+{\theta }^{2}\left(\frac{3{B}^{2}}{2\mu }+\frac{{B}_{0}^{2}}{2\mu }\right)\right)\tau =\left({p}_{0}+{\theta }^{2}p+{\theta }^{2}\left(\frac{3{B}_{0}^{2}}{2\mu }+\frac{{B}^{2}}{2\mu }\right)\right){\tau }_{0},\\ u={u}_{0}±\left(p+\frac{{B}^{2}}{2\mu }-{p}_{0}-\frac{{B}_{0}^{2}}{2\mu }\right)\sqrt{-\frac{\tau -{\tau }_{0}}{p+\frac{{B}^{2}}{2\mu }-{p}_{0}-\frac{{B}_{0}^{2}}{2\mu }}},\end{array}$ (9)

where ${\theta }^{2}=\frac{\gamma -1}{\gamma +1}$ and ${B}_{0}=\frac{k}{{\tau }_{0}}$ .

If $\left[q\right]\ne 0$ in (7), we obtain the combustion wave curve in the $\left(\tau ,p\right)$ plane

$\begin{array}{l}{D}_{u}\left(0\right):\left(\tau -{\theta }^{2}{\tau }_{0}\right)\left(p+{\theta }^{2}\left({p}_{0}+\frac{{B}_{0}^{2}}{2\mu }+\frac{3{B}^{2}}{2\mu }\right)\right)\\ =\left(1-{\theta }^{4}\right){\tau }_{0}{p}_{0}+\frac{{\theta }^{2}{\tau }_{0}}{\mu }\left[{B}_{0}^{2}\left(3-{\theta }^{2}\right)+{B}^{2}\left(1-3{\theta }^{2}\right)\right]+2{\theta }^{2}{q}_{0}.\end{array}$ (10)

Draw two straight lines from $\left({\tau }_{0},{p}_{0}\right)$ and they are tangent to the above curve. We call the tangent points A with $\tau <{\tau }_{0}$ and B with $\tau >{\tau }_{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 Figure 1).

Figure 1. The combustion wave curve in the (τ, p).

From the known Jouguet’s rule in  , there are at most three different kinds of wave series that can be linked to the state $\left(l\right)\equiv \left({\tau }_{l},{p}_{l},{u}_{l},{q}_{l}\right)$ :

1) ${S}_{u}\left(l\right)$ or ${R}_{u}\left(l\right)$ (containing no combustion waves),

2) $\left(i\right)+WDF\left(i\right)$ or $\left(i\right)+CJDF\left(i\right)+R\left(CJDF\left(i\right)\right)$ (containing no DT waves),

3) $SDT\left(l\right)$ or $CJDT\left(l\right)+R\left(CJDT\left(l\right)\right)$ (containing no DF waves),

where $i\equiv i\left(l\right)\equiv \left({u}_{i},{p}_{i},{\tau }_{i},q\right)$ is the state at $S\left(l\right)$ 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 $\left(\tau ,p\right)$ -plane we have

${R}_{u}\left(l\right):p{\tau }^{\gamma }={p}_{-}{\tau }_{-}^{\gamma },\left(0

$\begin{array}{l}{S}_{u}\left(l\right):\left(\tau -{\theta }^{2}{\tau }_{l}\right)\left(p+{\theta }^{2}\left({p}_{l}+\frac{{B}_{l}^{2}}{2\mu }+\frac{3{B}^{2}}{2\mu }\right)\right)\\ =\left(1-{\theta }^{4}\right){\tau }_{l}{p}_{l}+\frac{{\theta }^{2}{\tau }_{l}}{\mu }\left[{B}_{l}^{2}\left(3-{\theta }^{2}\right)+{B}^{2}\left(\left(1-3{\theta }^{2}\right)\right],\left(p>{p}_{l}\right),\end{array}$

$\begin{array}{l}SDT\left(l\right):\left(\tau -{\theta }^{2}{\tau }_{l}\right)\left(p+{\theta }^{2}\left({p}_{l}+\frac{{B}_{l}^{2}}{2\mu }+\frac{3{B}^{2}}{2\mu }\right)\right)\\ =\left(1-{\theta }^{4}\right){\tau }_{l}{p}_{l}+\frac{{\theta }^{2}{\tau }_{l}}{\mu }\left[{B}_{l}^{2}\left(3-{\theta }^{2}\right)+{B}^{2}\left(1-3{\theta }^{2}\right)\right]+2{\theta }^{2}{q}_{0},\left(p>{p}_{A}\right),\end{array}$

$\begin{array}{l}WDF\left(i\right):\left(\tau -{\theta }^{2}{\tau }_{i}\right)\left(p+{\theta }^{2}\left({p}_{i}+\frac{{B}_{i}^{2}}{2\mu }+\frac{3{B}^{2}}{2\mu }\right)\right)\\ =\left(1-{\theta }^{4}\right){\tau }_{i}{p}_{i}+\frac{{\theta }^{2}{\tau }_{i}}{\mu }\left[{B}_{i}^{2}\left(3-{\theta }^{2}\right)+{B}^{2}\left(1-3{\theta }^{2}\right)\right]+2{\theta }^{2}{q}_{0},\left({\left({p}_{D}\right)}_{i}

$R\left(CJDT\left(l\right)\right):p{\tau }^{\gamma }={p}_{A}{\tau }_{C}^{\gamma },\left(p<{p}_{A}\right),$

$R\left(CJDF\left(l\right)\right):p{\tau }^{\gamma }={\left({p}_{B}\right)}_{i}{\left({\tau }_{B}\right)}_{i}^{\gamma },\left(p<{\left({p}_{B}\right)}_{i}\right).$

Denote $W\left(l\right):={W}_{S}\left(l\right)\cup {W}_{D}\left(l\right)$ , where ${W}_{S}\left(l\right)$ denotes $\left({W}_{S}\left(l\right),{q}_{l}=0\right)$ or $\left({W}_{S}\left(l\right),{q}_{l}>0\right)$ or both of them, and ${W}_{D}\left(l\right)$ denotes ${W}_{DT}\left(l\right)\cup {W}_{DF}\left(l\right)$ , here ${W}_{DT}\left(l\right):=SDT\left(l\right)\cup CJDT\left(l\right)\cup R\left(CJDT\left(l\right)\right)$ and ${W}_{DF}\left(l\right):=WDF\left(i\right)\cup CJDF\left(i\right)\cup R\left(CJDT\left(i\right)\right)$ .

Now we study the combustion wave curves in the $\left(u,p\right)$ plane and construct the backward combustion wave curve $\stackrel{←}{W}\left(l\right)$ from the state $\left(l\right)=\left({\tau }_{l},{p}_{l},{u}_{l},{q}_{l}\right)$ .

From (9), for the backward wave ${\stackrel{←}{S}}_{\tau }\left({Q}_{0\tau }\right)$ we get

$u={u}_{0}-\sqrt{\left(p+\frac{{B}^{2}}{2\mu }-{p}_{0}-\frac{{B}_{0}^{2}}{2\mu }\right)\left({\tau }_{0}-\tau \right)},u<{u}_{0},$

and from (10), we know

$\tau =\frac{\left(1-{\theta }^{4}\right){\tau }_{0}{p}_{0}+\frac{{\theta }^{2}{\tau }_{0}}{\mu }\left[{B}_{0}^{2}\left(3-{\theta }^{2}\right)+{B}^{2}\left(1-3{\theta }^{2}\right)\right]+2{\theta }^{2}{q}_{0}}{p+{\theta }^{2}\left({p}_{0}+\frac{{B}_{0}^{2}}{2\mu }+\frac{3{B}^{2}}{2\mu }\right)}+{\theta }^{2}{\tau }_{0},$

and

$\tau -{\tau }_{0}=\frac{\left(1-{\theta }^{2}\right){\tau }_{0}\left({p}_{0}-p\right)+\frac{{\theta }^{2}{\tau }_{0}}{\mu }\left({B}_{0}^{2}-{B}^{2}\right)+2{\theta }^{2}{q}_{0}}{p+{\theta }^{2}\left({p}_{0}+\frac{{B}_{0}^{2}}{2\mu }+\frac{3{B}^{2}}{2\mu }\right)},$

thus we obtain the backward combustion wave curve in the $\left(u,p\right)$ plane (see Figure 2).

$\begin{array}{l}{\stackrel{←}{D}}_{\tau }\left(0\right):u\\ ={u}_{0}-\sqrt{\left(p+\frac{{B}^{2}}{2\mu }-{p}_{0}-\frac{{B}_{0}^{2}}{2\mu }\right)\frac{\left(1-{\theta }^{2}\right){\tau }_{0}\left(p-{p}_{0}\right)+\frac{{\theta }^{2}{\tau }_{0}}{\mu }\left({B}^{2}-{B}_{0}^{2}\right)-2{\theta }^{2}{q}_{0}}{p+{\theta }^{2}\left({p}_{0}+\frac{{B}_{0}^{2}}{2\mu }+\frac{3{B}^{2}}{2\mu }\right)}}.\end{array}$ (11)

Now denote the backward DF and DT wave curve by ${\stackrel{←}{W}}_{DF}\left(l\right)$ and ${\stackrel{←}{W}}_{DT}\left(l\right)$ , respectively, where

${\stackrel{←}{W}}_{DF}\left(l\right):=\stackrel{\to }{WDF}\left({i}_{s}\right)\cup \stackrel{\to }{CJDF}\left({i}_{s}\right)\cup \stackrel{←}{R}\left(\stackrel{\to }{CJDF}\left({i}_{s}\right)\right),$

${\stackrel{←}{W}}_{DT}\left(l\right):=\stackrel{\to }{SDT}\left(l\right)\cup \stackrel{\to }{CJDT}\left(l\right)\cup \stackrel{←}{R}\left(\stackrel{\to }{CJDT}\left(l\right)\right).$

Figure 2. The combustion wave curve in the plane (u, p).

Denote the backward wave curve $\stackrel{←}{W}\left(l\right)$ which can be linked to the state $\left(l\right)=\left({\tau }_{l},{p}_{l},{u}_{l},{q}_{l}\right)$ , then

$\stackrel{←}{W}\left(l\right):={\stackrel{←}{W}}_{S}\left(l\right)\cup {\stackrel{←}{W}}_{D}\left(l\right),$

${\stackrel{←}{W}}_{S}\left(l\right):=\left({\stackrel{←}{W}}_{S}\left(l\right),{q}_{l}=0\right)\text{or}\left({\stackrel{←}{W}}_{S}\left(l\right),{q}_{l}>0\right),$

${\stackrel{←}{W}}_{D}\left(l\right):={\stackrel{←}{W}}_{DF}\left(l\right)\cup {\stackrel{←}{W}}_{DT}\left(l\right).$

Similarly, we can construct the forward wave curve $\stackrel{\to }{W}\left(r\right)$ that can be linked to the state $\left(r\right)=\left({\tau }_{r},{p}_{r},{u}_{r},{q}_{r}\right)$ .

Since the image of J in $\left(\tau ,p,u\right)$ is a straight line which parallels to the τ-axis and the projection on the plane $\left(u,p\right)$ is a point, J is a plane curve in $\left(\tau ,p,u\right)$ and the projection on the plane $\left(u,p\right)$ 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 $\stackrel{←}{W}\left(l\right)$ and $\stackrel{\to }{W}\left(r\right)$ . 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 $\beta$ as small as possible, where $\beta$ is defined as oscillation frequency of $T\left(\xi \right)$ between the set $\left\{\xi \in {R}^{1}:T\left(\xi \right)\le {T}_{i}\right\}$ and the set $\left\{\xi \in {R}^{1}:T\left(\xi \right)>{T}_{i}\right\}$ ;

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. $\stackrel{←}{W}\left(l\right)={\stackrel{←}{W}}_{S}\left(l\right)\cup {\stackrel{←}{W}}_{DF}\left(l\right)\cup {\stackrel{←}{W}}_{DT}\left(l\right)$ , $\stackrel{\to }{W}\left(r\right)={\stackrel{\to }{W}}_{S}\left(r\right)$ . When there exists only one intersection point of $\stackrel{←}{W}\left(l\right)$ and $\stackrel{\to }{W}\left(r\right)$ , we obtain the unique solution is a detonation wave solution $\stackrel{\to }{DT}+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ if ${p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma }$ , or $\stackrel{\to }{DT}+J+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ if ${p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }$ .

When there are three intersection points of $\stackrel{←}{W}\left(l\right)$ and ${\stackrel{\to }{W}}_{S}\left(r\right)$ (see Figure 3), from the modified global entropy condition A, we discard the intersection point of ${\stackrel{\to }{W}}_{S}\left(r\right)$ and ${\stackrel{←}{W}}_{DT}\left(l\right)$ . Denote the intersection point of ${\stackrel{\to }{W}}_{S}\left(r\right)$ and ${\stackrel{←}{W}}_{S}\left(l\right)$ by ${\ast }_{S}$ and the intersection point of ${\stackrel{\to }{W}}_{S}\left(r\right)$ and ${\stackrel{←}{W}}_{DF}\left(l\right)$ by ${\ast }_{DF}$ . Denote the temperature at the point ${\ast }_{S}$ , ${\ast }_{DF}$ on ${\stackrel{\to }{W}}_{S}\left(r\right)$ by ${T}_{S}$ , ${T}_{DF}$ , respectively. The temperature at ${\ast }_{DF}$ on ${\stackrel{←}{W}}_{DF}\left(l\right)$ is greater than ${T}_{i}$ since the combustion process is exothermic.

Subcase 1.1. ${p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma }$ .

Figure 3. ql > 0, qr = 0 and there are there interactions.

When ${T}_{r}\le {T}_{i}$ , then $\beta \left({\ast }_{S}\right)=0$ , $\beta \left({\ast }_{DF}\right)=2$ , from the condition B, we select ${\ast }_{S}$ and obtain a non-combustion wave solution $\stackrel{←}{R}$ or $\stackrel{←}{S}+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ (Figure 4(a)).

When ${T}_{r}>{T}_{i}$ , then $\beta \left({\ast }_{S}\right)=1$ , $\beta \left({\ast }_{DF}\right)=1$ , from the condition C, we select ${\ast }_{DF}$ and obtain a combustion wave solution $\stackrel{\to }{DT}+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ (Figure 4(b)).

Subcase 1.2. ${p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }$ .

From the condition A, we discard the possible detonation DT wave solution and find that the possible Riemann solution is $\stackrel{←}{R}$ or $\stackrel{←}{S}+J+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ or $\stackrel{\to }{DT}+J+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ . According to the modified global entropy conditions we obtain the unique Riemann solution as follows (see Figure 5).

1) When ${T}_{r}>{T}_{i}$ , ${T}_{DF}>{T}_{i}$ , then $\beta \left({\ast }_{S}\right)=1$ , $\beta \left({\ast }_{DF}\right)=1$ , from the condition C, we select ${\ast }_{DF}$ and obtain a combustion wave solution containing a DF (Figure 5(b)).

2) When ${T}_{r}>{T}_{i}$ , ${T}_{DF}\le {T}_{i}$ ( $⇒{T}_{S}\le {T}_{i}$ ), then $\beta \left({\ast }_{S}\right)=1$ , $\beta \left({\ast }_{DF}\right)=3$ , from the condition B, we select ${\ast }_{S}$ and obtain a non-combustion wave solution (Figure 5(a)).

3) When ${T}_{r}\le {T}_{i}$ , ${T}_{S}\le {T}_{i}$ , then $\beta \left({\ast }_{S}\right)=0$ , $\beta \left({\ast }_{DF}\right)=2$ , from the condition B, we select ${\ast }_{S}$ and obtain a non-combustion wave solution (Figure 5(a)).

4) When ${T}_{r}\le {T}_{i}$ , ${T}_{S}>{T}_{i}$ ( $⇒{T}_{DF}>{T}_{i}$ ), then $\beta \left({\ast }_{S}\right)=2$ , $\beta \left({\ast }_{DF}\right)=2$ , from the condition C, we select ${\ast }_{DF}$ and obtain a combustion wave solution containing a DF (Figure 5(b)).

Case 2. ${q}_{l}>0,{q}_{r}=0$ and there are two intersection points of $\stackrel{←}{W}\left(l\right)$ and ${\stackrel{\to }{W}}_{S}\left(r\right)$ (see Figure 6).

Subcase 2.1. ${p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma }$ .

In this case, we select the point ${\ast }_{S}$ or ${\ast }_{DT}$ and obtain the possible solutions $\stackrel{←}{S}+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ or $\stackrel{\to }{DT}+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ . Now we select the unique Riemann solution as follows.

When ${T}_{r}>{T}_{i}$ , then $\beta \left({\ast }_{S}\right)=1$ , $\beta \left({\ast }_{DT}\right)=1$ , from the condition C, we select ${\ast }_{DT}$ and obtain a combustion wave solution $\stackrel{\to }{DT}+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ (Figure 7(a)).

When ${T}_{r}\le {T}_{i}$ , then $\beta \left({\ast }_{S}\right)=0$ , $\beta \left({\ast }_{DT}\right)=2$ , from the condition B, we select ${\ast }_{S}$ and obtain a non-combustion wave solution $\stackrel{←}{S}+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ (Figure 7(b)).

Subcase 2.2. ${p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }$ .

There are two possibilities: one is that there exists only one intersection point of $\stackrel{←}{W}\left(l\right)$ and $\stackrel{\to }{W}\left(r\right)$ and we obtain the unique Riemann solution is $\stackrel{\to }{DT}+\stackrel{>}{J}+\stackrel{\to }{R}$

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 $\stackrel{\to }{S}$ , the other one is that there are three possible solutions which are the noncombution wave solution $\stackrel{←}{S}$ or $\stackrel{←}{R}+J+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ , or the DF combustion wave solution $\stackrel{\to }{DT}+J+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ , or the DT combustion wave solution $\stackrel{\to }{DT}+\stackrel{<}{J}+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ . 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 ${T}_{r}>{T}_{i}$ , ${T}_{DF}>{T}_{i}$ , then $\beta \left({\ast }_{S}\right)=1$ , $\beta \left({\ast }_{DF}\right)=1$ , from the condition C, we select ${\ast }_{DF}$ and obtain a combustion wave solution containing a DF (Figure 8(b)).

2) When ${T}_{r}>{T}_{i}$ , ${T}_{DF}\le {T}_{i}$ ( $⇒{T}_{S}\le {T}_{i}$ ), then $\beta \left({\ast }_{S}\right)=1$ , $\beta \left({\ast }_{DF}\right)=3$ , from the condition B, we select ${\ast }_{S}$ 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 ${T}_{r}\le {T}_{i}$ , ${T}_{S}\le {T}_{i}$ , then $\beta \left({\ast }_{S}\right)=0$ , $\beta \left({\ast }_{DF}\right)=2$ , from the condition B, we select ${\ast }_{S}$ and obtain a non-combustion wave solution (Figure 8(a)).

4) When ${T}_{r}\le {T}_{i}$ , ${T}_{S}>{T}_{i}$ ( $⇒{T}_{DF}>{T}_{i}$ ), then $\beta \left({\ast }_{S}\right)=2$ , $\beta \left({\ast }_{DF}\right)=2$ , from the condition C, we select ${\ast }_{DF}$ and obtain a combustion wave solution containing a DF (Figure 8(b)).

Case 3. ${q}_{l}>0,{q}_{r}>0$ and the gas on the both sides are unburnt. In this case ,we know that $\stackrel{←}{W}\left(l\right)={\stackrel{←}{W}}_{S}\left(l\right)\cup {\stackrel{←}{W}}_{DF}\left(l\right)\cup {\stackrel{←}{W}}_{DT}\left(l\right)$ , $\stackrel{\to }{W}\left(r\right)={\stackrel{\to }{W}}_{S}\left(r\right)\cup {\stackrel{\to }{W}}_{DF}\left(r\right)\cup {\stackrel{\to }{W}}_{DT}\left(r\right)$ .

If the intersection point of $\stackrel{←}{W}\left(l\right)$ and $\stackrel{\to }{W}\left(r\right)$ is unique, the solution is $\stackrel{\to }{DT}+\stackrel{←}{DT}$ if ${p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma }$ , or $\stackrel{\to }{DT}+J+\stackrel{←}{DT}$ if ${p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }$ . Otherwise, there are two possible subcases: one is that there is an intersection point of ${\stackrel{←}{W}}_{S}\left(l\right)$ and ${\stackrel{\to }{W}}_{S}\left(r\right)$ , the other is that there is no intersection point of ${\stackrel{←}{W}}_{S}\left(l\right)$ and ${\stackrel{\to }{W}}_{S}\left(r\right)$ .

Case 3.1. In the former subcase (see Figure 9), we discuss it in the following two subcases.

Subcase 3.1.1. ${p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma }$ .

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 $\beta =0$ for (a), and it holds that $\beta =2$ for (b), (c) and (d). From the condition B, we select the intersection point 1 and obtain the unique non-combustion wave solution $\stackrel{←}{R}$ or $\stackrel{←}{S}+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ .

Subcase 3.1.2. ${p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }$ . 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 $\stackrel{←}{R}$ or $\stackrel{←}{S}+J+\stackrel{\to }{R}$ or $\stackrel{\to }{S}$ . The only difference is that here the contact discontinuity appears.

Case 3.2. In the latter subcase, there are only two possibilities: $\stackrel{←}{W}\left(l\right)$ intersects ${\stackrel{\to }{W}}_{DT}\left(r\right)$ only or $\stackrel{\to }{W}\left(+\right)$ intersects ${\stackrel{←}{W}}_{DT}\left(l\right)$ only. We just need to consider the former. If the intersection point is unique, the solution is $\stackrel{\to }{DT}+\stackrel{←}{DT}$ if ${p}_{-}{\tau }_{-}^{\gamma }={p}_{+}{\tau }_{+}^{\gamma }$ , or $\stackrel{\to }{DT}+J+\stackrel{←}{DT}$ if ${p}_{-}{\tau }_{-}^{\gamma }\ne {p}_{+}{\tau }_{+}^{\gamma }$ , otherwise, there are at most three intersection points (see Figure 11).

Subcase 3.2.1. ${p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }$ .

From the condition A, the intersection point of ${\stackrel{←}{W}}_{DT}\left(l\right)$ and ${\stackrel{\to }{W}}_{DT}\left(r\right)$ should be discarded. We denote the intersection point of ${\stackrel{\to }{W}}_{DT}\left(r\right)$ and ${\stackrel{←}{W}}_{S}\left(l\right)$ by ${}^{\ast }S$ and denote the intersection point of ${\stackrel{\to }{W}}_{DT}\left(+\right)$ and ${\stackrel{←}{W}}_{DF}\left(l\right)$ by ${}^{\ast }DF$ ,

Figure 9. There is an intersection point of ${\stackrel{←}{W}}_{S}\left(l\right)$ and ${\stackrel{\to }{W}}_{S}\left(r\right)$ .

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. $\stackrel{←}{W}\left(-\right)$ intersects ${\stackrel{\to }{W}}_{DT}\left(r\right)$ only.

respectively. We denote the temperature at the point ${}^{\ast }S$ , ${}^{\ast }DF$ on ${\stackrel{\to }{W}}_{DT}\left(r\right)$ by ${T}_{S}$ , ${T}_{DF}$ , respectively (see Figure 12).

Since ${T}_{S}>{T}_{i}$ we have ${T}_{DF}>{T}_{i}$ , then $\beta \left({}^{\ast }S\right)=2$ , $\beta \left({}^{\ast }DF\right)=2$ . From the condition C, we select ${}^{\ast }DF$ and obtain a combustion wave solution $\stackrel{\to }{DF}+\stackrel{←}{DT}$ .

Subcase 3.2.2. ${p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma }$ . 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 $\stackrel{\to }{DF}+\stackrel{←}{DT}$ . 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).

Figure 12. Solutions in Subcase 3.2.1.

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

$\left(B,\rho ,u\right)\left(x,0\right)=\left\{\begin{array}{ll}\left({B}_{l},{\rho }_{l},{u}_{l}\right),\hfill & -\infty (12)

for arbitrary ${x}_{1},{x}_{2}\in 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 $\stackrel{\to }{SDT}$ or $\stackrel{\to }{WDF}$ and a contact discontinuity: $\stackrel{>}{J}\stackrel{\to }{SDT}$ , $\stackrel{<}{J}\stackrel{\to }{SDT}$ , $\stackrel{>}{J}\stackrel{\to }{WDF}$ , $\stackrel{<}{J}\stackrel{\to }{WDF}$ ;

The interaction of a combustion wave $\stackrel{\to }{SDT}$ or $\stackrel{\to }{WDF}$ and a shock wave: $\stackrel{\to }{SDT}\stackrel{\to }{S}$ , $\stackrel{\to }{WDF}\stackrel{\to }{S}$ .

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

Case 1. $\stackrel{>}{J}\stackrel{\to }{SDT}$

In the $\left(x,t\right)$ plane, $\stackrel{>}{J}$ and $\stackrel{\to }{SDT}$ 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 $\left(u,p\right)$ plane. It is obvious that ${u}_{l}={u}_{m}>{u}_{r}$ , ${p}_{l}>{p}_{m}$ and ${\tau }_{l}>{\tau }_{m}$ .

From the analysis of the wave curves in the $\left(u,p\right)$ plane, we know that $\stackrel{←}{W}\left(l\right)={\stackrel{←}{W}}_{S}\left(l\right)\cup {\stackrel{←}{W}}_{DF}\left(l\right)\cup {\stackrel{←}{W}}_{DT}\left(l\right)$ , and $\stackrel{\to }{W}\left(r\right)={\stackrel{\to }{W}}_{S}\left(r\right)$ .

There are two possible cases: there is only one intersection point of ${\stackrel{←}{W}}_{DT}\left(l\right)$ and ${\stackrel{\to }{W}}_{S}\left(r\right)$ , or there are three intersection points of $\stackrel{←}{W}\left(l\right)$ and ${\stackrel{\to }{W}}_{S}\left(r\right)$ . From the arguments in the former section, we know the result is as follows:

When ${p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma }$ , from the modified global entropy conditions we obtain that

$\stackrel{>}{J}\stackrel{\to }{SDT}\to \left\{\begin{array}{l}\stackrel{\to }{SDT}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{←}{S}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{←}{R}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{\to }{DF}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\end{array}$ (13)

where “ $\to$ ” means the result of the wave interaction.

Figure 13. The interaction of $\stackrel{>}{J}$ and $\stackrel{\to }{SDT}$ .

When ${p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }$ , the results are given by

$\stackrel{>}{J}\stackrel{\to }{SDT}\to \left\{\begin{array}{l}\stackrel{\to }{SDT}+J+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{←}{S}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{←}{R}+J+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{\to }{DF}+J+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S}.\end{array}$ (14)

Case 2. $\stackrel{<}{J}\stackrel{\to }{SDT}$

In this case (Figure 14), we know ${u}_{l}={u}_{m}>{u}_{r}$ , and for the contact discontinuity it holds that ${p}_{l}<{p}_{m}$ and ${\tau }_{l}<{\tau }_{m}$ . Similar with the discussions as the above case, we know that when ${p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma }$ ,

$\stackrel{<}{J}\stackrel{\to }{SDT}\to \left\{\begin{array}{l}\stackrel{←}{S}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{←}{R}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{\to }{SDT}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{\to }{DF}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\end{array}$ (15)

when ${p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }$ , 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. $\stackrel{>}{J}\stackrel{\to }{WDF}$

For this case, we know that $\stackrel{>}{J}$ and $\stackrel{\to }{WDF}$ will encounter each other at a finite time in the $\left(x,t\right)$ 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 $\left(u,p\right)$ plane.

Figure 14. The interaction of $\stackrel{<}{J}$ and $\stackrel{\to }{SDT}$ .

It is easily shown that ${u}_{l}={u}_{m}$ and ${p}_{l}>{p}_{m}$ , ${\tau }_{l}>{\tau }_{m}$ . We proceed as follows (Figure 15). Based on the analysis in the $\left(u,p\right)$ plane, we notice specially that the temperature is just the ignition point on the wave front of the $\stackrel{\to }{WDF}$ , i.e. the state $\left(m\right)=\left({i}_{m}\right)$ where the symbol “ $\left({i}_{m}\right)$ ” means the ignition point corresponding to the state (m). Since there are three intersection points of $\stackrel{←}{W}\left(l\right)$ and ${\stackrel{\to }{W}}_{S}\left(r\right)$ , from the global entropy conditions we know that when ${p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma }$ ,

$\stackrel{>}{J}\stackrel{\to }{WDF}\to \left\{\begin{array}{l}\stackrel{←}{S}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{←}{R}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{\to }{DF}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\end{array}$ (16)

when ${p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }$ ,

$\stackrel{>}{J}\stackrel{\to }{WDF}\to \left\{\begin{array}{l}\stackrel{←}{S}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{←}{R}+J+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{\to }{DF}+J+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S}.\end{array}$ (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. $\stackrel{<}{J}\stackrel{\to }{WDF}$

In this case (Figure 16), we know ${u}_{l}={u}_{m}$ , and for the contact discontinuity we know ${p}_{l}<{p}_{m}$ and ${\tau }_{l}<{\tau }_{m}$ . There are two possibilities: one is that there is only one intersection point of ${\stackrel{←}{W}}_{DT}\left(l\right)$ and ${\stackrel{\to }{W}}_{S}\left(r\right)$ , the other one is that there are three intersection points of $\stackrel{←}{W}\left(l\right)$ and ${\stackrel{\to }{W}}_{S}\left(r\right)$ . Therefore we know the results are as follows:

When ${p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma }$ , the results are as follows

$\stackrel{<}{J}\stackrel{\to }{DF}\to \left\{\begin{array}{l}\stackrel{\to }{DT}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{←}{S}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{←}{R}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{\to }{DF}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\end{array}$ (18)

when ${p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }$ , the results are the same as the above except that J appears here.

Figure 15. The interaction of $\stackrel{>}{J}$ and $\stackrel{\to }{WDF}$ .

Figure 16. The interaction of $\stackrel{<}{J}$ and $\stackrel{\to }{WDF}$ .

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. $\stackrel{\to }{SDT}\stackrel{\to }{S}$

In this case (Figure 17), we know for $\stackrel{\to }{SDT}$ it holds that ${\lambda }_{m}<{\sigma }_{1}<{\lambda }_{l}$ , and for $\stackrel{←}{S}$ it holds that ${\lambda }_{r}<{\sigma }_{2}<{\lambda }_{m}$ , where ${\sigma }_{1}$ and ${\sigma }_{2}$ are respectively the velocity of $\stackrel{\to }{SDT}$ and $\stackrel{←}{S}$ . Therefore, $\stackrel{←}{S}$ will overtake $\stackrel{\to }{SDT}$ at a finite time and a new Riemann problem is formed. We still solve this problem in the $\left(u,p\right)$ plane. $\left(m\right)\in {\stackrel{←}{W}}_{DT}\left(l\right)$ , $\left(r\right)\in \stackrel{←}{S}\left(m\right)$ and the curves ${\stackrel{←}{W}}_{S}\left(l\right)\cup {\stackrel{←}{W}}_{DF}\left(l\right)$ are located below the curve ${\stackrel{\to }{W}}_{S}\left(m\right)$ according to the condition A. From Lemma 2.6., we know ${\stackrel{←}{W}}_{S}\left(l\right)\cup {\stackrel{←}{W}}_{DF}\left(l\right)$ are located below $\stackrel{\to }{R}\left(m\right)\left(u\le {u}_{m}\right)$ and the half straight line $p={p}_{m}\left(u>{u}_{m}\right)$ . Considering that $\stackrel{\to }{R}\left(m\right)$ and $\stackrel{\to }{R}\left(r\right)$ from Lemma 2.6., we know that ${\stackrel{\to }{W}}_{S}\left(r\right)$ is located above the curve $\stackrel{\to }{R}\left(m\right)\left(u\le {u}_{m}\right)$ and the half

Figure 17. The interaction of $\stackrel{\to }{SDT}$ and $\stackrel{←}{S}$ .

straight line $p={p}_{m}\left(u>{u}_{m}\right)$ . Therefore, ${\stackrel{\to }{W}}_{S}\left(r\right)$ does not intersect with ${\stackrel{←}{W}}_{S}\left(l\right)\cup {\stackrel{←}{W}}_{DF}\left(l\right)$ and the result of the interaction is expressed by

$\stackrel{\to }{SDT}\stackrel{\to }{S}\to \left\{\begin{array}{l}\stackrel{\to }{SDT}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},{p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma },\\ \stackrel{\to }{SDT}+J+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},{p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }.\end{array}$ (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. $\stackrel{\to }{WDF}\stackrel{\to }{S}$

In this case (Figure 18), similar with the discussions as the just above case, we know that $\stackrel{←}{S}$ will overtake $\stackrel{\to }{SDT}$ 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 $\left(m\right)\in {\stackrel{←}{W}}_{DF}\left(l\right)$ , $\left(r\right)\in \stackrel{←}{S}\left(m\right)$ and (r) is located above the curve ${\stackrel{←}{W}}_{DF}\left(l\right)$ , we know the temperature at this point is higher than the ignition temperature, therefore when ${p}_{l}{\tau }_{l}^{\gamma }={p}_{r}{\tau }_{r}^{\gamma }$ , the results is described by

$\stackrel{\to }{WDF}\stackrel{\to }{S}\to \left\{\begin{array}{l}\stackrel{\to }{WDF}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\\ \stackrel{\to }{DT}+\stackrel{\to }{R}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\stackrel{\to }{S},\end{array}$ (20)

when ${p}_{l}{\tau }_{l}^{\gamma }\ne {p}_{r}{\tau }_{r}^{\gamma }$ , 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 $\stackrel{\to }{WDF}$ and $\stackrel{←}{S}$ .

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).

Conflicts of Interest

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

Cite 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

References

1. 1. Chen, G.Q. and Wang, D.H. (2002) Global Solutions of Nonlinear Magnetohydrodynamics with Large Initial Data. Journal of Differential Equations, 182, 344-376. https://doi.org/10.1006/jdeq.2001.4111

2. 2. Cowling, T.G. (1957) Magnetoydrodynamics. Interscience, New York.

3. 3. Helliwell, J.B. (1962) Gas-Ionizing Shock and Combustion Waves in Magnetogasdynamics. Journal of Fluid Mechanics, 14, 405-419. https://doi.org/10.1017/S0022112062001329

4. 4. Hu, Y.B. and Sheng, W.C. (2013) The Riemann Problem of Conservation Laws in Magnetogasdynamics. Communications on Pure and Applied Analysis, 12, 755-769. https://doi.org/10.3934/cpaa.2013.12.755

5. 5. Liu, Y.J. and Sun, W.H. (2019) The Riemann Problem for the Simplified Combustion Model in Magnetogasdynamics, Submitted.

6. 6. Mareev, V.A. (1975) Self-Similar Magnetogasdynamic Flows Accompanied by detonation and Combustion Waves. Fluid Dynamics, 10, 86-93. https://doi.org/10.1007/BF01023785

7. 7. Pu, X.K. and Guo, B.L. (2013) Global Existence and Convergence Rates of Smooth Solutions for the Full Compressible MHD Equations. Zeitschrift für Angewandte Mathematik und Physik, 64, 519-538. https://doi.org/10.1007/s00033-012-0245-5

8. 8. Sekhar, T.R. and Sharma, V.D. (2010) Riemann Problem and Elementary Wave Interactions in Isentropic Magnetogasdynamics. Nonlinear Analysis: Real World Applications, 11, 619-636. https://doi.org/10.1016/j.nonrwa.2008.10.036

9. 9. Shen, C. (2011) The Limits of Riemann Solutions to the Isentropic Magnetogasdynamics. Applied Mathematics Letters, 24, 1124-1129. https://doi.org/10.1016/j.aml.2011.01.038

10. 10. Courant, R. and Friedrichs, K.O. (1948) Supersonic Flow and Shock Waves. Interscience, New York.

11. 11. Liu, Y.J. and Sheng, W.C. (2011) The Generalized Riemann Problem for Gas Dynamic Combustion. Applied Mathematics and Mechanics (English Series), 32, 1079-1090. https://doi.org/10.1007/s10483-011-1482-6

12. 12. Sun, M.N. (2017) Singular Solutions to the Riemann Problem for a Macroscopic Production MODEL. Zeitschrift fur Angewandte Mathematik und Mechanik, 97, 916-931. https://doi.org/10.1002/zamm.201600171

13. 13. Zhang, T. and Zheng, Y.X. (1989) Riemann Problem for Gasdynamic Combustion. Journal of Differential Equations, 77, 203-230. https://doi.org/10.1016/0022-0396(89)90142-3