1. Introduction
Unstart is an undesirable phenomenon that may occur in high-speed intakes during which the shock system gets disgorged and may exhibit self-sustained oscillations (Wagner et al. Reference Wagner, Yuceil, Valdivia, Clemens and Dolling2009; Chang et al. Reference Chang, Li, Xu, Bao and Yu2017; Im & Do Reference Im and Do2018). These oscillations, known as ‘buzz’, have been reported and studied widely in supersonic intakes by various researchers (Ferri & Nucci Reference Ferri and Nucci1951; Dailey Reference Dailey1955; Sterbentz & Davids Reference Sterbentz and Davids1955; Fisher, Neale & Brooks Reference Fisher, Neale and Brooks1970; Nagashima, Asanuma & Obokata Reference Nagashima, Asanuma and Obokata1972; Van Wie, Kwok & Walsh Reference Van Wie, Kwok and Walsh1996; Trapier, Duveau & Deck Reference Trapier, Duveau and Deck2006; Soltani & Sepahi-Younsi Reference Soltani and Sepahi-Younsi2016). The majority of experiments on intake unstart simulate an increase in backpressure by using mechanical blockage at the end of the isolator to throttle the flow, and the resulting flow system resembles a forward-facing cavity (shock oscillations ahead of a forward-facing cavity were studied by Engblom et al. Reference Engblom, Yuceil, Goldstein and Dolling1996); hence the quarter-wave resonator model ($f = a_0/4l$, where $l$ is the geometric length between the cowl lip and the throat) has been proposed as a good model for scaling the frequency of the flow oscillations (Hankey & Shang Reference Hankey and Shang1980). Furthermore, Newsome (Reference Newsome1984) has also reported a similar expression to estimate the frequency of unstart based on average duct Mach number ($M$), local speed of sound ($c$) and characteristic length ($l$), which is the distance between the cowl lip and the throat ($f= (2n+1)c(1-M^2)/4l$). The quarter-wave resonator model is a simple method and is able to estimate buzz frequency reasonably well in supersonic intakes. However, in hypersonic intakes, there is no consensus on the validity of this model. Tan, Shu & Zhi-Long (Reference Tan, Shu and Zhi-Long2009) have used the model given by Newsome (Reference Newsome1984) and have reported discrepancies with the experimentally observed frequencies in a hypersonic intake. On the other hand, recently, Berto et al. (Reference Berto, Benini, Wyatt and Quinn2020) and the authors (Devaraj et al. Reference Devaraj, Jutur, Rao, Jagadeesh and Anavardham2020) have reported intake unstart at hypersonic Mach numbers, and have observed that the frequencies of oscillation are in reasonable agreement with the predictions based on the quarter-wave resonator model. The geometric differences between supersonic and hypersonic intakes – i.e. the fact that in the former, the length of the forebody is an order of magnitude lower than the distance between cowl lip and throat, while in the latter the two quantities are comparable (as shown in figure 1) – entail different flow physics. Most previous researchers have used the geometric length as the characteristic length scale, even for hypersonic intakes. We feel that the lack of consensus is mainly due to this extrapolation of the length scale from the supersonic intakes, and that a length scale based on the flow physics can extend the validity of the quarter-wave resonator model to hypersonic intakes as well.
Predicting the mode of unstart is as important as characterizing its frequency. In supersonic intakes, the flow oscillations (buzz) observed during unstart can be distinguished based on the flow features that drive these oscillations. Ferri & Nucci (Reference Ferri and Nucci1951) have observed that the shock intersection gives rise to a vortex sheet, which separates the flow at the cowl surface upon impingement and initiates the buzz. In contrast, Dailey (Reference Dailey1955) observed the buzz to be initiated by the ramp-side separation caused by the shock–boundary-layer interaction (SBLI). Fisher et al. (Reference Fisher, Neale and Brooks1970) report both types of instabilities at the same frequency but with significantly different amplitudes. Most of the studies on unstart in hypersonic intakes report oscillatory behaviour driven by ramp-side separation (Tan et al. Reference Tan, Shu and Zhi-Long2009; Li et al. Reference Li, Gao, Jang and Yang2013). Zhang et al. (Reference Zhang, Tan, Chen, Yuan and Zhang2016) have reported both types of oscillatory behaviour in a hypersonic mixed-compression intake at different Mach numbers. While the literature has distinguished the oscillatory flows based on the driving flow features, there is no simple criterion on the geometrical design or the operating parameters that distinguish between the two modes of unstart. In considering several previous experimental studies, we noticed that the location of the external shock system with respect to the cowl lip (shown in figure 2), as observed in schlieren images, has a significant bearing on the flow features during unstart. This aspect has not been given sufficient attention in previous studies.
Motivated by these two important observations, in the present work we carry out experiments to study unstart induced by throttling in a hypersonic mixed-compression intake at different internal contraction ratio ($ICR$) conditions, which in turn correspond to different external compression shock locations with respect to the cowl lip. Our experiments aim to identify an appropriate length scale for frequency estimation at high throttling ratios (in non-reacting flows) and a criterion for distinguishing the mode of oscillatory unstart. These conditions mimic the high backpressure generated due to supersonic combustion and represent the worst-case scenario that can occur in real flight, where sustained oscillations may occur.
2. Methodology
The present work considers a mixed-compression intake consisting of external compression through two ramps at inclinations of $10.7^{\circ }$ and $25.6^{\circ }$ with the horizontal, and a constant-area isolator with a height ($h$) of 6.2 mm, as shown in figure 3. The sidewalls in the isolator and the cowl region are provided with BK7 glass to enable optical access for schlieren. The intake model is equipped with a movable cowl as well as a flap located at the end of the isolator to independently vary the internal contraction ratio ($ICR$) and throttling ratio ($TR$), which are defined in (2.1a,b), where $A_{is}$, $A_{th}$ and $A_{i}$ represent the cross-sectional area of the isolator, throat and area at the cowl lip, respectively:
The experiments are conducted at a Mach number of $6 \pm 0.1$ in a hypersonic wind tunnel at the Indian Institute of Science. The facility is of blow-down type with air stored at high pressure on one end, driving the flow to vacuum on the other, providing a steady test time of 5 s. The current experiments are carried out at a stagnation pressure of $7.55 \pm 0.04 \times 10^5$ Pa and a stagnation temperature of $423.80 \pm 9.47$ K. Accordingly, a flight-realistic Reynolds number of $8.73 \pm 0.19\times 10^6$ m$^{-1}$ has been achieved.
The experiments are conducted on the intake model by varying $TR$ during the steady test time, from a starting position of $TR = 0$ to the final $TR$, for different $ICR$ conditions. Three different $ICR$ values, 1.19, 1.28 and 1.37, are considered. A nominal condition corresponding to a case without throttling ($TR=0$) and a throttling condition at $TR = 0.56$ are considered at all three $ICR$ values. In addition, a high $TR$ of 0.68 is also considered for the $ICR$ values of 1.19 and 1.37 to investigate the effect of a higher $TR$ at different $ICR$ conditions. The intake model is instrumented with eight flush-mounted Kulite sensors, as shown in figure 3, capable of measuring unsteady pressure with an uncertainty of 2 %. The unsteady pressures from the Kulite sensors are acquired at a sampling rate of 100 kHz. Simultaneous schlieren visualization is also carried out for each experiment. A nano-pulsed laser is used as the light source (with 10 ns pulse duration) and is synchronized with the Photron SA4 high-speed camera to acquire images at 10 000 fps. Further details of the facility and instrumentation are given in Devaraj et al. (Reference Devaraj, Jutur, Rao, Jagadeesh and Anavardham2020). The dynamic content associated with the flow is analysed with the help of fast Fourier transform (FFT) of the unsteady pressure signals as well as dynamic mode decomposition (DMD) of the schlieren image sequence. The methodology used for the DMD analysis is given in Rao & Karthick (Reference Rao and Karthick2019).
3. Results
3.1. Started flows
The typical flow features observed during a started flow through the intake are presented in figure 4, which corresponds to the nominal condition of $TR = 0$ at different $ICR$ values. Two oblique shocks, OS1 and OS2, emanate from the external compression ramps. The flow turns back into the isolator through cowl shocks CS1 and CS2 and remains supersonic through the isolator region across all three $ICR$ conditions. The change in $ICR$ leads to variation in the strength of the cowl shocks CS1 and CS2, and accordingly, differences are observed in the isolator flowfield. Figure 4(d) shows the non-dimensional pressure profile along the isolator for different $ICR$ conditions. Peak pressures of $\approx$50 times the free-stream pressure in the isolator are observed at $ICR = 1.19$ and $ICR =1.28$, while the same at $ICR = 1.37$ is $\approx$60 times the free-stream pressure. It is important to note that a variation in $ICR$ leads to a change in impingement location of the external compression shock on the cowl surface. At $ICR = 1.37$ and $ICR = 1.28$, the external compression shock impinges on the cowl surface, while it misses the cowl lip for $ICR = 1.19$.
3.2. Unstarted flows
The flow features of the intake are modified significantly with the introduction of throttling. For the scope of the present work, the focus is on the final flow features after achieving the desired throttling. In the present study, unstarted flows are observed at $TR \geq 0.56$ for all three $ICR$ cases. The time-averaged schlieren images corresponding to $ICR$ cases of 1.37 and 1.28 at $TR = 0.56$ are shown in figure 5. A detached cowl shock (DCS) with a type V shock interaction, forming a supersonic jet SJ that impinges onto the ramp side and terminates very close to the cowl leading edge, can be observed. The movement of this shock system ahead of the cowl is not significant, and hence a time-averaged image is presented. The flow structure comprising the DCS and SJ is observed farther downstream at $ICR = 1.28$ than at $ICR = 1.37$. The FFT of the pressure signal and DMD of the schlieren visualizations do not show any dominant frequencies associated with the flow features (see figure 6a) for $TR = 0.56$ at $ICR = 1.37$, with a peak pressure of approximately 200 times the free-stream pressure. Similar behaviour is observed for the case of $TR = 0.56$ and $ICR = 1.28$, which is not shown here. However, at a higher $TR$ of 0.68, a dominant peak at 3686 Hz corresponding to a low-amplitude oscillatory unstarted flow is observed in both the FFT of the pressure signal and the DMD of the schlieren images. However, the flow features are not altered significantly compared to those of $TR = 0.56$.
In contrast, at545 $ICR = 1.19$, high-amplitude oscillatory unstarted flow is observed for $TR \geq 0.56$. Figure 7 presents the sequence of instantaneous schlieren images at $TR = 0.56$ for this $ICR$. The flow features in common with $ICR \geq 1.28$ are OS1, OS2, DCS and SJ. The two main differences observed at this $ICR$ compared to the other cases are that (a) the shock system is being pulled downstream closer to the cowl leading edge, and (b) the ramp-side separation bubble (SLE) is seen upstream of the point of shock interaction (at $x/h=14$). The frame corresponding to the most upstream location of the expelled shock system is taken as the reference time $t_{ref}$ for this case. At $t = t_{ref} + 0.4$ ms, the shock system has moved downstream, behind the sidewall, and SJ is seen near the cowl leading edge and enters the isolator. In this instance, a high-pressure region is created near the exit of the isolator, and an upstream motion of the shock system is initiated. The subsequent frames corresponding to $t = t_{ref} + 0.8$ and $t_{ref} + 1$ ms show the upstream movement of the flow, thereby completing the cycle. A plot of the pressure signal at sensor location B5, the corresponding FFT of the pressure signal and the DMD of the schlieren images are shown in figure 8. A dominant peak is observed at a frequency of 950 Hz at $TR = 0.56$. On the other hand, $TR = 0.68$ reveals a dominant peak at 1100 Hz with an increase in the amplitude of pressure oscillation. In both these cases, peak pressures of about 200 times the free-stream pressure are observed, with high-amplitude oscillation in comparison to the case $ICR = 1.37$ (figure 6). However, the flow features corresponding to this higher-$TR$ case, which were already reported in Devaraj et al. (Reference Devaraj, Jutur, Rao, Jagadeesh and Anavardham2020), are similar to those of $TR = 0.56$. At $ICR = 1.19$, the separation bubble present on the ramp side deflects SJ away from the ramp wall and allows it to enter the isolator. The entrainment observed at $ICR = 1.19$ is driven by the ramp-side separation, similarly to the observations made by Dailey (Reference Dailey1955); correspondingly, the amplitude of the oscillations is also found to be relatively high. The low-amplitude oscillations observed at $ICR = 1.37$ are characteristic of the mechanism proposed by Ferri & Nucci (Reference Ferri and Nucci1951). According to Ferri & Nucci (Reference Ferri and Nucci1951), the shear layer emanating from the point of the shock interaction is responsible for the instabilities during the unstart.
3.3. A criterion distinguishing between the two modes of unstart
Figure 9(a) shows a plot representing the higher-throttling-ratio conditions considered in the present work within the purview of Kantrowitz and isentropic limits of intake operation. The Mach number at the entrance of the internal compression section during the started operation ($M_i$) for various $ICR$ conditions is calculated based on two-dimensional Reynolds-averaged Navier–Stokes simulations using the open-source suite SU2 Falcon (Palacios et al. Reference Palacios2013). At $ICR = 1.19$, the forebody shock misses the cowl. An increase in $ICR$ corresponds to the impingement of the forebody shock onto the cowl surface, thereby reducing $M_i$ for the case of $ICR = 1.37$. Because of this reduction in $M_i$ for the same throttling ratio ($TR = 0.68$), the $ICR = 1.37$ falls below the isentropic limit, while the same for $ICR = 1.19$ is above the isentropic limit. Consequently, high- and low-amplitude oscillatory unstarted flows are observed at $ICR = 1.19$ and $ICR = 1.37$, respectively. Figure 9(b) presents a comparison of the non-dimensional location of shock impingement with free-stream Mach number from various studies, including the present one, where the effect of throttling has led to unstart. The $\varDelta /h$ values are obtained from the schlieren images reported in the respective studies, with $\varDelta$ denoting the distance of the shock impingement from the cowl lip in the direction normal to the free-stream flow during the started operation. When $\varDelta /h > 0$ (the shock misses the cowl), the boundary-layer separation at the ramp side is identified as the driving flow feature (Dailey criterion), and a high-amplitude oscillatory unstart is observed. High-amplitude oscillatory unstarts have been reported by Tan et al. (Reference Tan, Shu and Zhi-Long2009) and Li et al. (Reference Li, Gao, Jang and Yang2013), where the intake operates with $\varDelta /h > 0$. On the other hand, the shear layer emanating from the point of shock interaction is observed to be the driving flow feature (Ferri criterion) for the low-amplitude oscillatory unstart when $\varDelta /h < 0$ (the shock impinges on the cowl surface). From figure 9(b), it is evident that $\varDelta /h = 0$, which corresponds to shock-on-lip (SOL) condition during the started operation, is the boundary between the Dailey and Ferri criteria for hypersonic mixed-compression intakes operating with free-stream Mach numbers between 5 and 7. It is important to note that SOL here corresponds to that observed from the experimental schlieren images. The actual SOL criterion, therefore, can be useful in identifying the driving flow features and type of oscillations that can occur during unstart. The appendix gives additional supporting evidence to emphasize that the actual SOL is the distinguishing criterion.
3.4. Buzz frequency scaling
In our previous work (Devaraj et al. Reference Devaraj, Jutur, Rao, Jagadeesh and Anavardham2020), we have shown that the flow features in the intake differ significantly between static and dynamic flap experiments in lower-$TR$ cases. For example, the intake remains in the started condition in a dynamic flap experiment, while it is unstarted during a static flap experiment at $TR = 0.3$. Moreover, in a flight scenario, the combustion process is initiated after the establishment of the started flow through the intake. Hence, a dynamic throttling would closely mimic the practical scenario at lower-$TR$ conditions, where there is a difference between static and dynamic flap experiments. On the other hand, the oscillatory unstarted flows observed at higher $TR$ are independent of the mode of flap operation. In addition, Tan et al. (Reference Tan, Shu and Zhi-Long2009) and Li et al. (Reference Li, Gao, Jang and Yang2013) have shown that the frequency remains relatively unchanged beyond a particular throttling ratio. Hence, a scaling analysis for the limiting case of $TR \to 1$ will be applicable independent of the mode of flap operation.
It is evident from the FFT and DMD plots that the unstarted flow exhibits self-sustained oscillations at a well-defined dominant frequency at $TR = 0.68$. In these high-$TR$ conditions ($TR \to 1$), the blockage imposed by the flap gives rise to a configuration similar to that of a forward-facing cavity; accordingly, we propose a suitable length scale ($l^*$) that is defined as the extent of the subsonic region in the unstarted flow. In the present study, during low-amplitude oscillatory unstarted flow, the shear layer emanating from the point of shock interaction terminates near the cowl to form a virtual forward-facing cavity. Hence, the distance between the point of shock interaction and the impingement of the shear layer on to the cowl surface is an appropriate length scale (figure 10a(i)). On the other hand, during the high-amplitude oscillatory unstarted flow, the appropriate length scale would be the distance between the leading edge of the separation bubble and the throat, which is formed at the exit of the isolator (figure 10a(ii)). With this choice of the appropriate length scale in each of these cases, the frequency estimated using the quarter-wave resonator model matches the experimental values reasonably well. To further strengthen the validity of this scaling, table 1 compares the frequency values at the highest $TR$ from different experimental studies of mixed-compression intakes against theoretical estimates. For intakes that exhibit unstart with supersonic spillage, the separation region extends up to the leading edge of the forebody, with subsonic flow prevailing over the ramp side as well as the isolator (Tan et al. Reference Tan, Shu and Zhi-Long2009; Li et al. Reference Li, Gao, Jang and Yang2013). Hence the appropriate length scale in this case would be the axial distance between the throat and the leading edge of the intake (figure 10a(iii)). Figure 10(b) represents the experimentally observed frequency values normalized with theoretical estimates for different $TR$ ($TR > 0.5$) over a wide range of Mach numbers (ranging from $M_{\infty } = 2$ to $M_{\infty } = 6$), from both previous literature and the present study. It is interesting to note that the assumption of similarity to a forward-facing cavity leads to reasonably accurate frequency estimates (with a maximum difference of $\approx$20 %), even for $TR$ values as low as $TR = 0.55$. It may be noted that some researchers (Dailey Reference Dailey1955; Rodi, Emami & Trexler Reference Rodi, Emami and Trexler1996; Trapier et al. Reference Trapier, Duveau and Deck2006) have reported higher harmonics of the fundamental mode ($f = a_0/4l$) to be dominant in supersonic intakes, the reasons for which are unclear. Furthermore, when unstart is induced by combustion heat release (as in Im et al. Reference Im, Baccarella, McGann, Liu, Wermer and Do2016), chemical reactions may play a significant role, and identifying the exact length scale may be difficult. Detailed studies using experiments and high-fidelity numerical simulations are required to clarify the dynamics of shock oscillations in the context of chemically reacting flows.
4. Conclusion
Collating data corresponding to high-throttling-ratio cases from the current and previous studies pertaining to mixed-compression intakes, two modes of unstart can be identified. High-amplitude oscillatory unstarted flow is driven by the separation bubble on the ramp side (Dailey criterion). Low-amplitude oscillatory unstarted flow is driven by the shear layer emanating from the shock interactions (Ferri criterion). The two significant conclusions of this work are as follows: (a) the actual SOL (corresponding to the unthrottled condition) is the demarcating boundary between the two modes of oscillatory unstart observed upon high throttling; (b) the extent of the subsonic region (as observed from the experimental schlieren images) is the appropriate length scale to be used in the quarter-wave resonator model.
Supplementary movies
Supplementary movies are available at https://doi.org/10.1017/jfm.2021.230.
Acknowledgements
The authors would like to thank BrahMos Aerospace Private Limited, New Delhi, and an Indian Institute of Science Startup Grant for funding this research work. The authors would also like to thank Shivashankar and Shantha Kumar for their help in operating the 0.5 m Hypersonic Wind Tunnel, Murthy and Raju for their help in fabrication, Parthiban (Novatek Automation) for designing the mechanism of the flap and the cowl, Ramji V. for his help in setting up the apparatus for schlieren visualization and members of the Laboratory for Hypersonic and Shockwave Research at the Indian Institute of Science. The first author would like to acknowledge discussions with A. Bajpai and Dr T. Singh.
Declaration of interests
The authors report no conflict of interest.
Appendix
To further emphasize that the actual SOL is the demarcating criterion for the two modes of unstart, additional evidence corresponding to representative high-$TR$ cases at $ICR = 1.61$ and $ICR = 1.85$ is shown in figures 11 and 12, which correspond to $\varDelta /h < 0$. For these $ICR$ cases, the flow features in the unthrottled conditions, as shown in figure 11(a,c), are significantly different from those of the cases shown in § 3.1. However, the flow features upon throttling (figure 11b,d) are similar to those shown in figure 5. Moreover, the magnitude of pressure and the dynamic content as seen in figure 12 are similar to those presented in figure 6.