Development and verification of a high-speed compressible reactive flow solver in OpenFOAM

Abstract

This paper proposes a newly developed density-based compressible solver called DCRFoam on a basis of OpenFOAM platform for simulating detonation phenomena with detailed chemical kinetics. It employs AUSM+-up scheme to determine the convective fluxes and a four steps low storage Runge-Kutta time integration method for time integration. A total variation diminishing (TVD) and symmetrical flux limiter sgva is implemented for interpolation process. Also, the dynamic mesh adaptation method is applied to refine or coarsen the computational mesh dynamically since the resolution requirements of supersonic combustion are localized. The numerical solver has been validated systematically with many numerical and experimental results reported in the literature, including those from non-reacting cases (Sod’s shock tube problem, two-dimensional shock/cylinder interaction) to reacting flow cases (the ignition sequence in a shock tube, oblique detonation wave (ODW) initiating process and detonation propagating in a bifurcated channel). The results obtained by current schemes are compared with original central upwind Kurganov-Noelle-Petrova (KNP) scheme in official OpenFOAM. The corresponding result demonstrates the robustness and accuracy of the current solver.