3.1.1. Mixing process

The mixture generally has a non-uniform state with temperature, pressure and ER gradients (Boulal et al. Reference Boulal, Vidal and Zitoun2016) in practical detonation-based engines. Here, considering a more realistic engine, the mixing process is simulated under the supersonic air flow and transverse fuel injection before the mixture ignition. The mixing process time is set as t = 0.84 ms, and two groups of transverse fuel jets with high stagnate pressure are injected into the chamber, as presented in figure 1(b). As a result, an inhomogeneous hydrogen–air mixture is generated in the combustion chamber under the influence of the main supersonic flow. The snapshots of the mass fraction of H₂ during the mixing process under the main flow are given in figure 3. Thanks to the high dynamic pressure of the main flow, the fuel jet is deflected rapidly because of the low dynamic pressure ratio between the fuel and air. As a blockage effect is generated by the fuel jets, the penetration depth is gradually increasing. With the aid of Kelvin–Helmholtz (K-H) instability, vast fuel vortices propagate in the downstream chamber. Under the influence of the flow instability and diffusion effect, the uniformity of the hydrogen–air mixture is augmented. Since a high Br and low diffusion effect is formed in the 2-D simulation, a high mass fraction of H₂ in the upper and lower walls is observed.

Figure 3.

Variations of H₂ mass fraction during the mixing process under supersonic incoming air: (a) t = 0.22, (b) t = 0.26, (c) t = 0.34, (d) t = 0.50, (e) t = 0.70 and ( f) t = 0.82 ms.

After the end of the mixing process, the distributions of mixture fraction along the horizontal and transverse directions in the selected domain of x = 0–40 cm are plotted in figure 4. The mixture fraction is defined as follows:

(3.1)\begin{equation}f = \frac{{(1/2{M_H})({Y_H} - {Y_{H,2}}) - (1/{M_O})({Y_O} - {Y_{O,2}})}}{{(1/2{M_H}({Y_{H,1}} - {Y_{H,2}})) - (1/{M_O})({Y_{O,1}} - {Y_{O,2}})}},\end{equation}

where M denotes molecular weight, Y represents elements mass fraction and subscripts 1 and 2 denote fuel jet and incoming air, respectively. Note that the molecular and mass fraction of H in the fuel jet is 1, while the molecular and mass weight of O is 0.

Figure 4.

The mixture fraction distribution in the (a) horizontal and (b) transverse dictions at the end of the mixing process.

In figure 4(a), a highly fluctuating value is formed in the chamber between 0 and 1, while with the mixing propagating into the downstream chamber, the uniformity of the mixture is augmented. In the transverse direction, the mixture fraction still has a lot of fluctuations even in the downstream tube, and the mixture fraction value is in the range of 0.2–0.6. There is a more uniform mixture in the central domain of the chamber in contrast to the upper and lower walls. Part of the reason is that this is a 2-D simulation without a 3-D influence that has a larger diffusion effect and lower Br. Generally, the mixtures in the horizontal and transverse directions have high-concentration gradients and are far from the stoichiometric state.

3.1.2. Global flame acceleration characteristics

After the completion of the mixing process, the combustion evolution in the inhomogeneous mixture without fluid obstacles is depicted in figure 5, corresponding to case 3. At t = 1.09 ms, two regions with high temperatures are formed in the upper and lower walls, attributed to the hot spot ignition. The lower-temperature areas in figure 5(a) indicate the presence of fuel vortices. As seen in box A2, a substantial number of lower-temperature products are generated with the highly turbulent flow since there are a lot of fuel-rich mixtures due to the K-H and Rayleigh–Taylor (R-T) instabilities. Nevertheless, the combustion temperature significantly decreases (box A2), caused by the high H₂ concentration in the upstream chamber. As the flame propagates, the temperature increases from 1700 to 2400 K; see box B1, owing to the turbulent mixing reaction, while the temperature in the flame tip is still low (circle C in figure 5e). The temperature distribution in the combustion product is, therefore, extremely uneven. Next, with the flame front further propagating, it elongates steadily, and the combustion temperature is increased; see box D in figure 5( f), and therefore, the unburned mixture is preheated due to the increasing intensity of the LSW. However, the detonation transition fails because of the long-time FA and dramatically inhomogeneous mixture.

Figure 5.

Snapshots of the temperature contours during the FA in the inhomogeneous supersonic mixture without transverse jet obstacles (case 3): (a) t = 1.09, (b) t = 1.12, (c) t = 1.24, (d) t = 1.54, (e) t = 1.70 and ( f) t = 1.95875 ms.

Referring to the former literature (Peng et al. Reference Peng, Huang, Deiterding, Luan, Xing and You2018; Zhao et al. Reference Zhao, Deiterding, Liang, Cai and Wang2023a), to activate DDT, the jet obstacles are employed in the premixed mixture to stimulate DDT. Four groups of transverse jets are utilized (case 2). The temperature variations are given in figure 6. After the jet is injected into the chamber, a perturbed penetration and oblique shock wave are formed, preheating the unburned mixture. Furthermore, when the flame front propagates in this region, the transverse jet also provides a blockage effect to some degree as a solid obstruction. The flame–turbulence interaction is observed accordingly, resulting in a conspicuously increasing flame surface area (box A). Influenced by K-H and R-T instabilities, a corrugated flame front is generated, feeding back to the FA (Gamezo et al. Reference Gamezo, Ogawa and Oran2008).

Figure 6.

Time sequence of temperature fields showing the propagation of the flame–turbulence, flame–shock and shock–shock interactions during the FA and DDT in case 2. SW-shock wave, IHS-auto-ignited hot spot, DW-detonation wave: (a) t = 0.2450, (b) t = 0.2750, (c) t = 0.3200, (d) t = 0.3800, (e) t = 0.4550, ( f) t = 0.5000, (g) t = 0.5150, (h) t = 0.5300, (i) t = 0.5450 and ( j) t = 0.5600 ms.

Consequently, the LSW is strengthened, causing a preheated zone with high temperature at t = 0.5000 ms. The third group of the transverse jets is deflected by the main flow, and an auto-ignited hot spot is generated as manifested in figure 6( f–h) as a result of the intensive shock wave. Consequently, the onset of detonation occurs because the hot spot ignition evolves into detonation combustion at the upper boundary wall. This DDT mechanism is caused by the flame-kernel ignition through shock–boundary interaction, as reported in Goodwin's study, where a slightly higher Br was employed (Goodwin et al. Reference Goodwin, Houim and Oran2016). Such a DDT mechanism is also reported in Grogan & Ihme's (Reference Grogan and Ihme2015, Reference Grogan and Ihme2017) studies and it is found that the detonation transition though the flame-kernel ignition is due to the shock wave amplification by coherent energy release (SWACER) mechanism. Finally, the detonation wave propagates into the downstream chamber, as shown in figure 6(i, j).

