Convergence properties of detonation simulations
Abstract
We present a high-resolution convergence study of detonation initiated by a temperature gradient in a stoichiometric hydrogen-oxygen mixture using the PENCIL CODE and compare with a code that employs a fifth order weighted essentially non-oscillating (WENO) scheme. With Mach numbers reaching 10-30, a certain amount of shock viscosity is needed in the PENCIL CODE to remove or reduce numerical pressure oscillations on the grid scale at the position of the shock. Detonation is found to occur for intermediate values of the shock viscosity parameter. At fixed values of this parameter, the numerical error associated with those small wiggles in the pressure profile is found to decrease with decreasing mesh width ? like ? down to ?. With the WENO scheme, solutions are smooth at ?, but no detonation is obtained for ?. This is argued to be an artifact of a decoupling between pressure and reaction fronts.
- Publication:
-
Geophysical and Astrophysical Fluid Dynamics
- Pub Date:
- March 2020
- DOI:
- 10.1080/03091929.2019.1668382
- arXiv:
- arXiv:1902.03816
- Bibcode:
- 2020GApFD.114...58Q
- Keywords:
-
- Combustion;
- numerical methods;
- detonation;
- shock waves;
- chemical reaction;
- Physics - Fluid Dynamics;
- Physics - Computational Physics
- E-Print:
- 18 pages, 10 figures, 1 table, Geophys. Astrophys. Fluid Dynam., in press, special issue "Physics and Algorithms of the Pencil Code"