Supersonic Autoignitive Reaction Waves

• Supersonic autoignitive reaction waves are compressible reactive structures where finite-rate autoignition coupled with inlet velocity produces shock-free, steady propagation.
• The formulation bridges classical detonation and deflagration by integrating chemical kinetics in the preheat zone, allowing a continuous range of steady solutions.
• Reduced models, including reactive Burgers and Navier–Stokes equations, reveal key parameters such as activation energy and induction delay that determine flow regimes.

Searching arXiv for the cited papers and closely related work to ground the article in the current literature. Supersonic autoignitive reaction waves are compressible reactive-wave structures in which finite-rate autoignition, residence time, and wave propagation remain coupled under supersonic flow conditions. In the modern formulation, the front is governed by imposed inflow velocity and chemical autoignition kinetics rather than solely by diffusive flame propagation or by the shock-supported Chapman–Jouguet construction of classical detonation theory. In simplified analog models, the same class of phenomena appears during shock-induced ignition transients, where energy release amplifies forward-travelling compression waves and can generate internally supersonic reaction fronts. Across these formulations, the defining issues are the sensitivity of induction delay to compression, the relation between ignition length and hydrodynamic residence time, and the possibility of shock-free or weak-detonation-like supersonic propagation (Morii et al., 2023, Tang et al., 2012, Morii et al., 14 Jul 2025).

1. Concept and relation to deflagration and detonation

In a one-dimensional reactive-flow system with unburned premixed gas entering from the inlet boundary and burned gas exiting from the outlet boundary, the classical viewpoint is that steady-state solutions exist only when the inlet velocity matches either the velocity of a deflagration wave, determined by the burning-rate eigenvalue in the subsonic regime, or the velocity of a detonation wave, dictated by the Chapman–Jouguet condition in the supersonic regime. The general concept of the autoignitive reaction wave broadens that picture by including chemical reaction in the preheat zone, so that ignition occurs at a distance set approximately by the product of inlet velocity and ignition delay, $x_{ig}\approx u_{in}\tau_{ig}$, and steady solutions can arise over a continuous range of inlet velocities in both subsonic and supersonic bands (Morii et al., 2023).

Within this framework, a subsonic solution slightly above $S_L$ is described as a “fast deflagration” or autoignitive flame, whereas a supersonic solution slightly above $D_{CJ}$ is the “supersonic autoignitive wave.” The 2023 theory emphasizes that when the inlet velocity exceeds $D_{CJ}$, a distinct steady-state solution without any leading shock is possible, and that this wave is devoid of the typical shock wave commonly seen in detonation waves and lacks the detonation cell structure (Morii et al., 2023).

A later formulation defines an autoignitive reaction wave as a one-dimensional reaction front whose internal structure and propagation are governed solely by the imposed inflow velocity and the chemical autoignition kinetics, rather than by shock compression or flame-front diffusion. Under supersonic inlet conditions, fresh gas enters supersonically, autoignites along a continuous reaction zone, and exits still supersonic. On that basis, the wave satisfies the classical definition of a weak detonation, namely a flow that is supersonic both upstream and downstream of the reaction zone (Morii et al., 14 Jul 2025).

2. Governing equations and reduced descriptions

In full one-dimensional reactive flow, the governing equations comprise continuity, momentum, and total enthalpy or total energy, together with a chemical source term. In one formulation, the source term is written as $Q̇=\omega_f q$, with the fuel-consumption rate given by an Arrhenius law,

$\omega_f = \rho k_0 T^N \exp(-E_a/RT) Y_f .$

For a stationary wave written in the moving coordinate $\xi=x-Ut$, the governing equations reduce to

$$\frac{d(\rho u)}{d\xi}=0,\qquad \frac{d(\rho u^2+p)}{d\xi}=0,\qquad \frac{d(\rho u h + pu)}{d\xi}=\frac{\rho Q̇}{U},$$

with upstream and downstream equilibrium states matched across the wave (Morii et al., 2023).

A reduced but analytically tractable representation is Fickett’s reactive Burgers model in Lagrangian form,

$$\partial_t \rho + \partial_x p = 0,\qquad p=\tfrac12[\rho^2+\lambda_r Q].$$

In the two-step induction–reaction chain-branching model,