The transient temperature contour and the corresponding AMR distribution in case 2 are given in figure 7. The AMR works well in large temperature and pressure gradients, especially in the regions having shock wave, flame front and jet obstacles. These AMR regions also dynamically move as the flame tip and shock wave evolve. Hence, the current AMR threshold values are sound.

Figure 7.

The transient temperature contour and the corresponding AMR distribution: (a-i) t = 0.37 and (b-i) t = 0.40 ms.

Additionally, a comparison of the flame tip versus time in case 1 and case 2 is plotted in figure 8. The flame front position in the premixed mixture without a jet is larger than the case with the transverse jet in the incipient stage. This is due to the Br influence stemming from the transverse jet when a host of fluids is injected into the chamber, introducing a high-pressure region at the beginning of the flame front, which has a negative influence on FA at the initial time. This result is also reported in our previous study (Zhao et al. Reference Zhao, Liang, Deiterding, Cai and Wang2022a).

Figure 8.

Evolution of axial position of the flame front in the cases for premixed mixtures with and without transverse jet obstacles.

By contrast, the transverse jets provide vast perturbed vortices and a moderate Br, introducing plenty of flame–turbulence, flame–shock and shock–shock interactions and increasing the flow instabilities and the flame surface area, as discussed above. Subsequently, the deflagration flame evolves into detonation rapidly, as marked by circle A. Therefore, in the later stage, many transverse jets ejecting into the supersonic mixture benefit DDT when the flame front approaches the perturbed region. This is consistent with the previous research where one transverse jet obstacle is adopted in the static mixture (Frolov et al. Reference Frolov, Smetanyuk, Aksenov and Koval’2017; McGarry & Ahmed Reference McGarry and Ahmed2017; Peng et al. Reference Peng, Huang, Deiterding, Luan, Xing and You2018; Cheng et al. Reference Cheng, Zhang, Liu and Wang2021a). Furthermore, the present jet obstacles further introduce an oblique shock wave due to the supersonic flow, which is more reliable in increasing the intensity of the FA and LSW (Zhao et al. Reference Zhao, Deiterding, Liang, Cai and Wang2023a).

Consequently, inspired by case 2, four groups of head-on jets are also employed in the inhomogeneous supersonic mixture. To increase the uniformity of the mixture, the jet obstacles are composed of air. The FA and detonation transition snapshots are given in figure 9. The temperature variations in incipient time are similar in the above case. The combustion temperature is significantly as low as 1800 K; see box A1.

Figure 9.

Snapshots of the temperature fields during the FA and DDT under the influence of transverse jet obstacles (case 4). F-T – flame-turbulence, SW – shock wave, LSW – leading shock wave: (a) t = 1.06, (b) t = 1.105, (c) t = 1.165, (d) t = 1.255, (e) t = 1.285, ( f) t = 1.345, (g) t = 1.405, (h) t = 1.435, (i) t = 1.475, ( j) t = 1.5225, (k) t = 1.53125 and (l) t = 1.570 ms.

Nevertheless, as the air jets are injected into the chamber, as seen in figure 9(c), oblique shock waves form. An apparent flame–turbulence interaction is observed when the flame reaches the air jet regions (see figure 9d–f), increasing the flame surface area dramatically, and then the flame speed augments accordingly. Furthermore, the combustion temperature is increased to 2700 K (boxes B1 and B2) due to the increasing uniformity of the mixture when introducing a lot of perturbed jet flows. As a result, preheated zones are generated ahead of the flame tip because of the strong shock wave. Many shock waves focus on and collide with the oblique shock wave generated by the transverse jet. These features result in prominent preheated zones, as shown in figure 9(g–i). The elongated flame front (circle C) generates a dramatically increased flame surface area, resulting in a high energy release ratio (Gamezo et al. Reference Gamezo, Ogawa and Oran2008) and a significantly increased combustion temperature in contrast to figure 5 without jet obstacles. A stronger LSW is formed subsequently. The LSW preheats the turbulent jet, as figure 9(i–k) shows. Finally, a localized explosion point is formed against the lower wall, as figure 9(i) shows. Detailed information on the DDT and the subsequent detonation propagation will be discussed in §§ 3.3 and 3.4.

For the sake of revealing the failure of the DDT in the premixed and inhomogeneous mixtures without transverse jet obstacles, the positions of the flame front, LSW and the corresponding flame front propagation velocity are superimposed in figure 10(a,b), respectively. The FA in the premixed supersonic mixture can be classified into four stages. Stage I (0.2–0.33 ms): the flame front propagation velocity is decreased rapidly due to the influence of the supersonic flow. When it is close to the flow speed of the unburned mixture, the flame front goes into stage II (0.33–0.58 ms), and the flame propagation speed continues to increase as the flame surface area increases. In the next stage (0.58–0.68 ms), the flame propagation speed remains stable, which can be considered as a deflagration flame because its speed approaches half of the CJ speed. It is observed that the speed is around 1000 m s⁻¹, which is also close to the speed of sound in the combustion products. This means that DDT usually occurs when the flame tip propagation speed is close to the speed of sound in the product gas because there is no continuous transition from the deflagration combustion in the lower branch to the detonation model along the Hugoniot curve (Zhao et al. Reference Zhao, Liang, Deiterding, Cai and Wang2022a, Reference Zhao, Deiterding, Liang, Wang, Cai and Duell2023b).

Figure 10.

The positions of the flame front, LSW and the corresponding flame propagation velocity as a time sequence in the (a) premixed (case 1) and (b) non-uniform (case 3) supersonic mixture. The blue curve denotes the flame front propagation velocity.

Nevertheless, after that, due to the lack of flow instability and blockage effect within the smooth chamber and an increasing length between the flame tip and LSW (L_fs), the flame propagation speed decreases steadily. Hence, as reported in former DDT simulations (Oran et al. Reference Oran, Chamberlain and Pekalski2020; Gamezo et al. Reference Gamezo, Bachman and Oran2021), a constant L_fs, namely, the situation of a shock–flame complex (Gamezo et al. Reference Gamezo, Ogawa and Oran2008; Oran et al. Reference Oran, Chamberlain and Pekalski2020), is an essential condition to activate the detonation transition. In contrast, an increasing L_fs results in a failure of detonation initiation. As such, DDT is difficult to form in a deflagration flame within a smooth tube despite being in a premixed supersonic mixture without considering a perturbation configuration.

