Set-valued solutions for non-ideal detonation

The existence and structure of a steady-state gaseous detonation propagating in a packed bed of solid inert particles are analyzed in the one-dimensional approximation by taking into consideration frictional and heat losses between the gas and the particles. A new formulation of the governing equations is introduced that eliminates the difficulties with numerical integration across the sonic singularity in the reactive Euler equations. With the new algorithm, we find that when the sonic point disappears from the flow, there exists a one-parameter family of solutions parameterized by either pressure or temperature at the end of the reaction zone. These solutions (termed "set-valued" here) correspond to a continuous spectrum of the eigenvalue problem that determines the detonation velocity as a function of a loss factor.