1,439
Views
6
CrossRef citations to date
0
Altmetric
Articles

Analytical and numerical study of the expansion effect on the velocity deficit of rotating detonation waves

, , , &
Pages 761-774 | Received 17 Aug 2019, Accepted 09 Apr 2020, Published online: 07 May 2020

Figures & data

Figure 1. Schematic of an RDE.

Figure 1. Schematic of an RDE.

Figure 2. Typical RDE flow field in a two-dimensional simulation a-fresh propellants; b-reaction products; c-reaction products from the previous cycle.

Figure 2. Typical RDE flow field in a two-dimensional simulation a-fresh propellants; b-reaction products; c-reaction products from the previous cycle.

Figure 3. Ideal detonation-inert model in the wave-fixed reference frame

Figure 3. Ideal detonation-inert model in the wave-fixed reference frame

Figure 4. ZND model in the steady quasi-1D analysis.

Figure 4. ZND model in the steady quasi-1D analysis.

Figure 5. Profiles of the flow field calculated by the one-step kinetic model (a) pressure; (b) temperature; (c) expansion slope; (d) Mach number.

Figure 5. Profiles of the flow field calculated by the one-step kinetic model (a) pressure; (b) temperature; (c) expansion slope; (d) Mach number.

Figure 6. Comparison of p-v diagram between expansion and non-expansion detonation waves.

Figure 6. Comparison of p-v diagram between expansion and non-expansion detonation waves.

Figure 7. Normalised velocity of detonation Dex/Dnonex (colour and line labels) under the variation of the pre-detonation state, where Σ denotes (dσ/dx)/σ0 (a) with a constant temperature of 450 K; (b) with constant pressure 2 atm.

Figure 7. Normalised velocity of detonation Dex/Dnon−ex (colour and line labels) under the variation of the pre-detonation state, where Σ denotes (dσ/dx)/σ0 (a) with a constant temperature of 450 K; (b) with constant pressure 2 atm.

Figure 8. Comparison of the reaction rate when the pre-detonation temperature changes at 2 atm (a) reaction rate (b) extra work in reaction process in front of sonic point (In the p-v diagram, the extra work can be defined as the area between the leading shock process line and reaction-expansion process line shown in Fig. ).

Figure 8. Comparison of the reaction rate when the pre-detonation temperature changes at 2 atm (a) reaction rate (b) extra work in reaction process in front of sonic point (In the p-v diagram, the extra work can be defined as the area between the leading shock process line and reaction-expansion process line shown in Fig. 6).

Figure 9. Profiles of the flow field calculated by two-step kinetics model with pre-detonation (a) pressure; (b) temperature; (c) expansion slope; (d) Mach number.

Figure 9. Profiles of the flow field calculated by two-step kinetics model with pre-detonation (a) pressure; (b) temperature; (c) expansion slope; (d) Mach number.

Figure 10. Normalised velocity of detonation and equilibrium reaction progress parameter (colour and line labels) under the pre-detonation state variation, where Σ denotes (dσ/dx)/σ0 (a) the normalised deficit of velocities Dex/Dnonex with a constant temperature of 450 K; (b) normalised deficit of velocities Dex/Dnonex with constant pressure of 2 atm; (c) equilibrium reaction progress parameter λ with a constant temperature of 450 K; (d) equilibrium reaction progress parameter λ with a constant pressure of 2 atm.

Figure 10. Normalised velocity of detonation and equilibrium reaction progress parameter (colour and line labels) under the pre-detonation state variation, where Σ denotes (dσ/dx)/σ0 (a) the normalised deficit of velocities Dex/Dnon−ex with a constant temperature of 450 K; (b) normalised deficit of velocities Dex/Dnon−ex with constant pressure of 2 atm; (c) equilibrium reaction progress parameter λ with a constant temperature of 450 K; (d) equilibrium reaction progress parameter λ with a constant pressure of 2 atm.

Figure 11. Definition of geometric parameters in the RDE flow field.

Figure 11. Definition of geometric parameters in the RDE flow field.

Figure 12. Comparison between the numerical and expansion model for the deficit of the detonation velocity when the inlet stagnation state changes (a) with a constant total temperature of 450 K; (b) with a constant total pressure of 5 atm.

Figure 12. Comparison between the numerical and expansion model for the deficit of the detonation velocity when the inlet stagnation state changes (a) with a constant total temperature of 450 K; (b) with a constant total pressure of 5 atm.