In figure 10(b), flame propagation in the inhomogeneous mixture can be classified into three stages without stage VI in the inhomogeneous mixture. Stage I is the ignition time, and it has a high flame propagation speed but reduces abruptly owing to the extremely fuel-rich mixture. Under the effect of the main flow, the flame–turbulence interaction increases again in stage II. In stage II, the propagation speed continues to increase as time evolves, whereas the required time is increased from 0.3 ms to 0.6 ms compared with the FA ratio in the premixed mixture as plotted in figure 10(a); hence, the increasing ratio of the flame propagation speed is reduced and has more fluctuation. After that, in stage III, the flame speed approaches half of the CJ value, whereas its speed slightly decreases as L_fs increases. Hence, the detonation transition is still a failure.

3.2.1. Premixed mixture

In the incipient stage of the FA in the premixed supersonic mixture for case 1, a typical flame front structure and the corresponding temperature variation at t = 0.3100 ms are given in figure 11. A breached flame tip is formed due to the two hot spot ignitions in the upper and lower walls, which subsequently merge in the middle region. Analysis of species mass reveals that the HO₂ radical forms early in the front of the flame structure, followed by O, H and OH, initiating a chain-branching chemical reaction. In contrast, the peak value of HO₂ is a lower order than the peak value of H, O and OH radicals. The peak mass fraction of HO₂ reaches 1.5 × 10⁻⁵. As it is depleted rapidly, the concentrations of H, O and OH are increased significantly, accomplished by an increasing temperature from 306 to 2000 K; hence, the chain-branching chemical is continued through the consumption of HO₂ and H₂O₂. The peak value of H, O and OH radicals reach 0.0039, 0.0217 and 0.0087 with the mass fraction of H₂ and O₂ decreasing through the branching reactions of H + O₂ = O + OH and O + H₂ = H + OH. As a result, the active reaction zone is around a range between 0.525 and 0.60, which has a single layer. This flame front structure is in agreement with that of a premixed flame, as discussed by Law (Reference Law2010). In the middle region of the flame front structure, due to the flow instability, some residual mixture exists in the downstream product; see circle B1 in figure 11(a 2). Thus, the flame front structure is more complicated, and the reaction front has two layers. As a result, there is a hump in the variations of species mass fraction and temperature curves.

Figure 11.

Variations of temperature and species mass fraction across the flame front (a) near the wall and (b) middle region at the incipient FA for the premixed mixture in case 1 (t = 0.3100 ms).

As the flame front further propagates into a stable state with maximum flame propagation speed, typical flame structures in the premixed mixture, including the species mass information at different flame positions at t = 0.6800 ms, are plotted in figure 12. A lot of wrinkled flames appear in the flame front. Figure 12(a,b) gives the mass fractions at y/Ly = 0.125 and 0.625, corresponding to lines B1 and B2, respectively. The species curves show the same result: the HO₂ forms early in the front of the flame structure, followed by H and O, resulting in main chain-branching reactions and mixture burning.

Figure 12.

Variations of temperature and species mass fraction across the flame front in (a) lower (0.25/2.0) and (b) upper (1.25/2.0) regions at the incipient FA in case 1 at t = 0.6800 ms.

Furthermore, as the mass fraction of HO₂ reaches a peak value of 3.2 × 10⁻⁴, it is consumed sharply, with the mass augmentations of OH, O and H to peak values of O = 0.018, H = 0.0037 and OH = 0.009. Hence, the chain-branching reactions keep going, forming a gradual increase in H₂O product and an abrupt decrease in H₂ and O₂. Thus, the energy release ratio increases, increasing the combustion temperature from 428 to 2200 K. Due to the much higher pressure and temperature in the preheated mixture, the flame thickness of 1.2 × 10⁻³ m is much lower compared with the initial stage.

In contrast, the grid resolution can still achieve 48 points per flame thickness, which is sufficient to resolve the combustion characteristics. At the flame position y/Ly = 0.625, the species mass fraction is similar to that at y/Ly = 0.125. The HO₂ mass fraction appears early at the flame tip, followed by the H, O and OH radicals. With the rapid consumption of HO₂, H₂ and O₂, the OH and product of H₂O increase, leading to an increase in the combustion temperature.

3.2.2. Inhomogeneous supersonic mixture

Concerning the FA in the inhomogeneous supersonic mixture in case 3, the flame structure information at t = 0.1540 ms, around the two finger flame tips (y/Ly = 0.125, 0.75), is shown in figure 13. In the lower flame structure, the mass fraction is hugely uneven. In the flame tip, the HO₂ radical appears first, followed by H, O and OH, with the mass fraction of HO₂ augments, resulting in the increased temperature from 327 to 1420 K. But the increased ratio of the temperature is relatively lower than that in the premixed mixture. Much unburned H₂ is still left behind with a high mass fraction of H₂ = 0.01. As such, the second combustion layer continues in the downstream product so that another flame structure layer with a hump shape appears downstream. With the HO₂ increasing again, accompanied by the increases of H, O and OH radicals, the mass fractions of H₂ and O₂ are subsequently consumed rapidly, further increasing the temperature to 2010 K. The combustion process is, therefore, more complicated in contrast to the premixed mixture. Hence, the temperature contour near the lower wall exhibits an extreme maldistribution.

Figure 13.

Variations of temperature and species mass fraction across the flame front in (a) lower (0.25/2.0) and (b) upper (1.5/2.0) regions during FA for the inhomogeneous mixture in case 3 (t = 0.1540 ms).

Line A2 has just one layer of the flame front structure. The species mass fraction shows that the peak values of H and O are 0.0072 and 0.0018, respectively, so that the H mass fraction is almost two times that of the premixed mixture due to the fuel-rich mixture. Thereby, a significantly lower peak value of OH mass fraction of 0.00078 appears, contributed by the lack of O₂. Hence, this results in an extreme fuel-rich reaction, leading to a lower increased temperature from 325 to 1380 K. Consequently, many unburned fuels are left behind in the downstream chamber with a high mass fraction, H₂ = 0.014. Due to the fuel-rich combustion, a large flame thickness is generated when compared with the premixed mixture.

Figure 14 further gives more detailed information on the temperature, hydrogen, ER and flame front structure with species information along the y-direction position at t = 1.54 ms. To picture the variation of temperature, hydrogen, ER and the radical element distribution, 3-D plots are presented on the right side, where the x-axis is in a range of from 30.0 cm to 36.0 cm, while discrete y-axis coordinates are selected at 0.25, 0.5, 0.75, 1.25. 1.5 and 1.75 cm. The corresponding specific values are plotted in the z-axis direction of the 3-D plot.

