A numerical stability investigation of strong ZND detonations for Majda's model
arXiv:1011.1561
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
References in corpus (8)
- An efficient shooting algorithm for Evans function calculations in large systems
- Stability of strong ideal-gas shock layers
- A stability index for detonation waves in Majda's model for reacting flow
- On the One-dimensional Stability of Viscous Strong Detonation Waves
- One-dimensional stability of parallel shock layers in isentropic magnetohydrodynamic
- Numerical error analysis for Evans function computations: a numerical gap lemma, centered-coordinate methods, and the unreasonable effectiveness of continuous orthogonalization
- Existence and stability of viscoelastic shock profiles
- Efficient numerical stability analysis of detonation waves in ZND