Implementation of an OpenFOAM solver for shock and detonation simulation at high pressure

The high pressure condition widely used in industry makes it necessary to evaluate the impact of real gas behaviour in numerical simulation, especially for shock and detonation problems. Since most open source solvers are based on ideal or perfect gas assumption, we have developed an OpenFOAM solver for Real gas Shock and Detonation simulation, referred to as RSDFoam. We have considered the nonidealities in the equation of state, thermodynamic functions and mass action law brought about by the interaction between the fluid particles and by their finite volume. Connection between blastFoam and Cantera was built to utilize the CFD capabilities from the former and the efficient chemistry solver and real gas models in Cantera. The new solver has been thoroughly validated against analytical solutions, previous numerical simulation results and experimental data. For non-reactive flow, shock tube and oblique shock problems were solved. For reactive cases, constant volume reactor, steady detonation speed, reaction zone structure, and cellular structure were numerically studied. The satisfactory agreement demonstrates the accuracy and robustness of the solver for real gas based shock and detonation simulation.