Figure 14.

The (a) temperature, (b) equivalation ratio (ER) and flame (c–g) front structure distribution at t = 0.15400 ms for the inhomogeneous supersonic mixture (case 3) along different horizontal positions in the y-direction.

A high fluctuation of the temperature distribution is observed after the flame front, especially in the regions near the lower and upper walls. A gradually increased temperature is noted near the upper and lower walls. Hence, the combustion reaction is still continued in the downstream product. Although there is a much larger flame front surface area than in the premixed mixture, the temperature at the front is still significantly lower. In Zheng et al.'s (Reference Zheng, Kaplan, Houim and Oran2019) study, the FA occurred in an ER gradient normal to the flame propagation direction, and they also found that there is a larger surface area in the inhomogeneous mixture compared with the homogenous mixture, but this does not result in a high energy release ratio. This is due to the large surface area mainly located in the fuel-lean regions. Therefore, the temperature highly depends on the ER distribution.

As shown in the H₂ and ER maps, the upper wall has a much thicker fuel layer, and the ER distribution is highly uneven. A high fuel-rich condition exists in the upper half-domain, and the ER value is around 2.5. Yet, a much lower fuel concentration is observed near the middle domain (ER is approximately 0.8), resulting in two colour distributions in the combustion product. The ER map of the combustion product includes grey (ER = 0) and brown $(\textrm{ER} = \infty )$, representing that the fuel and oxidizer have been consumed, respectively. The H, O and OH mass fraction distributions also yield a high fluctuation in the downstream combustion. There is a high mass fraction of H in the upper combustion product, while a high O mass fraction is formed in the lower combustion product. Such a feature further leads to a local reaction in the downstream tube due to a further mixing process between the H₂ and O₂, especially in the connecting region between grey and brown; see circles C1 and C2 in figure 14(c). The HO₂ concentration mainly occurs in the flame front structure, with a maximum of 2.0 × 10⁻⁵. It is an essential species for initiating induction reaction, and a much thicker induction reaction zone is formed in the fuel-rich layer. Consequently, a much lower temperature occurs in the upper flame front. In general, the temperature contour distribution is particularly uneven.

To sum up, the flame front structure in the inhomogeneous supersonic mixture is more complicated, characterized by multiple reaction layers compared with the premixed mixture. Additionally, the flame thickness in the inhomogeneous mixture is much higher than that in the premixed mixture. The presence of extreme fuel-lean and fuel-rich mixtures with a perturbed distribution results in a significantly lower combustion temperature with a highly uneven distribution, generating a lower FA ratio. However, the combustion temperature is further increased due to the continuous reaction in the downstream flame resulting from the turbulent mixing process.

3.3.1. Mechanism of onset of detonation

Here, more instantaneous snapshots of the temperature and pressure variations surrounding the explosion point are given in figure 15, corresponding to case 4. A localized explosion point takes place, and it propagates into the downstream lower wall, as seen in the series of circles A. Meanwhile, a high-pressure region is generated in the lower wall accompanied by the detonation wave (see circles B), forming a transverse wave shedding from the detonation front.

Figure 15.

Snapshots of the temperature and pressure contours showing the propagation of the explosion point for case 4: (a1) t = 1.56375, (b1) t = 1.56625, (c1) t = 1.56875 and (d1) t = 1.57500 ms.

To thoroughly picture the mechanism of detonation transition, figure 16 shows the transient variations of temperature, pressure and OH before forming the localized explosion point with a time step of 0.00125 ms. A remarkable shock wave is generated at the lower wall and propagates in the downstream mixture (figure 16a 2), forming an intense Mach stem subsequently (figure 16b 2), where the maximum pressure reaches 3.0 MPa. In addition, the fierce shock wave also results in a high combustion temperature (circle A1). As a result, a localized explosion point forms in the lower wall (figure 16c 1), forming a high explosion pressure and OH concentration (circles B1 and C1). The onset of detonation is triggered, propagating in the downstream mixture. Hence, the strong Mach stem is responsible for the current DDT mechanism, which generates a high-pressure and high-temperature spot sufficient for the detonation transition.

Figure 16.

Contours of the temperature, pressure and OH variation illustrating the activation of the onset of DDT in case 4. SW – shock wave, MS – Mach stem, EP – explosion point: (a1) t = 1.56125, (b1) t = 1.56250 and (c1) t = 1.56375 ms.

Figure 17 presents a schematic depiction elucidating the progression of the shock wave and flame front evolution observed during the DDT initiation, as depicted in figure 16. Preceding the advancement of the flame front, a robust shock wave propagates downstream, with the reflected shock wave rebounding off the lower wall, thereby engendering a Mach stem configuration and thereby inducing a reactive gradient. Following this, a high-temperature explosion point manifests on the upper wall, subsequently evolving into an explosion point. Thus, the inception of the Mach stem structure emerges as a pivotal procedural facet in this phenomenon.

Figure 17.

Physical illustration of the shock wave (SW), Mach stem (MS) and flame front evolution during the onset of detonation. DW – detonation wave.

The DDT mechanism is caused by the gradient mechanism (Oran et al. Reference Oran, Chamberlain and Pekalski2020). This mechanism is also reported in Gamezo et al.'s (Reference Gamezo, Ogawa and Oran2008) and Goodwin et al.'s (Reference Goodwin, Houim and Oran2016) studies, where a high Br was employed in an obstacle-laden chamber, resulting in an intense Mach stem after undergoing a series of FA.

The current DDT mechanism can be further quantitively analysed by the pressure and temperature records, as shown in figure 18, where the red curves represent the temperature, and the other lines denote the corresponding pressure variations along the x direction along the lower wall. Note that a rapidly increased pressure value from 1.17 to 2.85 MPa (circle A) is formed thanks to the Mach stem; see figure 16(b 2), and the temperature of the unburned mixture is preheated to 820 K, resulting in a reactive gradient. Under these conditions, the DDT is activated with a localized explosion point with high pressure and temperature values.

Figure 18.

Variations of temperature and the corresponding pressure during the onset of detonation in case 4 along the x-direction in the lower wall, the red lines detonate the temperature. The transient time sequences are t = 1.5575 ms, 1.55875 ms, 1.5600 ms, 1.56125 ms, 1.56250 ms, 1.56375 ms and 1.56625 ms.

3.3.2. Mechanism of detonation reinitiation