$$\partial_t \lambda_i = -H(\lambda_i)\exp[(\rho-1)/\epsilon],\qquad \partial_t \lambda_r = K(1-H(\lambda_i))(1-\lambda_r)^\nu ,$$

where $S_L$0, $S_L$1, and $S_L$2 is the Heaviside trigger when $S_L$3. With the immediate post-shock state used as reference, $S_L$4, the key nondimensional parameter is

$S_L$5

which measures the sensitivity of induction delay to temperature times the relative duration of energy release (Tang et al., 2012).

For self-sustained supersonic reaction waves in the extended Fickett analogue, the same model class is written as

$S_L$6

with nondimensionalization by the Chapman–Jouguet speed $S_L$7 and the reference induction-zone thickness. The sensitivity parameter $S_L$8 then acts as a non-dimensional activation energy controlling the exponential dependence of induction delay on shock strength (Radulescu et al., 2011).

3. Autoignition–compression feedback and front acceleration

In the Fickett shock-ignition problem, the mechanism is exposed directly by the characteristic structure. Along the forward characteristic family $S_L$9,

$D_{CJ}$0

while along a particle path,

$D_{CJ}$1

Energy release along particle paths amplifies forward-travelling pressure waves. Those waves then pre-compress the medium in the induction layer ahead of the reaction zone and modify the induction delays of successive particles. The induction delay, in turn, controls how long the forward-travelling waves remain inside the reaction zone and therefore modulates their further amplification. The paper’s closed-form asymptotic analysis gives the reaction-front trajectory

$D_{CJ}$2

and the initial acceleration

$D_{CJ}$3

so the reaction zone accelerates in proportion to the activation energy, the induction-to-reaction time ratio, and the heat release (Tang et al., 2012).

The same compression–autoignition coupling appears in higher-fidelity simulations. In the supersonic lifted hydrogen flame studied by Cheng et al., the first strong shock intersection along the jet axis at around $D_{CJ}$4 compresses and heats locally mixed pockets of $D_{CJ}$5 and $D_{CJ}$6, producing a reaction induction zone marked by isolated $D_{CJ}$7 pockets. After a short induction delay, an autoignition initiation zone forms, with rapid temperature rise and $D_{CJ}$8 production, and the lifted flame base appears at $D_{CJ}$9. In that case, pressure work is the dominant compressibility-induced heating mechanism: $D_{CJ}$0 near shocks, whereas $D_{CJ}$1 in high-shear regions and is typically less than $D_{CJ}$2 of $D_{CJ}$3 (Huang et al., 2021).

Huang and Zhang’s reflected-shock-tube simulations show that multiple autoignition hot spots can coexist. They identify a wall-surface re-compression hot spot, a reflected-shock/rarefaction-wave interaction hot spot, and a bulk-induction hot spot in the compressed mixture. The associated reaction fronts are all supersonic relative to local acoustics, with reported values $D_{CJ}$4, $D_{CJ}$5, and $D_{CJ}$6. The third hot spot launches a right-running reaction wave that catches the reflected shock, and the resulting coupled structure becomes a Chapman–Jouguet detonation at $D_{CJ}$7 (Huang et al., 2021).

4. Regime criteria, steady solutions, and weak-detonation interpretation

The Fickett asymptotics provide a local supersonic-front criterion directly. With acoustic speed ahead of the reaction zone $D_{CJ}$8, the front speed

$D_{CJ}$9

satisfies $Q̇=\omega_f q$0, so the front is locally supersonic. The same exposition distinguishes between a merely supersonic instantaneous front speed and a coherent running-up internal wave. It states that a truly “running-up” supersonic autoignitive wave forms only if the initial acceleration is large enough for coherent $Q̇=\omega_f q$1 amplification to steepen into an internal shock; empirically, the onset occurs for $Q̇=\omega_f q$2, and the threshold can be written

$Q̇=\omega_f q$3

The same regime map associates $Q̇=\omega_f q$4 with a “slow flame,” $Q̇=\omega_f q$5 with an intermediate regime, and $Q̇=\omega_f q$6 with a sustained supersonic “autoignitive detonation” regime (Tang et al., 2012).

In the one-dimensional steady theory of premixed inflow, the no-shock supersonic criterion is simply

$Q̇=\omega_f q$7

