Simulations of pulsating one-dimensional detonations with true fifth order accuracy

  • Authors:

  • Andrew K. Henrick;Tariq D. Aslam;Joseph M. Powers

  • Affiliations:

  • Department of Aerospace and Mechanical Engineering, University of Notre Dame, Notre Dame, IN 46556, United States;Dynamic Experimentation Division, Los Alamos National Laboratory, Group DX-2, MS P952, Los Alamos, NM 87545, United States;Department of Aerospace and Mechanical Engineering, University of Notre Dame, Notre Dame, IN 46556, United States

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2006

Quantified Score

Hi-index 31.45

Visualization

Abstract

A novel, highly accurate numerical scheme based on shock-fitting coupled with fifth order spatial and temporal discretizations is applied to a classical unsteady detonation problem to generate solutions with unprecedented accuracy. The one-dimensional reactive Euler equations for a calorically perfect mixture of ideal gases whose reaction is described by single-step irreversible Arrhenius kinetics are solved in a series of calculations in which the activation energy is varied. In contrast with nearly all known simulations of this problem, which converge at a rate no greater than first order as the spatial and temporal grid is refined, the present method is shown to converge at a rate consistent with the fifth order accuracy of the spatial and temporal discretization schemes. This high accuracy enables more precise verification of known results and prediction of heretofore unknown phenomena. To five significant figures, the scheme faithfully recovers the stability boundary, growth rates, and wave-numbers predicted by an independent linear stability theory in the stable and weakly unstable regime. As the activation energy is increased, a series of period-doubling events are predicted, and the system undergoes a transition to chaos. Consistent with general theories of non-linear dynamics, the bifurcation points are seen to converge at a rate for which the Feigenbaum constant is 4.66+/-0.09, in close agreement with the true value of 4.669201.... As activation energy is increased further, domains are identified in which the system undergoes a transition from a chaotic state back to one whose limit cycles are characterized by a small number of non-linear oscillatory modes. This result is consistent with behavior of other non-linear dynamical systems, but not typically considered in detonation dynamics. The period and average detonation velocity are calculated for a variety of asymptotically stable limit cycles. The average velocity for such pulsating detonations is found to be slightly greater than the Chapman-Jouguet velocity.