The next flame propagation is further presented in figure 19. The detonation wave does not survive and is instead extinguished. As illustrated, the detonation wave is decoupled from a flame and a shock wave, as shown in the lower flame front and box A in the temperature and pressure contours, respectively. A gradually reduced temperature and pressure form. As such, a detonation quenching occurs. This is attributed to the extreme fuel-lean mixture in the downstream chamber. As shown in the series of boxes A in the ER map, a significantly lower ER value of around 0.15 exists ahead of the flame front. Such a lower ER distribution is caused by the air transverse jet obstacles located in x = 36 cm, which introduce a large amount of air, as presented in figure 9( j). As a result, sustained detonation propagation highly depends on a suitable ER.

Figure 19.

Snapshots of temperature, pressure and corresponding ER in the detonation extinction for case 4. SW – shock wave: (a1) t = 1.58125, (b1) t = 1.58375 and (c1) t = 1.59375 ms.

A donation reinitiation is noted after the flame further propagates into the downstream chamber; see figure 20. First, a remarkable Mach stem is generated against the lower wall; see circle B1. It propagates into the region having a high ER value (circle C2), located in the lower boundary wall compared with the middle part in circle C1. Then, a localized explosion appears (figure 20b 1), propagating in the lower wall region with a high ER value; see the series of arrows in the temperature contours. Finally, more detonation combustion is formed at t = 1.64125 ms where more areas with a high ER mixture exist. Therefore, sustaining detonation highly relies on a continuously suitable mixture. Although there is a quite small and elongated region of the suitable and sensitive mixture, the detonation wave propagation can still survive.

Figure 20.

The detonation reinitiation rendered by temperature, pressure and ER evolution in the inhomogeneous supersonic mixture (case 4): (a1) t = 1.61875, (b1) t = 1.62625, (c1) t = 1.63125 and (d1) t = 1.64125 ms.

To elucidate the flow structure comprehensively, a schematic representation delineating the evolution of the shock wave, LSW and flame front antecedent to detonation initiation is provided in figure 21. Herein, the blue-coloured region demarcated along the lower periphery signifies a composition conducive to sensitivity. Initially, an incident shock wave propagates downstream adjacent to the bottom wall. Subsequent to this propagation, a Mach stem is engendered, subsequently precipitating the localization of an explosion point. This critical juncture serves as the impetus for detonation combustion, which subsequently propagates throughout the sensitive and elongated region.

Figure 21.

Schematic illustration of the flow structure before the onset of detonation. SW – shock wave, SM – sensitive mixture, LSW – leading shock wave, MS – Mach stem, LEP – localized explosion point.

The temperature and pressure records during the detonation reinitiation along the x direction from t = 1.62125 to t = 1.63375 ms are also superimposed in figure 22, where the red curves denote the temperature, and the other lines represent the pressure. Increasing pressure and temperature have occurred; see arrow A, where the peak pressure and temperature reach T = 1504 K and p = 1.95 MPa at t = 1.6225 ms, respectively. Hence, under the high pressure and temperature conditions post-shock, the onset of detonation is triggered at t = 1.62375 ms. Next, the detonation propagates downstream (figure 20d 1). The DDT mechanism is, therefore, the same as the first localized explosion point, as discussed in figures 16 and 17, which is caused by the gradient mechanism because of the intense Mach stem. Such a DDT mechanism is consistent with the result reported by Lee, Knystautas & Yoshikawa's (Reference Lee, Knystautas and Yoshikawa1980) experiment study, referred to as the SWACER mechanism.

Figure 22.

Temperature variations and the corresponding pressure during the detonation reinitiation in case 4 along the x direction in the lower wall. The transient time sequences are t = 1.62125 ms, 1.6225 ms, 1.62375 ms, 1.625 ms, 1.62625 ms, 1.62875 ms, 1.63125 ms and 1.63375 ms.

In summary, the inhomogeneous supersonic mixture exhibits more dynamic characteristics, including the detonation transition, quenching and reinitiation, which primarily stem from the extremely inhomogeneous ER distributions. The results indicate that the detonation cannot be sustained when the ER falls below 0.15. The DDT mechanism is ascribed to the Zel’dovich gradient mechanism (Zel'Dovich et al. Reference Zel'Dovich, Librovich, Makhviladze and Sivashinsky1970) thanks to a strong Mach stem forming in the boundary wall. As such, in the smooth tube, the DDT mechanism resembles previous numerical and experimental findings, where a robust Mach stem was generated in a solid-laden chamber with a high Br in the FA (Boeck et al. Reference Boeck, Hasslberger and Sattelmayer2014; Goodwin et al. Reference Goodwin, Houim and Oran2016) and DDT scenarios (Vollmer et al. Reference Vollmer, Ettner and Sattelmayer2012; Zheng et al. Reference Zheng, Kaplan, Houim and Oran2019) occurred in the hydrogen–air or methane–air mixtures. Moreover, the detonation reinitiation mechanism has the same DDT mechanism. Still, the formation of the DDT requires a highly sensitive mixture region so that a discretized low ER mixture can be utilized to block the detonation propagation. Yet, the detonation reinitiation can be activated later when the intense shock wave propagates in a highly sensitive mixture.

3.4.1. Global propagation characteristics

The initial detonation propagation in the supersonic non-uniform mixture is illustrated in figure 23. A smooth flame front scenario is generated owing to the overdriven detonation and low increased temperature ratio between the detonation wave and post-LSW. A sudden lower-temperature region in the leading flame front is formed; see box A1; however, an increased temperature subsequently forms when the detonation wave impacts the upper wall, generating a strong Mach stem. A curved flame appears in box A2.

Figure 23.

Snapshots of the detonation propagation with temperature evolution in the inhomogeneous supersonic mixture (case 4). TW – transverse wave, SL – slip line, TP – triple point, SW – shock wave: (a) t = 1.65125, (b) t = 1.6625, (c) t = 1.67125 and (d) t = 1.67750, (e) t = 1.7, ( f) t = 1.7075, (g) t = 1.715 and (h) t = 1.7225 ms.