When chemical reaction in the preheat zone is retained, ignition occurs at $Q̇=\omega_f q$8, so steady-state solutions are not restricted to the isolated eigenvalues $Q̇=\omega_f q$9 and $\omega_f = \rho k_0 T^N \exp(-E_a/RT) Y_f .$0; instead, a continuous family of steady solutions arises for inlet velocities slightly above each one (Morii et al., 2023).

The 2025 analysis recasts that supersonic branch in classical detonation language. Neglecting diffusive and viscous fluxes under sufficiently supersonic inlet conditions reduces the governing equations to

$\omega_f = \rho k_0 T^N \exp(-E_a/RT) Y_f .$1

which are identical to the Rayleigh-flow equations if $\omega_f = \rho k_0 T^N \exp(-E_a/RT) Y_f .$2. On this basis, the paper proves the mathematical equivalence between supersonic autoignitive reaction waves and classical weak detonations. Its universal realization conditions are a fluid-dynamic requirement, $\omega_f = \rho k_0 T^N \exp(-E_a/RT) Y_f .$3, and a chemical-kinetic requirement, $\omega_f = \rho k_0 T^N \exp(-E_a/RT) Y_f .$4 (Morii et al., 14 Jul 2025).

A recurrent point of interpretation is whether a shock-free supersonic autoignitive wave should be separated from detonation theory. The 2023 formulation stresses the absence of a leading shock and of detonation cells, whereas the 2025 formulation identifies the same supersonic branch with weak detonation through Rayleigh-flow equivalence. These descriptions are not mutually exclusive: one emphasizes observable structure, the other thermodynamic classification (Morii et al., 2023, Morii et al., 14 Jul 2025).

5. Stability, pulsation, and transition phenomena

Self-sustained supersonic reaction waves need not remain steady. In the extended Fickett detonation analogue, stable and pulsating supersonic waves are both predicted, and increasing reaction-rate sensitivity drives the system from steady propagation to stable limit cycles and then to chaos through the classical Feigenbaum route. Linear analysis gives a critical sensitivity $\omega_f = \rho k_0 T^N \exp(-E_a/RT) Y_f .$5, with period doubling at $\omega_f = \rho k_0 T^N \exp(-E_a/RT) Y_f .$6, $\omega_f = \rho k_0 T^N \exp(-E_a/RT) Y_f .$7, and $\omega_f = \rho k_0 T^N \exp(-E_a/RT) Y_f .$8, and the exposition states that the dynamics become fully chaotic beyond $\omega_f = \rho k_0 T^N \exp(-E_a/RT) Y_f .$9. The physical mechanism is again the coherence between internal wave motion and energy release: stronger compression shortens induction delay, shifts the exothermic zone closer to the shock, and allows forward-travelling waves to spend more time in phase with heat release, whereas weakened compression lengthens the induction zone and breaks that coherence (Radulescu et al., 2011).

The ignition-transient study in Fickett’s model also reports both subsonic and supersonic internal flame-propagation regimes, consistent with experiment and previous reactive Euler models. That result connects the transient ignition problem to the broader family of self-accelerating internal waves in detonation development (Tang et al., 2012).

In wedge-stabilized oblique-detonation-wave configurations, stability is controlled by finite ignition delay and boundary-layer augmentation of the leading oblique shock. The ignition criterion is written as an induction-length constraint,

$\xi=x-Ut$0

with $\xi=x-Ut$1 set by the oblique-shock relation and $\xi=x-Ut$2 determined from the shocked state. Bachman and Goodwin show that at $\xi=x-Ut$3, $\xi=x-Ut$4 gives a burning boundary layer without detonation, $\xi=x-Ut$5 gives a prompt stable oblique detonation wave, and $\xi=x-Ut$6 produces a receding wave followed by redetonation near the leading edge. The proposed cycle repeats indefinitely with period $\xi=x-Ut$7, corresponding to $\xi=x-Ut$8, and is attributed to the coupling of boundary-layer growth, separation, back-pressure rise, and renewed over-compression at the leading edge (Bachman et al., 2020).

6. Canonical configurations, numerical evidence, and applications

