rhoCentralRfFoam: An OpenFOAM solver for high speed chemically active flows – Simulation of planar detonations

Abstract

A new solver developed within the framework of OpenFOAM 2.3.0, called rhoCentralRfFoam which can be interpreted like an evolution of rhoCentralFoam, is presented. Its use, performing numerical simulations on initiation and propagation of planar detonation waves in combustible mixtures H2–Air and H2–O2–Ar, is described. Unsteady one dimensional (1D) Euler equations coupled with sources to take into account chemical activity, are numerically solved using the Kurganov, Noelle and Petrova second order scheme in a domain discretized with finite volumes. The computational code can work with any number of species and its corresponding reactions, but here it was tested with 13 chemically active species (one species inert), and 33 elementary reactions. A gaseous igniter which acts like a shock-tube driver, and powerful enough to generate a strong shock capable of triggering exothermic chemical reactions in fuel mixtures, is used to start planar detonations. The following main aspects of planar detonations are here, treated: induction time of combustible mixtures cited above and required mesh resolutions; convergence of overdriven detonations to Chapman–Jouguet states; detonation structure (ZND model); and the use of reflected shocks to determine induction times experimentally. The rhoCentralRfFoam code was verified comparing numerical results and it was validated, through analytical results and experimental data.