Besides that, a remarkable transverse wave occurs in the downstream combustion product. A host of vortices appear in the slip line due to the Richtmyer–Meshkov (R-M) instability (Mahmoudi et al. Reference Mahmoudi, Karimi, Deiterding and Emami2014; Mahmoudi & Mazaheri Reference Mahmoudi and Mazaheri2015); see figure 23(d). A triple point is also generated. A curved flame with a much lower temperature appears again (box A3), but an increased temperature is formed later when it propagates in the downstream region (box A4). Multiple triple points (box A5) and a partially smooth flame front (box A6) with lower temperature form in the leading flame tip; see figure 23(g,h). However, the combustion temperature in the downstream product is much higher than the initial overdriven detonation propagation. This is reasonable because the uniformity of the upstream mixture is much higher. Moreover, the consumption of the unburned pocket is further continued in the downstream flame through the diffusion effect and turbulent mixing generated by the K-H and R-M instabilities, as suggested by previous investigations (Austin Reference Austin2003; Radulescu et al. Reference Radulescu, Sharpe, Lee, Kiyanda, Higgins and Hanson2005; Radulescu et al. Reference Radulescu, Sharpe, Law and Lee2007; Mahmoudi & Mazaheri Reference Mahmoudi and Mazaheri2015; Han et al. Reference Han, Wang and Law2019). In general, dramatically complex dynamic behaviours of the detonation propagation, including the curved flame and cellular instability in the inhomogeneous supersonic mixture, occur.

For quantitively evaluating the detonation propagation, the flame front position at the embryonic stage and the corresponding averaged flame front propagation velocity with time evolution near the lower wall (y/Ly = 0.1/2.0), middle region (y/Ly = 1.0/2.0) and upper wall (y/Ly = 1.9/2.0) are plotted in figure 24(a,b), respectively. The detonation flame front first appears at the lower wall, whereas when it propagates in the upper domain and impacts the upper wall, the axial position of the flame front moves quickly, as presented in the temperature contours during the stage from t = 1.65625 to 1.68625 ms. However, the flame propagation in the middle region is more complicated under the influence of the fuel-lean mixture. Notably, there is a deviant velocity oscillation, especially in the lower and middle areas of the flame front. The average velocity in the upper region is gradually decreasing owing to the decrease of detonation overdriven degree. In the middle region flame, the flame propagation velocity decreases rapidly from 2700 to 2200 m s⁻¹ at the developing stage but rapidly increases to 2700 m s⁻¹ during the period of t = 0.1675–1.6825 ms and then gradually decreases to 2200 m s⁻¹. The propagation velocity in the lower region has a similar trend.

Figure 24.

The (a) axial position of the flame front and the (b) corresponding averaged flame front velocity at different y-direction positions near the lower wall (y/Ly = 0.1/2.0), middle region (y/Ly = 1.0/2.0) and upper wall (y/Ly = 1.9/2.0) in case 4.

After that, a more regular oscillation of the averaged velocity occurs with high amplitude, but the average velocity in three different positions is close to 2300 m s⁻¹. Hence, the inhomogeneous mixture increases the instability of the propagation speed. In Han et al.'s (Reference Han, Wang and Law2019) study, a fluctuating average detonation velocity with a different amplitude was also observed in a concentration gradient with a transverse and linear distribution. They reported a much higher velocity oscillation in the fuel-rich region compared with the fuel-lean region. In the current study, an oscillation with high amplitude is also observed in the lower wall mixture, which has a high ER, but a small amplitude appears in the middle region with low ER.

Furthermore, by subtracting the flow velocity of the mixture (around 500 m s⁻¹), the detonation propagation velocity is approximately 1800 m s⁻¹. The theoretical CJ velocity for detonation propagation in the uniform mixture is 1974.95 m s⁻¹, which is calculated by Cantera when using Burke et al.'s (Reference Burke, Chaos, Ju, Dryer and Klippenstein2012) chemical reaction mechanism. Consequently, there is around an 8.8 % deficit in the average propagation velocity compared with the theoretical CJ speed. This velocity deficit is consistent with previous research reported by Kessler et al. (Reference Kessler, Gamezo and Oran2012) and Han et al. (Reference Han, Wang and Law2019). The flame front is, therefore, inclined remarkably under the influence of velocity deficit with different amplitudes due to the inhomogeneous mixture.

3.4.2. Flame front characteristics in detonation quenching and reinitiation

To elucidate the macroscopic phenomena mentioned above, such as the curved flame front and velocity oscillation, figure 25 provides more detailed information on the flame propagation associated with the temperature, pressure, ER and combustion product during the initial time. The combustion ER map also consists of two colours, grey and brown, representing that the fuel and oxidizer are depleted, respectively, as discussed before. At t = 1.65625 ms, a lower ER value (0.15) exists ahead of the flame front in the upper region (box C1), but a more sensitive mixture is generated in box C2 so that there is a transverse concentration gradient. Moreover, a concentration gradient is also formed in the longitudinal direction. Hence, the ER distributions become highly complex, resulting in more propagation dynamic characteristics.

Figure 25.

Snapshots of over-driven detonation propagation rendered by (first column) temperature, (second column) pressure, (third column) ER and (fourth column) product showing the decupled detonation front (case 4). DW – detonation wave, MS – Mach stem, TW – transverse wave, DC – decoupled flame, DR – detonation reinitiation: (a1) t = 1.65625, (a2) t = 1.65875, (a3) t = 1.66250, (a4) t = 1.66500, (a5) t = 1.66750 and (a6) t = 1.67125 ms.

First, an intense Mach stem is formed after the overdriven detonation wave collides with the upper wall, introducing a high-pressure region and intense detonation combustion. The transverse wave is observed subsequently; see figure 25(b 2). When the detonation wave propagates in the lower ER region, a partial flame front is decoupled with a flame and a shock wave, rapidly forming a curved flame and lower-temperature region (box A1). This is also called partial detonation quenching. The current result is in accordance with previous literature when the detonation propagates in a lower ER region from a lean to a leaner gradient. In Boulal et al. (Reference Boulal, Vidal and Zitoun2016) experimental investigation, detonation quenching occurs in an inhomogeneous mixture with a lean-to-leaner distribution. In their research, a reduced number of transverse waves is observed and even the detonation quenching occurs in the fuel-lean mixture, where the ER value is around 0.2. Hence, this value is close to the current research. It is noted that the detonation quenching does not mean combustion extinction. A detonation reinitiation occurs later. As reported by Zheng et al.'s (Reference Zheng, Kaplan, Houim and Oran2019) study, the detonation wave was rapidly decoupled into a shock and a flame when the flame front travelled through an extremely inhomogeneous mixture. Han et al. (Reference Han, Wang and Law2019) also numerically reported the decoupled flame front. They suggested the LSW is decoupled from the chemical reaction layer due to the quite fuel-lean mixture with a lower OH mass fraction. However, in their works, the inhomogeneous mixture has a linear transverse concentration gradient. A more complicated mixture is not considered.