The one-dimensional methane–air simulations used to introduce the general autoignitive-wave concept illustrate how inlet velocity selects among flashback, steady thin reaction zones, and shock-free supersonic propagation. The calculations use the full 1D compressible Navier–Stokes equations with one-step chemistry, an in-house solver “myRNS,” $\xi=x-Ut$9, $$\frac{d(\rho u)}{d\xi}=0,\qquad \frac{d(\rho u^2+p)}{d\xi}=0,\qquad \frac{d(\rho u h + pu)}{d\xi}=\frac{\rho Q̇}{U},$$0, ignition delay $$\frac{d(\rho u)}{d\xi}=0,\qquad \frac{d(\rho u^2+p)}{d\xi}=0,\qquad \frac{d(\rho u h + pu)}{d\xi}=\frac{\rho Q̇}{U},$$1, $$\frac{d(\rho u)}{d\xi}=0,\qquad \frac{d(\rho u^2+p)}{d\xi}=0,\qquad \frac{d(\rho u h + pu)}{d\xi}=\frac{\rho Q̇}{U},$$2, and $$\frac{d(\rho u)}{d\xi}=0,\qquad \frac{d(\rho u^2+p)}{d\xi}=0,\qquad \frac{d(\rho u h + pu)}{d\xi}=\frac{\rho Q̇}{U},$$3 (Morii et al., 2023).

Inlet condition	$$\frac{d(\rho u)}{d\xi}=0,\qquad \frac{d(\rho u^2+p)}{d\xi}=0,\qquad \frac{d(\rho u h + pu)}{d\xi}=\frac{\rho Q̇}{U},$$4	Reported behavior
$$\frac{d(\rho u)}{d\xi}=0,\qquad \frac{d(\rho u^2+p)}{d\xi}=0,\qquad \frac{d(\rho u h + pu)}{d\xi}=\frac{\rho Q̇}{U},$$5 $$\frac{d(\rho u)}{d\xi}=0,\qquad \frac{d(\rho u^2+p)}{d\xi}=0,\qquad \frac{d(\rho u h + pu)}{d\xi}=\frac{\rho Q̇}{U},$$6	$$\frac{d(\rho u)}{d\xi}=0,\qquad \frac{d(\rho u^2+p)}{d\xi}=0,\qquad \frac{d(\rho u h + pu)}{d\xi}=\frac{\rho Q̇}{U},$$7	flashback
$$\frac{d(\rho u)}{d\xi}=0,\qquad \frac{d(\rho u^2+p)}{d\xi}=0,\qquad \frac{d(\rho u h + pu)}{d\xi}=\frac{\rho Q̇}{U},$$8 $$\frac{d(\rho u)}{d\xi}=0,\qquad \frac{d(\rho u^2+p)}{d\xi}=0,\qquad \frac{d(\rho u h + pu)}{d\xi}=\frac{\rho Q̇}{U},$$9	$$\partial_t \rho + \partial_x p = 0,\qquad p=\tfrac12[\rho^2+\lambda_r Q].$$0	steady thin reaction zone, no shock
$$\partial_t \rho + \partial_x p = 0,\qquad p=\tfrac12[\rho^2+\lambda_r Q].$$1 $$\partial_t \rho + \partial_x p = 0,\qquad p=\tfrac12[\rho^2+\lambda_r Q].$$2	$$\partial_t \rho + \partial_x p = 0,\qquad p=\tfrac12[\rho^2+\lambda_r Q].$$3	flashback
$$\partial_t \rho + \partial_x p = 0,\qquad p=\tfrac12[\rho^2+\lambda_r Q].$$4 $$\partial_t \rho + \partial_x p = 0,\qquad p=\tfrac12[\rho^2+\lambda_r Q].$$5	$$\partial_t \rho + \partial_x p = 0,\qquad p=\tfrac12[\rho^2+\lambda_r Q].$$6	steady supersonic autoignitive wave, no shock

