Here we analyze properties of an equation that we previously proposed to model the dynamics of unstable detonation waves [A. R. Kasimov, L. M. Faria, and R. R. Rosales, Model for shock wave chaos, Phys. Rev. Lett., 110 (2013), 104104]. The equation is utIt describes a detonation shock at x = 0 with the reaction zone in x < 0. We investigate the nature of the steady-state solutions of this nonlocal hyperbolic balance law, the linear stability of these solutions, and the nonlinear dynamics. We establish the existence of instability followed by a cascade of period-doubling bifurcations leading to chaos.
Introduction.A detonation is a shock wave that propagates in a reactive medium where exothermic chemical reactions are ignited as a result of the heating by the shock compression. The energy released in these reactions, in turn, feeds back to the shock in the form of compression waves and thus sustains the shock motion. The dynamics of such shock-reaction coupling is highly nonlinear due to the sensitivity of the chemical reactions to temperature, making the problem significantly more challenging than shock dynamics in nonreactive media. A steady planar detonation wave is rarely observed in experiments. Complex time-dependent and multidimensional structures tend to develop [14,26]. Numerical simulations of the equations of reactive gas dynamics are able to reproduce at a qualitative level the complex structures observed in experiments (see, e.g., [40,1,30]). However, obtaining physical insights into the basic mechanisms of the instability requires simplified modeling and remains challenging.In one dimension, the instabilities of the reactive shock wave manifest themselves in the form of a "galloping detonation" [15,14], wherein the shock speed oscillates around its steady value. It has been shown through extensive numerical experiments that as the activation energy, E, a parameter in the equations measuring the temperature sensitivity of the chemical reactions, is varied, the shock speed transitions from a constant to an oscillatory function. Further increase of E leads to a period-doubling bifurcation cascade, which ultimately results in the shock moving at a chaotic speed [29,20]. The mechanism for such instabilities is still not completely understood.In this paper, we show that the model introduced in [23], which consists of a single nonlocal partial differential equation (PDE), is capable of reproducing the complexity observed in one-dimensional simulations of reactive Euler equations. The model possesses traveling wave solutions precisely analogous to the ZND theory (named after Zel'dovich [44], von Neumann [41], and Döring [6], who independently developed the