In the current decoupled detonation flame front associated with partial detonation quenching, an accumulation of increased pressure is yielded in this front, as presented by box B1 in figure 25(b 4). When the intense shock wave propagates into these regions with a higher ER value (boxes C3 and C4), the detonation reinitiation is triggered in the partially flame front; see figure 25(b 5). The detonation reinitiation was also reported in the previous experimental study, where the detonation decay could reinitiate in the non-uniform mixture with 14.0 %–19.0 % hydrogen (Kuznetsov et al. Reference Kuznetsov, Alekseev, Dorofeev, Matsukov and Boccio1998). This value is close to the current result. Then, a more uniform flame front is observed since the flame front propagates in a region having high and uniform ER regions. Such a result also contributes to an uneven distribution of the combustion product of H₂O, and a gradually increasing penetration region with a lower H₂O mass fraction is generated, as shown in circles D1–D4. Accompanying the formation of detonation reinitiation, a high mass fraction of H₂O is yielded again in the flame front; see circles D5 and D6. Therefore, the complicated dynamic characteristics, including the curved flame front, partially decoupled flame front and detonation reinitiation, are combined results in the high and low sensitive mixtures along the parallel and transverse directions and transverse wave.

Figure 26 depicts a schematic representation detailing the intricate evolution of detonation wave bending, quenching and subsequent reinitiation. Here, the blue region signifies combustion products, while the blue line traces the path of the detonation wave. Additionally, the red, green and yellow areas correspond to high, low and optimal ERs, respectively. When the detonation wave advances into the low ER region, the combustion expansion fails to sustain a self-sustaining detonation, owing to fuel scarcity. Consequently, this results in partial bending and decoupling of the detonation wave front. However, upon traversing into the optimal ER region represented by yellow, the detonation wave reignites, thereby restoring its flat-fronted propagation.

Figure 26.

Physical illustration of the curved detonation wave, partially decoupled detonation and detonation re-initiation. DW – detonation wave, DQ – detonation quenching, DR – detonation reinitiation.

To draw the influence of an inhomogeneous mixture on the detonation propagation speed, figure 27 gives the instantaneous flame propagation velocity near the (red line) lower and (blue line) upper walls at different periods, corresponding to figure 25(a 2,a 3), respectively. All flame front propagation velocities have a high oscillation within 1750–3500 m s⁻¹ due to the inherent unstable feature of the detonation wave (Austin Reference Austin2003). A rather higher amplitude of the flame propagation velocity in the upper wall is formed compared with the lower flame propagation velocity. This contributes to the intense Mach stem, which increases the intensity of the detonation wave.

Figure 27.

Transient detonation propagation velocity near (red curve) lower wall (y/Ly = 0.1/2.0) and (blue curve) upper wall (y/Ly = 1.9/2.0) in case 4 at different periods (a) t = 1.65875–1.66045 ms, (b) t = 1.6625–1.6642 ms.

3.4.3. Detonation front structure and instability

To compare the differences between the decoupled flame and detonation wave at t = 1.6625 ms, the flame front structure is superimposed in figure 28(c,d), and for lines B (y/Ly = 1.4/2.0) and C (y/Ly = 1.9/2.0) in figure 28(a). A very fuel-lean mixture is formed with a value of H₂ = 0.007. After the shock wave, the H, O and OH increase rapidly, followed by HO₂ and H₂O₂ radicals, while the peak values of H and O radicals are pretty small at H = 0.00033 and O = 0.0122. Subsequently, with the consumption of H₂ and HO₂, less energy release is formed, resulting in a lower increased temperature from 1054 to 1430 K. Hence, the main contribution of the increasing temperature is formed by the decoupled LSW. To thoroughly analyse the connection between the decoupled flame front's shock wave and flame structure, figure 29 gives the temperature, pressure and OH variations along line B. As demonstrated, after the intense post-shock period, the temperature of the fresh mixture is increased rapidly; then a lower increasing ratio appears within the chemical reaction. Here, a rather lower peak mass fraction of OH = 0.0056 appears. Moreover, the distance between the LSW and the reaction front is around 0.4 mm so that the LSW is decoupled from the reaction layer due to the extreme fuel-lean mixture. As a result, the simulation of the extreme fuel-lean and fuel-rich mixtures needs a detailed hydrogen–air chemistry reaction mechanism because of the complex dynamic behaviours, especially for the detonation quenching and partially decoupled flame front.

Figure 28.

The (a) temperature contour, (b) OH contour, and the corresponding mass fraction of the flame front structure near the (c) decoupled detonation flame (line B, y/Ly = 1.4/2.0) and (d) upper wall (line C, y/Ly = 1.9/2.0) at t = 1.6625 ms in case 4.

Figure 29.

The temperature, pressure and OH variations along line B for the decoupled detonation front at t = 1.6625 ms (see figure 28a).

In addition to the detonation front structure along line C in figure 28(d), thanks to the intense shock wave, a rapid increase of mass fraction appears in the O and H radicals, followed by HO₂ and H₂O₂. The peak values of H, O and OH reach 0.00264, 0.04376 and 0.0206, respectively, much higher than the deflagration flame, as discussed in the premixed mixture, especially in the OH radical. Moreover, compared with the deflagration flame, the H and O radicals appear early in the flame front instead of the HO₂ and H₂O₂. As such, the flame front structure is significantly different due to the early LSW. Moreover, the peak value of HO₂ approaches 4.5 × 10⁻⁵ through the chain-breaching reaction of H + O₂(+M) = HO₂(+M), which is almost two times that of the deflagration flame front. This is due to the HO₂ radical being associated with pressure, which is more vital for the detonation flame with high pressure. After that, with the abrupt consumption of HO₂ and H₂O₂, the mass fraction of OH and H₂O radicals is increased accordingly with the main chain-branching reactions through HO₂ + H = OH + OH and HO₂ + OH = H₂O + O₂, resulting in a high energy release ratio, and feeding back to a dramatic increase in temperature that reaches 3200 K.

