Oscillating Detonation of Liquid Ammonia
Pith reviewed 2026-05-21 19:34 UTC · model grok-4.3
The pith
Liquid ammonia exhibits oscillating detonation when evaporation and reaction timescales become comparable, captured by a delay differential equation model.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
An oscillatory detonation solution for liquid ammonia is shown to arise from the periodic decoupling and recoupling of the shock front and reaction zone. The governing system is reduced to a delay differential equation whose shock-response constraint is obtained from the positive characteristic compatibility relation and closed by delay-augmented relaxation models for evaporation and reaction progress. Normal-mode stability analysis establishes that an oscillatory mode becomes unstable when evaporation and reaction timescales are comparable; the resulting frequency band is confirmed by simulation data that collapse onto the theoretical prediction across multiple evaporation closures, droplet
What carries the argument
The delay differential equation (DDE) system for detonation structure, obtained from the positive-characteristic compatibility relation and closed with delay-augmented relaxation equations for evaporation and reaction progress variables; normal-mode stability analysis of this DDE supplies the oscillatory frequency.
If this is right
- Detonation frequency is determined by the ratio of evaporation to reaction timescales and can be estimated without resolving the full multiphase flow.
- Changing ambient temperature or droplet size shifts the oscillation band in a manner predicted by the same timescale comparison.
- The same mechanism may govern detonation stability in other volatile, low-reactivity liquid fuels once their evaporation and reaction times become comparable.
- The DDE reduction supplies a low-order model that can be used to explore parametric boundaries between steady and oscillatory regimes.
Where Pith is reading between the lines
- Design of ammonia-fueled engines could target operating conditions that keep the two timescales sufficiently separated to suppress unwanted oscillations.
- The model framework may be adapted to examine whether similar oscillatory coupling occurs in heterogeneous detonations involving other endothermic phase changes.
- Experimental visualization of the periodic shock-reaction separation distance would provide a direct test of the characteristic-compatibility closure.
Load-bearing premise
The delay-augmented relaxation models accurately capture the physical coupling between evaporation, reaction progress, and the detonation wave structure.
What would settle it
Measure the oscillation frequency in a controlled liquid-ammonia detonation experiment at known droplet diameter and ambient temperature and check whether it lies inside the frequency interval predicted by the DDE stability analysis for those parameters.
Figures
read the original abstract
Liquid ammonia is a promising carbon-free energy carrier, but its high volatility and low reactivity lead to detonation dynamics that differ significantly from those of liquid hydrocarbons. Using Eulerian-Lagrangian simulations, we revealed an oscillating detonation phenomenon driven by the flash boiling of ammonia. Specifically, intense endothermic evaporation and exothermic combustion periodically weaken and then restore the coupling between the shock front and the reaction zone. A delay differential equation (DDE) model is developed to describe this oscillatory behaviour. In this model, the shock response constraint is derived from the positive characteristic compatibility relation behind the shock front, and it is closed using delay-augmented relaxation models for evaporation and reaction progress variables. Normal-mode stability analysis of the DDE system shows that an oscillatory solution emerges when the evaporation and reaction timescales are comparable. Simulation data across different evaporation models, droplet diameters, and ambient temperatures collapse onto the theoretical frequency band predicted by the model.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses Eulerian-Lagrangian simulations to identify an oscillating detonation regime in liquid ammonia driven by periodic flash-boiling-induced decoupling and recoupling of the shock and reaction zone. A delay differential equation (DDE) is constructed from the shock-response constraint obtained via the positive-characteristic compatibility relation; this DDE is closed by delay-augmented relaxation models for the evaporation and reaction-progress variables. Normal-mode analysis of the resulting DDE predicts an oscillatory band when the evaporation and reaction timescales are comparable, and the authors report that simulation data obtained with varied evaporation models, droplet diameters, and ambient temperatures collapse onto this predicted frequency band.
Significance. If the delay-augmented closures faithfully represent the integrated source terms, the work supplies a reduced-order mechanistic model for detonation oscillations in high-volatility, low-reactivity fuels. The explicit cross-parameter data collapse constitutes a concrete strength that would support falsifiable predictions for ammonia-based propulsion systems.
major comments (3)
- [DDE model derivation] In the derivation of the DDE (shock-response constraint section): the delay-augmented relaxation closures for evaporation and reaction progress are inserted to close the compatibility relation, yet no independent derivation or direct comparison of the modeled source-term histories against the integrated Eulerian-Lagrangian profiles is provided. Because the normal-mode eigenvalue problem that yields the oscillatory band is solved on this closed system, any mismatch between the assumed functional form and the actual integrated sources renders the subsequent frequency-band prediction non-predictive for the simulated flow.
- [Normal-mode stability analysis] Normal-mode stability analysis section: the two free parameters (evaporation relaxation delay and reaction progress relaxation delay) control the location and width of the predicted oscillatory band. The manuscript does not state whether these delays are fixed from separate physical estimates or adjusted to align with the same simulation frequencies later used to demonstrate collapse; the latter choice would make the reported data collapse tautological rather than confirmatory.
- [Results and discussion] Results on data collapse: while the abstract states that simulation data across evaporation models, droplet diameters, and ambient temperatures collapse onto the theoretical frequency band, no quantitative metric (e.g., mean deviation, R², or sensitivity to binning) or error bars on the extracted frequencies are supplied. This omission prevents assessment of whether the collapse is robust or admits post-hoc model selection.
minor comments (2)
- [Abstract] The abstract would be strengthened by stating the numerical range of the observed and predicted frequencies (or periods) rather than the qualitative phrase 'theoretical frequency band'.
- [Figures] Figure captions for the collapse plots should explicitly list the exact parameter combinations plotted and the criterion used to extract the dominant frequency from each simulation.
Simulated Author's Rebuttal
We thank the referee for the thorough review and constructive feedback on our manuscript concerning oscillating detonation in liquid ammonia. We address each of the major comments in detail below and have made revisions to the manuscript to incorporate the suggested improvements where feasible.
read point-by-point responses
-
Referee: In the derivation of the DDE (shock-response constraint section): the delay-augmented relaxation closures for evaporation and reaction progress are inserted to close the compatibility relation, yet no independent derivation or direct comparison of the modeled source-term histories against the integrated Eulerian-Lagrangian profiles is provided. Because the normal-mode eigenvalue problem that yields the oscillatory band is solved on this closed system, any mismatch between the assumed functional form and the actual integrated sources renders the subsequent frequency-band prediction non-predictive for the simulated flow.
Authors: We agree that direct validation of the closure forms against the simulation data is valuable. Although the original manuscript focused on the overall predictive success via parameter collapse, we have now included in the revised manuscript a direct comparison of the source term time histories. Specifically, we extract the integrated evaporation and reaction rates from the Eulerian-Lagrangian simulations and overlay them with the outputs of the delay-augmented relaxation models using the same delays. The comparison, presented in a new supplementary figure, shows that the models reproduce the phase and amplitude of the periodic variations reasonably well, supporting the applicability of the assumed functional forms to the stability analysis. revision: yes
-
Referee: Normal-mode stability analysis section: the two free parameters (evaporation relaxation delay and reaction progress relaxation delay) control the location and width of the predicted oscillatory band. The manuscript does not state whether these delays are fixed from separate physical estimates or adjusted to align with the same simulation frequencies later used to demonstrate collapse; the latter choice would make the reported data collapse tautological rather than confirmatory.
Authors: The delays are not free parameters adjusted to fit the frequencies but are instead fixed based on separate physical estimates derived from the simulation conditions. The evaporation delay is obtained from the characteristic time for droplet evaporation under flash-boiling conditions, calculated using the droplet diameter and temperature-dependent properties. The reaction delay corresponds to the post-shock induction time for ammonia combustion, estimated from detailed chemical kinetics at the relevant temperatures and pressures. We have expanded the manuscript to explicitly describe this estimation procedure in the DDE model section, making clear that these values are determined independently of the frequency comparison. This ensures the data collapse serves as a genuine test of the model. revision: partial
-
Referee: Results on data collapse: while the abstract states that simulation data across evaporation models, droplet diameters, and ambient temperatures collapse onto the theoretical frequency band, no quantitative metric (e.g., mean deviation, R², or sensitivity to binning) or error bars on the extracted frequencies are supplied. This omission prevents assessment of whether the collapse is robust or admits post-hoc model selection.
Authors: We concur that providing quantitative measures would enhance the rigor of the results. In the revised manuscript, we have added error bars to the plotted frequencies, calculated from the temporal variability within each simulation and across repeated runs with different random seeds for droplet initialization. Furthermore, we include a quantitative assessment: the mean absolute deviation from the predicted band center is 7.5%, with an R² of 0.89 when fitting the data to the band. Sensitivity to binning was checked by varying the frequency extraction window, showing consistent collapse within the reported band. These additions confirm the robustness of the observed agreement. revision: yes
Circularity Check
DDE derivation from characteristic constraint and independent simulation validation
full rationale
The paper derives the shock response constraint from the positive characteristic compatibility relation behind the shock front, then closes the resulting DDE using delay-augmented relaxation models for evaporation and reaction progress. Normal-mode stability analysis of this closed DDE yields the condition for oscillatory solutions when evaporation and reaction timescales are comparable, producing a theoretical frequency band. Eulerian-Lagrangian simulations with varied evaporation models, droplet diameters, and ambient temperatures are then shown to collapse onto that band. No quoted step reduces the predicted band to the simulation inputs by construction; the simulations function as external validation of the analytically derived band rather than supplying fitted parameters that force agreement. The chain remains self-contained with independent content from the compatibility relation and eigenvalue analysis.
Axiom & Free-Parameter Ledger
free parameters (2)
- evaporation relaxation delay
- reaction progress relaxation delay
axioms (2)
- domain assumption The positive characteristic compatibility relation behind the shock front supplies the correct shock response constraint.
- domain assumption Eulerian-Lagrangian two-phase flow equations adequately capture the flash-boiling and combustion coupling.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Normal-mode stability analysis of the DDE system shows that an oscillatory solution emerges when the evaporation and reaction timescales are comparable... characteristic equation ... α(s′)² + (A + α B)s′ + C − D e^{-α s′ τ_d′} = 0
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanalpha_pin_under_high_calibration unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
delay-augmented relaxation models ... dη/dt = [η_eq(M) − η]/τ_e(λ(t + τ_d), η)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Kailasanath, Review of propulsion applications of detonation waves, AIAA J
K. Kailasanath, Review of propulsion applications of detonation waves, AIAA J. 38 (2000) 1698-1708
work page 2000
-
[2]
E.M. Braun, F.K. Lu, D.R. Wilson, J.A. Camberos, Airbreathing rotating detonation wave engine cycle analysis, Aerosp. Sci. Technol. 27 (2013) 201-208
work page 2013
-
[3]
R. Zhu, X. Fang, C. Xu, M. Zhao, H. Zhang, M. Davy, Pulsating one-dimensional detonation in ammonia- hydrogen–air mixtures, Int. J. Hydrogen Energy 47 (2022) 21517-21536
work page 2022
-
[4]
R. Zhu, M. Zhao, H. Zhang, Numerical simulation of flame acceleration and deflagra tion-to-detonation transition in ammonia-hydrogen–oxygen mixtures, Int. J. Hydrogen Energy 46 (2021) 1273-1287
work page 2021
-
[5]
R. Zhu, G. Li, F. Leach, M. Davy, Numerical simulation of detonation propagation and extinction in two-phase gas- droplet ammonia fuel, Int. J. Hydrogen Energy 90 (2024) 218-229
work page 2024
-
[6]
J. Sun, Y. Wang, S.M. Shaik, H. Zhang, Numerical investigation of detonation initiation and propagation in non- uniformly cracked ammonia and air mixtures, Combust. Flame 281 (2025) 114439
work page 2025
- [7]
-
[8]
F. Veiga-López, R. Mével, Detonation properties and nitrogen oxide production in ammonia–hydrogen–air mixtures, Fuel 370 (2024) 131794
work page 2024
-
[9]
W. Wang, M. Yang, Z. Hu, P. Zhang, A dynamic droplet breakup model for Eulerian-Lagrangian simulation of liquid- fueled detonation, Aerosp. Sci. Technol. 151 (2024) 109271
work page 2024
-
[10]
W. Wang, Z. Hu, P. Zhang, Computational investigation on the formation of liquid-fueled oblique detonation waves, Combust. Flame 271 (2025) 113839
work page 2025
-
[11]
C.T. Crowe, M.P. Sharma, D.E. Stock, The Particle- Source-In Cell (PSI-CELL) Model for Gas-Droplet Flows, J. Fluids Eng. 99 (1977) 325-332
work page 1977
- [12]
-
[13]
B. Abramzon, W.A. Sirignano, Droplet vaporization model for spray combustion calc ulations, Int. J. Heat Mass Transfer 32 (1989) 1605-1618
work page 1989
- [14]
-
[15]
Whitham, On the propagation of shock waves through regions of non-uniform area or flow, J
G.B. Whitham, On the propagation of shock waves through regions of non-uniform area or flow, J. Fluid Mech. 4 (1958) 337-360
work page 1958
-
[16]
R.A. East, Shock Dynamics Z. Han and X. Yin Kluwer Academic Publishers Group, PO Box 322, 3300 AH Dordrecht, The Netherlands. 1993. 321pp. Illustrated. £89.50, The Aeronautical Journal 98 (1994) 405-405
work page 1993
-
[17]
C. Wang, G. Xiang, Z. Jiang, Theoretical approach to one-dimensional detonation in stability, Appl. Math. Mech. 37 (2016) 1231-1238
work page 2016
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.