For the $$\partial_t \rho + \partial_x p = 0,\qquad p=\tfrac12[\rho^2+\lambda_r Q].$$7 case, the reported steady thin reaction zone sits at $$\partial_t \rho + \partial_x p = 0,\qquad p=\tfrac12[\rho^2+\lambda_r Q].$$8, the temperature rises from $$\partial_t \rho + \partial_x p = 0,\qquad p=\tfrac12[\rho^2+\lambda_r Q].$$9 to equilibrium $$\partial_t \lambda_i = -H(\lambda_i)\exp[(\rho-1)/\epsilon],\qquad \partial_t \lambda_r = K(1-H(\lambda_i))(1-\lambda_r)^\nu ,$$0, and the pressure rises from $$\partial_t \lambda_i = -H(\lambda_i)\exp[(\rho-1)/\epsilon],\qquad \partial_t \lambda_r = K(1-H(\lambda_i))(1-\lambda_r)^\nu ,$$1 to $$\partial_t \lambda_i = -H(\lambda_i)\exp[(\rho-1)/\epsilon],\qquad \partial_t \lambda_r = K(1-H(\lambda_i))(1-\lambda_r)^\nu ,$$2. For the $$\partial_t \lambda_i = -H(\lambda_i)\exp[(\rho-1)/\epsilon],\qquad \partial_t \lambda_r = K(1-H(\lambda_i))(1-\lambda_r)^\nu ,$$3 case, the shock-free supersonic autoignitive wave appears at $$\partial_t \lambda_i = -H(\lambda_i)\exp[(\rho-1)/\epsilon],\qquad \partial_t \lambda_r = K(1-H(\lambda_i))(1-\lambda_r)^\nu ,$$4, with equilibrium temperature $$\partial_t \lambda_i = -H(\lambda_i)\exp[(\rho-1)/\epsilon],\qquad \partial_t \lambda_r = K(1-H(\lambda_i))(1-\lambda_r)^\nu ,$$5, pressure $$\partial_t \lambda_i = -H(\lambda_i)\exp[(\rho-1)/\epsilon],\qquad \partial_t \lambda_r = K(1-H(\lambda_i))(1-\lambda_r)^\nu ,$$6, and reaction-zone thickness of order $$\partial_t \lambda_i = -H(\lambda_i)\exp[(\rho-1)/\epsilon],\qquad \partial_t \lambda_r = K(1-H(\lambda_i))(1-\lambda_r)^\nu ,$$7; the deflagration thickness is reported as $$\partial_t \lambda_i = -H(\lambda_i)\exp[(\rho-1)/\epsilon],\qquad \partial_t \lambda_r = K(1-H(\lambda_i))(1-\lambda_r)^\nu ,$$8 (Morii et al., 2023).

Large-eddy simulation has extended the concept to multidimensional turbulent combustors. In the sparse-Lagrangian MMC-LES study of the supersonic lifted hydrogen flame, mean and rms temperature, velocity, species mole fractions, and mixture fraction agreed with laser-Raman data within $$\partial_t \lambda_i = -H(\lambda_i)\exp[(\rho-1)/\epsilon],\qquad \partial_t \lambda_r = K(1-H(\lambda_i))(1-\lambda_r)^\nu ,$$9, and the mean lifted flame base was located at $S_L$00. The associated chemical explosive mode analysis showed that temperature dominates the explosive tendency in the central fuel jet, the $S_L$01 radical dominates on the fuel-rich side, and $S_L$02 dominates on the fuel-lean side (Huang et al., 2021).

Zhu et al. proposed a dynamic combustion model for supersonic turbulent combustion in which the sub-grid PDF is represented by two delta functions and the filtered source term is

$S_L$03

The model is designed to recover both the mixing-limited PaSR limit and the homogeneous-autoignition limit. In the strut–cavity flame-holder test case, the dynamic-model field $S_L$04 rises to $S_L$05 in the wake and cavity shear layers, and the peak local heat release in the wake increases from $S_L$06 to $S_L$07, a reported $S_L$08 increase that the paper interprets as the footprint of a supersonic autoignitive reaction wave embedded in a supersonic background flow (Zhu et al., 2023).

The practical implications drawn in the literature are similarly broad. The shock-free supersonic branch has been proposed as a combustor mode with high thermal efficiency from supersonic combustion without strong shocks, reduced wall-pressure oscillations and mechanical stresses compared with detonation-driven modes, elimination of detonation cell tearing, and a continuous range of operating points for $S_L$09 (Morii et al., 2023). The weak-detonation interpretation extends the same concept to scramjets, pulse-detonation engines, and astrophysical settings such as Type Ia supernovae, with the common realization conditions $S_L$10 and finite autoignition delay over a sufficient residence length (Morii et al., 14 Jul 2025).