J. Fluid Mech. (2022), vol. 949, A37, doi:10.1017/jfm.2022.769 Transitional hypersonic flow over slender cone/flare geometries Cameron S. Butler1 and Stuart J. Laurence1,† 1Department of Aerospace Engineering, University of Maryland – College Park, College Park, MD 20742, USA (Received 15 December 2021; revised 23 June 2022; accepted 16 August 2022) Experiments are performed in a Mach-6 shock tunnel to examine the laminar-to-turbulent transition process associated with a sudden increase in surface angle on a slender body. A cone/flare geometry with a 5◦ frustum and compression angles ranging from 5◦ to 15◦ allow a range of mean flow configurations, spanning an attached shock-wave/boundary-layer interaction to a fully separated one; the unit Reynolds number of the flow is also varied to modify the state of incoming second-mode boundary-layer disturbances. Ultra-high-speed schlieren visualizations provide a global picture of the flow development, supplemented by high-frequency surface pressure measurements. For the 5◦ compression, the unsteady flow field is dominated by the second-mode waves, whose breakdown to turbulence is generally accelerated (compared with the straight-cone configuration) by encountering the angle change. As the compression angle is increased to induce separation, lower-frequency disturbances appear along the separated shear layer that exhibit much larger amplification rates than the incoming second-mode waves; the latter effectively freeze in amplitude downstream of the separation point before rapidly breaking down upon reattachment. The shear-layer disturbances become dominant at the largest compression angle tested. Radiation of disturbance energy to the external flow is consistently observed: this generally occurs along mean flow features (flare, separation or reattachment shocks) for the second-mode disturbances and spontaneously for the shear-layer waves. The combined application of spectral proper orthogonal decomposition and a global bispectral analysis allows the identification of important unsteady flow structures and the association of these with prominent nonlinear interactions in the various configurations. Key words: compressible boundary layers, hypersonic flow, transition to turbulence † Email address for correspondence: stuartl@umd.edu © The Author(s), 2022. Published by Cambridge University Press. This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/ licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited. 949 A37-1 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss mailto:stuartl@umd.edu http://creativecommons.org/licenses/by/4.0 http://creativecommons.org/licenses/by/4.0 http://crossmark.crossref.org/dialog?doi=https://doi.org/10.1017/jfm.2022.769&domain=pdf https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence 1. Introduction The design of practical hypersonic vehicles is constrained by the extreme thermo- mechanical surface loads which occur when travelling within the stratosphere at high Mach number. This is further complicated by the dramatic increase in surface heat flux and skin friction which accompanies the laminar-turbulent transition of the vehicle boundary layer. Due to incomplete understanding of the hypersonic transition process and inadequate modelling capabilities, these issues are presently avoided by over-designed thermal protection systems; however, these come with significant weight penalties, placing undesirable constraints on the vehicle and mission. Safe and efficient vehicle design thus requires advances in our understanding of the underlying physics of hypersonic boundary-layer transition. In the low-disturbance environments typical of hypersonic flight (and sufficiently quiet ground-test facilities), boundary-layer transition on slender, smooth bodies is characterized by the linear growth of instabilities within the boundary layer, leading up to nonlinear modal interactions and breakdown (Fedorov 2011). For axisymmetric or near-axisymmetric geometries, the dominant instability mechanism above Mach numbers of approximately four is the Mack or second mode. Following its discovery (Mack 1975), the second mode has received substantial attention, both experimentally (for example, Demetriades 1974; Stetson & Kimmel 1992; Laurence et al. 2012; Laurence, Wagner & Hannemann 2014, 2016; Casper et al. 2016; Kennedy et al. 2018; Craig et al. 2019) and theoretically/computationally (for example, Fedorov & Tumin 2011; Sivasubramanian & Fasel 2014; Unnikrishnan & Gaitonde 2020), with studies tending to focus on smoothly varying surface geometries. The outer mould line of true flight vehicles does not always vary smoothly in the streamwise direction however, and may exhibit a sudden increase in angle, e.g. at a control surface or intake. Such abrupt angle changes at supersonic conditions will introduce a shock wave and accompanying shock-wave/boundary-layer interaction (SWBLI). These SWBLIs can give rise to new instability mechanisms, particularly in the separated case – see, for example, Roghelia et al. (2017), Guiho, Alizard & Robinet (2016), Sidharth et al. (2018) – but would also be expected to affect disturbances propagating from the upstream boundary layer, potentially leading to complex unsteady interactions. In the present work, we focus on the nominally two-dimensional interactions produced by a sudden increase in angle of a sharp, slender cone, i.e. resulting in a cone/flare geometry. Much of the prior work on two-dimensional SWBLIs has sought to characterize transitional effects on flow topology (e.g. separation length and unsteadiness) and thermo-mechanical surface loading (Heffner, Chpoun & Lengrand 1993; Benay et al. 2006; Running et al. 2018). An attempt to measure the disturbances generated within a SWBLI was made by Benitez et al. (2020), who employed focused laser differential interferometry (FLDI) to interrogate the separated boundary layer on a cone-cylinder-flare model. Low-frequency (50–170 kHz) travelling waves were identified downstream of the compression corner under quiet flow conditions, but could not be located within the separation region. Point-like measurement techniques such as FLDI, however, necessarily preclude a global view of instability development. Only a limited number of studies have elucidated transition dynamics or the impact of the SWBLI on pre-existing disturbances. For example, Balakumar, Zhao & Atkins (2002) used both linear stability theory and direct numerical simulation (DNS) to study two-dimensional, fixed-frequency disturbances in a flat-plate boundary layer encountering a 5.5◦ compression corner at Mach 5.4. Linear stability theory revealed the existence of multiple unstable modes within narrow regions of the separation bubble, while 949 A37-2 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow DNS showed the second-mode waves to be of neutral stability while traversing the separated shear layer but to grow exponentially upstream of separation and downstream of reattachment. The unstable low-frequency mode within the separation region was found to have a frequency of 38 % of the dominant second-mode frequency. This low-frequency mode was also shown to emanate energy from the separated boundary layer at a point just above the corner. Novikov, Egorov & Fedorov (2016) carried out DNS of three-dimensional, broad-spectrum wavepackets on this same configuration. Both oblique-wave- and second-mode-dominated wavepackets were examined: the latter were found to be neutrally stable within the upstream part of the separation region before amplifying downstream. Strong forcing resulted in significant streamwise stretching of the wavepacket tail and the formation of a turbulent spot downstream of reattachment. Recently, a full transition scenario for a laminar boundary layer encountering an axisymmetric 15◦ compression ramp at Mach 5 was computed by Lugrin et al. (2020) using ‘quasi-DNS’. White-noise forcing was used to excite a range of convective instabilities, and spectral proper orthogonal decomposition (SPOD) was then applied to the unsteady results to identify key flow structures. These authors observed a transition process dominated by streamwise streaks resulting from the nonlinear interaction of oblique first-mode waves. The shear layer and reattachment regions triggered linear amplification of these streaks, ultimately leading to breakdown. Finally, the present authors (Butler & Laurence 2021b) experimentally examined the transitional Mach-6 flow over a slender, 5◦ cone/flare configuration with a 10◦ angle increase, sufficient to generate a limited separation bubble at the conditions tested. Both second-mode disturbances and lower-frequency shear-layer instability waves along the boundary of the separated region were observed with an ultra-high-speed schlieren visualization technique and high-frequency pressure measurements. The second-mode waves were seen to radiate energy along the separation shock when incident upon it; for lower Reynolds numbers, the separation region inhibited second-mode growth, while at higher Reynolds numbers, rapid breakdown occurred near reattachment. In the present work we extend this earlier analysis to a range of flare angles, with compression angles from 5◦ to 15◦ examined. This choice of angles encompasses a range of mean flow fields from fully attached through to substantially separated, while limiting the deflection to within a range that one might realistically expect to encounter on a practical hypersonic vehicle. The global diagnostic and analysis techniques introduced in our earlier work (Butler & Laurence 2021a,b) are employed to examine both the behaviour of incoming second-mode disturbances as they traverse the SWBLI and the instability mechanisms intrinsic to the SWBLI itself. In § 2 we describe the experimental facility, test article and diagnostics used. In § 3 we provide a primarily qualitative overview of the observed fluid phenomena, then in § 4 a spectral analysis derived from the schlieren measurements is presented. Results from a modal-reduction technique, SPOD, are described in § 5, which are then used to inform a bispectral analysis of the nonlinear interactions present in § 6. Conclusions are drawn in § 7. 2. Experimental methodology 2.1. Facility overview All experiments were performed in HyperTERP, a small-scale reflected shock tunnel operated by the University of Maryland. A contoured Mach-5.95 nozzle with a 22 cm exit diameter was employed, exhausting into a 30.5 cm diameter free-jet test section equipped with 15.2 cm diameter windows. A more detailed description of the facility is 949 A37-3 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence 5 10 15 Time (ms) 0 0.5 1.0 1.5 2.0 P 0 (M Pa ) Re33 Re45 Re52 Figure 1. Sample reservoir pressure traces at each condition. Condition Rem (106 m−1) P0 (bar) T0 (K) U∞ (m s−1) P∞ (Pa) T∞ (K) ρ∞ (kg m−3) δc (mm) Re33 3.33 10.0 890 1263 630 110 0.020 1.81 Re45 4.49 13.5 890 1263 851 110 0.027 1.56 Re52 5.24 15.8 890 1263 996 110 0.032 1.44 Table 1. Typical HyperTERP test conditions for this investigation. given in Butler & Laurence (2019). The total specific enthalpy was held approximately constant at 0.89 MJ kg−1, corresponding to a temperature of 890 K, with unit Reynolds numbers, Rem = ρ∞U∞/μ∞, of 3.33, 4.49 and 5.24 × 106 m−1 achieved by varying the reservoir pressure (viscosity is calculated here using Sutherland’s law). Reference reservoir (subscript 0) and corresponding free stream properties (subscript ∞) for each condition are detailed in table 1; the methodology used to calculate these conditions and an overview of the free stream characterization are given in Butler & Laurence (2021b). Typical reservoir pressure traces for all conditions are provided in figure 1. In each case, the steady test time is approximately 6 ms, during which the unsteadiness in pressure is of the order of 2 % (standard deviation). Note that runs at condition Re52 are often affected by spikes in pressure (e.g. at t ≈ 6 ms in figure 1) which must dissipate before data reduction. Shot-to-shot variation in the mean pressure is of the order of 1.3 %, while systematic uncertainty (combined calibration and nonlinearity) is estimated as 1.6 %. Shock-speed estimates are accurate to ±5 m s−1, contributing 0.4 % uncertainty in the reservoir temperature. 2.2. Test article and instrumentation The test article for this study was a cone/flare model comprising a 5◦ half-angle, stainless steel frustum of 410 mm length and interchangeable, 76.2 mm long, Delrin flares. Flare half-angles of 5◦, 10◦, 15◦ and 20◦ were employed, corresponding to a straight continuation of the frustum surface and compression angles of +5◦, +10◦ and +15◦. The model with the +10◦ configuration is pictured in figure 2. The nose radius was measured to be 0.10 mm using a SmartScope optical gauge. The straight-cone configuration was included to provide comparisons against an undisturbed boundary layer. 949 A37-4 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow 410 76,2 Blank Blank U1 U2 U3 U4 D1 D2 D3 D4 D5 Figure 2. Model schematic for the +10◦ configuration showing the sensor layout; all dimensions are in millimetres. Station U1 U2 U3 U4 D1 D2 D3 D4 D5 s (mm) 326 353 381 399 415 422 434 447 460 Table 2. Surface coordinates of PCB stations. A single streamwise row of PCB 132B38 pressure transducers was installed along the upper surface of the cone, held in place using clear nail polish (Ort & Dosch 2019). The sensor locations are indicated in figure 2, and corresponding distances along the cone surface from the nosetip are provided in table 2. These sensors were sampled at 2 MHz with a 600 kHz low-pass filter to remove aliasing effects. The factory-supplied calibration was used in converting the voltage signals to pressure, though this does not account for the varying frequency response of the sensors (Ort & Dosch 2019). It should additionally be noted that Ort & Dosch (2019) observed parasitic resonances at frequencies above 300 kHz for this sensor model and corresponding peaks were observed in many of the spectra in the present experiments; results near this frequency should thus be treated with caution. To analyse the pressure-disturbance data, average power spectral densities (PSDs) were computed over the steady test duration using Welch’s method with Hann windows of length 128. Disturbance N factors, i.e. the spatial integral of the amplification rate, could be computed for each frequency according to �N( f , xi) = 1 2 ln ( PSD( f , xi) PSD( f , x0) ) , (2.1) where PSD( f , xi) refers to the power of frequency f at streamwise location xi and x0 refers to the most upstream station. Similarly, the maximum second-mode N factor could be computed by instead considering only the peak PSD within the frequency range of the second mode at each station. The model was initially installed at approximately zero incidence (pitch and yaw) with the aid of a laser level by aligning the model seam and PCB array with the horizontal and vertical centrelines of the nozzle. To further refine the model pitch, measurements were also performed with an additional PCB sensor installed on the underside of the model opposite station U1, as second-mode frequencies are known to be highly sensitive to angle-of-attack variations. From these measurements, it was determined that the pitch angle was less than 0.2◦, with any residual offset corresponding to the upper ray of the model (where measurements were performed) lying on the flow-leeward side. Although no comparable measurements were made regarding model yaw, small yaw 949 A37-5 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence variations are expected to have a negligible effect on the flow behaviour in the region of interest. 2.3. Calibrated schlieren High-speed flow field visualization, obtained using a standard Z-type schlieren with a horizontal knife edge, serves as the primary means of flow interrogation throughout this work. Light pulses of 20 ns duration were provided by a Cavilux HF laser and a Phantom v2512 camera was used to capture the images at frame rates between 440 kHz and 822 kHz (the lower frame rate being used for the larger flare angle), allowing resolution of spectral content up to 411 kHz. The frame rates were chosen to be higher than the Nyquist sampling rate of the dominant boundary-layer disturbances (typically second-mode waves). The field of view ranged from 512 × 32 pixels at 822 kHz to 640 × 64 pixels at 440 kHz, with the camera rotated to maximize the region of flow visualized. The magnification of the optical set-up resulted in a length scale of 0.115 mm pixel−1 for the +15◦ configuration and 0.139 mm pixel−1 for all others. The theoretical, undisturbed boundary-layer thickness at the corner junction, δc, ranged from 1.4–1.8 mm depending on the condition (particular values are provided in table 1). This gives a visualized boundary-layer resolution of at least 10 pixels, though it should be noted that this number will be lower far upstream of the corner and along the compression flares due to the reduced boundary-layer thicknesses there. The pixel intensities in each image were converted to integrated density gradients using the calibration procedure described by Hargather & Settles (2012) and Kennedy (2019). A plano-convex, spherical lens (focal length 10 m, diameter 25.4 mm) with a known deflection angle profile was placed into the field of view and a reference image captured. The intensity profile along the lens face was then mapped to the known deflection profile, assuming zero deflection to correspond to 92 % of the background pixel intensity (to account for the 8 % absorption specific to the calibration lens). For a weak lens, the deflection angle, ε, is given as a function of radial distance from the lens centre, r, by r/f = tan ε ≈ ε. (2.2) The deflection profile is then mapped to the density gradient of the test gas according to ε = κL n∞ ∂ρ ∂y , (2.3) where L is the integration path length of the light, κ is the Gladstone–Dale constant, ρ refers to the gas density, y is the direction normal to the knife edge and n∞ is the index of refraction of air at laboratory conditions. Note that this formulation assumes the density gradient profile to be constant over the integration path; if this is not the case, the term L∂ρ/∂y should be replaced by a corresponding integral. The ultra-high frame rate employed throughout this work allows us to perform spectral analysis on the schlieren data without assuming a propagation speed for the disturbances (as was necessary, for example, in Kennedy et al. 2018). To this end, after converting each image over the steady test duration into a map of density gradient using (2.3), PSD curves were computed for each pixel in the field of view using Welch’s method with Hann windows of length 64. This process then allows us to visualize the spatial distribution of frequency content anywhere within the field of view. The number of frames included in these calculations varied from 2400 to 6000 (75 to 187 realizations), depending on the frame rate of the test. The one exception to this is the condition Re45 test of the 949 A37-6 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow U3 D1 D3 D5 350 360 371 381 391 402 412 423 433 444 454 X (mm) Figure 3. Reference-subtracted image sequence captured for the +0◦ configuration at condition Re33; the inter-image spacing, �t, is 9.7 µs. +15◦ configuration. The latter half of this test was contaminated by anomalous turbulence and, thus, only 1300 images were kept. For reference, the exact number of frames used for each test is reported later in this article in table 4. 3. General behaviour 3.1. Straight (+0◦) continuation We begin by examining the behaviour of the undisturbed boundary layer in the +0◦ configuration through the time-resolved schlieren sequences presented in figures 3 (condition Re33) and 4 (condition Re52). Exemplary schlieren sequences from combinations of flow condition and flare angle not provided here in the main article are included for reference in the online supplementary material (available at https://doi.org/ 10.1017/jfm.2022.769), together with animations of selected combinations. The images have been rotated to align the abscissa with the frustum and X refers to the distance along the cone surface from the nose tip. At the top of each sequence is shown an average flow-on image to highlight the mean boundary-layer profile. Subsequent images are reference subtracted using a 40 image sliding average (corresponding to time intervals of between 49 and 91 µs, depending on the frame rate) to emphasize transient features. Red triangles in the mean flow image denote the locations of PCB sensors, while red bars are used to bracket regions of interest and approximate the propagation of the disturbances as determined visually. Mean PSDs for select PCB sensors at each condition are presented in figure 5. In the first reference-subtracted image of figure 3, distinct, rope-like waves are seen upstream of U3, revealing the presence of a second-mode wavepacket, with an additional wavepacket visible in the vicinity of D3. The upstream wavepacket then appears to intensify as it propagates along the cone surface, and the last three images of the sequence show distortions to its regular, periodic structure, which we interpret as an initial sign of breakdown, i.e. the precursor of transition. While this schlieren sequence is intended to be representative of the overall observed behaviour, the transition process in such 949 A37-7 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 https://doi.org/10.1017/jfm.2022.769 https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence U3 D1 D3 359 369 380 390 401 411 422 432 442 453 X (mm) Figure 4. Reference-subtracted image sequence captured for the +0◦ configuration at condition Re52; �t = 10.2 µs. 100 10–2 10–4 100 10–2 10–4 100 10–2 10–4 200 400 600 200 400 600 200 400 600 U1 U3 D1 D3 D5 Frequency (kHz) Frequency (kHz) Frequency (kHz) P S D ( P a2 H z– 1 ) (a) (b) (c) Figure 5. Power spectral densities at select PCB stations for the +0◦ configuration at (a) condition Re33, (b) condition Re45 and (c) condition Re52. conventional facilities is stochastic by nature and there is occasionally significant variation in wavepacket behaviour, with some breaking down over the frustum and others traversing the entire cone surface without showing obvious distortions to their structure. The PCB spectra shed additional light on these qualitative trends (it should be emphasized that these spectra are representative of the average behaviour, rather than the isolated phenomena revealed in the images). Between U1 and U3, the peak second-mode frequency decreases from approximately 210 kHz to 190 kHz while the maximum N factor increases by 0.6. Additional amplification is observed up to D3, where the maximum N factor has increased by an additional 0.6. The deterioration of the wavepacket features downstream aligns with the amplitude saturation and spectral broadening observed in the spectra for sensors D3 and D5. Note that the 300 kHz peak at station D5 is likely an example of anomalous sensor resonance and not attributable to flow features. 949 A37-8 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow Wavepackets at condition Re52 almost universally exhibit transitional behaviour, which we interpret here to correspond to a clear loss of the regular, periodic structures observed upstream and the onset of a more chaotic, disordered appearance; this is exemplified in figure 4. The highlighted wavepacket in this sequence breaks down as it reaches X = 400 mm, resulting in what appears to be a young turbulent spot over much of the extension. The PCB spectra indicate that mean saturation occurs between U1 and U3, as the fundamental peak (now at 280 kHz) experiences no growth while surrounding frequencies amplify by nearly an order of magnitude. This spectral broadening continues along the length of the model, though the second-mode peak is still apparent at D5, indicating that (on average) the boundary layer is not yet in a fully developed turbulent state. The sensor resonance in D3 and D5 now overlaps with the second-mode peak however, making it difficult to discern the true amplitude of the disturbances (again, the strong resonance peak near 300 kHz should be ignored as unphysical). The behaviour of the boundary-layer disturbances at condition Re45 generally lies between what has been presented for conditions Re33 and Re52, with a decrease in second-mode amplitude between sensors D3 and D5 indicating saturation early on the cone extension. 3.2. The +5◦ flare When compared with the straight extension, the +5◦ compression flare serves largely to promote disturbance growth and transition. Figure 6 presents an image sequence captured at condition Re45 with the +5◦ configuration, again concentrating on an incoming second-mode wavepacket. The field of view for this test is focused primarily on the flare to better capture the downstream development of the boundary-layer disturbances. First, from the mean image, we note that the boundary layer appears fully attached at the corner. When the wavepacket encounters the SWBLI, transient flow structures are radiated away along the expected location of the corner shock (though the shock itself is too weak to be clearly visible): this is indicated by the red arrow in the fourth and fifth images of the sequence. This radiation, apparent for nearly all wavepackets, has a periodicity similar to the second-mode structures and typically appears to emanate from the tail of a wavepacket (this is most clearly seen in movie 1, provided in the online supplementary material). Similar radiation of disturbance energy along weak, mean flow discontinuities appeared in computations performed by Sawaya et al. (2018) of second-mode waves interacting with two-dimensional wall deformations. Downstream of the corner the wavepackets retain their rope-like appearance, but the structure angle of the disturbances decreases such that they appear more parallel to the model surface. In the schlieren sequence, the tail of the wavepacket shows signs of breakdown around X = 440 mm. This occurs as the head of the packet, which has maintained its periodic structure, is leaving the field of view. Amplification of the second mode is seen in the PCB spectra of figure 7(b) until station D3, similar to the straight extension at this condition (Re45). Importantly, however, the broadband amplification observed between sensors D1 and D3 exceeds that of the +0◦ configuration, implying an overall acceleration of the transition process relative to the undisturbed boundary layer. We also note that the PCB spectra along the flare, particularly D3 and D5, demonstrate a shift in the second-mode peak to higher frequencies consistent with the reduced boundary-layer thickness. The behaviour just described for the Re45 condition is also generally representative of condition Re33, though the boundary layer is less transitional along the flare at the lower Rem. This is evident from the reduced spectral broadening seen in the PCB spectra in figure 7(a). At condition Re52 (rightmost plot of figure 7), the PCB spectra along the flare (D1, D3 and D5) show significant amplification in frequencies above and below the 949 A37-9 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence 412 423 433 444 454 464 X (mm) D1 D3 D5 Figure 6. Reference-subtracted image sequence captured for the +5◦ configuration at condition Re45; �t = 7.3 µs. 100 10–2 10–4 100 10–2 10–4 100 10–2 10–4 200 400 600 Frequency (kHz) 200 400 600 Frequency (kHz) 200 400 600 Frequency (kHz) U1 U3 D1 D3 D5 P S D ( P a2 H z– 1 ) (a) (b) (c) Figure 7. Power spectral densities for select PCB stations at the +5◦ configuration at (a) condition Re33, (b) condition Re45 and (c) condition Re52. second-mode peak. It is unclear to what extent the spike in second-mode power at D3 and D5 is caused by sensor resonance, but in any case, the signal has become largely broadband along the flare, indicative of the onset of a young turbulent boundary layer. Indeed, the full test video shows this to be typically the case, with many wavepackets already breaking down upstream of the flare. When contrasted with the straight-cone case, the results here indicate that the attached SWBLI here accelerates boundary-layer transition for unit Reynolds numbers above that of condition Re33. 949 A37-10 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow Re33 Re45 Re52 +10◦Xsep (mm) 391 396 403 +10◦Xatt (mm) 432 428 424 +15◦Xsep (mm) <370 <370 394 +15◦Xatt (mm) 440 433 427 Table 3. Approximate separation and reattachment locations for the +10◦ and +15◦ configurations. 3.3. The +10◦ flare The development of second-mode wavepackets encountering the +10◦ flare was reviewed in detail by Butler & Laurence (2021a). The adverse pressure gradient imposed by this compression is sufficient to separate the boundary layer upstream of the flare for all Reynolds numbers considered. Approximate separation and reattachment locations were identified as points at which the pseudo-streamline profiles reported in § 4 exhibited a sudden change in slope; reasonable agreement (to within 2 mm) was obtained between these results and preliminary numerical simulations (Butler 2021). These locations are represented in each image sequence by vertical, dashed lines and are given in table 3; the uncertainty in each location is estimated as ±3 mm. Figure 8 depicts the typical transitional behaviour of wavepackets which encounter the separated region at condition Re45; a corresponding movie, movie 2, is provided in the online supplementary material. In the mean image, the boundary layer separates between sensor U3 and the corner junction, and reattaches in the vicinity of D3. At the separation point, the instability waves lift off the cone surface and propagate largely within the separated shear layer. The incoming wavepacket can be seen to undergo substantial morphing downstream of separation, with some features propagating near the wall within the recirculation zone. As in the +5◦ case, periodic structures are radiated away when the wavepacket encounters the SWBLI (highlighted by red arrows in the second, third and fourth images of the sequence); however, this now occurs along the expected location of the separation shock. Downstream of reattachment the packet’s structure has become distorted and retained little of its ‘rope-like’ appearance. Such behaviour is not fully representative of condition Re45 however, as the state of the incoming wavepacket may substantially alter its development through the SWBLI. In the first image of the sequence in figure 8, for example, a wavepacket is visible downstream of D3 that has largely retained its structure. The latter behaviour is more representative of wavepackets seen at condition Re33. Beyond the clear second-mode signature upstream of separation (sensors U1 and U3) in the middle image of figure 9, there is evidence of a harmonic developing in the U3 spectra at 500 kHz (a phenomenon also seen at Re33 around 440 kHz). Sensor D1 shows a marked decrease in second-mode power, which may reflect the tendency for the disturbances to lift off the surface; note also, however, the increase in low-frequency content relative to sensor U3, peaking at 85 kHz. A similar peak is seen around 73 kHz at condition Re33 and we have demonstrated previously that these lower-frequency spectral features correspond to distinct, shear-layer disturbances (Butler & Laurence 2021a). The increasingly broadband nature of the spectra of sensors D3 and D5 match the transitional behaviour generally observed in images. At condition Re52, the separation region has shrunk significantly and the instantaneous flow structures on the flare are generally distorted and chaotic in appearance. This is also reflected in the PCB spectra, where even sensor D1 exhibits primarily broadband content. 949 A37-11 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence U3 D1 D3 386 396 406 417 427 438 X (mm) Figure 8. Reference-subtracted image sequence captured for the +10◦ flare at condition Re45; �t = 9.7 µs. The dashed vertical lines indicate the approximate separation and reattachment locations. 100 10–2 10–4 100 10–2 10–4 100 10–2 10–4 P S D ( P a2 H z– 1 ) 200 400 600 Frequency (kHz) 200 400 600 Frequency (kHz) 200 400 600 Frequency (kHz) U1 U3 D1 D3 D5 (a) (b) (c) Figure 9. Representative PSDs at select PCB stations for the +10◦ configuration at (a) condition Re33, (b) condition Re45 and (c) condition Re52. 3.4. The +15◦ flare As evident from the mean flow-on image in figure 10 (captured at condition Re45), the +15◦ compression flare produces a large separation bubble which at conditions Re33 and 949 A37-12 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow U3 D1 D3 380 388 396 405 413 422 430 438 X (mm) Figure 10. Reference-subtracted image sequence captured for the +15◦ configuration at condition Re45; �t = 5.5 µs. The dashed vertical line indicates the approximate reattachment location. Re45 extends far enough upstream to fully envelop sensor U3. Approximate separation and reattachment locations for this configuration are given in table 3. These were determined in the same manner as previously, though the separation location could not be determined at conditions Re33 and Re45 as the boundary layer appeared to separate upstream of the field of view in both cases. For conditions Re33 and Re45, second-mode wavepackets appear only intermittently within the field of view for this flare angle (note that this also made it difficult to define a pseudo-streamline, which will lead to additional uncertainty in the reattachment location) and the dominant transient phenomenon instead appears to be shear-layer disturbances that originate over the separation bubble. These disturbances manifest themselves as long, wavy structures primarily aligned with the direction of propagation: an example may be seen in the instantaneous Re45 image sequence of figure 10 (a corresponding movie, movie 3, can be found in the online supplementary material), which is also representative 949 A37-13 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence 10–2 100 102 10–2 100 102 U1 U3 D1 D3 D5 0 50 100 150 200 0 50 100 150 200 P S D ( P a2 H z– 1 ) Frequency (kHz) Frequency (kHz) (a) (b) Figure 11. Representative PSDs at select PCB stations for the +15◦ configuration at (a) condition Re33 and (b) condition Re45. of the behaviour at Re33. These disturbances appear to undergo substantial intensification once the shear layer begins to recompress from X = 412 mm; then, upon approaching reattachment, they emanate ribbon-like features into the shock layer, as annotated by red arrows in the first three images of the sequence. Near reattachment, the disturbances take on a thin, layered structure that persists downstream on the flare. The final image of figure 10 shows a further ribbon of energy detaching from the shear layer as another shear disturbance approaches the flare. The PCB spectra for conditions Re33 and Re45 are presented in figure 11. Note that excessive sensor resonance contaminated the PCB data at high frequencies for most experiments with this flare angle; thus, the plotted frequency range is limited to 200 kHz and no spectra are presented for Re52. The dominant features for both conditions are peaks present at stations U3 and D1 within the separation bubble which, as we shall see, correspond to the shear-layer disturbances. At condition Re45, this disturbance appears at 90 kHz in the U3 spectrum and proceeds to amplify and shift to lower frequencies downstream, peaking at 50 kHz at station D1. This peak quickly becomes indistinct further along the flare, with rapid spectral broadening occurring between D3 and D5 downstream of reattachment. Similar behaviour is observed at Re33, though the U3 and D1 peaks now occur at 60 kHz and 40 kHz, respectively. Although the appearance of second-mode waves at conditions Re33 and Re45 was rare for this flare angle, several such instances were observed, a typical example of which is presented in figure 7 of the supplementary material. At condition Re52 however, second-mode waves re-emerge as the primary disturbance, with shear-layer waves appearing only intermittently. Several particularly intense bursts of turbulence were observed in this experiment which caused the separation bubble to collapse and reform in a transient process. The online movie 4 depicts such an event, and the image sequence of figure 8 of the supplementary material shows a wavepacket passing through the SWBLI during the recovery process. We may also compute the root-mean-square (r.m.s.) values of the mean-subtracted PCB signals to study the impact of each configuration on the surface pressure fluctuations. The r.m.s. fluctuations over the frequency range 15–280 kHz are shown for the +0◦, 949 A37-14 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow 0 0.1 0.2 0.3 0.4 0.5 0.6 0 0.1 0.2 0.3 0.4 0.5 0.6 0 0.1 0.2 0.3 0.4 0.5 0.6 +0° +5° +10° +15° 1.1 1.2 1.3 1.4 1.5 1.6 1.8 2.0 1.8 2.0 2.2 2.4 Rex (×106) p̃/ p c Rex (×106) Rex (×106) (a) (b) (c) Figure 12. Normalized r.m.s. values of surface pressure fluctuations (15–280 kHz) computed for (a–c) conditions Re33, Re45 and Re52. Filled symbols indicate sensors that were determined to lie beneath a separated-flow region. The dashed line in each case indicates the location of the compression corner: upstream sensors are U1 to U4, while downstream sensors are D1 to D5. +5◦ and +10◦ configurations at each condition in figure 12, normalized by the surface pressure computed from the Taylor–MacColl solution. Filled symbols correspond to sensors determined to be between the separation and reattachment locations in table 3. Note that sensor D2 is only shown for the +10◦ case and that sensor D4 appears to give anomalously low readings in some cases. First, we note that the normalized fluctuation magnitude grows as Rem is increased, which is to be expected based on the overall boundary-layer states noted earlier. Somewhat surprisingly though, we see that all configurations show nearly identical development in fluctuation levels at condition Re33. At condition Re45, the compression flares cause substantial amplification of the surface pressure fluctuation relative to the undisturbed case, though this growth is generally not observed until station D3 (Rex = 1.95 × 106). The slight decrease in p̃ between U3 and U4 for the +10◦ configuration is likely caused by the onset of separation between this sensor pair: we have noted that the second-mode disturbances propagate primarily along the shear layer, and, thus, the separated-flow region may act as a kind of ‘buffer zone’ between the disturbances and the cone surface. For condition Re52, increased fluctuation levels for +5◦ and +10◦ are observed even at station D1; this is to be expected based on the rapid onset of turbulence on the flare observed in the image sequences. These fluctuations peak at station D3 (Rex = 2.27 × 106) on the flare and decrease downstream. Note that, for the +5◦ configuration, the mean inviscid pressure would increase by a factor of 1.9 across the flare shock, relaxing to 1.8 further downstream (corresponding factors for +10◦ are 3.4 and 3.1, and, for +15◦, they are 5.5 and 4.8). Thus, comparing the compression fluctuation levels to those of the +0◦ case, the increase in pressure fluctuations for all flare angles is somewhat smaller than that in the mean pressure. 4. Schlieren spectral analysis 4.1. High-frequency behaviour In figure 13 we present spatial distributions of the integrated disturbance power for each configuration at condition Re33 over the frequency band 170–270 kHz (corresponding to the dominant second-mode frequencies determined from the PCB measurements). 949 A37-15 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence 5 0 5 0 5 0 5 0 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 –9.5 –10 –9.5 –10 –9.5 –10 –9.5 –10 X (mm) Y (m m ) Y (m m ) Y (m m ) Y (m m ) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) Figure 13. Spatial contours of average PSD from 170–270 kHz at condition Re33 for each flare configuration. Separation/reattachment locations are indicated by vertical lines. Note that the data for the +15◦ case was obtained at only 440 kHz, meaning that frequency content between 220–270 kHz would be aliased down to 170–220 kHz. The dashed line in each image traces the wall-normal location of maximum PSD strength within this second-mode frequency range and can be interpreted as a pseudo-streamline along which the disturbances tend to propagate (and was used earlier to determine approximate separation and reattachment locations). The full-range spectra corresponding to each pixel along these pseudo-streamlines are presented in the left column of figure 14. The right column of figure 14 illustrates the change in N factor, which provides a better picture of the local growth rate of disturbances. These N factors are in terms of fluctuations in the magnitude of the density gradient recorded by the schlieren apparatus (with the reference location being the upmost point of the corresponding visualization). If we are to extend these results to the density fluctuations that would typically be of more interest, we must invoke a parallel-flow assumption (Kennedy et al. 2018); this assumption will become questionable across and immediately downstream of the corner and within regions of flow separation, so caution should be exercised in interpreting these results in this way. Note that the abscissa in each of the N-factor plots is the dimensionless stability Reynolds number, R, given by R = √ RemX, (4.1) and the ordinate is also given in terms of the dimensionless frequency F′ = Fδc/Ue, where Ue is the edge velocity calculated from the Taylor–MacColl cone solution. The topmost contour in figure 13 illustrates the behaviour of the straight-cone second-mode content, which saturates and then gradually diminishes in intensity downstream of X = 435 mm. The streamline spectra in figure 14(a) support this interpretation, as evidence of spectral broadening is observed even upstream of the extension. The growth/decay characteristics of the second mode become more interesting when the boundary layer interacts with a compression corner. For the +5◦ compression 949 A37-16 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow 343 368 393 418 442 1100 1150 1200 F re q u en cy ( k H z) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) Fδ c/ U e 0 50 –10 –9.8 –9.6 –9.4 –9.2 150 200 100 250 300 405 419 434 448 462 F re q u en cy ( k H z) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) 0 50 –10.2 –10.0 –9.8 –9.6 –9.4 –9.2 –9.0 –8.8 150 200 100 250 300 0 50 150 200 100 250 300 1160 1180 1200 1220 1240 0 50 150 200 100 250 300 0 0.1 –0.1 –0.2 0.2 0.3 0.4 0.5 0.6 �N �N 0 0.07 0.14 0.22 0.29 0.36 0.43 Fδ c/ U e 0 0.07 0.14 0.22 0.29 0.36 0.43 0 0.2 –0.2 0.4 0.6 0.8 376 390 404 418 432 F re q u en cy ( k H z) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) 0 50 –10.2 –10.0 –9.8 –9.6 –9.4 –9.2 –9.0 –8.8 150 200 100 250 300 1120 1140 1160 1180 1200 0 50 150 200 100 250 300 �N Fδ c/ U e 0 0.07 0.14 0.22 0.29 0.36 0.43 0 0.2 –0.2 0.4 0.6 0.8 R 376 390 405 420 434 x (mm) F re q u en cy ( k H z) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) 0 50 –10.0 –9.5 –9.0 –8.5 –8.0 150 200 100 250 300 1120 1140 1160 1180 1200 0 50 150 200 100 250 300 �N Fδ c/ U e 0 0.07 0.14 0.22 0.29 0.36 0.43 0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 (a) (b) (c) (d) Figure 14. Power spectral densities (left panels) and N-factor distributions (right panels) computed along pseudo-streamlines at condition Re33 for the (a) +0◦, (b) +5◦, (c) +10◦ and (d) +15◦ configurations. The dashed and solid vertical lines indicate corner and (where relevant) separation/reattachment locations. (second contour of figure 13), the wavepackets amplify substantially on the flare, peaking at around X = 450 mm (R = 1220). The pseudo-streamline spectra downstream of this point broaden (figure 14b), correlating well with the wavepacket development observed in the instantaneous flow images and the spectrum of PCB D5. Substantial amplification of 949 A37-17 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence content below 50 kHz is observed in the N-factor plot (figure 14b) and there is evidence of additional disturbance development at 75 kHz. The N-factor plot also demonstrates second-mode amplification at higher frequencies along the flare (F′ � 0.3), matching the PCB observations. For the +10◦ flare (third contour of figure 13), the second-mode energy is seen to amplify along the frustum until the separation point at approximately X = 391 mm, after which it undergoes a rapid decay; this is perhaps linked to the radiation phenomenon noted earlier, though other effects may also be at work. The amplitude of the second-mode fluctuations appears to freeze within the separated shear layer before growing substantially downstream of reattachment alongside low-frequency content near 50 kHz (seen in figure 14c). The pronounced growth in frequencies below 100 kHz which occurs from X = 415–429 mm (R = 1175–1200) is completely absent on the straight cone, suggesting at least one new instability mechanism generated by the separated shear layer. Indeed, the N factors show two distinct bands of growth within the separation bubble: one at 70–80 kHz that peaks within the separation region and decays slightly as the boundary layer reattaches, and another at 30–40 kHz that continues to amplify downstream. The bottom contour of figure 13 (+15◦ flare) shows consistent growth of high-frequency content along the shear layer, in contrast to the +10◦ results. It is worth noting, however, that the streamline spectra in figure 14(d) show no distinct peak in the second-mode range, meaning that this behaviour is caused by broadband amplification rather than modal growth. This is consistent with the relative lack of distinct wavepackets observed in the schlieren images. Instead, low-frequency disturbances can be seen developing as far upstream as X = 390 mm (R = 1140) in the pseudo-streamline spectra, dropping from 65 kHz at onset to 40 kHz at the corner. Amplification of higher frequencies does not begin until downstream of X = 432 mm (R = 1200). Similarly to the +10◦ case, the N-factor contour plot for the +15◦ compression is dominated by two bands of growth, now at approximately 15 kHz and 45 kHz. It is worth noting that the N factors for +15◦ reach much greater values than any other configuration, indicating that the amplification rates for the shear-layer instabilities are higher than those for the second mode. Many of the observations regarding the behaviour of the second mode for condition Re33 also hold at condition Re45. Spatial contours of the second-mode strength (now integrated within the range 200–300 kHz) are depicted in figure 15, and spectra along the pseudo-streamlines are given in figure 16. The straight-cone case again shows consistent amplification of second-mode content leading up to saturation. While the magnitude of the fluctuations has increased compared with Re33, the disturbances do not attain as large a change in N factor due to the earlier onset of saturation. The spatial contour for the +5◦ Re45 configuration (second contour in figure 15) shows high-frequency content reaching off-wall distances greater than the upstream boundary-layer thickness downstream of X ≈ 440 mm. This may be attributed to intermittent turbulent behaviour of the flare boundary layer and is consistent with the dramatic broadband excitation along the pseudo-streamline in figure 16(b); these observations reinforce the conclusion that the SWBLI promotes transition at this condition. Around this same point, the N-factor spectrum shows significant amplification of content below 30 kHz and within bands around 80 kHz and 150 kHz. These bands may potentially correspond to new disturbances, though it should be noted that content at 150 kHz is particularly weak at the upstream edge of the field of view (where the reference power is considered). For the +10◦ case, the second-mode disturbances freeze in amplitude within the upstream portion of the separation bubble (third contour of figure 15), but begin to amplify (weakly) within the downstream portion as they approach reattachment, as in the computations of Novikov et al. (2016). Dramatic spectral broadening is observed 949 A37-18 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow 5 0 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 –9.5 –9.5 –9.0 –9.5 –9.0 –9.8 –9.6 –9.4 –9.2 –10 X (mm) Y (m m ) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) 5 0 Y (m m ) 5 0 Y (m m ) 5 0 Y (m m ) Figure 15. Spatial contours of average PSD from 200–300 kHz at condition Re45 for each flare configuration. Separation/reattachment locations are indicated by vertical lines. in the vicinity of reattachment near X = 430 mm, or R = 1390 (figure 16c). Although it is difficult to discern in the PSD pseudo-streamline spectrum, the N-factor spectrum for this test shows evidence of weakly amplifying disturbances at 80 kHz within the separation bubble from R = 1300–1330, which we have previously shown to correspond to intermittent shear-layer disturbances (Butler & Laurence 2021a). In the +15◦ results we again see a clear shift away from second-mode-dominated behaviour. The bottom-most spatial contour of figure 15 shows amplification of high-frequency content through the entire separation bubble, with the most significant growth again occurring just upstream of reattachment. There also appears to be elevated content along the expected location of the reattachment shock, which was not observed at +10◦ but is likely related to the previously discussed radiation phenomenon. As in the corresponding Re33 case however, the pseudo-streamline spectra in figure 16(d) (both raw and N-factor plots) indicate no significant second-mode peak. Instead, the dominant feature is again the shear-layer instability, which is seen from X ≈ 387 mm (R ≈ 1320) and reduces in frequency from approximately 85 kHz to 50 kHz at the corner. Substantial low-frequency amplification occurs downstream of the corner near the reattachment zone. Switching focus to condition Re52, in figure 17 the spatial PSD contours now correspond to 230–330 kHz and the intermittently transitional nature of the incoming boundary layer has a significant effect in all configurations. For the +0◦ case, this manifests itself as high-frequency content lying significantly above the pseudo-streamline downstream of X ≈ 403 mm (R ≈ 1450), where the disturbance strength now peaks. This also corresponds to the location beyond which significant spectral broadening is seen in the pseudo-streamline spectra of figure 18(a). Just downstream of the +5◦ compression, the high-frequency content amplifies briefly but rapidly before apparently saturating (second contour of figure 17). This is accompanied by nearly instantaneous spectral broadening in figure 18(b), with the power of lower frequencies rising to meet or exceed the second-mode power. The high-frequency content in the +10◦ case grows steadily along the frustum and decays slightly downstream of separation from X ≈ 403–413 mm in the third contour 949 A37-19 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence 360 378 396 414 431 1280 1320 1360 1400 F re q u en cy ( k H z) Fδ c/ U e 0 50 –10.0 –10.2 –9.8 –9.6 –9.4 –8.8 –9.2 –9.0 150 200 100 250 300 405 419 434 448 462 F re q u en cy ( k H z) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) 0 50 –10.0 –9.8 –9.6 –9.4 –9.2 –9.0 –8.8 150 200 100 250 300 0 50 150 200 100 250 300 1360 1380 1400 1420 136013401320 1380 1400 1420 14601440 0 50 150 200 100 250 300 0 0.1 –0.1 –0.2 0.2 0.3 0.4 0.5 0.6 �N �N 0 0.06 0.12 0.19 0.25 0.31 0.37 0 0.2 –0.2 0.4 0.6 1.0 1.2 0.8 �N 0 0.2 –0.2 0.4 0.6 1.0 0.8 380 395 409 423 437 F re q u en cy ( k H z) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) 0 50 –10 –9.8 –9.6 –9.4 –9.2 –9.0 –8.6 –8.8 150 200 100 250 300 0 50 150 200 100 250 300 R 376 390 405 420 434 x (mm) F re q u en cy ( k H z) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) 0 50 –10.0 –9.5 –9.0 150 200 100 250 300 1300 1320 1340 1360 1380 1400 0 50 150 200 100 250 300 �N Fδ c/ U e 0 0.06 0.12 0.19 0.25 0.31 0.37 Fδ c/ U e 0 0.06 0.12 0.19 0.25 0.31 0.37 Fδ c/ U e 0 0.06 0.12 0.19 0.25 0.31 0.37 0 0.5 1.0 1.5 (a) (b) (c) (d) Figure 16. Same as figure 14, but for condition Re45. of figure 17. As the flow reattaches, this content amplifies and spreads out to cover an off-wall distance significantly greater than the upstream boundary-layer thickness, again implying transition. This behaviour is mirrored by the +15◦ configuration, where the turbulent state of the flare boundary layer is even more obvious. Both separated 949 A37-20 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow 5 0 370 380 390 400 410 420 430 440 450 370 380 390 400 410 420 430 440 450 370 380 390 400 410 420 430 440 450 370 380 390 400 410 420 430 440 450 –8.5 –9.4 –9.2 –9.0 –8.8 –9.4 –9.2 –9.0 –8.8 –9.4 –9.2 –9.0 –8.8 –9.0 X (mm) Y (m m ) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) 5 0 Y (m m ) 5 0 Y (m m ) 5 0 Y (m m ) Figure 17. Spatial contours of average PSD from 230–330 kHz at condition Re52 for each flare configuration. Separation/reattachment locations are indicated by vertical lines. configurations display effectively instantaneous spectral broadening at the corner in figures 18(c) and 18(d). The low-frequency signature of the shear-layer mode has largely vanished from the streamline spectra for the +15◦ extension even upstream of this point. 4.2. Low-frequency behaviour The consistent generation of low-frequency content, particularly in the separated cases, justifies a spatial analysis with the frequency range chosen to bracket the potential shear-layer instability observed thus far in the PCB and schlieren spectra. Figure 19 shows the integrated power within this range, from 40 kHz to 80 kHz, for each configuration at condition Re33. Mild growth occurs downstream of the +5◦ corner in the vicinity of the pseudo-streamline, likely associated with the low-frequency N-factor peak shown in figure 14. This is dwarfed however by the growth observed in the separated regions of the +10◦ and +15◦ configurations. For the +10◦ extension, much of the amplification occurs from X ≈ 410–430 mm as the separation bubble is compressing towards reattachment. As mentioned previously, both pseudo-streamline contours (figure 14c) show growth of a distinct 75 kHz band within this region, whereas the amplification seen along the flare is likely due to the 40 kHz peak seen in the N-factor plot. These observations largely hold for the +15◦ case, but the disturbance amplification can be discerned as far upstream as X ≈ 390 mm and reaches a much greater magnitude. There is also evidence of energy emanating away from the shear layer at approximately X = 415 mm, which is consistent with the instantaneous behaviour noted in the reference-subtracted sequences. Figure 20 shows that the +5◦ compression causes greater low-frequency (now 50–90 kHz) amplification at condition Re45 than at condition Re33. The spatial amplification seen along the flare for this configuration most likely corresponds to the N-factor peak at around 80 kHz in figure 16(b), as both features begin to dissipate at approximately X = 455 mm (R = 1430). The +10◦ compression again shows 949 A37-21 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence (c) (a) (b) (d) 375 389 404 418 432 x (mm) 383 395 407 418 430 x (mm) F re q u en cy ( k H z) 0 50 150 200 100 250 300 350 0 50 150 200 100 250 300 350 –10.0 –8.4 –8.6 –8.8 –9.0 –9.2 –9.4 –9.6 –9.8 lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) –9.5 –9.0 –8.0 –8.5 352 373 395 416 437 349 371 392 413 434 F re q u en cy ( k H z) 0 50 150 200 100 250 300 350 0 50 150 200 100 250 300 350–8.6 –8.8 –9.0 –9.2 –9.4 –8.6 –8.8 –9.0 –9.2 –9.4 lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) Figure 18. Power spectral densities along pseudo-streamlines at condition Re52 for the (a) +0◦, (b) +5◦, (c) +10◦ and (d) +15◦ configurations. The dashed and solid vertical lines indicate corner and (where relevant) separation/reattachment locations. 5 0 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 –8.5 –9.5 –9.0 –9.0 –9.0 –9.5 –9.0 –9.5 X (mm) Y (m m ) 5 0 Y (m m ) 5 0 Y (m m ) 5 0 Y (m m ) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) Figure 19. Spatial contours of average PSD from 40–80 kHz at condition Re33 for each flare configuration. Separation/reattachment locations are indicated by vertical lines. 949 A37-22 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow 5 0 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 375 385 395 405 415 425 435 445 455 –8.6 –9.2 –9.0 –8.8 –8.6 –9.2 –9.0 –8.8 –9.2 –9.4 –9.0 –9.5 –9.0 –8.8 X (mm) Y (m m ) 5 0 Y (m m ) 5 0 Y (m m ) 5 0 Y (m m ) lo g 1 0 PS D (( ∂ ρ 2 /∂ y) /H z) Figure 20. Spatial contours of average PSD from 50–90 kHz at condition Re45 for each flare configuration. evidence of the shear-layer instability developing over the separation bubble, particularly within the downstream portion from X ≈ 415–425 mm; this was obscured in the raw pseudo-streamline spectra (figure 16c), partly because the pseudo-streamline is tuned to the second-mode disturbances and does not follow the maximum amplification of low-frequency structures (which we see in figure 20 occurs instead above the pseudo-streamline). This content amplifies further downstream of reattachment, but the pseudo-streamline spectra suggest this amplification to be caused primarily by spectral broadening and transition. The shear-layer instability for the +15◦ configuration in the lowermost contour again develops far upstream of the corner. There is a brief but notable drop in PSD directly above the corner at X = 412 mm, where significant energy is seen emanating from the shear layer. This drop is not captured in figure 16(d) because the maximum low-frequency fluctuations again occur further from the wall than the pseudo-streamline. The growth and slight decay of the shear disturbances illustrated here correlate well with the production terms for the baroclinic instability identified by Dwivedi et al. (2019) for a two-dimensional compression ramp, suggesting a similar mechanism may be at work. Another particularly strong ray of energy originates from the reattachment zone in the spatial contour, at which point the pseudo-streamline spectrum in figure 16(d) shows significant spectral broadening. At condition Re52 (not shown), the low-frequency content experiences transition-related amplification in all configurations. The separated cases no longer show significant evidence of the shear-layer instability, though apparently unrelated low-frequency content does amplify within the separation bubbles. 5. Spectral proper orthogonal decomposition The structure and evolution of disturbances may also be analysed through the SPOD methodology of Towne, Schmidt & Colonius (2018). This technique provides a set of orthogonal modes which oscillate at distinct frequencies and evolve in space and time. These modes may generally be thought of as coherent flow structures, ranked by the total 949 A37-23 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence Extension Condition Nf (frames) Lhann (frames) Realizations Re33 3720 128 56 +0◦ Re45 3800 128 58 Re52 3690 128 56 Re33 4900 256 37 +5◦ Re45 6000 256 45 Re52 4500 256 34 Re33 5000 256 38 +10◦ Re45 5920 256 45 Re52 5000 256 38 Re33 2400 128 36 +15◦ Re45 1300 64 39 Re52 3000 128 45 Table 4. Parameters for SPOD computations. energy they contain. Resolving these modes may allow us to identify the flow structures responsible for the amplification peaks observed in § 4; also, as will be shown in § 6, this technique is particularly useful when coupled with bispectral analysis, as it provides a method to probe the physical structure of frequencies which dominate nonlinear growth mechanisms. In the present implementation, Hann windows of length 64, 128 or 256 were used (depending on the frame rate of the test) with 50 % overlap to group the images of each sequence, resulting in a minimum of 34 flow realizations for a given wind-tunnel run. The specific parameters for each SPOD computation are given in table 4. 5.1. Attached-flow cases As has been established, the boundary layer on the straight-cone configuration is largely dominated by the second-mode instability: this is confirmed by the SPOD analysis. The energies of the SPOD modes computed for each condition are given in figure 21, with the highest-ranked mode indicated by a bold line in each case. Throughout this section, only the highest-ranked mode will be considered for analysis, as it represents the greatest energy content and appears to capture the dominant instability mechanisms (in the left plot of figure 21, for example, the second-mode energy near 200 kHz is much more prominent in the highest-ranked mode than the others). Note that the second-mode peak appears reduced at Re52: this is due to turbulent content which possesses significant flow energy but is less structurally coherent. Eigenvalue contours for select SPOD modes for condition Re33 on the straight cone are given in figure 22. Because they are computed from schlieren images, these eigenvalues correspond to spatially coherent oscillations in the vertical density gradient, ∂ρ/∂y. Note that the colour map for each SPOD mode has been scaled independently to provide the best qualitative visualization of that mode’s coherent flow structures (as is the case throughout). The 184 kHz contours (figure 22b) show the typical rope-like structure of the second-mode disturbances, amplifying from the upstream end of the field of view until dissipating from around X = 400 mm. While these second-mode disturbances are dominant, lower-frequency content is also present. Figure 22(a) shows coherent structures at 34 kHz, with a shallower angle than the second-mode features, amplifying as far upstream as X = 380 mm. It is possible that these structures correspond to the first mode, 949 A37-24 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow Frequency (kHz) Frequency (kHz) Frequency (kHz) E n er g y ( A rb .) 100 200 300 100 200 300 100 200 300 10–4 10–3 10–2 10–1 100 (a) (b) (c) Figure 21. The SPOD mode energy for the +0◦ extension at condition Re33 (left panel), condition Re45 (middle panel) and condition Re52 (right panel). 34 kHz 184 kHz 281 kHz 340 360 380 400 420 440 460 340 360 380 400 420 440 460 340 360 380 400 420 440 460 X (mm) Y (m m ) 0 2.5 Y (m m ) 0 2.5 Y (m m ) 0 2.5 (a) (b) (c) Figure 22. The SPOD mode shapes for the +0◦ configuration at condition Re33 for frequencies of (a) 34 kHz, (b) 184 kHz and (c) 281 kHz. as similar features were observed and speculated to be first-mode waves in Casper et al. (2016) and Kennedy et al. (2022). The large spike in SPOD energy from 250–300 kHz in figure 21 is exemplified by the 280 kHz mode (figure 22c), which corresponds to uniform flickering; this is thus a noise mode from the schlieren apparatus and can be safely ignored. The enhanced growth of lower-frequency content at condition Re52 is demonstrated by the first three contours of figure 23; such features are also generally representative of the low-frequency content at condition Re45. The 32 kHz mode (figure 23a) closely resembles the corresponding one at Re33, while the 59 kHz structures exhibit a similar angle and occupy generally the same region of the boundary layer, suggesting a possible first harmonic. The structures of the 144 kHz mode (figure 23c), meanwhile, closely resemble those of the dominant second mode but at half the frequency and, thus, may represent a subharmonic. The second-mode structures are now observed at higher frequencies and overpower the flickering noise, with the 294 kHz mode (figure 23d) demonstrating strong energy even at the upstream edge of the field of view. All modal structures begin to dissipate or lose coherence from X ≈ 430 mm: this is attributed to the transitional 949 A37-25 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence X (mm) Y (m m ) 0 2.5 Y (m m ) 0 2.5 Y (m m ) 0 2.5 Y (m m ) 0 2.5 (b) (a) (c) (d) 294 kHz 350 360 370 380 390 400 410 420 430 440 450 350 360 370 380 390 400 410 420 430 440 450 350 360 370 380 390 400 410 420 430 440 450 350 360 370 380 390 400 410 420 430 440 450 144 kHz 59 kHz 32 kHz Figure 23. The SPOD mode shapes for the +0◦ configuration at condition Re52 for frequencies of (a) 32 kHz, (b) 59 kHz, (c) 144 kHz and (d) 294 kHz. nature of the boundary layer at such high Rem, which degrades the coherency of the flow structures. The SPOD energy spectra for the +5◦ compression (not shown here) still demonstrate distinct peaks surrounding the second-mode fundamental, but with additional regions of lower-frequency elevation. At condition Re33, for instance, the highest-ranked mode contains substantially more energy than the second highest at around 70–80 kHz; this is also the case near 150 kHz at condition Re45. The mode shapes computed for condition Re33 reveal complex growth of disturbances along the flare, as seen in the contours of figure 24. Some of this growth appears to correspond to modulation of the structures observed in the straight-cone case. For example, the 32 kHz mode (figure 24a) resembles that seen in figure 22(a), tentatively identified as being generated by first-mode waves. Second-mode content at 199 kHz (figure 24f ) can be seen passing through the corner and intensifying until approximately X = 450 mm, matching the saturation point seen in the pseudo-streamline spectrum of figure 14(b). Rope-like structures can also be seen developing just downstream of the corner at 180 kHz (figure 24e) – the growth of such apparent lower-frequency second-mode disturbances is somewhat surprising given the reduced boundary-layer thickness on the flare. Other lower-frequency structures however do not have a clear analogue in the +0◦ case. For example, coherent structures at 77 kHz (figure 24b) appear downstream of the corner and seemingly correspond to a significant N-factor peak along the pseudo-streamline (figure 14b). Another mode shape with substantial energy in this region (125 kHz) is shown in figure 24(c), though the relevant structures are less clearly periodic than those at the other frequencies shown. As will be discussed in § 6, this growth region correlates well with the appearance of nonlinear interactions, particularly for the 77 kHz and 125 kHz structures. An additional mode shape with nonlinear relevance is shown at 161 kHz (figure 24d); this amplifies downstream of the corner but then loses coherence towards the end of the flare. Low-frequency content is once more seen to develop downstream of the corner at condition Re45 in figure 25. The 35 kHz structures again resemble the possible first-mode waves identified in the straight-cone case; the 71 kHz and 100 kHz mode shapes (figures 25b and 25c) develop slightly further downstream but share structural similarities 949 A37-26 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow ( f ) (e) (d) (c) (b) (a) 32 kHz 77 kHz 125 kHz 161 kHz 180 kHz 199 kHz X (mm) Y (m m ) 2.5 0 5.0 Y (m m ) 2.5 0 5.0 Y (m m ) 2.5 0 5.0 Y (m m ) 2.5 0 5.0 Y (m m ) 2.5 0 5.0 Y (m m ) 2.5 0 5.0 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 Figure 24. The SPOD mode shapes for the +5◦ configuration at condition Re33 for frequencies of (a) 32 kHz, (b) 77 kHz, (c) 125 kHz, (d) 161 kHz, (e) 180 kHz and ( f ) 199 kHz. with the 35 kHz mode, and, thus, may represent first and second harmonics. This would reinforce the observation that the attached SWBLI modulates existing disturbances, in this case promoting harmonic development. Particular note should also be taken of the 151 kHz mode, as this corresponds to a significant N-factor peak for the test (figure 16b). As was noted earlier, this elevated N factor may be largely due to the relative weakness of the signal at this frequency upstream (which is reflected by the absence of structures in the SPOD contours); nevertheless, on the flare we see well-defined, coherent structures at this frequency. The SPOD mode shapes within the second-mode band, e.g. at 244 kHz (figure 25e), show evidence of dissipation due to the onset of turbulence downstream of X = 450 mm, and energy can be seen radiating away from the boundary layer along the corner shock. High-frequency content is observed at 382 kHz (figure 25f ) from around X = 420 mm, likely corresponding to transitional structures. The rapid breakdown induced by the compression corner at condition Re52 is evident in the structure of the SPOD modes shown in figure 26. Low-frequency content appears abruptly at the corner (61 kHz and 136 kHz modes, figures 26a and 26b); these structures have a slightly irregular appearance indicative of contamination by intermittent turbulence. The lack of coherent structures along the frustum at 136 kHz is somewhat surprising given the spectral peak exhibited by the U3 PCB sensor in figure 7. The higher-frequency mode shapes associated with the second mode generally show a rapid breakdown in structure 949 A37-27 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence ( f ) (e) (d) (c) (b) (a) 35 kHz 100 kHz 151 kHz 244 kHz 382 kHz 71 kHz Y (m m ) 2.5 0 5.0 Y (m m ) 2.5 0 5.0 Y (m m ) 2.5 0 5.0 Y (m m ) 2.5 0 5.0 Y (m m ) 2.5 0 5.0 Y (m m ) 2.5 0 5.0 X (mm) 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 Figure 25. The SPOD mode shapes for the +5◦ configuration at condition Re45 for frequencies of (a) 35 kHz, (b) 71 kHz, (c) 100 kHz, (d) 151 kHz, (e) 244 kHz and ( f ) 382 kHz. at the flare, which is to be expected based on the immediate spectral broadening and disappearance of the second-mode peak along the pseudo-streamline. Nevertheless, the 270 kHz mode, which corresponds closely to the U3 peak second-mode frequency at this condition, persists for a short distance downstream of the corner; this indicates that some wavepackets traverse the SWBLI without immediately breaking down. 5.2. Separated-flow cases A discussion of the major SPOD mode shapes for the +10◦ extension at the three conditions has been provided in Butler & Laurence (2021b). Here we give a brief overview, focusing on those shapes corresponding to frequencies that exhibit significant nonlinear interactions in § 6. Three mode shapes for condition Re33 are shown in figure 27. The second-mode waves at 215 kHz exhibit the expected rope-like structure, weakening only slightly as they traverse the corner separation, with some evidence of energy passing into the separation bubble. The 42 kHz mode shape consists of long-wavelength features approximately aligned with the shear layer; these amplify significantly immediately downstream of separation and maintain their amplitude through reattachment. This behaviour is consistent with that of the 40 kHz N-factor peak which develops along the pseudo-streamline (figure 14c). Finally, the mode shape representative of the previously 949 A37-28 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow (c) (b) (a) 61 kHz 0 2.5 136 kHz 0 2.5 270 kHz 350 360 370 380 390 400 410 420 430 440 450 350 360 370 380 390 400 410 420 430 440 450 350 360 370 380 390 400 410 420 430 440 450 0 2.5 Y (m m ) Y (m m ) Y (m m ) X (mm) Figure 26. The SPOD mode shapes for the +5◦ configuration at condition Re52 for frequencies of (a) 61 kHz, (b) 136 kHz and (c) 270 kHz. (c) (b) (a) 380 390 400 410 420 430 440 380 390 400 410 420 430 440 380 390 400 410 420 430 440 42 kHz 0 2.5 5.0 71 kHz 0 2.5 5.0 215 kHz 0 2.5 5.0 Y (m m ) Y (m m ) Y (m m ) X (mm) Figure 27. The SPOD mode shapes for the +10◦ configuration at condition Re33 at (a) 42 kHz, (b) 71 kHz and (c) 215 kHz. Separation and reattachment locations are indicated by the dashed vertical lines. noted shear-layer disturbance at 71 kHz is shown in figure 27(b). Figure 28 displays two mode shapes for condition Re45. The second-mode structures at 257 kHz maintain their amplitude downstream of separation, then take on a modified appearance downstream of reattachment. The 61 kHz mode amplifies within the downstream portion of the separation bubble and appears significantly augmented by the onset of transition downstream of reattachment. As evidenced by the streamline spectra of figures 14(d) and 16(d), the large separation bubble created by the +15◦ configuration for conditions Re33 and Re45 continues the trend of shifting the dominant modal content to lower frequencies, with significant content under 70 kHz but little-to-no prominent high-frequency content. This is reflected in the SPOD energy spectra, where the dominant peak at these conditions now corresponds to the shear-layer disturbances at around 40–50 kHz. The highest-rank modes were also elevated slightly around 80–90 kHz, corresponding to the frequency range of the first harmonic of the shear disturbances. 949 A37-29 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence (b) (a) 380 390 400 410 420 430 440 380 390 400 410 420 430 440 61 kHz 0 2.5 5.0 257 kHz 0 2.5 5.0 Y (m m ) Y (m m ) X (mm) Figure 28. The SPOD mode shapes for the +10◦ configuration at condition Re45 at (a) 61 kHz and (b) 257 kHz. (a) Y (m m ) (b) Y (m m ) (c) Y (m m ) (d) Y (m m ) 380 390 400 410 420 430 440 380 390 400 410 420 430 440 380 390 400 410 420 430 440 380 390 400 410 420 430 440 X (mm) 14 kHz 0 2.5 5.0 7.5 41 kHz 0 2.5 5.0 7.5 58 kHz 0 2.5 5.0 7.5 83 kHz 0 2.5 5.0 7.5 Figure 29. The SPOD mode shapes for the +15◦ configuration at condition Re33 for frequencies of (a) 14 kHz, (b) 41 kHz, (c) 58 kHz and (d) 83 kHz. The reattachment location is indicated by the dashed vertical line. The SPOD mode shapes presented in figure 29 depict the structures of the prominent low-frequency disturbances at condition Re33. The 41 kHz mode (figure 29b) closely resembles the shear-layer disturbance identified for the +10◦ configuration (figure 27b) and appears strongest just upstream of reattachment, from X = 420–430 mm. The 58 kHz mode (figure 29c) develops further upstream within the separation region and radiates energy away from the shear layer. Significant energy is also seen in the 83 kHz mode (figure 29d) as the separation region compresses, appearing as longitudinal streaks emanating from the boundary layer at X = 420 mm. It is notable that this corresponds approximately to the harmonic frequency of the shear-layer disturbances. 949 A37-30 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow Very-low-frequency (14 kHz) structures aligned with the pseudo-streamline are observed amplifying rapidly in the immediate vicinity of reattachment in figure 29(a). These latter features likely correspond to the lower band of the N-factor peak in figure 14(d). At condition Re45, the +15◦ configuration is still largely characterized by low-frequency modal structures, which are illustrated by figure 30. The first SPOD contour, computed at 26 kHz, shows features concentrated around the downstream end of the separation bubble and corresponds approximately to the subharmonic of the shear mode; this will be addressed further in § 6. The mode shape at 52 kHz (figure 30b) matches the expected structure of the shear-layer disturbance and again amplifies starting directly above the cone/flare junction, while the 60 kHz mode (figure 30c) shows content developing further upstream with energy radiating away from the shear layer. The angle of these radiating features (which remains approximately constant) is suggestive of frozen Mach radiation, which would indicate that the generating disturbances are moving supersonically with respect to the external flow. To confirm this, in figure 31 we present disturbance propagation speeds computed based on the cross-correlation methodology described by Butler & Laurence (2021b). The cross-correlation was performed along a pseudo-streamline now constructed using the maximum intensity of the shear-layer disturbances, with the intensity time-series bandpass filtered from 40 kHz to 120 kHz. On this plot we also show the edge velocity, Ue, on the straight-cone section, together with the quantity Ue(1 − 1/Me), which is the relative sonic line for disturbances. Note that the actual flow velocity above the separated region will be slightly smaller than Ue, since the fluid in this region has been processed by the separation shock. In any case, the disturbance speed upstream of the corner (Up ≈ 620 ms−1) lies well below the sonic line, indicating that the disturbances here are indeed locally supersonic. Modal content at 103 kHz (figure 30d) is composed almost entirely of ribbon-like waves emanating from the shear layer ahead of reattachment. The angle of these waves is now close to that of the shear layer itself, suggesting more a subsonic ‘flapping’ motion of the generating disturbances rather than the frozen radiation of the lower-frequency modes. This is consistent with the propagation speed seen in this region in figure 31, which is now much closer to the relative sonic line (also remembering that Ue is a slight overestimate). It is notable that the relevant frequency corresponds closely to the first harmonic of the shear-layer disturbances. The potential nonlinear development of these structures will be studied further in the next section. High-frequency (223 kHz) features are observed in figure 30(e) within the severely thinned boundary layer downstream of reattachment. For condition Re52, the mode shapes were generally dominated by turbulent structures, especially on the flare. Second-mode structures were seen upstream of separation, but there was little evidence of shear-layer disturbances. 6. Bispectral analysis The high frame rates employed throughout this work also facilitate the use of higher-order spectral techniques such as bispectral analysis; unlike previous applications of this technique however, which were limited to measurements by individual sensors, the present schlieren technique allows a global picture to be derived. Computation of the normalized bispectrum, or bicoherence, allows one to identify the primary nonlinear growth mechanisms at work. More precisely, the bicoherence provides a measure of the degree of quadratic phase coupling present between three frequencies f1, f2 and f1 + f2. The bispectrum for a frequency triplet is defined by B( f1, f2) = E[X( f1)X( f2)X∗( f1 + f2)], (6.1) 949 A37-31 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence (a) Y (m m ) (b) Y (m m ) (c) Y (m m ) (d) Y (m m ) 380 390 400 410 420 430 440 380 390 400 410 420 430 440 380 390 400 410 420 430 440 380 390 400 410 420 430 440 380 390 400 410 420 430 440 X (mm) 0 2.5 5.0 7.5 0 2.5 5.0 7.5 0 2.5 5.0 7.5 0 2.5 5.0 7.5 (e) Y (m m ) 0 2.5 5.0 7.5 26 kHz 52 kHz 60 kHz 103 kHz 223 kHz Figure 30. The SPOD mode shapes for the +15◦ configuration at condition Re45 for frequencies of (a) 26 kHz, (b) 52 kHz, (c) 60 kHz, (d) 103 kHz and ( f ) 223 kHz. 380 390 400 410 420 430 440 X (mm) U p (m s – 1 ) Ue (1 – 1/Me) Ue 800 600 1000 1200 400 Figure 31. Calculated propagation speed of the shear disturbances for the +15◦ configuration at condition Re45. Here Ue and Me are the boundary-layer edge velocity and Mach number on the straight cone, and the location of the corner is indicated by the dashed vertical line. 949 A37-32 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow +0◦ +5◦ +10◦ +15◦ Condition Re33 Re45 Re33 Re45 Re33 Re45 Re33 Re45 Confidence Level 0.026 0.025 0.020 0.016 0.019 0.016 0.040 0.075 Table 5. The 95 % confidence levels for zero squared bicoherence at each experimental condition. where E is the expectation operator and X( f ) refers to the Fourier transform component at frequency f for the time series of interest. The squared bicoherence is computed as |b( f1, f2)|2 = |B( f1, f2)|2 S( f1)S( f2)S( f1 + f2) , (6.2) which utilizes the normalization method of Brillinger (1965), S( f ) being the PSD at frequency f . There is some disagreement in the literature as to the proper normalization method for bispectral analysis, with some researchers (e.g. Kimmel & Kendall 1991; Craig et al. 2019) opting for the definition given by Kim & Powers (1979). This latter method of normalization however has received scrutiny from Hinich & Wolinsky (2005) for ‘artificially’ bounding the result between 0 and 1 and in some cases can destroy evidence of nonlinearity within experimental datasets. In the context of the present work, the Kim and Powers normalization has been found to ruin the symmetries inherent to the bispectrum and reduce the prominence of expected resonant interactions. There also exists precedent for the use of the Brillinger normalization for boundary-layer transition studies (Chokani 1999; Chokani et al. 2005). Welch’s method is employed for the bispectrum computations using windows of width 64. There is some uncertainty in identifying statistically significant levels of bicoherence. Elgar & Guza (1988) found the 95 % confidence level for zero squared bicoherence to be approximately 3/m, where m is the number of segments. Using this estimator, we arrive at the values in table 5 for the minimum statistically significant squared bicoherence values at each condition. All interactions highlighted for discussion in this work exceed these minimum thresholds. In what follows, we examine nonlinear interactions both in localized regions of the flow field and from a global perspective, with a particular focus on linking interactions to mode shapes observed in the SPOD analysis of the previous section. For brevity, we will use the notation [f1,f2,f1+f2] to refer to a specific bispectral interaction, with the omitted units understood to be kHz. This discussion will concentrate on conditions Re33 and Re45, as the interactions were less prominent and meaningful in the generally transitional-turbulent flow within the region of flow visualization at condition Re52. 6.1. Attached-flow cases We begin by examining the bicoherence spectra observed in the +0◦ configuration to discern the undisturbed nonlinear behaviour of the wavepackets. Representative bicoherence spectra for condition Re33 are shown in figure 32 for three distinct streamwise stations, annotated with red lines in the top mean flow image. Each of these regions corresponds to 21 pixels along the pseudo-streamline for which the bicoherence spectra have been averaged. The dashed line connecting the vertices [Fs/4, Fs/4] and [Fs/2, 0], where Fs is the sampling frequency, separates the inner triangle (IT) from the outer triangle (OT). The OT (located to the right of this line) is distinct in that it is non-redundant but relies upon aliased information, i.e. for any frequency pair [f1,f2] within the OT, 949 A37-33 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence X (mm) f1 (kHz) f 2 ( k H z) 350 100 200 300 f1 (kHz) 100 200 300 f1 (kHz) b2 100 200 300 360 371 381 391 402 412 423 433 444 454 0 100 200 0.05 0 0.10 (a) (b) Figure 32. (a) Averaged image for the +0◦ configuration at condition Re33 with red lines denoting regions for computing the average bicoherence; (b) bicoherence spectra for each region with the contours ordered streamwise from left to right. f1 + f2 will lie above the Nyquist frequency, Fs/2. Information within this region thus requires special scrutiny. At the most upstream station, the dominant nonlinear interaction is fundamental resonance at 210 kHz (i.e. [210,210,420]), as expected based on prior literature. This resonance shifts to lower frequencies downstream to match the fundamental second-mode peak frequency. At the second station, illustrated by the middle contour plot of figure 32, the bicoherence spectrum exhibits elevated levels around [165,30,195]. Such peaks which occur in the vicinity of the second mode and involve low frequencies are typically attributed to sideband interactions, though it is potentially significant that 30 kHz roughly corresponds to the low-frequency (possibly first-mode) structures observed in the SPOD analysis. The final contour of figure 32 demonstrates that these resonant and sideband interactions persist far along the extension, where the instantaneous images have shown intermittent breakdown of wavepackets. We can better visualize the spatial development of these individual interactions as in figure 33, where the r.m.s. bicoherence has been computed within two select regions of the frequency plane (annotated in figure 32 with dashed boxes) for each pixel within the field of view. Figure 33(b) corresponds to fundamental resonance: we see the interaction is generally concentrated around the pseudo-streamline and amplifies along the main cone body. The energy signature of this interaction begins to disperse over the extension as intermittent breakdown occurs, but maintains its presence until the end of the field of view. Note that the downstream part of this contour should be treated with caution, as the resonant peak shifts to frequencies too low to be captured by the considered region far downstream. Figure 33(a) illustrates the spatial distribution of the low-frequency coupling, which appears to be concentrated along the boundary-layer edge; however, it does appear along the pseudo-streamline from X = 415–440 mm, matching the peak intensity of the 34 kHz SPOD structures in figure 22(a). The compression-corner configurations generally lead to the development of new nonlinear exchange mechanisms. Figure 34 demonstrates this point for the +5◦ configuration at condition Re33 (note that the frame rate here has been increased to 822 kHz, which facilitates better resolution of interactions). The most upstream bispectral plot verifies the presence of fundamental resonance near 200 kHz as well as low-frequency coupling in the immediate vicinity of the corner. From the spatial contour of figure 35(e), 949 A37-34 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow 340 360 380 400 420 440 460 340 360 380 400 420 440 460 0 0.05 (b) (a) 0 2 4 0 0.05 b2 0 2 4 X (mm) Y (m m ) Y (m m ) Figure 33. Spatial bicoherence contours on the straight cone at condition Re33 for frequency triplets of (a) [165,30,195] and (b) [190,190,380]. X (mm) (b) (a) f1 (kHz) f 2 ( k H z) 100 200 300 400 f1 (kHz) 100 200 300 400 f1 (kHz) 100 200 300 400 b2 0 100 300 200 0.05 0 0.10 412 423 433 444 454 464 Figure 34. As in figure 32 but for the +5◦ configuration at condition Re33. we see that this resonance enters from upstream of the field of view concentrated along the pseudo-streamline and persists downstream of the corner until approximately X = 430 mm, at which point the interaction dissipates and moves closer to the wall. Indeed, the second streamwise station in figure 34 demonstrates the continued presence of 200 kHz resonance, but it is absent at the third station. Multiple additional resonant interactions develop at the second station of figure 34 around frequencies of 70–80 kHz and 150–160 kHz. In addition to the interaction between these two frequency bands ([160,70,230]), there is a large degree of low-frequency coupling in the vicinity of [160,10,170]. As seen in figures 35(a)–35(c), these interactions peak sharply from X = 420–430 mm. This region corresponds closely to where the 77 kHz and 161 kHz SPOD mode shapes presented in figure 24 first develop, with the 161 kHz structure appearing slightly further upstream. The interaction at [160,70,230] can thus best be explained as deconstructive interference between second-mode wavepackets at 230 kHz and the low-frequency content to produce additional low-frequency structures. Notably, the interaction at [160,10,130] exhibits substantial energy from just upstream of the corner and above the pseudo-streamline; this part of the contour likely corresponds to the corner shock, as mean flow interactions are captured by this bi-frequency region. The resonant interaction at 120 kHz that peaks at the third station in figure 34 has the same general spatial distribution (figure 35d) as the preceding interactions, but shifted downstream. This interaction can be linked to the SPOD structure at 125 kHz (figure 24c), 949 A37-35 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 C.S. Butler and S.J. Laurence X (mm) (a) (b) (c) (d) (e) Y (m m ) Y (m m ) Y (m m ) Y (m m ) Y (m m ) 0 5 0 5 0 5 0 0.05 0.10 0 5 0 0.05 0.10 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 0 5 0 0.05 0.10 0 0.1 0.2 0 0.05 0.10 b2 Figure 35. Spatial bicoherence contours for the +5◦ flare at condition Re33 for frequency triplets of (a) [75,75,150], (b) [160,10,170], (c) [160,70,230], (d) [120,120,240] and (e) [200,200,400]. which peaks in amplitude near X = 440 mm and shows a broadly similar distribution of energy. Spatial bicoherence distributions for several of the notable interactions on the +5◦ flare at condition Re45 are given in figure 36. The first three contours correspond to resonant interactions at 35 kHz, 70 kHz and 230 kHz. The dominant feature for all these interactions is a band of substantially elevated signals lying well above the pseudo-streamline on the flare. This band originates from X = 415–435 mm depending on the frequency and much of this content, especially downstream, likely corresponds to turbulent features. Nevertheless, we can deduce other regions of nonlinear significance in these contours. Resonance of the second-mode fundamental at 230 kHz (figure 36c) is seen along the pseudo-streamline at the beginning of the field of view but rapidly drops off downstream of the corner. The 70 kHz resonance develops slightly upstream of the corner and is also clearly concentrated along the pseudo-streamline until approximately 420 mm. Both the 35 kHz resonance and the [135,10,145] interaction (figure 36d) spike in amplitude dramatically from 410–420 mm, with much of the signal concentrated above the pseudo-streamline; this behaviour of the low-frequency resonance is consistent with the emergence of the SPOD mode in figure 25(b). The energy concentration of the [135,10,145] contour correlates with a locally well-defined SPOD structure at 140 kHz not shown in the preceding discussion and likely corresponds to a mean flow interaction with the corner shock. Synthesizing the results from both conditions, the +5◦ compression causes a localized spike in nonlinear interactions early on the flare, which may facilitate energy transfer and alter the transition process significantly. 949 A37-36 ht tp s: // do i.o rg /1 0. 10 17 /jf m .2 02 2. 76 9 Pu bl is he d on lin e by C am br id ge U ni ve rs ity P re ss https://doi.org/10.1017/jfm.2022.769 Transitional hypersonic cone/flare flow X (mm) (a) (b) (c) (d) Y (m m ) Y (m m ) Y (m m ) Y (m m ) 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 410 420 430 440 450 460 470 0 5 0 5 0 5 0 5 b2 0 0.05 0.10 0 0.05 0.10 0 0.05 0.10 0 0.05 0.10 Figure 36. Spatial bicoherence contours for the +5◦ flare at condition Re45 for frequency triplets of (a) [35,35,70], (b) [70,70,140], (c) [230,230,460] and (d) [135,10,145]. X (mm) (a) (b) (c) Y (m m ) Y (m m ) Y (m m ) b2 0 5 0 0.05 0.10 0 5 0 0.05 0.10 380 390 400 410 420 430 440 380 390 400 410 420 430 440 380 390 400 410 420 430 440 0 5 0 0.05 0.10 Figure 37. Spatial bicoherence contours for the +10◦ flare at condition Re33 for frequency triplets of (a) [45,45,90], (b) [220,20,240] and (c) [220,220,440]. 6.2. Separated-flow cases There was a general decrease in nonlinear interactions for the +10◦ configuration compared with the +5◦ configuration at condition Re33; figure 37 shows the spatial distribution of three of the more prominent interactions. Fundamental resonance at [220,220,440] and low-frequency coupling of the second-mode waves at [220,20,240] are captured upstream of the separation point but quickly die off within the separated region, wh