A numerical stability investigation of strong ZND detonations for Majda's model
Abstract
We carry out a systematic numerical stability analysis of ZND detonations of Majda's model with Arrhenius-type ignition function, a simplified model for reacting flow, as heat release and activation energy are varied. Our purpose is, first, to answer a question of Majda whether oscillatory instabilities can occur for high activation energies as in the full reacting Euler equations, and, second, to test the efficiency of various versions of a numerical eigenvalue-finding scheme suggested by Humpherys and Zumbrun against the standard method of Lee and Stewart. Our results suggest that instabilities do not occur for Majda's model with Arrhenius-type ignition function, nor with a modified Arrhenius-type ignition function suggested by Lyng--Zumbrun, even in the high-activation energy limit. We find that the algorithm of Humpherys--Zumbrun is in the context of Majda's model - times faster than the one described in the classical work of Lee and Stewart and - times faster than an optimized version of the Lee--Stewart algorithm using an adaptive-mesh ODE solver
Keywords
Cite
@article{arxiv.1011.1561,
title = {A numerical stability investigation of strong ZND detonations for Majda's model},
author = {Blake Barker and Kevin zumbrun},
journal= {arXiv preprint arXiv:1011.1561},
year = {2010}
}