An adaptive numerical scheme for high-speed reactive flow on overlapping grids

  • Authors:
  • William D. Henshaw;Donald W. Schwendeman

  • Affiliations:
  • Centre for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore, CA;Department of Mathematical Sciences, Rensselaer Polytechnic Institute, Troy, NY

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

Quantified Score

Hi-index 31.52

Visualization

Abstract

We describe a method for the numerical solution of high-speed reactive flow in complex geometries using overlapping grids and block-structured adaptive mesh refinement. We consider flows described by the reactive Euler equations with an ideal equation of state and various stiff reaction models. These equations are solved using a second-order accurate Godunov method for the convective fluxes and a Runge-Kutta time-stepping scheme for the source term modeling the chemical reactions. We describe an extension of the adaptive mesh refinement approach to curvilinear overlapping grids. Numerical results are presented showing the evolution to detonation in a quarter-plane provoked by a temperature gradient and the propagation of an overdriven detonation in an expanding channel. The first problem, which considers a one-step Arrhenius reaction model, is used primarily to validate the numerical method, while the second problem, which considers a three-step chain-branching reaction model, is used to illustrate mechanisms of detonation failure and rebirth for the channel geometry.