Furthermore, we analysed the local detonation combustion and product. The profiles of the mass fraction of H₂ and the local Ma number along lines A, B and C at t = 1.6625 ms are given in figure 30(a,b), respectively. Residual hydrogen exists in the downstream flame; see circles A and B, especially in line A, and the H₂ mass fraction reaches 0.012, implying that the H₂ is not consumed fully behind the detonation front. Several previous studies have also found these results in homogenous and inhomogeneous mixtures (Austin Reference Austin2003; Radulescu et al. Reference Radulescu, Sharpe, Lee, Kiyanda, Higgins and Hanson2005, Reference Radulescu, Sharpe, Law and Lee2007; Mahmoudi & Mazaheri Reference Mahmoudi and Mazaheri2015; Han et al. Reference Han, Wang and Law2019). Subbotin (Reference Subbotin1975) and Gamezo, Desbordes & Oran (Reference Gamezo, Desbordes and Oran1999) suggested that the unburned material is highly related to the detonation instability, and the unburned material appears in a highly unstable detonation propagation. In Sharpe's (Reference Sharpe2001) work, they reported that the formation of an unburned material is generated by the departure of the slip line and secondary triple point. Mahmoudi & Mazaheri (Reference Mahmoudi and Mazaheri2015) also found an unreacted mixture when using a high mesh resolution simulation considering the viscosity. More recently, the unburned material was observed in an inhomogeneous mixture (Han et al. Reference Han, Wang and Law2019), and Han et al. (Reference Han, Wang and Law2019) revealed that more unburned pockets after the detonation front are generated when steepening the mixture concentration gradient. In the present research, due to the concentration gradient in both the longitudinal and horizontal directions, the residual hydrogen appears more random and mainly depends on the ER value. At t = 1.6625 ms, the ER value is much higher ahead of the lower detonation front, as noted in figure 25(c 2) near the bottom wall. As a result, a high bulk of residual hydrogen occurs when the denotation wave propagates in these regions.

Figure 30.

Variations of macroscopic variables of (a) H₂ and (b) Mach number along lines A, B and C in figure 28(a).

By insight into the Mach number profile, the average Mach number ahead of the flame front is close to 1.5. A steep increase in Mach number is due to the LSW. After the combustion process, the pressure expands with a decrease in pressure and density. Regarding the Mach number, the ideal ZND solution quickly approaches the local speed of sound, namely, Ma = 1, while in the current inhomogeneous supersonic mixture, it needs more time to reach the sonic point. Such a feature is attributed to the intense flow instabilities generated by the non-equilibrium chemical reaction in the downstream chamber, where a continuous reaction of the unburned pocket appears under the influences of the R-M and K-H instabilities. In Radulescu et al.'s (Reference Radulescu, Sharpe, Law and Lee2007) study, the delayed equilibrium reaction was also reported when comparing the detonation structure and the variations of the sonic surface. By employing the same method, Han et al. (Reference Han, Wang and Law2019) also found that the non-equilibrium reaction is elongated when steeping the concentration gradient. Hence, the present result is qualitatively consistent with the previous study.

To reveal the detonation instability, figure 31 gives the next detonation propagation scenarios from x = 66 to x = 73 cm, accompanied by the pressure, ER and combustion product of H₂O. First, a curved flame with a lower temperature appears due to the fuel-lean mixture (box C1 in figure 31c 1), while a detonation reinitiation is observed when this flame front propagates into a high ER mixture. A remarkable triple point is generated and propagates from the lower wall to the upper wall (circles A1–A3). At t = 1.715 ms, a complicated flame front is generated. A detonation wave is decoupled again (see DC2 in figure 31a 3) to a flame and a shock wave since it propagates in the lower ER region (circle C2), followed by a lower combustion product (circle D2), while a multiple-head detonation wave is generated in the upper flame front; see boxes B1 and B2. Yet, the cell distribution is irregular, and therefore, an irregular detonation wave is formed. At t = 1.7225 ms, the detonation wave propagates downstream, detonation reinitiation in circle D3 and a smooth flame front with a large detonation cellular size are observed (box B3). As a result, the inhomogeneous promotes the detonation instability with an irregular wave. As the flame front travels into the fuel-lean mixture, an increasing induced reaction zone forms in circle A5.

Figure 31.

Snapshots of detonation propagation rendered by (first row) temperature, (second row) pressure, (third row) ER and (fourth row) product showing the partially decupled detonation front and irregular detonation cellular structure in case 4. DC – decoupled flame, TP – triple point, D-R – detonation re-initiation: (a1) t = 1.7, (a2) t = 1.7075, (a3) t = 1.715, (a4) t = 1.7225 and (a5) t = 1.7300 ms.

Figures 32 and 33 further provide detailed detonation wave structures and corresponding density schlieren enlargements for part of the time in figure 31. When the decoupled shock propagates to an appropriate ER area, detonation reinitiation occurs, and part of the detonation wave appears in a low-temperature area downstream of combustion, as shown in the green area in figure 32(a). At the same time, the structure of the detonation wave is complex. The theoretical detonation cell size in the premixed mixture is around 9.0 mm when using the Cantera tool calculation. However, an irregular detonation propagation is generated in the computational domain, where the detonation cell size is small in the upper domain. In comparison, the detonation cell size in the lower half-domain is quite large. This is mainly due to a large transverse concentration gradient in the mixture. In general, the detonation cell size is small compared with the theoretical size. The numerical results of the previous study (Hu et al. Reference Hu, Zhang, Khoo and Jiang2005) also obtained a relatively smaller cell size than the experiment, which is explained by the higher heat loss at the wall in the experiment (Fickett & Davis Reference Fickett and Davis2000). Due to the use of adiabatic wall BCs and the limitations of 2-D simulation, it is still reasonable to obtain a small cell size in the numerical simulation.

Figure 32.

Flame front structure and detonation propagation instability in inhomogeneous supersonic mixture rendered by temperature contour (a) and destiny schlieren (b) at t = 1.7075 ms. TW – transverse wave, TP – triple point, SL – slipping line, MS – Mach stem.

Figure 33.

Flame front structure and partial decouple flame front rendered by temperature contour (a) and destiny schlieren (b) at t = 1.7150 ms. TW – transverse wave, TP – triple point, SW – shock wave, MS – Mach stem.

The temperature distribution of the detonation wave is uneven, and the temperature in the lower half-domain is relatively high. This is mainly because the ER distribution in the longitudinal direction is particularly uneven, leading to different triple point structure distributions and an uneven temperature distribution. The density schlieren in figure 32(b) shows that multiple transverse waves are generated in a smaller range downstream of the detonation front. In contrast, high-intensity transverse waves extend to a long area downstream in the entire combustion chamber.

Figure 33 further shows the detonation extinction and complexity of the detonation wave structure. At this moment, part of the detonation wave is extinguished due to a lack of fuel near the lower wall, resulting in a low-temperature area on the flame front and downstream. But the decoupled shock wave still propagates forward. Meanwhile, multiple triple points and many transverse waves are generated in the upper part. In contrast, the stronger transverse waves in the combustion products still propagate up and down in the combustion chamber. It can, therefore, be concluded from these two moments that, in an inhomogeneous supersonic mixture, the instability and complexity of detonation wave propagation are formed due to inhomogeneous ER distribution.