One-dimensional hydrogen-oxygen detonation problem #7
Description
I attempted to reproduce the one-dimensional hydrogen-oxygen detonation problem described in the manual (stoichiometric H₂/O₂ mixture, initial temperature of 295 K, initial pressure of 0.2 bar). The detonation is initiated by applying an initial condition of 3000 K and 10 bar over a small region near the closed end. I used the Evans 7-species, 8-reaction model, as in the reference paper, but the simulation kept diverging. The time step was set to 1e-12 μs, with a grid size of 5 μm. I applied filtering in the x-direction twice (once every 4 time steps and once every 20 time steps). The boundary condition is set as non-reflective. To stabilize the simulation, 20 grid points are used in the y and z directions. The numerical schemes used are WENO5 and CD6. Switching to Li's 9-species, 19-reaction mechanism also resulted in divergence. If anyone have any suggestions for improving this case?