Abstract:As described in the classic works of Lee--Stewart and Short--Stewart, the numerical evaluation of linear stability of planar detonation waves is a computationally intensive problem of considerable interest in applications. Reexamining this problem from a modern numerical Evans function point of view, we derive a new algorithm for their stability analysis, related to a much older method of Erpenbeck, that, while equally simple and easy to implement as the standard method introduced by Lee--Stewart, appears to be potentially faster and more stable.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: how to perform the numerical stability analysis of ZND (Zeldovich–von Neumann–Döring) detonation waves more efficiently and stably. Specifically, the author re - examines the classical Lee - Stewart algorithm and derives a new algorithm to evaluate the stability of ZND detonation waves from the perspective of the modern numerical Evans function.
### Background and Problem Description
1. **Computational Complexity**:
- The linear stability analysis of ZND detonation waves is a highly computationally intensive problem. In particular, when physical parameters change, the stability of the detonation wave will change, resulting in phenomena such as pulsation, rotation, and cellular - like solutions.
- Due to the number of physical parameters (for example, there are four parameters for polyatomic gases) and the difficulty of a single calculation, this problem has proven to be very numerically intensive.
2. **Limitations of Existing Methods**:
- In their classic 1990 paper, Lee and Stewart introduced a standard algorithm that for the first time accurately calculated the stability boundary of one - dimensional detonation waves. Although this method is straightforward and easy to implement, a comprehensive survey of all regions of the parameter space is still computationally very expensive.
- Although computing power has increased significantly in the past few decades, making previously difficult calculations now feasible, it is still meaningful to look for more efficient algorithms.
### Advantages of the New Algorithm
1. **Direction Selection**:
- The new algorithm integrates from \(x = -\infty\) to \(x = 0\), rather than from \(x = 0\) to \(x = -\infty\). The advantage of this is that we are seeking a generalized eigenfunction that decays exponentially at \(x = -\infty\). Therefore, in the forward direction (\(-\infty\to 0\)), the desired solution will grow exponentially, while the error mode will be exponentially suppressed.
- In contrast, when integrating in the reverse direction (\(0\to -\infty\)), the desired solution will decay exponentially, while the error mode will be exponentially amplified, which is numerically undesirable.
2. **Computational Efficiency**:
- Through this improvement, the new algorithm is numerically more stable and faster, thus reducing the computation time and making interactive numerical exploration possible.
- Subsequent work has shown that, compared with the existing standard methods, the new algorithm performs excellently in practical applications, especially in cases of high activation energy or the "square - wave" limit.
### Conclusion
This paper proposes a new numerical algorithm for more efficient and stable analysis of the stability of ZND detonation waves. By changing the integration direction and using the Evans function method, the new algorithm is more numerically stable and has higher computational efficiency, providing a powerful tool for in - depth understanding of detonation phenomena.