Convergence properties of detonation simulations

• 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.

Ämnesord

NATURVETENSKAP  -- Fysik -- Astronomi, astrofysik och kosmologi (hsv//swe)

NATURAL SCIENCES  -- Physical Sciences -- Astronomy, Astrophysics and Cosmology (hsv//eng)

NATURVETENSKAP  -- Fysik (hsv//swe)

NATURAL SCIENCES  -- Physical Sciences (hsv//eng)

Nyckelord

Combustion

numerical methods

detonation

shock waves

chemical reaction

Publikations- och innehållstyp

ref (ämneskategori)

art (ämneskategori)